Moholso: A MATLAB program to determine crustal thickness by an isostatic and a global gravitational model
2012 (English)In: Computers & Geosciences, ISSN 0098-3004, E-ISSN 1873-7803, Vol. 44, 177-183 p.Article in journal (Refereed) Published
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 Veiling 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 (Moholso.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.
Place, publisher, year, edition, pages
2012. Vol. 44, 177-183 p.
Inversion method, Isostasy, MATLAB, Moho depth, Mohorovicic discontinuity, Vening Meinesz-Moritz hypothesis
IdentifiersURN: urn:nbn:se:kth:diva-99399DOI: 10.1016/j.cageo.2011.10.012ISI: 000306034100019ScopusID: 2-s2.0-84861910399OAI: oai:DiVA.org:kth-99399DiVA: diva2:542294
QC 201207312012-07-312012-07-302012-07-31Bibliographically approved