Search

On the basis of Space-Wise Least Square method, three numerical methods including Cholesky decomposition, pre-conditioned conjugate gradient and Open Multi-Processing parallel algorithm are applied into the determination of gravity field with satellite gravity gradiometry data. The results show that, Cholesky decomposition method has been unable to meet the requirements of computation efficiency when the computer hardware is limited. Pre-conditioned conjugate gradient method can improve the computation efficiency of huge matrix inversion, but it also brings a certain loss of precision. The application of Open Multi-Processing parallel algorithm could achieve a good compromise between accuracy and computation efficiency. Key words: satellite gravity gradiometry; Cholesky decomposition; pre-conditioned conjugate gradient; open multi-processing parallel algorithm; data processing

The global marine gravity anomaly model is predominantly recovered from along-track radar altimeter data. Despite significant advances in gravity anomaly recovery, the improvement of the gravity anomaly model remains constrained by the absence of cross-track geoid gradients and the reduction in radar altimeter data, especially in coastal and high-latitude regions. ICESat-2 laser altimetry, with a three-pair laser beam configuration, a small footprint, and a near-polar orbit, facilitates the determination of cross-track geoid gradients and provides valid observations in certain regions. We present an ICESat-2 altimeter data processing strategy that includes the determination of cross-track geoid gradients and the combination of along-track and cross-track geoid gradients. Utilizing these methods, we developed a new global marine gravity model, SDUST2022GRA, from radar and laser altimeter data. Different weight determination methods were applied to each type of altimeter datum. The precision and spatial resolution of SDUST2022GRA were assessed against published altimeter-derived global gravity anomaly models (DTU17, V32.1, NSOAS22) and shipborne gravity measurements. SDUST2022GRA achieved a global precision of 4.43 mGal , representing an improvement of approximately 0.22 mGal over existing altimeter-derived models. In local coastal and high-latitude regions, SDUST2022GRA showed an enhancement of 0.16–0.24 mGal compared to the other models. The spatial resolution of SDUST2022GRA is approximately 20 km in certain regions, which is slightly superior to the other models. The percentage contribution of ICESat-2 to the improvement of the gravity anomaly model is 4.3 % in low- to mid-latitude regions by comparing SDUST2022GRA with ICESat-2 to SDUST2021GRA without ICESat-2, and this is increasing in coastal regions. These assessments suggest that SDUST2022GRA is a reliable global marine gravity anomaly model. The SDUST2022GRA data are freely available at https://doi.org/10.5281/zenodo.8337387 (Li et al., 2023). [ABSTRACT FROM AUTHOR]

The Philippine Sea, located at the edge of the northwestern Pacific Ocean, possesses complex seabed topography. Developing a high-precision bathymetric model for this region is of paramount importance, as it provides fundamental geoinformation essential for Earth observation and marine scientific research, including plate motion, ocean circulation, and hydrological characteristics. The gravity–geologic method (GGM), based on marine gravity anomalies, serves as an effective bathymetric prediction technique. To further strengthen the prediction accuracy of conventional GGM, we introduce the improved GGM (IGGM). The IGGM considers the effects of regional seafloor topography by employing weighted averaging to more accurately estimate the short-wavelength gravity component, along with refining the subsequent modeling of long-wavelength gravity component. In this paper, we focus on seafloor topography modeling in the Philippine Sea based on the IGGM, combining shipborne bathymetric data with the Scripps Institution of Oceanography (SIO) V32.1 gravity anomaly. To reduce computational complexity, the optimal parameter values required for IGGM are first calculated before the overall regional calculation, and then, based on the terrain characteristics and distribution of sounding data, we selected four representative local sea areas as the research objects to construct the corresponding bathymetric models using GGM and IGGM. The analysis indicates that the precision of the IGGM models in four regions is improved to varying degrees, and the optimal calculation radius is 2 ′. Based on the above finding, a high-precision 1′×1′ bathymetric model of the Philippine Sea (5–35° N, 120–150° E), known as the BAT_PS model, is constructed using IGGM. Results demonstrate that the BAT_PS model exhibits a higher overall precision compared to the General Bathymetric Chart of the Oceans (GEBCO), topo_25.1, and DTU18 models at single-beam shipborne bathymetric points. [ABSTRACT FROM AUTHOR]

