The construction of the high-resolution Moho depth model is significant for studying the characteristics of the complex tectonic movement (seafloor spreading, plate subduction phenomena) in Papua New Guinea. We calculate the region’s Moho relief and lithosphere thinning factor using the XGM 2019e gravity field model and nonlinear fast gravity inversion method under the GEMMA Moho depth model’s constraint considering the influence of lithosphere thermal gravity anomaly. The calculation result shows that the Moho depth is between 6—34 km, forming two large depressions in Woodlark Basin (WB) and Solomon Sea Plate (SSP) with deep scattered islands. In addition, the findings suggest that Significant differences exist in the shape and tectonic movement intensity of the North and South oceanic crust at the WB. Nevertheless, the lithosphere extends evenly in Manus Basin (MB). WB collided with the Solomon Islands at a higher angle than the SSP subducted under Bismarck Sea Plate (BSP); strong earthquakes may frequently occur on both sides and in deeper positions at West New Britain Trench in the future.

Three-dimensional gravity inversion technology involves inferring the underground density structure based on observed gravity anomaly data. In addition to gravity inversion based on physics-driven methods, deep learning, as a purely data-driven technique, is increasingly gaining attention in geophysical inversion problems. However, purely data-driven methods rely on the implicit relationships within the data during the inversion process, which results in a lack of clear physical significance. This study proposes a three-dimensional gravity inversion method that integrates physical equations with deep learning. Based on the U-Net architecture, the gravity forward equation is incorporated as a physical constraint term, and a composite loss function—comprising three-dimensional mean squared error, a depth-weighting function, and three-dimensional intersection-over-union loss—is constructed to enhance inversion accuracy. Numerical experiments indicate that this method outperforms traditional algorithms in terms of density recovery accuracy and boundary clarity. When applied to gravity anomaly data from the Tangshan earthquake region in China, this method successfully inverted the three-dimensional subsurface density structure, revealing a high-density anomaly beneath the seismic source area, which provides important evidence for understanding the regional earthquake generation mechanism.

The tectonic setting of Taiwan and its surrounding regions is characterized by the complex interaction between the northwest-oriented Ryukyu subduction zone and the east-oriented Manila subduction zone. Within this subduction framework, the elastic thickness of the lithosphere (Te) serves as a critical parameter for elucidating the mechanical behavior of the area. In this study, we employed the admittance–correlation method to estimate Te values across Taiwan and adjacent territories. The findings indicate that sedimentary loading results in an overestimation of the maximum Te by approximately 50 km; after adjustment, the Te values range from 0 to 60 km throughout the study area. On Taiwan, Te values predominantly lie between 20 and 30 km, decreasing to 10–20 km near the margins adjacent to the Ryukyu and Manila subduction fronts. The Philippine Sea Plate exhibits comparatively higher Te values, ranging from 40 to 65 km. The spatial distribution of Te broadly corresponds with major tectonic subdivisions. Statistical analyses reveal a weak negative correlation between Te and surface heat flow (r = −0.44) and a weak positive correlation with shear-wave velocity anomalies at a depth of 100 km (r = 0.22), suggesting that the thermal structure exerts only a moderate influence on lithospheric strength in this region. Nonetheless, within oceanic crustal domains, the relationship between Te and oceanic crustal age largely adheres to models of crustal cooling and lithospheric thickening, consistent with isotherm depths of approximately 200–400 °C. Additionally, dynamic topography associated with slab subduction may locally diminish Te by up to 25 km. Cross-sectional profiles through northern Taiwan and the Philippine Sea block reveal pronounced coupling between subduction geometry and Te distribution. The observed spatial patterns of Te reflect the mechanical imprint of prolonged tectonic evolution, with the orientation of Te gradients generally aligned with the direction of maximum principal compressive stress. Collectively, these results suggest that subduction geometry and tectonic processes are important factors influencing the spatial variability and evolutionary trajectory of lithospheric strength in Taiwan and its environs.

The existence of aftershocks in an earthquake sequence can impact the analysis of the mainshock. In this study, we present a method for deleting an aftershock sequence based on the spatial relationship between earthquakes and faults. This method improves the performance of space window selection in the classical K-K method by eliminating aftershocks with an ideal fault buffer zone. The determination of fault buffer zones is based on a trial-and-error analysis of 69,714 earthquake records from the China Seismic Network Center (CENC) collected between 1980 and 2020. We selected 20 typical big earthquakes ($M_L$7.0–8.0 or $~Ms$6.6–8.0; for earthquakes above magnitude $Ms$7 or $M_{L}7.2$, $M_L$ is approximately equal to $Ms$) as the mainshocks to establish the fault buffer zones. We also propose an empirical formula to determine the distance of the fault buffer zone by counting the aftershock deletion effect at different buffer distances. Compared with the classical K-K method, our method considers the correlation between the spatial distribution of aftershocks and faults, eliminates earthquake groups that are not related to the mainshock, greatly reduces the spatial range of aftershocks, improves the performance of deleting aftershocks of different magnitudes, and provides a new rule and reference for aftershock deletion.