MohoIso: A MATLAB program to determine crustal thickness by an isostatic and a global gravitational model
- 作者：Mohammad Bagherbandi mohbag@kth.se
- 关键词：Inversion method ; Isostasy ; MATLAB ; Moho depth ; Mohorovi膷i膷 discontinuity ; Vening Meinesz&ndash ; Moritz hypothesis
- 刊名：Computers and Geosciences
- 出版年：2012
- 期刊代码：16_00983004
- 类别：ear
- 出版时间：July, 2012
- 卷：44
- 期：Complete
- 页码：177-183
- 文件大小：829 K
摘要
This paper focuses on the modeling of the boundary between Earth's crust and upper mantle using a gravimetric-isostatic model. Here a MATLAB code is presented based on the gravimetric-isostatic model i.e. the Vening Meinesz-Moritz model. Inverse problems in isostasy consist in making the isostatic anomalies to be zero under a certain isostatic hypothesis. The Vening Meinesz-Moritz problem is to determine the Moho depth such that the compensating attraction totally compensates the Bouguer gravity anomaly on the Earth's surface, implying that the isostatic anomaly vanishes on the Earth's surface. The main idea is easy but the theoretical analysis is somewhat difficult. Here a practical method to recover the Moho depth from the gravity data is used in the MATLAB code (MohoIso.m) based on the Vening Meinesz-Moritz method. The code has been designed based on different sub-codes. The body of the main code works according to the vectorization technique, because this technique causes that the speed of code increases. One of the important possible limitations for the code is over-flow and under-flow for higher degrees in the fully normalized associated Legendre function. This problem occurs in the subroutine applied in this study, it limits the numerical study up to degrees 1800-2000.