A Method for Full-Depth Sound Speed Profile Reconstruction Based on Average Sound Speed Extrapolation

: The speed of sound in seawater plays a crucial role in determining the accuracy of multi-beam bathymetric measurements. In deep-sea multibeam measurements, the challenge of inadequate longitudinal coverage of sound speed profiles arises from variations in seafloor topography, meteorological conditions, measurement equipment, and operational efficiency, resulting in diminished measurement precision. Building upon the EOF (Empirical Orthogonal Function), a method employed to analyze spatiotemporal data such as sound speeds, this paper addresses the limitations of the EOF method caused by the shallowest sampling depth of the sound speed profile samples. It proposes two methods for EOF reconstruction of measured sound speed profiles extended to full water depth by splicing measured sound speed profiles at non-full water depths with historical average sound speed profiles of the surveyed sea area. Specially, Method 2 introduces the latest metaheuristic optimization algorithm, CPO (Crested Porcupine Optimizer), which exhibited superior performance on multiple standard test functions in 2024. The study reconstructs randomly sampled measured sound speed profiles using the two proposed methods and commonly employed substitution and splicing methods, followed by a comparative analysis of the experimental outcomes. At a sampling depth of 200 m, Method 2 demonstrates performance superior to other methods, with RMSE, MAE, MAPE, and R 2 values of 0.9511 m/s, 0.8492 m/s, 0.0566%, and 0.9963, respectively. Method 1 yields corresponding values of 0.9594 m/s, 0.8492 m/s, 0.0568%, and 0.9962, respectively. Despite its slightly inferior performance compared with Method 2, it offers substantial advantages over the substitution and splicing methods. Varying the sampling depth of measured sound speed profiles reveals that Methods 1 and 2 exhibit inferior reconstruction performance in shallow water compared with the substitution and splicing methods. Nevertheless, when the sampling depth surpasses the depth range of initial spatial modes with abrupt variations, both methods achieve notably higher reconstruction accuracy compared with the substitution and splicing methods, reaching a stabilized state. Sound ray tracing reveals that the reconstructed sound speed profiles from both methods meet the stringent accuracy standards for bathymetric measurements, achieving an effective beam ratio of 100%. The proposed methods not only provide rapid reconstruction of sound speed profiles, thereby improving the efficiency of multibeam bathymetric surveys, but also provide references for the reasonable determination of sampling depths of sound speed profiles to ensure reconstruction accuracy, demonstrating practical application value.


Introduction
The precision of multibeam seafloor topography measurements is contingent upon numerous factors, among which the speed of sound in seawater emerges as pivotal.Principally, it impacts the computation of seafloor measurement point coordinates by shaping sound ray trajectories, and the accuracy and dependability of measurement outcomes hinge directly upon the correctness of the sound speed profile [1][2][3][4][5].Conventionally, the sound speed profile manifests as a function of depth-dependent sound velocity, mirroring the longitudinal variations in sound speed structure within the local marine domain, chiefly dictated by marine environmental factors like seawater temperature, salinity, and pressure [6,7].Due to the influence of factors such as solar radiation, precipitation, ocean currents, tides, internal waves, and others on the marine environment, the sound speed profile exhibits complex spatiotemporal variations [8].Therefore, it is necessary to timely acquire accurate sound speed profiles of the measured sea area for real-time or post-survey correcting of multibeam survey data.
Currently, there are two main methods for measuring seawater sound speed: direct methods and indirect methods [9,10].Direct methods involve using equipment such as AML Minos SVP or SVP-plus sound velocity profilers to directly obtain sound speed profiles in seawater through fixed-point measurements.Indirect methods utilize devices like CTD, XCTD, XBT, MVP, etc., to obtain temperature and salinity information in seawater through survey or fixed-point measurements, calculating the sound speed in seawater using suitable empirical formulas.When conducting deep-sea bathymetric surveys, whether using direct methods or indirect methods, the following situations are often encountered, leading to the inability to obtain a complete sound speed profile for the entire sea depth.
(1) Complex variations in seafloor topography in the surveyed sea area.This includes two situations.First, during the navigation progress of a survey vessel from a dock to deep sea through shallow waters.The water depth gradually increases, leading to the collected sound speed profile not achieving full coverage of the water's depth.Second, complexity in the topographical changes within the surveyed sea area.If the selection of sound profile stations is unreasonable (such as choosing points in shallow waters), it results in an insufficient vertical coverage range of the sound speed profile, leading to incomplete post-survey sound speed corrections.
(2) Insufficient effective detection depth of equipment.The effective detection depth of XBT is approximately 800 m, while that of XCTD is around 1100 m.For deep-sea topography measurements in water depths of over a thousand meters, the detection depth of such survey equipment falls far short of the requirements for comprehensive coverage of an area's topography measurements.
(3) Operational inefficiency.The "Hydrographic Survey Specification of China" stipulates that sound speed acquisition operations should ideally be conducted daily.When measuring sound speed using CTD or sound velocity profilers to obtain a sound speed profile of the entire water depth in areas in which the sea is 4000 m deep, a conservative estimate indicates that it would take 3-4 h [11], from stopping the vessel to deploying and then recovering the equipment, to obtain a profile.Operational efficiency is thus low.
(4) Meteorological conditions.During operations, interruptions due to high winds and waves at sea are common, leading to the inability to obtain a complete sound speed profile of the entire sea depth.
The National Oceanic and Atmospheric Administration (NOAA) of the United States clearly stipulates that during the measurement process, the actual non-full water depth sound speed profile must be vertically extended to the maximum depth of the actual field water depth [12].Field practices have found that using non-full water depth sound speed profiles for speed corrections can cause the depth measurement strips to exhibit undulations, impacting the accurate representation of seafloor topography [13].Therefore, reconstructing sound speed profiles of the full water depth using measured profiles is of significant importance for enhancing depth measurement accuracy and operational efficiency.
In reconstructing sound speed profiles, using mathematical means for fitting and extrapolation is a common and efficient approach, especially suitable for nearshore or offshore areas.This method is particularly applicable because, compared to shallow water areas such as estuaries and coastlines where sound speed variations are irregular, the sound speed variations in nearshore or offshore waters exhibit certain regularities and can be to some extent predicted [14].Wu Tongyuan et al. [15] employed a polynomial model to fit and extrapolate the measured sound speed profiles.The results showed that, among various extrapolations based on fourth-order fitting, third-order polynomial extrapolation still maintained high accuracy at relatively long distances.This method is highly effective when dealing with constant sound speed gradients or minor gradient changes.However, it may lead to significant errors when encountering sudden features such as a sound channel axis.To represent the sound speed structures of different water layers such as the mixed layer, thermocline, deep sound channel layer, and deep isothermal layer, structured parameter models of sound speed profiles are commonly used.Munk [16] proposed a sound speed profile model for mid-latitude deep-sea sound channels using an exponential expression based on the layering characteristics of deep-sea density structures.Davis et al. [17] proposed the Generalized Digital Environmental Model (GDEM), which describes temperature and salinity profiles in seawater using more complex function combinations.Building upon analytical function models such as Munk and GDEM, Zhang Xu et al. [18] proposed a 9-parameter sound speed profile model, LSSPM (Layered Sound Speed Profile Model), which is more suitable for the Chinese nearshore environment.This model significantly reduces the number of parameters while integrating the advantages of both Munk and GDEM models, thereby improving the efficiency of describing sound speed structures.Although parameterized models facilitate the expression of sound speed profiles, their idealized representation cannot fully describe complex marine environments.
EOF is a widely used method for representing sound speed profiles.It can accurately represent the complex variations of sound speed profiles using the first few eigenvectors obtained by EOF decomposition and a small number of parameters.Research by LeBlanc et al. [19] shows that EOF is the most effective basis function for describing sound speed profiles.Park et al. [20] demonstrated that only a few representative empirical orthogonal functions are needed to represent the main differences in the sound field of a marine area.Shen Yuanhai et al. [21] demonstrated the feasibility of using EOF to represent shallow water sound speed profiles using data from the South China Sea.Sun Wenchuan et al. [22] demonstrated the effectiveness of the EOF method in representing sound speed profiles.He Li et al. [23] represented sound speed profiles in the East China Sea using a small number of empirical orthogonal functions.Zhang Zhenmai et al. [24] reconstructed the sound speed profile of the entire ocean depth using finite-depth measured sound speed data combined with the EOF method.Zhang Wei et al. [25] extrapolated incomplete temperature and salinity profiles, calculated sound speed using empirical formulas, selected three sound speed points with significant changes, solved EOF coefficients, and reconstructed the entire ocean depth sound speed profile using the first three EOFs.Cheng Fang et al. [26] used a small amount of full-depth sound speed profile sampling data in the surveyed area as a basis and the EOF method to extrapolate the sound speed profiles that were not sampled to full depth.Although the EOF method can simplify the representation of sound speed profiles, it is not suitable for extrapolating non-full-depth sound speed profiles because this method is limited by the minimum sampling depth of the sample data.
Inspired by the research results of Yang Fan et al. [27] and based on the reconstruction of sound speed profiles using EOF decomposition, this paper addresses the practical problem of the EOF method's inability to extend measured non-full-depth sound speed profiles to full depth.It considers splicing the full-depth historical average sound speed profile of an area with a measured sound speed profile of that area to extend the measured profile to full depth.On this basis, two methods are proposed to reconstruct the measured sound speed profile.One is to combine the extended measured sound speed profile with the sample data of historical sound speed profiles to form a new sound speed matrix, perform EOF analysis on this matrix, and reconstruct the extended measured sound speed profile.The other is to directly perform an EOF analysis on the matrix of historical sound speed profiles and use the rapid optimization ability of the CPO algorithm to find the optimal reconstruction coefficients for the extended measured sound speed profile.For convenience, these two methods are respectively referred to as Method 1 and Method 2. This article is organized as follows: In Section 1, the significance of sound speed profile reconstruction is briefly introduced, along with an overview of existing methods and the associated challenges.Section 2 provides a detailed explanation of EOF, highlighting the novel approach proposed in this study.The experimental setup, including the data utilized, the experimental process, and the evaluation metrics for the results, is presented in Section 3. Section 4 presents a comprehensive analysis of the experimental results, discussing the findings in detail.Finally, in Section 5, the relevant work discussed in the article is summarized and potential avenues for future research are proposed.

