150 likes | 157 Views
This study explores methods for combining various gravity-related datasets to compute accurate geoid models in Texas. Preliminary results show promising outcomes, but further work is needed to refine parameters and validate the results against independent datasets.
E N D
Data fusion for geoid computation – numerical tests in Texas area(preliminary results) Yan Ming Wang & Xiaopeng Li National Geodetic Survey, NOAA, USA International Symposium on Gravity, Geoid and Height Systems October 9-12, 2012Venice, Italy
Outline Objective of study Combination methods used Preliminary results Future work
Objective of study Many types of gravity related data sets available today: - satellite gravity models (long wavelength) - Airborne gravity (medium wavelength) - Surface gravity (short wavelength) - High resolution digital elevation and density models (Ultra short wavelength) - Geopotential numbers from leveling (short wavelength) - Deflections of the vertical (short wavelength) Goal: To combine data in an optimal way (old topic, but a new challenge: can it reach 1 cm geoid accuracy?)
Few words on the use of geopotential numbers C (1) The ellipsoidal height is NOT known at most historical leveling benchmarks, so that the disturbing potential can not be computed directly for all benchmarks To use the geopotential number C in determination of the gravity field, we compute initial disturbing potential as T0 = W0 – C – (U0 – γ0•h0) ≈ - C + γ0•h0 andh0 = H* + ζ = H + N
Few words on the use of geopotential numbers C (2) Gravity anomaly can be computed from the geopotential numbers
Combination methods(1) Over determined boundary values problems Least squares collocation Observation equation: l ... observation L … linear operator n … observation error The least squares solution:
Combination methods(2) Use of harmonic (Eigen) functions Task: using observations l to determine coefficients . Solution in matrix form:
Geopotential numbers The gravity anomaly dg can be computed from geopotential numbers as D(T0 )/Dh = g – γ0 = dg + (γ0 - γQ) γQnormal gravity on the telluroid.
Conclusions Surface data show rich high frequency of gravity field, but it shows also systematical difference from the airborne data. Airborne data is heavily smoother to remove high frequency dynamic errors. Preliminary results show the combination results are promising.
Future Work Using truncated Stokes kernel to let satellite model controls the long wavelength Refine parameters used in both methods Include the geopotential numbers in the computations – if it helps to improve the solution Compare the results against independent data sets, such GSVS11 GPS/leveling data (1cm relative geoid accuracy along 300km line)