Difference between pages "JSG T.26" and "JSG T.31"

From Icctwiki
(Difference between pages)
Jump to: navigation, search
(Created page with "<big>'''JSG 0.13: Integral equations of potential theory for continuation and transformation of classical and new gravitational observables'''</big> Chair:''Michal Šprlák (...")
 
(Created page with "<big>'''JSG 0.18: High resolution harmonic analysis and synthesis of potential fields'''</big> Chair: ''Sten Claessens (Australia)''<br> Affiliation:''Comm. 2 and GGOS'' __T...")
 
Line 1: Line 1:
<big>'''JSG 0.13: Integral equations of potential theory for continuation and transformation of classical and new gravitational observables'''</big>
+
<big>'''JSG 0.18: High resolution harmonic analysis and synthesis of potential fields'''</big>
  
Chair:''Michal Šprlák (Czech Republic)''<br>
+
Chair: ''Sten Claessens (Australia)''<br>
Affiliation:''Commission 2 and GGOS''
+
Affiliation:''Comm. 2 and GGOS''
  
 
__TOC__
 
__TOC__
  
===Introduction===
+
===Terms of Reference===
  
The description of the Earth's gravitational field and its temporal variations belongs to fundamental pillars of modern geodesy. The accurate knowledge of the global gravitational field is important in many applications including precise positioning, metrology, geophysics, geodynamics, oceanography, hydrology, cryospheric and other geosciences. Various observation techniques for collecting gravitational data have been invented based on terrestrial, marine, airborne and more recently, satellite sensors. On the other hand, different parametrization methods of the gravitational field were established in geodesy, however, with many unobservable parameters. For this reason, the geodetic science has traditionally been formulating various gravitational parameter transformations, including those based on solving boundary/initial value problems of potential theory, through Fredholm's integral equations.
+
The gravitational fields of the Earth and other celestial bodies in the Solar System are customarily represented by a series of spherical harmonic coefficients. The models made up of these harmonic coefficients are used widely in a large range of applications within geodesy. In addition, spherical harmonics are now used in many other areas of science such as geomagnetism, particle physics, planetary geophysics, biochemistry and computer graphics, but one of the first applications of spherical harmonics was related to the gravitational potential, and geodesists are still at the forefront of research into spherical harmonics. This holds true especially when it comes to the extension of spherical harmonic series to ever higher degree and order (d/o).
 +
The maximum d/o of spherical harmonic series of the Earth’s gravitational potential has risen steadily over the past decades. The highest d/o models currently listed by the International Centre for Global Earth Models (ICGEM) have a maximum d/o of 2190. In recent years, spherical harmonic models of the topography and topographic potential to d/o 10,800 have been computed, and with ever-increasing computational prowess, expansions to even higher d/o are feasible. For comparison, the current highest-resolution global gravity model has a resolution of 7.2” in the space domain, which is roughly equivalent to d/o 90,000 in the frequency domain, while the highest-resolution global Digital Elevation Model has a resolution of 5 m, equivalent to d/o ~4,000,000.
  
Traditionally, Stokes’s, Vening-Meinesz’s and Hotine’s integrals have been of interest in geodesy as they accommodated geodetic applications. In recent history, new geodetic integral transformations were formulated. This effort was mainly initiated by new gravitational observables that became gradually available to geodesists with the advent of precise GNSS (Global Navigation Satellite Systems) positioning, satellite altimetry and aerial gravimetry/gradiometry. The family of integral transformations has enormously been extended with satellite-to-satellite tracking and satellite gradiometric data available from recent gravity-dedicated satellite missions.
+
The increasing maximum d/o of harmonic models has posed and continues to pose both theoretical and practical challenges for the geodetic community. For example, the computation of associated Legendre functions of the first kind, which are required for spherical harmonic analysis and synthesis, is traditionally subject to numerical instabilities and underflow/overflow problems. Much progress has been made on this issue by selection of suitable recurrence relations, summation strategies, and use of extended range arithmetic, but further improvements to efficiency may still be achieved.  
  
Besides numerous efforts in developing integral equations to cover new observables in geodesy, many aspects of integral equations remain challenging. This study group aims for systematic treatment of integral transformation in geodesy, as many formulations have been performed by making use of various approaches. Many solutions are based on spherical approximation that cannot be justified for globally distributed satellite data and with respect to requirements of various data users requiring gravitational data to be distributed the reference ellipsoid or at constant geodetic altitude. On the other hand, the integral equations in spherical approximation possess symmetric properties that allow for studying their spatial and spectral properties; they also motivate for adopting a generalized notation. New numerically efficient, stable and accurate methods for upward/downward continuation, comparison, validation, transformation, combination and/or for interpretation of gravitational data are also of high interest with increasing availability of large amounts of new data.
+
There are further separate challenges in ultra-high d/o harmonic analysis (the forward harmonic transform) and synthesis (the inverse harmonic transform). Many methods for the forward harmonic transform exist, typically separated into least-squares and quadrature methods, and further comparison between the two at high d/o, including studying the influence of aliasing, is of interest. The inverse harmonic transform, including synthesis of a large variety of quantities, has received much interest in recent years. In moving towards higher d/o series, highly efficient algorithms for synthesis on irregular surfaces and/or in scattered point locations, are of utmost importance.
 +
 
 +