EOF Decomposition and Reconstruction of Sound Speed Profiles
The EOF analysis, also known as Principal Component Analysis (PCA) or Eigenvector Analysis, was initially introduced by the statistician Karl Pearson [28] in 1901.It is a method for analyzing the structure of a matrix of data and extracting its primary features.In the 1950s, Lorenz introduced the EOF method to the field of atmospheric science.This method can decompose meteorological field elements that vary over time into orthogonal spatial modes and time coefficients, effectively extracting the primary patterns of variation and their time series from the data.It is an important tool in the study of large-scale circulation, ocean dynamics, climate change, and other fields.The fundamental process of EOF decomposition and reconstruction of sound speed profiles is illustrated in Figure 1.
convenience, these two methods are respectively referred to as Method 1 and Method 2.
This article is organized as follows: In Section 1, the significance of sound speed profile reconstruction is briefly introduced, along with an overview of existing methods and the associated challenges.Section 2 provides a detailed explanation of EOF, highlighting the novel approach proposed in this study.The experimental setup, including the data utilized, the experimental process, and the evaluation metrics for the results, is presented in Section 3. Section 4 presents a comprehensive analysis of the experimental results, discussing the findings in detail.Finally, in Section 5, the relevant work discussed in the article is summarized and potential avenues for future research are proposed.

EOF Decomposition and Reconstruction of Sound Speed Profiles
The EOF analysis, also known as Principal Component Analysis (PCA) or Eigenvector Analysis, was initially introduced by the statistician Karl Pearson [28] in 1901.It is a method for analyzing the structure of a matrix of data and extracting its primary features.In the 1950s, Lorenz introduced the EOF method to the field of atmospheric science.This method can decompose meteorological field elements that vary over time into orthogonal spatial modes and time coefficients, effectively extracting the primary patterns of variation and their time series from the data.It is an important tool in the study of large-scale circulation, ocean dynamics, climate change, and other fields.The fundamental process of EOF decomposition and reconstruction of sound speed profiles is illustrated in Figure 1.Specifically, to decompose the sound speed profiles using the EOF method, first, the N sound speed profiles collected at different times in the measured sea area need to be interpolated in the depth direction into M standard layers, resulting in the standardized Here, ij c represents the sound speed value corresponding to the j water layer of the i sound speed profile.Each row comprises the sound speed values of N sound speed profile samples at the same water layer, and each column represents the sound speed values of the same sound speed profile at different water layers.By averaging each row of the standardized sound speed matrix, the historical average sound speed profile is obtained, where i c represents the average sound speed value of different sound speed profiles at the same water layer.The average sound speed profile is expanded into N columns by replicating it column-wise, resulting in the Specifically, to decompose the sound speed profiles using the EOF method, first, the N sound speed profiles collected at different times in the measured sea area need to be interpolated in the depth direction into M standard layers, resulting in the standardized sound speed matrix C M×N : Here, c ij represents the sound speed value corresponding to the j water layer of the i sound speed profile.Each row comprises the sound speed values of N sound speed profile samples at the same water layer, and each column represents the sound speed values of the same sound speed profile at different water layers.By averaging each row of the standardized sound speed matrix, the historical average sound speed profile T is obtained, where c i represents the average sound speed value of different sound speed profiles at the same water layer.The average sound speed profile is expanded into N columns by replicating it column-wise, resulting in the average sound speed matrix C M×N .Subtracting the average sound speed matrix from the standardized sound speed matrix, the perturbation matrix (anomalies matrix) ∆C M×N is obtained: The covariance matrix R M×M of the perturbation matrix is computed as: The essence of EOF analysis lies in analyzing the covariance matrix of the data, obtaining the eigenvalues and corresponding eigenvectors of the covariance matrix through eigen decomposition.Each eigenvector represents a spatial mode, and its corresponding eigenvalue indicates the importance of that mode to the entire dataset.Perform eigen decomposition on R M×M : Here, D M×M represents the matrix of eigenvalues of the covariance matrix.After arranging the eigenvalues in descending order, the following is obtained: F M×M represents the matrix of eigenvectors corresponding to the eigenvalue matrix.After synchronous adjustment with the eigenvalue matrix, it can be expressed as follows: Here, f k (k = 1, 2, . . ., M) represents the k order spatial mode.After performing EOF decomposition on the historical sound speed profiles in the measurement area, any sound speed profile within the region can be reconstructed using the leading spatial modes.
Here, c(z) represents the reconstructed sound speed profile, c(z) represents the average sound speed profile, z represents the depth of each standard layer, f i (z) represents the i order spatial mode, α i represents the corresponding time coefficient of the i order spatial mode, and K represents the order of the selected spatial modes.
By moving c(z) from Equation (7) to the left side and expressing it in matrix form, ∆C M×N is then obtained as follows: Thus, the matrix of time coefficients A K×N is obtained as follows: Here, T represents the time coefficient corresponding to the i-th sound speed profile.After obtaining the time coefficients, any sound speed profile within the measurement area can be reconstructed according to Equation (7).

