Details
| Original language | English |
|---|---|
| Pages (from-to) | 207-226 |
| Number of pages | 20 |
| Journal | Studia geophysica et geodaetica |
| Volume | 58 |
| Issue number | 2 |
| Publication status | Published - 10 Feb 2014 |
Abstract
A two-step free-positioned point mass method is used for regional gravity field modeling together with the remove-compute-restore (RCR) technique. The Quasi-Newton algorithm (L-BFGS-B) is implemented to solve the nonlinear problem with bound constraints in the first step, while in the second step the magnitudes of the point masses are re-adjusted with known positions in the least-squares sense. In order to reach a good representation of the gravity field, a number of parameter sets have to be defined carefully before the computations. The effects of four important parameter sets (depth limits, number of point masses, original/reduced basis functions and optimization directions) are investigated for regional gravity field modeling based on two numerical test cases with synthetic and real data. The results show that the selection of the initial depth and depth limits is of most importance. The number of point masses for obtaining a good fit is affected by the data distribution, while a dependency on the data variability (signal variation) is negligible. Long-wavelength errors in the predicted height anomalies can be reduced significantly by using reduced basis functions, and the radial-direction optimization proves to be stable and reliable for regular and irregular data scenarios. If the parameter sets are defined properly, the solutions are similar to the ones computed by least-squares collocation (LSC), but require fewer unknowns than LSC.
Keywords
- depth limits, free-positioned point masses, regional gravity field modeling
ASJC Scopus subject areas
- Earth and Planetary Sciences(all)
- Geophysics
- Earth and Planetary Sciences(all)
- Geochemistry and Petrology
Cite this
- Standard
- Harvard
- Apa
- Vancouver
- BibTeX
- RIS
In: Studia geophysica et geodaetica, Vol. 58, No. 2, 10.02.2014, p. 207-226.
Research output: Contribution to journal › Article › Research › peer review
}
TY - JOUR
T1 - Regional gravity field modeling using free-positioned point masses
AU - Lin, Miao
AU - Denker, Heiner
AU - Müller, Jürgen
N1 - Funding information:: The first author was financially supported by the China Scholarship Council (CSC).
PY - 2014/2/10
Y1 - 2014/2/10
N2 - A two-step free-positioned point mass method is used for regional gravity field modeling together with the remove-compute-restore (RCR) technique. The Quasi-Newton algorithm (L-BFGS-B) is implemented to solve the nonlinear problem with bound constraints in the first step, while in the second step the magnitudes of the point masses are re-adjusted with known positions in the least-squares sense. In order to reach a good representation of the gravity field, a number of parameter sets have to be defined carefully before the computations. The effects of four important parameter sets (depth limits, number of point masses, original/reduced basis functions and optimization directions) are investigated for regional gravity field modeling based on two numerical test cases with synthetic and real data. The results show that the selection of the initial depth and depth limits is of most importance. The number of point masses for obtaining a good fit is affected by the data distribution, while a dependency on the data variability (signal variation) is negligible. Long-wavelength errors in the predicted height anomalies can be reduced significantly by using reduced basis functions, and the radial-direction optimization proves to be stable and reliable for regular and irregular data scenarios. If the parameter sets are defined properly, the solutions are similar to the ones computed by least-squares collocation (LSC), but require fewer unknowns than LSC.
AB - A two-step free-positioned point mass method is used for regional gravity field modeling together with the remove-compute-restore (RCR) technique. The Quasi-Newton algorithm (L-BFGS-B) is implemented to solve the nonlinear problem with bound constraints in the first step, while in the second step the magnitudes of the point masses are re-adjusted with known positions in the least-squares sense. In order to reach a good representation of the gravity field, a number of parameter sets have to be defined carefully before the computations. The effects of four important parameter sets (depth limits, number of point masses, original/reduced basis functions and optimization directions) are investigated for regional gravity field modeling based on two numerical test cases with synthetic and real data. The results show that the selection of the initial depth and depth limits is of most importance. The number of point masses for obtaining a good fit is affected by the data distribution, while a dependency on the data variability (signal variation) is negligible. Long-wavelength errors in the predicted height anomalies can be reduced significantly by using reduced basis functions, and the radial-direction optimization proves to be stable and reliable for regular and irregular data scenarios. If the parameter sets are defined properly, the solutions are similar to the ones computed by least-squares collocation (LSC), but require fewer unknowns than LSC.
KW - depth limits
KW - free-positioned point masses
KW - regional gravity field modeling
UR - http://www.scopus.com/inward/record.url?scp=84898546247&partnerID=8YFLogxK
U2 - 10.1007/s11200-013-1145-7
DO - 10.1007/s11200-013-1145-7
M3 - Article
AN - SCOPUS:84898546247
VL - 58
SP - 207
EP - 226
JO - Studia geophysica et geodaetica
JF - Studia geophysica et geodaetica
SN - 0039-3169
IS - 2
ER -