Another question that has occupied geodesists for many decades is whether there is a substantial benefit to the use of oblate ellipsoidal (or spheroidal) harmonics instead of spherical harmonics. The limitations of the spherical harmonic series for use on or near the Earth’s surface are becoming more and more apparent as the maximum d/o of the harmonic series increase. There are still open questions about the divergence effect and the amplification of the omission error in spherical and spheroidal harmonic series inside the Brillouin surface.
 +
 
 +
The Hotine-Jekeli transformation between spherical and spheroidal harmonic coefficients has proven very useful, in particular for spherical harmonic analysis of data on a reference ellipsoid. It has recently been improved upon and extended, while alternatives using surface spherical harmonics have also been proposed, but the performance of the transformations at very high d/o may be improved further. Direct use of spheroidal harmonic series requires (ratios of) associated Legendre functions of the second kind, and their stable and efficient computation is also of ongoing interest.
  
 
===Objectives===
 
===Objectives===
  
* To consider different types of gravitational data, i.e., terrestrial, aerial and satellite, available today and to formulate their mathematical relation to the gravitational potential.
+
The objectives of this study group are to:
* To study mathematical properties of differential operators in spherical and Jacobi ellipsoidal coordinates, which relate various functionals of the gravitational potential.
+
* Create and compare stable and efficient methods for computation of ultra-high degree and order associated Legendre functions of the first and second kind (or ratios thereof), plus its derivatives and integrals.
* To complete the family of integral equations relating various types of current and foreseen gravitational data and to derive corresponding spherical and ellipsoidal Green’s functions.
+
* Study the divergence effect of ultra-high degree spherical and spheroidal harmonic series inside the Brillouin sphere/spheroid.
* To study accurate and numerically stable methods for upward/downward continuation of gravitational field parameters.
+
* Verify the numerical performance of transformations between spherical and spheroidal harmonic coefficients to ultra-high degree and order.
* To investigate optimal combination techniques of heterogeneous gravitational field observables for gravitational field modelling at all scales.
+
* Compare least-squares and quadrature approaches to very high-degree and order spherical and spheroidal harmonic analysis.
* To investigate conditionality as well as spatial and spectral properties of linear operators based on discretized integral equations.
+
* Study efficient methods for ultra-high degree and order harmonic analysis (the forward harmonic transform) for a variety of data types and boundary surfaces.
* To classify integral transformations and to propose suitable generalized notation for a variety of classical and new integral equations in geodesy.
+
* Study efficient methods for ultra-high degree and order harmonic synthesis (the inverse harmonic transform) of point values and area means of all potential quantities of interest on regular and irregular surfaces.
  
===Program of Activities===
+
===Program of activities===
  
* Presenting research results at major international geodetic and geophysical conferences, meetings and workshops.
+
* Providing a platform for increased cooperation between group members, facilitating and encouraging exchange of ideas and research results.
* Organizing a session at the forthcoming Hotine-Marussi Symposium 2017.
+
* Creating and updating a bibliographic list of relevant publications from both the geodetic community as well as other disciplines for the perusal of group members.
* Cooperating with related IAG Commissions and GGOS.
+
* Organizing working meetings at international symposia and presenting research results in the appropriate sessions.
* Monitoring activities of JGS members as well as other scientists related to the scope of JGS activities.
 
* Providing bibliographic list of relevant publications from different disciplines in the area of JSG interest.
 
  
===Members===
+
===Membership===
  
'' '''Michal Šprlák (Czech Republic), chair''' <br /> Alireza Ardalan (Iran) <br /> Mehdi Eshagh (Sweden) <br /> Will Featherstone (Australia) <br /> Ismael Foroughi (Canada) <br /> Petr Holota (Czech Republic) <br /> Juraj Janák (Slovakia) <br /> Otakar Nesvadba (Czech Republic) <br /> Pavel Novák (Czech Republic) <br /> Martin Pitoňák (Czech Republic) <br /> Robert Tenzer (China) <br /> Guyla Tóth (Hungary) <br />''
+
'' '''Sten Claessens (Australia), chair''' <br /> Hussein Abd-Elmotaal (Egypt) <br /> Oleh Abrykosov (Germany) <br /> Blažej Bucha (Slovakia) <br /> Toshio Fukushima (Japan) <br /> Thomas Grombein (Germany) <br /> Christian Gruber (Germany) <br /> Eliška Hamáčková (Czech Republic) <br /> Christian Hirt (Germany) <br /> Christopher Jekeli (USA) <br /> Otakar Nesvadba (Czech Republic) <br /> Moritz Rexer (Germany) <br /> Josef Sebera (Czech Republic) <br /> Kurt Seitz (Germany) <br />''