Variance Contribution Ratio of Spatial Modes
In EOF analysis, the variance contribution ratio is a crucial concept employed to quantify the contribution of each order of spatial modes to the entire sound speed field, with the calculation formula as follows: Here, σ i represents the variance contribution ratio of the i order spatial mode, λ i denotes the eigenvalue corresponding to the i order spatial mode (in descending order), and M is the total number of eigenvalues.
The cumulative variance contribution ratio of the first K order spatial modes is given by: When reconstructing sound speed profiles, concerning the selection of spatial mode orders, Davis [29] pointed out that if the cumulative variance contribution ratio of the first K order spatial modes exceeds or equals 95%, then it is considered that the main information of the sound speed field in that region can be precisely delineated by these K order spatial modes.

Method 1
Illustrated in Figure 2, the fundamental process of Method 1 for a designated measurement sea area unfolds as follows: x FOR PEER REVIEW 7 of 25 (3) Compute the average value of the matrix of standard sound speed profiles rowwise to derive the historical average sound speed profile.
(4) Determine the disparity between the sound speed at the concluding depth of the sound speed profile to be reconstructed and the sound speed at the depth corresponding to the identical depth of the average sound speed profile.Add this difference to the average sound speed profile from this depth to the last depth corresponding to it, then truncate this part of the average sound speed profile and concatenate it with the measured sound speed profile to be reconstructed to obtain an extrapolated measured sound speed profile.
(5) Combine the extrapolated measured sound speed profile with the original matrix of standard sound speed profiles to form a new matrix of sound speed profiles.In accordance with Section 2.1, perform EOF decomposition on this matrix to obtain spatial modes and temporal coefficients.Taking 95% as the minimum threshold, calculate the cumulative variance contribution ratio of the first K order spatial modes to determine the minimum number of orders of spatial modes necessary to reconstruct the sound speed profile, and then reconstruct the measured sound speed profile according to Equation (7).

CPO
CPO is a novel metaheuristic intelligent optimization algorithm proposed by Abdel-Basset et al. in 2024 [32].This algorithm emulates a plethora of defense behaviors exhibited by the crested porcupine, rendering it especially conducive for tackling large-scale optimization problems.The defense mechanism of the crested porcupine is categorized into four escalating types: visual, auditory, olfactory, and physical attacks.Within the framework of the CPO algorithm, the search space is depicted as illustrated in Figure 3, where regions A, B, C, and D symbolize four distinct defense stages.Region A represents the first defense stage, where the crested porcupine raises and flaps its quills upon recognizing the  (1) Choose historical full-depth sound speed profiles obtained at various stations and times within the measurement sea area (if temperature-salinity-depth (pressure) data are available at these stations, suitable empirical formulas should be employed to convert them into sound speed values).The choice of historical sound speed profiles limits the precision of sound speed profile EOF reconstruction.If there are significant differences in the properties of sound speed variations within the measurement area, it may lead to a decrease in the accuracy of sound speed profile EOF reconstruction, failing to accurately reflect the real sound speed variation.At present, there is no specific discussion of the spatiotemporal range concerning the selection of historical sound speed profiles.Although a narrower spatiotemporal range yields superior reconstruction outcomes, for extensive deep-sea terrain measurements, the reconstruction efficiency will be significantly diminished.In pertinent literature, numerous researchers confine the spatial scope within an interval of 5 • × 5 • , and a temporal scope within quarterly periods.Of course, the spatiotemporal range for selecting historical sound speed profiles should also be reasonably determined according to the actual hydrological conditions of the measurement sea area.If drastic changes occur in the hydrological elements within the measurement area, the spatial range can be narrowed down to 3 • × 3 • , and the temporal range limited to monthly intervals, concurrently augmenting the number of historical sound speed profile samples in the measurement area.
(2) Employ the Akima interpolation method [30,31] to interpolate the historical sound speed profiles chosen in step (1) into standard sound speed profiles with uniform intervals based on the specified depth interval, thereby assembling a matrix of standard sound speed profiles.To reconstruct the measured non-full-depth sound speed profiles, employ the same method to interpolate them to their maximum measured depth.
(3) Compute the average value of the matrix of standard sound speed profiles row-wise to derive the historical average sound speed profile.
(4) Determine the disparity between the sound speed at the concluding depth of the sound speed profile to be reconstructed and the sound speed at the depth corresponding to the identical depth of the average sound speed profile.Add this difference to the average sound speed profile from this depth to the last depth corresponding to it, then truncate this part of the average sound speed profile and concatenate it with the measured sound speed profile to be reconstructed to obtain an extrapolated measured sound speed profile.
(5) Combine the extrapolated measured sound speed profile with the original matrix of standard sound speed profiles to form a new matrix of sound speed profiles.In accordance with Section 2.1, perform EOF decomposition on this matrix to obtain spatial modes and temporal coefficients.Taking 95% as the minimum threshold, calculate the cumulative variance contribution ratio of the first K order spatial modes to determine the minimum number of orders of spatial modes necessary to reconstruct the sound speed profile, and then reconstruct the measured sound speed profile according to Equation (7).

Method 2 2.3.1. CPO
CPO is a novel metaheuristic intelligent optimization algorithm proposed by Abdel-Basset et al. in 2024 [32].This algorithm emulates a plethora of defense behaviors exhibited by the crested porcupine, rendering it especially conducive for tackling large-scale optimization problems.The defense mechanism of the crested porcupine is categorized into four escalating types: visual, auditory, olfactory, and physical attacks.Within the framework of the CPO algorithm, the search space is depicted as illustrated in Figure 3, where regions A, B, C, and D symbolize four distinct defense stages.Region A represents the first defense stage, where the crested porcupine raises and flaps its quills upon recognizing the predator, implementing the first defense strategy, namely visual defense.Region B represents the second defense stage, indicating that the predator is not deterred by the first defense mechanism and is moving toward the prey.At this point, the crested porcupine starts making noise to threaten the predator, implementing the second defense mechanism, namely auditory defense.Region C represents the third defense stage, where if the predator is not deterred by the first and second defense mechanisms and continues to approach the crested porcupine, it begins to secrete odors and spread them around to dissuade the predator from getting closer, implementing the third defense mechanism, namely olfactory defense.Region D represents the fourth defense stage, where the crested porcupine starts to physically attack the predator.The initial two defense mechanisms are affiliated with the exploratory behavior of CPO, whereas the subsequent two defense mechanisms are aligned with the exploitative behavior of CPO.The mathematical simulation of each defense mechanism is detailed in [32] and will not be reiterated in this paper.
It is noteworthy that the algorithm introduces a highly innovative strategy termed as the "population cyclic reduction technique", wherein solely threatened crested porcupines trigger their defense mechanisms rather than the entire population.Test results suggest that the algorithm demonstrates markedly superior performance in contrast to other commonly employed optimization algorithms.
J. Mar.Sci.Eng.2024, 12, x FOR PEER REVIEW 8 of 25 affiliated with the exploratory behavior of CPO, whereas the subsequent two defense mechanisms are aligned with the exploitative behavior of CPO.The mathematical simulation of each defense mechanism is detailed in [32] and will not be reiterated in this paper.
It is noteworthy that the algorithm introduces a highly innovative strategy termed as the "population cyclic reduction technique", wherein solely threatened crested porcupines trigger their defense mechanisms rather than the entire population.Test results suggest that the algorithm demonstrates markedly superior performance in contrast to other commonly employed optimization algorithms.