On 6 February 2023, the MW 7.8 Pazarcik and the MW 7.5 Elbistan earthquakes occurred in southeastern Turkey, close to the Syrian border, causing many deaths and a great deal of property destruction. The Pazarcik earthquake mainly damaged the East Anatolian Fault Zone (EAFZ). The Elbistan earthquake mainly damaged the Cardak fault (CF) and the Doğanşehir fault (DF). In this study, Sentinel-1A ascending (ASC) and descending (DES) orbit image data and pixel offset tracking (POT) were used to derive surface deformation fields in the range and azimuth directions induced by the Pazarcik and Elbistan earthquakes (hereinafter referred to as the Turkey double earthquakes). Utilizing GPS coordinate sequence data, we computed the three-dimensional surface deformation resulting from the Turkey double earthquakes. The surface deformation InSAR and GPS results were combined to invert the coseismic slip distribution of the EAFZ, CF, and DF using a layered earth model. The results show that the coseismic ruptures of the Turkey double earthquakes were dominated by left-lateral strike-slips. The maximum slip was 7.76 m on the EAFZ and about 8.2 m on the CF. Both the earthquakes ruptured the surface. The Coulomb failure stress (CFS) was computed based on the fault slip distribution and the geometric parameters of all the active faults within 300 km of the MW 7.8 Pazarcik earthquake's epicenter. The CFS change resulting from the Pazarcik earthquake suggests that the subsequent Elbistan earthquake was triggered by the Pazarcik earthquake. The Antakya fault experienced an increase in CFS of 8.4 bars during this double-earthquake event. Therefore, the MW 6.3 Uzunbağ earthquake on 20 February 2023 was jointly influenced by the Turkey double earthquakes. Through stress analysis of all the active faults within 300 km of the MW 7.8 Pazarcik earthquake's epicenter, the Ecemis segment, Camliyayla fault, Aadag fault, Ayvali fault, and Pula segment were all found to be under stress loading. Particularly, the Ayvali fault and Pula segment exhibited conspicuous stress loading, signaling a higher risk of future seismic activity. [ABSTRACT FROM AUTHOR]

The Philippine Sea, located at the edge of the Northwest Pacific Ocean, possesses complex seabed topography. Developing a high-precision bathymetric model for this region is of paramount importance as it provides fundamental geoinformation essential for Earth observation and marine scientific research, including plate motion, ocean circulation, and hydrological characteristics. The gravity-geologic method (GGM), based on marine gravity anomalies, serves as an effective bathymetric prediction technique. To further strengthen the prediction accuracy of conventional GGM, we introduce the improved GGM (IGGM). The IGGM considers the effects of regional seafloor topography by employing weighted averaging to more accurately estimate the short-wavelength gravity component, along with refining the subsequent modeling of long-wavelength gravity component. In this paper, we focus on seafloor topography modeling in the Philippine Sea based on the IGGM, combining shipborne bathymetric data with the SIO V32.1 gravity anomaly. To reduce computational complexity, the optimal parameter values required for IGGM are first calculated before the overall regional calculation, and then, based on the terrain characteristics and distribution of sounding data, we selected four representative local sea areas as the research objects to construct the corresponding bathymetric models using GGM and IGGM. The analysis indicates that the precisions of the IGGM models in four regions are improved to varying degrees, and the optimal calculation radius is 2′. Based on the above finding, a high-precision 1'×1' bathymetric model of the Philippine Sea (5° N–35° N, 120° E–150° E), known as the BAT_PS model, is constructed using IGGM. Results demonstrate that the BAT_PS model exhibits a higher overall precision compared to GEBCO, topo_25.1, and DTU18 models at single-beam shipborne bathymetric points. [ABSTRACT FROM AUTHOR]