Sound speed profiles (SSPs) must be detected simultaneously to perform a multibeam depth survey. Accurate real-time sound speed profile (SSP) acquisition remains a critical challenge in deep-sea multibeam bathymetry due to the limitations regarding direct measurements under harsh operational conditions. To address the issue, we propose a joint inversion framework integrating World Ocean Atlas 2023 (WOA23) temperature–salinity model data, historical in situ SSPs, and surface sound speed measurements. By constructing a high-resolution regional sound speed field through WOA23 and historical SSP fusion, this method effectively mitigates spatiotemporal heterogeneity and seasonal variability. The artificial lemming algorithm (ALA) is introduced to optimize the inversion of empirical orthogonal function (EOF) coefficients, enhancing global search efficiency while avoiding local optimization. An experimental validation in the northwest Pacific Ocean demonstrated that the proposed method has a better performance than that of conventional substitution, interpolation, and WOA23-only approaches. The results indicate that the mean absolute error (MAE), root mean square error (RMSE), and maximum error (ME) of SSP reconstruction are reduced by 41.5%, 46.0%, and 49.4%, respectively. When the reconstructed SSPs are applied to multibeam bathymetric correction, depth errors are further reduced to 0.193 m (MAE), 0.213 m (RMSE), and 0.394 m (ME), effectively suppressing the “smiley face” distortion caused by sound speed gradient anomalies. The dynamic selection of the first six EOF modes balances computational efficiency and reconstruction fidelity. This study provides a robust solution for real-time SSP estimation in data-scarce deep-sea environments, particularly for underwater autonomous vehicles. This method effectively mitigates the seabed distortion caused by missing real-time SSPs, significantly enhancing the accuracy and efficiency of deep-sea multibeam surveys.

Multibeam echosounder systems (MBESs) are widely used in unmanned surface vehicle swarms (USVs) to perform various marine bathymetry surveys because of their excellent performance. To address the challenges of systematic error superposition and edge beam error propagation in multibeam bathymetry surveying, this study proposes a novel error adjustment method integrating crossover error density clustering and beam incident angle (BIA) compensation. Firstly, a bathymetry error detection model was developed based on adaptive Density-Based Spatial Clustering of Applications with Noise (DBSCAN). By optimizing the neighborhood radius and minimum sample threshold through analyzing sliding-window curvature, the method achieved the automatic identification of outliers, reducing crossover discrepancies from ±150 m to ±50 m in the deep sea at a depth of approximately 5000 m. Secondly, an asymmetric quadratic surface correction model was established by incorporating the BIA as a key parameter. A dynamic weight matrix ω = 1/(1 + 0.5θ²) was introduced to suppress edge beam errors, combined with Tikhonov regularization to resolve ill-posed matrix issues. Experimental validation in the Western Pacific demonstrated that the RMSE of crossover points decreased by about 30.4% and the MAE was reduced by 57.3%. The proposed method effectively corrects residual systematic errors while maintaining topographic authenticity, providing a reference for improving the quality of multibeam bathymetric data obtained via USVs and enhancing measurement efficiency.