Basic Process
The objective of employing the CPO algorithm is to leverage its rapid convergence advantage in order to swiftly obtain optimal EOF reconstruction coefficients.As depicted in Figure 4, the basic process for Method 2, applied to the designated measurement marine area, is outlined as follows: (1)-( 4) follow the steps outlined in Section 2.2 (1)-( 4).
(5) Conduct EOF decomposition on the standard sound speed profile matrix to extract spatial modes and temporal coefficients.Compute the cumulative variance contribution rate of the first K spatial modes, with a minimum threshold of 95%, to determine the minimum number of spatial modes K involved in the reconstruction of measured sound speed profiles.This value of K serves as one of the input parameters for the CPO algorithm, termed as the problem dimension.Simultaneously, identify the maximum and minimum values of the temporal coefficients corresponding to each order of the K order spatial modes.Arrange the maximum values of the temporal coefficients in ascending order of order to form a row vector, constituting the second input parameter of the CPO algorithm, denoted as the upper bound of the search space.Similarly, arrange the minimum values of the temporal coefficients in ascending order of order to form a row vector, constituting the third input parameter of the CPO algorithm, termed as the lower bound of the search space.
(6) Analogous to other metaheuristic population-based algorithms, the CPO algorithm necessitates configuring the population size, maximum number of iterations, and fitness function.
(7) Employ the optimal temporal coefficients generated by the algorithm, in conjunction with the spatial modes derived from EOF decomposition, and employ Equation (7) to reconstruct the measured sound speed profiles.

Basic Process
The objective of employing the CPO algorithm is to leverage its rapid convergence advantage in order to swiftly obtain optimal EOF reconstruction coefficients.As depicted in Figure 4, the basic process for Method 2, applied to the designated measurement marine area, is outlined as follows:

Experimental Data
The experimental data utilized in this study were sourced from the global Argo realtime ocean observing network observations, accessed from the website of the Chinese (1)-( 4) follow the steps outlined in Section 2.2 (1)-( 4).
(5) Conduct EOF decomposition on the standard sound speed profile matrix to extract spatial modes and temporal coefficients.Compute the cumulative variance contribution rate of the first K spatial modes, with a minimum threshold of 95%, to determine the minimum number of spatial modes K involved in the reconstruction of measured sound speed profiles.This value of K serves as one of the input parameters for the CPO algorithm, termed as the problem dimension.Simultaneously, identify the maximum and minimum values of the temporal coefficients corresponding to each order of the K order spatial modes.Arrange the maximum values of the temporal coefficients in ascending order of order to form a row vector, constituting the second input parameter of the CPO algorithm, denoted as the upper bound of the search space.Similarly, arrange the minimum values of the temporal coefficients in ascending order of order to form a row vector, constituting the third input parameter of the CPO algorithm, termed as the lower bound of the search space.
(6) Analogous to other metaheuristic population-based algorithms, the CPO algorithm necessitates configuring the population size, maximum number of iterations, and fitness function.
(7) Employ the optimal temporal coefficients generated by the algorithm, in conjunction with the spatial modes derived from EOF decomposition, and employ Equation ( 7) to reconstruct the measured sound speed profiles.

