Incorporating topographic effects in surface wave tomography based on shortest-path ray tracing

Nanqiao Du,Tingwei Yang,Tao Xu,Qinya Liu
DOI: https://doi.org/10.1093/gji/ggae115
IF: 3.352
2024-03-07
Geophysical Journal International
Abstract:SUMMARY We propose a new method to obtain 3-D shear wave velocity structures for regions with undulating topography based on surface wave dispersion data. In our method, we assume that surface waves propagate along a curved free surface, and the dispersion effects of these waves can be modelled as a frequency-dependent ray tracing problem. We use the shortest-path method to address off-great-circle propagation arising from inhomogeneous elastic parameters and topographic variations. We then apply our method to both synthetic and real-data inversions and demonstrate that ignoring topographic effects may significantly distort the inverted images. Finally, we analysed the accuracy of our method and provided a rule-of-thumb principle to quantitatively assess the need to account for topographic effects in surface wave tomography for a region.
geochemistry & geophysics
What problem does this paper attempt to address?
### Problems the paper attempts to solve This paper aims to solve the problem of how to accurately consider the topographic effect when performing surface - wave tomography in areas with undulating terrain. Specifically, the paper proposes a new method to simulate the propagation path of surface waves under complex topographic conditions by the Shortest - Path Method (SPM), and based on this method, perform a three - dimensional shear - wave velocity structure inversion of surface - wave dispersion data. This method can more accurately handle the influence of topographic changes on the propagation of surface waves, thereby improving the imaging accuracy. ### Background and motivation 1. **Importance of surface - wave tomography**: - Surface - wave tomography is an important tool for studying the dynamic processes in the shallow and deep layers of the earth, especially using surface - wave data recorded from ambient noise signals, which can be obtained through dense seismic arrays. - Successful applications include imaging studies at regional, continental, and global scales. 2. **Limitations of existing methods**: - Existing surface - wave tomography methods usually assume that the topographic scale is much smaller than the wavelength and penetration depth, and thus ignore the topographic effect. - However, with the emergence of short - period nodal arrays and distributed acoustic sensors, the wavelength of short - period surface waves may be comparable to the topographic scale, especially in areas with significant topographic changes (such as the Qinghai - Tibet Plateau, onshore/offshore, etc.), and the topographic effect cannot be ignored. 3. **Influence of topographic effect**: - Ignoring the topographic effect may lead to significant distortion of the inversion image. - For example, Köhler et al. (2012) showed that a 2 - km topographic change along a 200 - km profile would cause a 0.7% underestimation of the 3 - second - period phase velocity. ### Method overview 1. **Surface - wave WKBJ theory**: - Adopt the surface - wave WKBJ theory, assuming that the surface wave propagates along a curved surface with a rapidly changing phase term and a slowly changing amplitude. - Derive the local dispersion relation and the corresponding Hamilton - Jacobi equation. 2. **Shortest - Path Method (SPM)**: - SPM is a ray - tracing method based on graph theory and can flexibly handle irregular surfaces or discontinuities. - By discretizing the free - form surface, use quadrilateral elements to adapt to topographic changes, and add auxiliary nodes at the edges of each element to increase the ray coverage. - Use the Dijkstra algorithm to systematically advance the wavefront and update the travel time of nodes. - Automatically determine the ray path from the source point to any node. 3. **Fréchet kernel**: - Calculate the gradient of the travel time with respect to the current dispersion velocity to obtain the discretized Fréchet kernel. - Represent the relationship between the perturbation of the travel time and the perturbation of the ray velocity in matrix form. 4. **Three - dimensional direct inversion**: - Adopt the direct inversion method to relate the one - dimensional phase/group velocity change to the perturbation of elastic parameters. - Combine with the Thomson - Haskell propagation matrix algorithm to obtain the three - dimensional Fréchet kernel. - Solve the linear equations by adding a regularization term and update the current model. ### Experimental verification 1. **Synthetic tests**: - Conduct synthetic tests on a simple one - dimensional layered model to verify the effectiveness of the method. - Simulate the travel - time field under Gaussian hill terrain and observe the influence of the terrain on the ray path. - Conduct a checkerboard model test to verify the processing effect of the topographic effect. 2. **Application to actual data**: - Apply the new method to the actual short - period seismic data in the Qingchengzi ore field on the Liaodong Peninsula in Northeast China. - The inversion results show that there are widely distributed low - velocity bodies in the central area, which are interpreted as low - velocity granite intrusions in the Paleoproterozoic Liaohe Group. - There are high - velocity bodies in the eastern and western parts of the Baiyun gold belt respectively, which are interpreted as high - velocity features in the marble - limestone. ### Conclusion The method proposed in this paper can effectively handle the influence of topographic changes on the propagation of surface waves and improve the accuracy of surface - wave tomography. This is of great significance for geological exploration, resource assessment, and other fields.