Revision as of 10:01, 29 April 2016

JSG 0.18: High resolution harmonic analysis and synthesis of potential fields

Chair: Sten Claessens (Australia)
Affiliation:Comm. 2 and GGOS

Terms of Reference

The gravitational fields of the Earth and other celestial bodies in the Solar System are customarily represented by a series of spherical harmonic coefficients. The models made up of these harmonic coefficients are used widely in a large range of applications within geodesy. In addition, spherical harmonics are now used in many other areas of science such as geomagnetism, particle physics, planetary geophysics, biochemistry and computer graphics, but one of the first applications of spherical harmonics was related to the gravitational potential, and geodesists are still at the forefront of research into spherical harmonics. This holds true especially when it comes to the extension of spherical harmonic series to ever higher degree and order (d/o). The maximum d/o of spherical harmonic series of the Earth’s gravitational potential has risen steadily over the past decades. The highest d/o models currently listed by the International Centre for Global Earth Models (ICGEM) have a maximum d/o of 2190. In recent years, spherical harmonic models of the topography and topographic potential to d/o 10,800 have been computed, and with ever-increasing computational prowess, expansions to even higher d/o are feasible. For comparison, the current highest-resolution global gravity model has a resolution of 7.2” in the space domain, which is roughly equivalent to d/o 90,000 in the frequency domain, while the highest-resolution global Digital Elevation Model has a resolution of 5 m, equivalent to d/o ~4,000,000.

The increasing maximum d/o of harmonic models has posed and continues to pose both theoretical and practical challenges for the geodetic community. For example, the computation of associated Legendre functions of the first kind, which are required for spherical harmonic analysis and synthesis, is traditionally subject to numerical instabilities and underflow/overflow problems. Much progress has been made on this issue by selection of suitable recurrence relations, summation strategies, and use of extended range arithmetic, but further improvements to efficiency may still be achieved.

There are further separate challenges in ultra-high d/o harmonic analysis (the forward harmonic transform) and synthesis (the inverse harmonic transform). Many methods for the forward harmonic transform exist, typically separated into least-squares and quadrature methods, and further comparison between the two at high d/o, including studying the influence of aliasing, is of interest. The inverse harmonic transform, including synthesis of a large variety of quantities, has received much interest in recent years. In moving towards higher d/o series, highly efficient algorithms for synthesis on irregular surfaces and/or in scattered point locations, are of utmost importance.

Another question that has occupied geodesists for many decades is whether there is a substantial benefit to the use of oblate ellipsoidal (or spheroidal) harmonics instead of spherical harmonics. The limitations of the spherical harmonic series for use on or near the Earth’s surface are becoming more and more apparent as the maximum d/o of the harmonic series increase. There are still open questions about the divergence effect and the amplification of the omission error in spherical and spheroidal harmonic series inside the Brillouin surface.

The Hotine-Jekeli transformation between spherical and spheroidal harmonic coefficients has proven very useful, in particular for spherical harmonic analysis of data on a reference ellipsoid. It has recently been improved upon and extended, while alternatives using surface spherical harmonics have also been proposed, but the performance of the transformations at very high d/o may be improved further. Direct use of spheroidal harmonic series requires (ratios of) associated Legendre functions of the second kind, and their stable and efficient computation is also of ongoing interest.

Objectives

The objectives of this study group are to:

  • Create and compare stable and efficient methods for computation of ultra-high degree and order associated Legendre functions of the first and second kind (or ratios thereof), plus its derivatives and integrals.
  • Study the divergence effect of ultra-high degree spherical and spheroidal harmonic series inside the Brillouin sphere/spheroid.
  • Verify the numerical performance of transformations between spherical and spheroidal harmonic coefficients to ultra-high degree and order.
  • Compare least-squares and quadrature approaches to very high-degree and order spherical and spheroidal harmonic analysis.
  • Study efficient methods for ultra-high degree and order harmonic analysis (the forward harmonic transform) for a variety of data types and boundary surfaces.
  • Study efficient methods for ultra-high degree and order harmonic synthesis (the inverse harmonic transform) of point values and area means of all potential quantities of interest on regular and irregular surfaces.

Program of activities

  • Providing a platform for increased cooperation between group members, facilitating and encouraging exchange of ideas and research results.
  • Creating and updating a bibliographic list of relevant publications from both the geodetic community as well as other disciplines for the perusal of group members.
  • Organizing working meetings at international symposia and presenting research results in the appropriate sessions.

Membership

Sten Claessens (Australia), chair
Hussein Abd-Elmotaal (Egypt)
Oleh Abrykosov (Germany)
Blažej Bucha (Slovakia)
Toshio Fukushima (Japan)
Thomas Grombein (Germany)
Christian Gruber (Germany)
Eliška Hamáčková (Czech Republic)
Christian Hirt (Germany)
Christopher Jekeli (USA)
Otakar Nesvadba (Czech Republic)
Moritz Rexer (Germany)
Josef Sebera (Czech Republic)
Kurt Seitz (Germany)