Experimental Data
The experimental data utilized in this study were sourced from the global Argo real-time ocean observing network observations, accessed from the website of the Chinese Argo Real-Time Data Center (http://www.argo.org.cn,accessed on 26 March 2024).A defined region in the Pacific Ocean (latitude 23-26 • N, longitude 135-138 • E) was selected as the study area.A total of 43 station measurements of temperature and salinity profiles collected during February from 2002 to 2008 were chosen as the experimental samples.One station was randomly selected as the validation station.The distribution of stations is illustrated in Figure 5.

Experimental Data
The experimental data utilized in this study were sourced time ocean observing network observations, accessed from th Argo Real-Time Data Center (http://www.argo.org.cn,accessed fined region in the Pacific Ocean (latitude 23-26° N, longitude 1 the study area.A total of 43 station measurements of tempera collected during February from 2002 to 2008 were chosen as t One station was randomly selected as the validation station.T is illustrated in Figure 5.As the Argo data provide temperature and salinity information at specific oceanic points, it is necessary to employ suitable empirical formulas to convert them into sound speed data.Following the comparative analysis conducted by Zhou Fengnian et al. [33], the Mackenzi formula and the EM stratified simplified sound speed formula were found to yield satisfactory computational accuracy across the entire depth range of 1-12,000 m.Because the Mackenzie formula has a relatively simple structure and remains unchanged across various depth layers, it was ultimately selected for sound speed calculation in this study.The formula is as follows: Here, C (m/s) represents sound speed, T ( • C) denotes temperature, S (ppt) signifies salinity, D (m) indicates depth.
Analysis of the sound speed profile data from various stations reveals that sampling depths typically hover around 1000 m.Due to the irregular sampling intervals of Argo, interpolating the sound speed profiles to standard layers becomes necessary for EOF analysis.Hence, this study employs the Akima interpolation method to interpolate each sound speed profile from 0 m to 1000 m with a 1 m depth interval.The interpolated sound speed profiles are illustrated in Figure 6.It is apparent from the figure that the sound speed values within each profile fall within the range of 1480-1540 m/s, indicating a consistent variation in sound speed.The structure of these sound speed profiles adheres to the characteristic features of typical deep-sea sound speed profiles.

24, 12, x FOR PEER REVIEW 10 of 25
As the Argo data provide temperature and salinity information at specific oceanic points, it is necessary to employ suitable empirical formulas to convert them into sound speed data.Following the comparative analysis conducted by Zhou Fengnian et al. [33], the Mackenzi formula and the EM stratified simplified sound speed formula were found to yield satisfactory computational accuracy across the entire depth range of 1-12,000 m.Because the Mackenzie formula has a relatively simple structure and remains unchanged across various depth layers, it was ultimately selected for sound speed calculation in this study.The formula is as follows:

Experimental Process
Assuming a measured depth of 200 m for the sound speed profile, the portions of the tested profile exceeding this depth are truncated, retaining only the data up to 200 m as the target for reconstruction.The experimental process, depicted in Figure 7, involves the application of four distinct methods to reconstruct the aforementioned sound speed profile.Subsequently, an evaluation and analysis of the reconstructed results are conducted using specified metrics.Following this, simulation calculations are performed to trace ray paths based on the reconstructed sound speed profiles obtained from each method, followed by an analysis of the outcomes.
lowed by an analysis of the outcomes.
It is pertinent to note that in practical multibeam operations, when confronted with issues such as untimely sound speed profile measurements or abrupt changes in seabed topography resulting in insufficient longitudinal coverage, operators often resort to contingency measures such as the nearest substitution method or nearest splicing method.This involves substituting the most recently measured or geographically closest sound speed profile for the current one, or cutting and splicing parts of them with the current profile to extend it.In areas with relatively stable sound speed variations, this practice not only simplifies the operation but also ensures the accuracy of sound ray tracing while extending the sound speed profile.However, in regions with complex sound speed variations, these methods may introduce significant measurement errors due to the inability to ascertain whether the sound speed variations between two points are consistent.In order to comprehensively compare the advantages and disadvantages of different methods, this study also includes these two methods as part of the experiments.

Reconstruction of the Measured Sound Speed Profile Using Four Methods
The nearest substitution method involves substituting the profile with the one measured closest in time or location to the current sound speed profile in use.Similarly, the nearest splicing method entails splicing the profile with the one measured closest in time or location to the current sound speed profile, using the splicing method consistent with the procedure outlined in Section 2.2, step (4), which is based on average sound speed splicing.The relevant parameters of the tested sound speed profile and the nearest (either closest in time or location) sound speed profiles are presented in Table 1.It is pertinent to note that in practical multibeam operations, when confronted with issues such as untimely sound speed profile measurements or abrupt changes in seabed topography resulting in insufficient longitudinal coverage, operators often resort to contingency measures such as the nearest substitution method or nearest splicing method.This involves substituting the most recently measured or geographically closest sound speed profile for the current one, or cutting and splicing parts of them with the current profile to extend it.In areas with relatively stable sound speed variations, this practice not only simplifies the operation but also ensures the accuracy of sound ray tracing while extending the sound speed profile.However, in regions with complex sound speed variations, these methods may introduce significant measurement errors due to the inability to ascertain whether the sound speed variations between two points are consistent.In order to comprehensively compare the advantages and disadvantages of different methods, this study also includes these two methods as part of the experiments.

Reconstruction of the Measured Sound Speed Profile Using Four Methods
The nearest substitution method involves substituting the profile with the one measured closest in time or location to the current sound speed profile in use.Similarly, the nearest splicing method entails splicing the profile with the one measured closest in time or location to the current sound speed profile, using the splicing method consistent with the procedure outlined in Section 2.2, step (4), which is based on average sound speed splicing.The relevant parameters of the tested sound speed profile and the nearest (either closest in time or location) sound speed profiles are presented in Table 1.The overall processes of the substitution and splicing methods are straightforward and require no further elaboration.Here, the paper only focuses on providing specific explanations for the proposed methods, namely Method 1 and Method 2.

1.
Method 1 (1) EOF decomposition.Combine the 42 historical sound speed profiles (excluding the tested sound speed profile) interpolated in Section 3.1 into a standard sound speed profile matrix and calculate the average value of each row to obtain the average sound speed profile.In accordance with Section 2.2, concatenate the average sound speed profile with the measured sound speed profile to obtain the sound speed profile for the entire water depth.Combine the extrapolated sound speed profile with the historical sound speed profile matrix to form a new standard sound speed matrix.In accordance with Section 2.1.1,perform EOF decomposition on this matrix to obtain the first seven order spatial modes as shown in Figure 8. (1) EOF decomposition.Combine the 42 historical sound speed profiles (excluding the tested sound speed profile) interpolated in Section 3.1 into a standard sound speed profile matrix and calculate the average value of each row to obtain the average sound speed profile.In accordance with Section 2.2, concatenate the average sound speed profile with the measured sound speed profile to obtain the sound speed profile for the entire water depth.Combine the extrapolated sound speed profile with the historical sound speed profile matrix to form a new standard sound speed matrix.In accordance with Section 2.1.1,perform EOF decomposition on this matrix to obtain the first seven order spatial modes as shown in Figure 8. Analysis of Figure 8 reveals that the first seven order spatial modes (except the 5th mode) tend toward zero in depths greater than 800 m, indicating that there are basically no disturbances in the sound speed profiles in water layers deeper than 800 m, and the sound speed variations tend to stabilize.The first, third, fourth, and sixth spatial modes all exhibit extreme values (either maximum or minimum) in the depth range of 100-200 m, indicating significant variations in the sound speed profiles in this depth interval.Further analysis, in conjunction with Figure 6, reveals that compared with other depth layers, there are indeed significant differences in the sound speed profiles in this water layer.
(2) Determination of spatial mode order.Calculate the variance contribution ratio and cumulative variance contribution ratio of the first K spatial modes.The variation of the variance contribution ratio and cumulative variance contribution ratio of the first 20 spatial modes is shown in Figure 9, and the specific values of the variance contribution ratio and cumulative variance contribution ratio of the first seven spatial modes are shown in Table 2. Analysis of Figure 8 reveals that the first seven order spatial modes (except the 5th mode) tend toward zero in depths greater than 800 m, indicating that there are basically no disturbances in the sound speed profiles in water layers deeper than 800 m, and the sound speed variations tend to stabilize.The first, third, fourth, and sixth spatial modes all exhibit extreme values (either maximum or minimum) in the depth range of 100-200 m, indicating significant variations in the sound speed profiles in this depth interval.Further analysis, in conjunction with Figure 6, reveals that compared with other depth layers, there are indeed significant differences in the sound speed profiles in this water layer.
(2) Determination of spatial mode order.Calculate the variance contribution ratio and cumulative variance contribution ratio of the first K spatial modes.The variation of the variance contribution ratio and cumulative variance contribution ratio of the first 20 spatial modes is shown in Figure 9, and the specific values of the variance contribution ratio and cumulative variance contribution ratio of the first seven spatial modes are shown in Table 2.

25% in cumulative varianc
There is a sharp decrease in the variance contribution ratio of the thi lowed by a gradual decrease in the variance contribution ratio of a modes.Therefore, it can be inferred that the main information of th surveyed area is concentrated in the first two spatial modes.The cum  From Table 2, it can be observed that the variance contribution ratio of the first two spatial modes is the highest, reaching 82.25% in cumulative variance contribution ratio.There is a sharp decrease in the variance contribution ratio of the third spatial mode, followed by a gradual decrease in the variance contribution ratio of all subsequent spatial modes.Therefore, it can be inferred that the main information of the sound field in the surveyed area is concentrated in the first two spatial modes.The cumulative variance contribution rate of the first seven spatial modes is 95.92%, which is greater than 95%.Therefore, it is decided to reconstruct the measured sound speed profiles using the first seven spatial modes.
(3) Reconstruction of measured sound speed profiles.Reconstruct the measured sound speed profiles according to Equation (7).

2.
Method 2 (1) EOF decomposition.Combine the 42 historical sound speed profiles (excluding the tested sound speed profile) interpolated in Section 3.1 into a standard sound speed profile matrix.The average value of each row is then calculated to obtain the average sound velocity profile.In accordance with Section 2.2, the average sound velocity profile is concatenated with the measured sound velocity profile to obtain the sound velocity profile over the entire water depth.Following Section 2.1.1,EOF decomposition is performed on this matrix to obtain the first seven spatial modes as shown in Figure 10.(2) Determination of spatial mode order.The variance contribution ratio and cumulative variance contribution ratio of the first K spatial modes are calculated.The variation of the variance contribution ratio and cumulative variance contribution ratio for the first 20 spatial modes is shown in Figure 11, and the specific values of the variance contribution ratio and cumulative variance contribution ratio for the first seven spatial modes are shown in Table 3.  (2) Determination of spatial mode order.The variance contribution ratio and cumulative variance contribution ratio of the first K spatial modes are calculated.The variation of the variance contribution ratio and cumulative variance contribution ratio for the first 20 spatial modes is shown in Figure 11, and the specific values of the variance contribution ratio and cumulative variance contribution ratio for the first seven spatial modes are shown in Table 3. (2) Determination of spatial mode order.The variance contribution ratio and cumulative variance contribution ratio of the first K spatial modes are calculated.The variation of the variance contribution ratio and cumulative variance contribution ratio for the first 20 spatial modes is shown in Figure 11, and the specific values of the variance contribution ratio and cumulative variance contribution ratio for the first seven spatial modes are shown in Table 3.    2 and 3 reveal that the trends in the variance contribution rate and cumulative variance contribution rate of each mode obtained by EOF decomposition using Method 1 and Method 2 are consistent, with only minor differences in numerical values.As the cumulative variance contribution rate of the first seven spatial modes is 95.90%, which exceeds 95%, the first seven spatial modes are chosen to reconstruct the measured sound velocity profile.
(3) Determination of upper and lower bounds for the CPO algorithm.The upper and lower bounds of the temporal coefficients corresponding to the first seven spatial modes are shown in Table 4.The upper and lower bounds are arranged into 7-dimensional row vectors in ascending order of mode number, forming the upper and lower bounds of the algorithm search space.
(4) Initialization of algorithm parameters.According to Section 2.3.2, the population size is set to 40, the maximum number of iterations is set to 500, and the problem dimension is set to 7, representing the seven reconstruction coefficients to be obtained.The fitness function is set as the RMSE between the reconstructed sound velocity profile and the extrapolated measured sound velocity profile.After parameter initialization, the optimization process begins.
(5) Reconstruction of measured sound velocity profile.Following Equation ( 7), the measured sound velocity profile is reconstructed.

Constant Gradient Sound Ray Tracing and Accuracy Evaluation
In addition to various metric evaluations for the reconstructed sound velocity profile, it is imperative to validate its ability to meet the requirements for sound velocity correction.This underscores the significance of reconstructing the sound velocity profile, as it serves the purpose of multibeam bathymetry, with bathymetric accuracy being the ultimate criterion for assessment.There are several methods for sound velocity correction, such as the equivalent sound velocity profile method [34,35], average sound velocity method [36], and ray tracing method [37].However, the prevailing multibeam bathymetric systems and post-processing software currently employ the constant gradient ray tracing method for the position calculation of sounding points.Despite its computational complexity, the constant gradient ray tracing method is widely adopted due to its high accuracy in sound velocity correction.The principles and computational methods of this method are detailed in reference [38] and will not be reiterated in this paper.
Regarding the maximum allowable error in bathymetry, the 6th edition of the International Hydrographic Organization (IHO) Standards for Hydrographic Surveys [39] stipulates that for second-order depth measurements (in waters deeper than 200 m), the maximum error of sounding points at a 95% confidence level should not exceed 20 m plus 10% of the depth.The formula for calculating the depth error limit is as follows: Here, TVU max (m) represents the maximum allowable depth error, d(m) denotes the depth, a and b are coefficients.For second-order depth measurements, a = 1.0 (m), b = 0.023 (dimensionless coefficient).
The 2022 edition of the Chinese National Standard for Hydrographic Surveys specifies that for fourth-order surveys (in waters deeper than 200 m) or maps with scales smaller than 1:50,000, the maximum allowable error of sounding points at a 95% confidence level should not exceed 20 m, and the maximum allowable error in depth should not exceed 2% of the depth.
Both standards provide overall regulations for the positional and depth errors of sounding points but do not further refine the limits of errors caused by the sound velocity profile.The National Oceanic and Atmospheric Administration (NOAA) of the United States established in 1999 that a deviation exceeding 0.25% of the depth due to differences between the used sound velocity profile and the actual profile constitutes an over-limit condition.However, this regulation only restricts the depth error caused by the sound velocity profile and does not provide clarification on the positional errors it induces.Therefore, in the accuracy evaluation conducted in this study, NOAA's regulation is adopted for depth errors.Regarding positional errors, considering that the Chinese National Standard for Hydrographic Surveys is more stringent than the IHO Standards for Hydrographic Surveys and to ensure measurement accuracy, it is decided to adhere to the regulations of this standard.
We first simulate the generation of 131 beams, with a beam angle of 1 • and an incidence angle ranging from −60 • to 60 • .We then utilize the constant gradient sound ray tracing method to conduct ray tracing at various depths using reconstructed sound velocity profiles obtained from four different methods, and compare these results with the ray tracing outcomes using the tested sound velocity profile.This aims to assess the efficacy of the sound velocity corrections at different depths as furnished by the reconstructed sound velocity profiles.

Evaluation Metrics
To comprehensively evaluate the reconstruction performance of each method, RMSE (Root Mean Squared Error), MAE (Mean Absolute Error), MAPE (Mean Absolute Percentage Error), and R 2 (R-squared) are selected as evaluation metrics.RMSE measures the magnitude of the reconstruction error, i.e., the difference between the reconstructed values and the true values.A smaller RMSE indicates better consistency between the reconstructed and true values.This metric is sensitive to large errors and can effectively reflect the model's ability to handle outliers or extreme values.MAE quantifies the average magnitude of differences between reconstructed and true values.Although it is less sensitive than RMSE, it provides a more intuitive result, with smaller values indicating better model performance.MAPE measures the relative error between reconstructed and true values, reflecting the extent to which reconstructed values deviate from true values.Smaller MAPE values indicate better reconstruction performance of the model.R 2 measures the model's fitting ability to the data.It is a dimensionless metric ranging from 0 to 1, with values closer to 1 indicating better model fitting performance.The formulas for calculation are as follows: Here, y i represents the sound velocity value at the i depth in the measured sound velocity profile, ŷi represents the sound velocity value at the i depth in the reconstructed sound velocity profile, y represents the average sound velocity in the measured sound velocity profile, and n represents the number of sound velocity points.

Analysis of Reconstruction Results with Fixed Sampling Depth of the Measured Sound Speed Profile
During the experiment, the sampling depth of the measured sound speed profile was set to 200 m, and four methods were employed to reconstruct it.Both Method 1 and Method 2 as proposed in this paper are based on the principle of EOF.The reconstruction coefficients obtained by the two methods are shown in Table 5, indicating that the coefficients obtained by both methods are generally consistent.The reconstructed sound speed profiles and corresponding metrics for each method are shown in Figure 12 and Table 6, respectively.From Table 6, it can be observed that for the measured sound speed profile with a sampling depth of 200 m used in this experiment, Method 2 achieves RMSE, MAE, MAPE, and R2 for the full-depth (1000 m) reconstructed sound speed profile of 0.9511 m/s, 0.8492 m/s, 0.0566%, and 0.9963, respectively.These metrics surpass those of other methods, sig-  From Table 6, it can be observed that for the measured sound speed profile with a sampling depth of 200 m used in this experiment, Method 2 achieves RMSE, MAE, MAPE, and R2 for the full-depth (1000 m) reconstructed sound speed profile of 0.9511 m/s, 0.8492 m/s, 0.0566%, and 0.9963, respectively.These metrics surpass those of other methods, signifying that Method 2 can more precisely depict the sound speed profile across the entire water depth, leveraging the measured sound speed profile at a specific location.While Method 1's performance across various metrics is marginally inferior to that of Method 2, the disparities are negligible.In contrast to the nearest substitution and nearest splicing approaches, Method 1 maintains evident superiority.The nearest substitution and nearest splicing approaches do not exhibit absolute superiority or inferiority between them.For each approach individually, utilizing the sound speed profile closest in time yields superior reconstruction effects compared to employing the sound speed profile nearest in location.Importantly, the RMSE of the reconstructed sound speed profile, acquired by splicing the nearest sound speed profile in location with the measured one, stands at 2.0551 m/s, markedly exceeding the corresponding value for other reconstructed profiles.This occurs because during the splicing process, the difference between the sound speed at the end of the measured sound speed profile and the sound speed at the same depth of the spliced sound speed profile needs to be calculated and added to all sound speeds of the spliced sound speed profile that are deeper than the sampling depth.As shown in Figure 12, the sound speed profile nearest in location is closer to the tested sound speed profile originally, but after adding the difference, the profile shifts to the left, away from the tested sound speed profile, leading to increased error.Therefore, whether nearest in time or nearest in location, the splicing approach may confront these issues, potentially yielding errors surpassing those of direct substitution.Only when the sound speed trends of two points align or the disparity in sound speed remains minimal can the splicing approach yield superior outcomes.Nevertheless, predicting sound speed fluctuations between two points proves unfeasible in practice.Thus, in regard to stability, the method advanced in this study also holds merit.
The sound speed error curves of different depths of each reconstructed sound speed profile are shown in Figure 13, and the corresponding statistical parameters are shown in Table 7.The sound speed error curves of different depths of each reconstructed sound speed profile are shown in Figure 13, and the corresponding statistical parameters are shown in Table 7.   14 illustrates that the sound speed errors of both Method 1 and Method 2 range from 0-2 m/s, with more than half of the errors falling below 1 m/s.When splicing the sound speed profile nearest in time with the measured sound speed profile, although 69% of sound speed errors are less than 1 m/s, 3% of sound speed errors are distributed between 2-3 m/s, resulting in inferior reconstruction results compared to Method 1 and Method 2. When substituting the measured sound speed profile with the sound speed profile nearest in time, 72% of sound speed errors are less than 1 m/s, which is the highest proportion among all reconstructed sound speed profiles.However, this method generates a wide range of sound speed errors, ranging from 0-4 m/s, leading to unsatisfactory reconstruction results.The method of splicing the sound speed profile nearest in location performs the worst, as evident from the figure, characterized by the broadest range of error distribution.According to Table 7, the maximum sound speed error reaches 4.2695 m/s. Figure 14 illustrates that the sound speed errors of both Method 1 and Method 2 range from 0-2 m/s, with more than half of the errors falling below 1 m/s.When splicing the sound speed profile nearest in time with the measured sound speed profile, although 69% of sound speed errors are less than 1 m/s, 3% of sound speed errors are distributed between 2-3 m/s, resulting in inferior reconstruction results compared to Method 1 and Method 2. When substituting the measured sound speed profile with the sound speed profile nearest in time, 72% of sound speed errors are less than 1 m/s, which is the highest proportion among all reconstructed sound speed profiles.However, this method generates a wide range of sound speed errors, ranging from 0-4 m/s, leading to unsatisfactory reconstruction results.The method of splicing the sound speed profile nearest in location performs the worst, as evident from the figure, characterized by the broadest range of error distribution.According to

Analysis of Reconstruction Results with Variable Sampling Depth of the Measured Sound Speed Profile
The analysis in Section 4.1 is conducted based on a measured sound speed profile with a sampling depth of 200 m.In order to comprehensively evaluate the performance of each method under different sampling depths of the measured sound speed profile, sampling depths are set from 0 m to 1000 m at intervals of 1 m.RMSE is employed as the evaluation metric.The performance of each method under varying sampling depths is depicted in Figure 15.
Figure 14.Distribution of sound speed errors for each reconstructed sound speed profile.The x-axis represents the range of sound speed errors, and the y-axis represents the ratio of sound speed points, whose errors fall within a certain range, to the total number of sound speed points in the sound speed profile.

Analysis of Reconstruction Results with Variable Sampling Depth of the Measured Sound Speed Profile
The analysis in Section 4.1 is conducted based on a measured sound speed profile with a sampling depth of 200 m.In order to comprehensively evaluate the performance of each method under different sampling depths of the measured sound speed profile, sampling depths are set from 0 m to 1000 m at intervals of 1 m.RMSE is employed as the evaluation metric.The performance of each method under varying sampling depths is depicted in Figure 15.Examination of Figure 15 reveals that for sampling depths of the measured sound speed profile below 152 m, both Method 1 and Method 2 demonstrate suboptimal reconstruction efficacy relative to the nearest substitution and nearest splicing methodologies.Method 1 exceeds an RMSE of 3.5 m/s, showcasing the least favorable performance among all methodologies.Conversely, Method 2 surpasses Method 1, with its efficacy steadily ameliorating as the sampling depth extends past 100 m.Upon reaching a depth of 147 m, both Method 1 and Method 2 commence demonstrating nearly indistinguishable efficacy, sustaining consistent performance with increasing sampling depth.Upon reaching a depth of 162 m, both methodologies surpass the substitution and splicing approaches, attaining an RMSE of under 1 m/s in sound speed profile reconstruction.This performance gradually stabilizes, providing a reference for determining the sampling depth of the measured sound speed profile.Combining this analysis with Figure 8, it is observed that the significant variations in the spatial modes occur predominantly within depths ranging from 100 m to 200 m.As the sampling depth of the measured sound speed profile surpasses this layer, the reconstruction performance of both Method 1 and Method 2 improves.Therefore, it is recommended that to attain satisfactory reconstruction outcomes, the sampling depth of the measured sound speed profile ideally should surpass the depths harboring substantial spatial mode variations, which are typically situated within shallow water layers.

Analysis of Sound Ray Tracing Results
To enhance computational efficiency, sound ray tracing was initiated at a depth of 50 m, with intervals of 25 m until the maximum water depth (1000 m) was reached.The results of sound ray tracing for the four methods indicate that except for the substitution method (location nearest), which exhibited some instances of edge beam measurements exceeding the limit, the methods all met the accuracy requirements proposed in Section 3.2.2 and the effective beam ratio reached 100%.Regarding the effective beam ratio, Zhu Xiaochen [40] pointed out in his doctoral dissertation that it represents the ratio of beams meeting the depth tolerance to the total number of beams, quantifying the effectiveness of sound speed correction for multibeam measurement data.It is worth noting that, in Table 6, the metrics of the substitution method (location nearest) are actually superior to those of the splicing method (location nearest).However, the results of sound ray tracing for the substitution method (location nearest) are inferior to those of the splicing method (location nearest).The reason is that, in the shallower part of the splicing depth, the sound speed profile nearest in location differs significantly from the actual sound speed profile (as shown in Figure 12).During sound ray tracing, errors accumulate layer by layer.In contrast, the splicing method matches the actual sound speed profile in the shallower depths, resulting in no errors during sound ray tracing in this layer.Although the difference between the splicing method and the actual sound speed profile is greater in the deeper part of the splicing depth compared to the substitution method, the absence of error accumulation in the shallower part results in higher overall tracking accuracy for the splicing method.
Figure 16 illustrates the sound ray tracing errors for Methods 1 and 2, including depth errors and position errors of depth measurement points.The trends of sound ray tracing errors for each method are similar.Due to space limitations, this paper specifically analyzes only the sound ray tracing errors of Method 2.
12, x FOR PEER REVIEW 22 of 25 are to the sound ray tracing results of the actual sound speed profile.When the depth is constant, as the incident angle of the beam increases, the depth error gradually increases, and the sound ray tracing results of edge beams are significantly inferior to those of central beams.
Figure 16d illustrates the position errors of sound ray tracing for Method 2. Compared to depth errors, position errors do not exhibit a singular trend of increasing or decreasing with depth.Instead, they fluctuate, indicating that at certain depths, the sound ray tracing results are better.However, for any depth, the position error is far below the accuracy indicators proposed in Section 3.2.2.

Conclusions
During the execution of a deep-sea multibeam bathymetric survey, the acquisition of representative sound speed profiles for the surveyed sea areas emerges is a pivotal task,  Figure 16c depicts the depth errors of sound ray tracing for Method 2. Here, "depth error" refers to the absolute difference between the sound speed correction mutual difference of the actual sound speed profile and the reconstructed sound speed profile and 0.25% of the water depth value.If this value is less than 0, it indicates that the sound speed correction results of the reconstructed sound speed profile meet the measurement accuracy requirements.Overall, as the depth increases, the mutual difference in depth becomes increasingly smaller than 0.25% of the water depth, indicating that the deeper the water, the closer the sound ray tracing results using the reconstructed sound speed profile are to the sound ray tracing results of the actual sound speed profile.When the depth is constant, as the incident angle of the beam increases, the depth error gradually increases, and the sound ray tracing results of edge beams are significantly inferior to those of central beams.
Figure 16d illustrates the position errors of sound ray tracing for Method 2. Compared to depth errors, position errors do not exhibit a singular trend of increasing or decreasing with depth.Instead, they fluctuate, indicating that at certain depths, the sound ray tracing results are better.However, for any depth, the position error is far below the accuracy indicators proposed in Section 3.2.2.

Conclusions
During the execution of a deep-sea multibeam bathymetric survey, the acquisition of representative sound speed profiles for the surveyed sea areas emerges is a pivotal task, indispensable for guaranteeing the quality of the collected data.Nonetheless, the acquisition of full-depth sound speed profiles presents formidable challenges, owing to factors encompassing the effective detection depth of equipment, meteorological conditions, variations in seabed topography, and operational efficiency considerations.The variation of sound speed in seawater, though intricate, still exhibits discernible regularities on specific spatiotemporal scales within the remote ocean.To tackle this complexity, an approach that involves leveraging a substantial dataset consisting of historical full-depth sound speed profiles obtained from maritime regions is proposed.Through the application of EOF analysis, the aim is to identify the underlying patterns of variability, which can facilitate the reconstruction of partially-depth measured sound speed profiles.To mitigate the limitation of EOF reconstruction depth posed by the shallowest historical sound speed profiles, the measured sound speed profile is integrated with the historical mean sound speed profile to create comprehensive full-depth representations.Building upon this framework, two methodologies are presented for reconstructing the measured sound speed profile.
In the experiments utilizing Argo data, the sampling depth of the measured sound speed profile is set at 200 m.Method 1 and Method 2 are compared with commonly used methods in practical measurements, such as the substitution method and the splicing methods.The experimental results indicate that Method 2 outperforms other methods, yielding a reconstructed sound speed profile with superior performance.This is evidenced by lower RMSE, MAE, MAPE, and higher R2 values compared to the alternatives.While Method 1 shows slightly inferior performance to Method 2, the discrepancies between the two methods are negligible, and both notably surpass the substitution and splicing methods.It is important to note that the efficacy of the substitution and splicing methods heavily relies on the hydrological conditions of the surveyed sea area.
Further experimentation involving varying depths of the measured sound speed profile reveals that Method 1 and Method 2 exhibit decreased effectiveness compared to the substitution method and the splicing methods when the measured sound speed profile is shallow.However, the experimental results demonstrate that that once the depth of the measured sound speed profile surpasses the depths where the primary spatial modes undergo significant changes, the performance of Method 1 and Method 2 improves significantly.It is recommended that measurement operators take this depth into consideration when determining the measurement depth of the sound speed profile.Compared to the current measurement depth, which typically extends to at least the depth where the sound of 24 channel axis appears (approximately 1000-1500 m), this optimized depth is considerably reduced, thereby enhancing operational efficiency significantly.
Finally, sound ray tracing of reconstructed sound speed profiles is conducted at different depths.Apart from the substitution method (location nearest) which exhibits some instances of edge beam measurements exceeding the limit, all the other methods meet the measurement accuracy requirements, with an effective beam ratio of 100%.
In future research, we will focus on improving the performance of the proposed methods in reconstructing sound speed profiles when the sampling depth is shallow and explore other methods with superior performance for sound speed profile reconstruction.

Figure 1 .
Figure 1.Process of EOF decomposition and reconstruction of sound speed profiles.

Figure 3 .
Figure 3. Visualization of CPO search space.

Figure 3 .
Figure 3. Visualization of CPO search space.

Figure 5 .Figure 5 .
Figure 5.The distribution of the sound profile stations, where the re validation station, the yellow solid circle denotes the station closest station, and the green solid circle indicates the station closest in measu station, with other points representing regular stations.
(m/s) represents sound speed, T (°C) denotes temperature, S (ppt) sig- nifies salinity, D (m) indicates depth.Analysis of the sound speed profile data from various stations reveals that sampling depths typically hover around 1000 m.Due to the irregular sampling intervals of Argo, interpolating the sound speed profiles to standard layers becomes necessary for EOF analysis.Hence, this study employs the Akima interpolation method to interpolate each sound speed profile from 0 m to 1000 m with a 1 m depth interval.The interpolated sound speed profiles are illustrated in Figure6.It is apparent from the figure that the sound speed values within each profile fall within the range of 1480-1540 m/s, indicating a consistent variation in sound speed.The structure of these sound speed profiles adheres to the characteristic features of typical deep-sea sound speed profiles.

Figure 6 .
Figure 6.Sound speed profile curves for each station, where the red solid line represents the sound speed profile of the validation station, the green solid line represents the sound speed profile closest in measurement time to the validation profile, the yellow solid line represents the sound speed profile closest in measurement location to the validation profile, and the light blue solid lines represent the sound speed profiles of the remaining stations.

3. 2 .Figure 6 .
Figure 6.Sound speed profile curves for each station, where the red solid line represents the sound speed profile of the validation station, the green solid line represents the sound speed profile closest in measurement time to the validation profile, the yellow solid line represents the sound speed profile closest in measurement location to the validation profile, and the light blue solid lines represent the sound speed profiles of the remaining stations.

Figure 8 .
Figure 8.The first 7 order spatial modes obtained by EOF decomposition of the new standard sound speed matrix.

Figure 8 .
Figure 8.The first 7 order spatial modes obtained by EOF decomposition of the new standard sound speed matrix.

Figure 9 .
Figure 9. Variation of variance contribution ratio and cumulative variance c first 20 spatial modes.

Figure 9 .
Figure 9. Variation of variance contribution ratio and cumulative variance contribution ratio of the first 20 spatial modes.

J 25 Figure 10 .
Figure 10.The first 7 spatial modes obtained from EOF decomposition of the standard sound velocity matrix.

Figure 11 .
Figure 11.Variation of variance contribution ratio and cumulative variance contribution ratio of the first 20 spatial modes.

Figure 10 .
Figure 10.The first 7 spatial modes obtained from EOF decomposition of the standard sound velocity matrix.

J 25 Figure 10 .
Figure 10.The first 7 spatial modes obtained from EOF decomposition of the standard sound velocity matrix.

Figure 11 .
Figure 11.Variation of variance contribution ratio and cumulative variance contribution ratio of the first 20 spatial modes.

Figure 11 .
Figure 11.Variation of variance contribution ratio and cumulative variance contribution ratio of the first 20 spatial modes.

J 25 Figure 12 .
Figure 12.Reconstructed sound speed profiles obtained by four methods.

Figure 13 .
Figure 13.Sound speed error curves of each reconstructed sound speed profile.

Figure 13 .
Figure 13.Sound speed error curves of each reconstructed sound speed profile.

Method 2 Figure 14 .
Figure14.Distribution of sound speed errors for each reconstructed sound speed profile.The x-axis represents the range of sound speed errors, and the y-axis represents the ratio of sound speed points, whose errors fall within a certain range, to the total number of sound speed points in the sound speed profile.

Figure 15 .Figure 15 .
Figure 15.Performance of each method in reconstructing sound speed profiles with variable sampling depths.

Figure 16 .
Figure 16.Sound ray tracing results for Methods 1 and 2. (a) Depth error for Method 1.(b) Position error for Method 1. (c) Depth error for Method 2. (d) Position error for Method 2.

Figure 16 .
Figure 16.Sound ray tracing results for Methods 1 and 2. (a) Depth error for Method 1.(b) Position error for Method 1. (c) Depth error for Method 2. (d) Position error for Method 2. of24

Table 1 .
Relevant parameters of tested sound speed profile and nearest (time/location) sound speed profiles.

Table 1 .
Relevant parameters of tested sound speed profile and nearest (time/location) sound speed profiles.

Table 2 .
Variance contribution ratio and cumulative variance contribution ra modes.

Table 2 ,
it can be observed that the variance contribution spatial modes is the highest, reaching 82.

Table 2 .
Variance contribution ratio and cumulative variance contribution ratio of the first 7 spatial modes.

Table 3 .
Variance contribution ratio and cumulative variance contribution ratio of the first 7 spatial modes.

Table 4 .
Upper and lower bounds of temporal coefficients corresponding to the first 7 spatial modes.

Table 5 .
Reconstruction coefficients obtained by Method 1 and Method 2 for sound speed profile.
Figure 12.Reconstructed sound speed profiles obtained by four methods.

Table 6 .
RMSE, MAE, MAPE, and R 2 of each reconstructed sound speed profile.

Table 7 .
Statistical parameters of sound speed errors for each reconstructed sound speed profile.

Table 7 .
Statistical parameters of sound speed errors for each reconstructed sound speed profile.

Table 7 ,
the maximum sound speed error reaches 4.2695 m/s.
Figure14.Distribution of sound speed errors for each reconstructed sound speed profile.The x-axis represents the range of sound speed errors, and the y-axis represents the ratio of sound speed points, whose errors fall within a certain range, to the total number of sound speed points in the sound speed profile.4.2.Analysis of Reconstruction Results with Variable Sampling Depth of the Measured Sound Speed ProfileThe analysis in Section 4.1 is conducted based on a measured sound speed profile