Summary
Before inverting Moho topography, the traditional Parker-Oldenburg method requires the determination of two important hyperparameters, the average Moho depth and Moho density contrast. The selection of these two hyperparameters will directly affect the inversion results. In this paper, a new method for estimating hyperparameters is proposed which is used to improve the Parker-Oldenburg method. The new method is improved by using simulated annealing to accurately estimate the average Moho depth and Moho density contrast based on the relationship between Moho depths and corresponding gravity anomalies at seismic control points. Synthetic tests show that compared to the improved Bott's method and the trial and error method, our method reduces the error in Moho density contrast and average Moho depth by 0.83% and 1.81% respectively. In addition, compared with the trial and error method, our method greatly improves the computational efficiency. In a practical example, we apply this method to invert the Moho topography in the northern South China Sea. The inversion results show that the Moho topography in the northern South China Sea ranges from 8.2 to 33 km. The root mean squared error between our Moho topography and the seismic validation points is 0.94 km. Compared with the CRUST 1.0 model, our Moho topography is more accurate.