Reconstruction of Ocean Front Model Based on Sound Speed Clustering and Its Effectiveness in Ocean Acoustic Forecasting

: As a mesoscale phenomenon of the ocean, the ocean front can directly affect the structural characteristics of sound speed proﬁles and further affect the acoustic propagation characteristics of the sea area. In this paper, we use the fuzzy C-means (FCM) algorithm to cluster the surface sound speed in the sea area of the Kuroshio Extension (KE) and detect the frontal zone of Kuroshio Extension (KEF). At the same time, the sound speed proﬁle (SSP) is used instead of the temperature proﬁle to establish the model of the sound speed ﬁeld in the front area of the Kuroshio Extension and to improve the theoretical model of the ocean front. Compared with the actual ocean front calculated by reanalysis data, the root means square error (RSME) of the transmission loss (TL) calculated by the model is controlled below 6 dB, which proves the validity of the model. Finally, we propose the melt function in the model to forecast the depth change of the acoustic convergence area. Compared with the actual calculation result based on reanalysis data, the root means square error (RSME) of the depth forecasting after the frontal zone is 43.3 m. This reconstruction method does not rely on the high spatial resolution data of the whole sea depth and can be of referential signiﬁcance to acoustic detection in the ocean front environment. d ) 1 August 2019. The detection results of the frontal zone were similar to those of the FCM, especially in winter. However, the continuity of the judgment results of the temperature gradient method is very poor compared with those of FCM, and this discontinuity may become worse due to different judgment standards. Compared with the traditional method of calculating the temperature gradient to detect the frontal zone, FCM is obviously more reliable in our research. environment, and the SD = ( a ) 100 m, ( b ) 300 m, and ( c ) 500 m; the TL of the model, and the SD = ( d ) 100 m, ( e ) 300 m, and ( f ) 500 m; the TL of actual ocean front, and the SD = ( g ) 100 m, ( h ) 300 m, and ( i ) 500 m. Due to the existence of the ocean front, there was a signiﬁcant downward shift in the distribution of TL after its passing through the frontal zone in the convergence area, which was a very obvious phenomenon is in both the model and the actual ocean front environments. The above phenomenon can help us achieve the forecasting of the depth change of the acoustic convergence area under ocean front environments. actual results, which proves the effectiveness of this method. It proves that the melt function was applied to effectively achieve the depth forecast of the convergence area in the ocean front environment.


Introduction
As an important phenomenon in marine dynamics, the ocean front represents a narrow transition zone between two different water masses [1]. The frontal zone is located in the convergence area, where there is a strong vertical movement, and the exchange of momentum, heat, and water vapor is very active [2]. Therefore, the ocean front has become an important research field of physical oceanography. The position and strength of the ocean front change with time [3]. For the detection of the ocean front, the dynamic method or the gradient judgment method of hydrological elements including temperature can be used [4,5], in which a critical value should be given first. If the calculated horizontal temperature gradient is greater than this critical value, it is considered that there is an ocean front [6].
The influence of the ocean front on acoustic propagation is also of great concern to hydroacoustic scholars [7]. When the horizontal gradients of temperature and salinity change, the ocean front appears in this area, which affects the structure of sound speed profile (SSP) [8], then the acoustic propagation characteristics of the sea area [9], and finally, the underwater detection [10]. From the perspective of the acoustic effect of the ocean front, the traditional dynamic methods or gradient methods can no longer meet the needs of research since the acoustic effect of the ocean front is based on SSP, rather than a single hydrological element. In recent years, there has been a new research method of using the fuzzy cluster analysis of SSP for detecting the ocean front. Fuzzy cluster analysis refers to constructing a fuzzy matrix according to the properties of the research object itself and determining the cluster relationship according to certain membership, aiming to perform clustering objectively and accurately [11]. Common fuzzy clustering algorithms include the k-means algorithm, fuzzy C-means (FCM) algorithm, and so on. The classification research of the SSP by cluster analysis was first used in the 21st century. Mandelberg [12] used the hierarchical clustering method to classify the SSP of the North Atlantic and the Northeast Pacific. Abiva et al. [13] used principal component analysis (PCA) and self-organizing map (SOM) to automatically cluster the SSP, which was applicable in terms of the sea area of the Strait of Gibraltar. In addition, Dubberley and Zingerelli [14] applied fuzzy clustering to oceanographic parameters related to acoustics and divided them into multiple categories, which proved the practicability of the method. Although the use of various cluster analysis methods to study the SSP has achieved some gratifying results, the necessity of using the high spatial resolution data of the whole sea depth may lead to some inconveniences in the actual ocean exploration [15].
Although many works have been carried out on the modeling of the ocean front [16][17][18], there is still little research conducted on the combination of theoretical model and actual ocean front. Since the previous models are still based on single parameters such as temperature or current speed [19], it is believed by some hydroacoustic scholars that they are no longer suitable for the research of the ocean front model of the sound speed field. Small J. et al. [20] by comparing the ocean front temperature model with the measured data, the effectiveness of the model in the sound field prediction was basically verified, and the root means square error (RMSE) calculated by the model was within 7 dB compared with the measured one.
In this paper, in order to reduce the dependence of ocean front modeling on the data of the whole sea area and the whole sea depth, we only cluster the surface sound speed of the sea with the fuzzy C-means (FCM) algorithm, to obtain the position and range of the frontal zone. At the same time, the theoretical model of the ocean front is improved by replacing the temperature profile with the SSP, and the model of the sound speed field of the Kuroshio Extension front (KEF) is established quickly and effectively. We compare the model with the actual sound speed distribution of the KEF calculated by reanalysis data and verify the validity of the model from the perspective of transmission loss (TL). Finally, it is found that the melt function in this model can forecast the depth of the acoustic convergence area in the KEF, which can be of reference to acoustic detection in the ocean front environment. Section 1 of this paper is the introduction. In Section 2 of this paper, the data and methods are introduced, including the physical oceanography and data introduction, the description of the ocean front characteristic model and the introduction of the technical route. Section 3 of this paper is results and discussion, including reconstruction results of ocean front sound speed field model and comparison results of TL calculations, and the effectiveness of forecasting the changes of depth in convergence areas. Section 4 of this paper is a summary and conclusion, and we summarize the whole research process and give relevant conclusions.

Physical Oceanography and Data Introduction
Kuroshio is the strongest warm current in the Northwest Pacific (NP) [21][22][23]. Due to the difference in intensity and position of Kuroshio in different seasons, the range and intensity of the KEF also show obvious seasonal differences [24,25]. Kuroshio and its extension have attracted much attention as ocean phenomena [26], but the changes in the underwater acoustic environment caused by Kuroshio and its extension on sound propagation have not been effectively analyzed. In this paper, our study area is selected within the range shown in Figure 1a. The KEF divides the sea area into three water masses, namely, the low temperature and low salinity water mass in the north of the frontal zone, the high temperature and high salinity water mass in the south of the frontal zone, and the frontal transitional water mass between them. the high temperature and high salinity water mass in the south of the frontal zone, and the frontal transitional water mass between them. The sea surface sound speed distribution of the study area on 1 January 2019. The KEF divides the sea area into three water masses, namely, the low temperature and low salinity water mass in the north of the frontal zone, the high temperature and high salinity water mass in the south of the frontal zone, and the frontal transitional water mass between them.
The data source used in this paper is the daily mean temperature, salinity, and depth data (reanalysis data) output by HYCOM [27,28] (the Hybrid Ocean Model) (since the ocean front in winter is significantly different from that in summer, we used the daily mean data in January, February, July, and August). The type is NetCDF (Network Common Data Format) with a horizontal resolution of 0.08° and vertical stratification of 40 layers. The speed of sound in seawater changes with temperature, salinity, and depth, and empirical Equations are usually used to express the relationship among them. A simplified empirical Equation was used to calculate the sound speed: where c refers to the speed of sound, T the temperature, S the salinity, and D the depth. This empirical Equation was established by Herman and Medwin [29]. As an experimental summary of a large number of sea sound speed measurement data, the empirical Equation has proven the accuracy of the Equation to be used directly. Since the sea depth is generally deep and the seabed terrain is relatively flat in the study area, we uniformly set the sea depth at 5000 m to eliminate the impact of the terrain on the sound propagation. Through the above settings, the 144° E was selected as the characteristic section (the direction is 38° N-36° N) to obtain the sound speed distribution of the KEF in different periods ( Figure 2). It was found that due to the difference in time, the position and intensity of the KEF also vary greatly. The sea surface sound speed distribution of the study area on 1 January 2019. The KEF divides the sea area into three water masses, namely, the low temperature and low salinity water mass in the north of the frontal zone, the high temperature and high salinity water mass in the south of the frontal zone, and the frontal transitional water mass between them.
The data source used in this paper is the daily mean temperature, salinity, and depth data (reanalysis data) output by HYCOM [27,28] (the Hybrid Ocean Model) (since the ocean front in winter is significantly different from that in summer, we used the daily mean data in January, February, July, and August). The type is NetCDF (Network Common Data Format) with a horizontal resolution of 0.08 • and vertical stratification of 40 layers. The speed of sound in seawater changes with temperature, salinity, and depth, and empirical Equations are usually used to express the relationship among them. A simplified empirical Equation was used to calculate the sound speed: where c refers to the speed of sound, T the temperature, S the salinity, and D the depth. This empirical Equation was established by Herman and Medwin [29]. As an experimental summary of a large number of sea sound speed measurement data, the empirical Equation has proven the accuracy of the Equation to be used directly. Since the sea depth is generally deep and the seabed terrain is relatively flat in the study area, we uniformly set the sea depth at 5000 m to eliminate the impact of the terrain on the sound propagation. Through the above settings, the 144 • E was selected as the characteristic section (the direction is 38 • N-36 • N) to obtain the sound speed distribution of the KEF in different periods ( Figure 2). It was found that due to the difference in time, the position and intensity of the KEF also vary greatly.

Description of Ocean Front Characteristic Model
For the study on the characteristic model of ocean front, we refer to the parameterized two-dimensional characteristic model of ocean front established by Carriere O. et al. [17]: where r and z are the horizontal and vertical coordinates, T 0 and T i are the temperature profiles, and m is the melt function. The ocean front temperature model proposed by Carriere O. is relatively simple in form. The input value is only the temperature profile, and a better ocean front modeling effect can be obtained. Moreover, the strength of the ocean front and the position of the frontal zone can also be achieved by adjusting the melt function. We improved the above model by replacing the temperature profile with the SSP to represent the two-dimensional parameterized characteristic model of ocean front based on the SSP: Appl. Sci. 2021, 11, 8461 4 of 22 where C 1 and C 2 are the sound speed profiles and m is the melt function (which is consistent with Equation (2)). The expression is as follows: where R is the distance between the two sound speed profiles, and the melt parameter a varies from −1.5 to 1.5. In this paper, since the distance between the north and south of the study area extends across two latitudes, R is defined as 200 km. Figure

Description of Ocean Front Characteristic Model
For the study on the characteristic model of ocean front, we refer to the parameter-  of the study area extends across two latitudes, R is defined as 200 km. Figure 3 shows the function curve of melt function m with the melt parameter a in the range of −1.5-1.5. Different melt parameters correspond to different frontal zone positions. Here, we define the initial position of the frontal zone as the position corresponding to when the value of the melt function m is 0.1, and the termination position of the frontal zone as the position corresponding to when the value of the melt function m is 0.9. According to Equation (3), in order to complete the establishment of the two-dimensional model of the sound speed field of the ocean front, we need to obtain the sound speed profiles of the water masses on both sides of the frontal zone and the melt function related to the distance between the two sound speed profiles. Since the KEF is relatively independent and the water masses on the north and south of the frontal zone are relatively stable, the sound speed profiles at (38° N, 144° E) and (36° N, 144° E) were selected to represent the sound speed profiles of the water masses on the north and south sides of the frontal zone (corresponding to R = 200 km mentioned above). Since different melt parameters a in the melt function correspond to different positions and widths of the frontal zone, the positions and widths of the frontal zone should be obtained first (clustering the surface sound speed of the study area with the fuzzy C-means (FCM) algorithm). Then the values of the corresponding melt parameters can be determined according to Figure 3 to quickly reconstruct the two-dimensional model of the sound speed field of the KEF. According to Equation (3), in order to complete the establishment of the two-dimensional model of the sound speed field of the ocean front, we need to obtain the sound speed profiles of the water masses on both sides of the frontal zone and the melt function related to the distance between the two sound speed profiles. Since the KEF is relatively independent and the water masses on the north and south of the frontal zone are relatively stable, the sound speed profiles at (38 • N, 144 • E) and (36 • N, 144 • E) were selected to represent the sound speed profiles of the water masses on the north and south sides of the frontal zone (corresponding to R = 200 km mentioned above). Since different melt parameters a in the melt function correspond to different positions and widths of the frontal zone, the positions and widths of the frontal zone should be obtained first (clustering the surface sound speed of the study area with the fuzzy C-means (FCM) algorithm). Then the values of the corresponding melt parameters can be determined according to Figure 3 to quickly reconstruct the two-dimensional model of the sound speed field of the KEF.
In order to reduce the dependence of ocean front modeling on the data of the whole sea area and the whole sea depth, only the surface sound speed of the study area was clustered with the fuzzy C-means (FCM) algorithm. FCM algorithm was proposed by Bezdek J.C. et al. in the 1980s [30,31]. FCM obtains the membership degree of each sample point to all cluster centers by optimizing the objective function, thereby determining the category of sample points to automatically classify the sample data [32]. In this paper, we use the surface sound speed as data points and cluster the surface sound speed through the FCM algorithm, which has achieved the effect of detecting the frontal zone. There are 26 * 26 = 676 data points in the study area surface, FCM divides 676 sound speed data x i (i = 1, 2,. . . , 676) into c groups, and identifies the cluster center of each group by minimizing the objective function of the dissimilarity index. The expression of the objective function is as follows: where Q is the membership matrix, V is the cluster center, r ik is the membership degree of the k sample point for i category, n is the number of the sample point (in this study, n = 676), c is the number of categories (considering the existence of ocean front, the study area is divided into three water masses, namely, high temperature and salinity water mass located in the south of the frontal zone, low temperature and salinity water mass located in the north of frontal zone south of the frontal zone, and transitional water mass of frontal zone between them, so we set c to 3), q is the weighted index (generally q = 2), and u k − V i is the distance between the k sample and the i category. The steps of FCM used in this study are as follows: (1) Initializing the membership matrix Q with a number between 0 and 1. This step can be completed through translation-standard deviation transformation and translationrange transformation.
The expression of translation-standard deviation transformation is as follows: The expression of translation-range transformation is as follows: After transformation, all x ij ∈ [0, 1] of the sound speed matrix can be obtained, so we complete the initialization of the membership matrix Q.
(2) Calculating the cluster center V(i) of the corresponding matrix Q(i).
(3) Updating the membership matrix Q(i) and comparing Then, the iteration can be stopped, and the surface sound speed data points are clustered.

Introduction of Technical Route
The whole research process is summarized, and the technical route is given as follows ( Figure 4). In order to reduce the dependence of ocean front modeling on the data of the whole sea area and the whole sea depth, we only used the sound speed data of sea surface and the SSP of (38 • N, 144 • E) and (36 • N, 144 • E) (calculated by Equation (1)). The range and position of the frontal zone can be obtained by using FCM for surface sound speed. Then, the appropriate melt parameter a and melt function m can be selected correspondingly according to the results of FCM, and the sound speed profiles of water masses on both sides can be used as the input of the model to obtain the two-dimensional model of the sound speed field of the ocean front. Compared with the previous research, our research requires less data in this model to quickly obtain the distribution of the sound speed of ocean front. Next, the BELLHOP ray model was used (BELLHOP is an underwater acoustic toolbox for the calculation of the sound field, and we can adjust the parameters through the environment file to get the corresponding sound rays information, multipath information, and so on) to calculate the TL and compare the results of the model with the actual results (calculated by reanalysis data) of the KEF to verify the effectiveness of our model. Finally, we focused on the change of the depth of the acoustic convergence area in the ocean front environment and realized the forecast of the depth of the convergence area by combining the depth of the convergence area in the model and the melt function.
ters through the environment file to get the corresponding sound rays information, multipath information, and so on) to calculate the TL and compare the results of the model with the actual results (calculated by reanalysis data) of the KEF to verify the effectiveness of our model. Finally, we focused on the change of the depth of the acoustic convergence area in the ocean front environment and realized the forecast of the depth of the convergence area by combining the depth of the convergence area in the model and the melt function. Figure 4. The technical route of the research process. We used FCM to cluster the surface sound speed and detect the frontal zone of KEF; At the same time, sound speed profiles were used to establish the model of the sound speed field through the melt function. Next, the model was compared with the actual ocean front calculated by reanalysis data, which proved its validity. Finally, we proposed the melt function in the model to forecast the depth change of the acoustic convergence area.

Results and Discussion
In this section, we cluster the surface sound speed of the sea with the fuzzy C-means (FCM) algorithm to obtain the position and range of the frontal zone. At the same time, the theoretical model of the ocean front is improved by replacing the temperature profile with the SSP, and the model of the sound speed field of the Kuroshio Extension front Figure 4. The technical route of the research process. We used FCM to cluster the surface sound speed and detect the frontal zone of KEF; At the same time, sound speed profiles were used to establish the model of the sound speed field through the melt function. Next, the model was compared with the actual ocean front calculated by reanalysis data, which proved its validity. Finally, we proposed the melt function in the model to forecast the depth change of the acoustic convergence area.

Results and Discussion
In this section, we cluster the surface sound speed of the sea with the fuzzy C-means (FCM) algorithm to obtain the position and range of the frontal zone. At the same time, the theoretical model of the ocean front is improved by replacing the temperature profile with the SSP, and the model of the sound speed field of the Kuroshio Extension front (KEF) is established quickly and effectively (Section 3.1). We compare the model with the actual sound speed distribution of the KEF calculated by reanalysis data and verify the validity of the model from the perspective of transmission loss (TL) (Section 3.2). Finally, it is found that the melt function in this model can forecast the depth of the acoustic convergence area in the KEF (Section 3.3).

Reconstruction Results of Ocean Front Sound Speed Field Model
We used FCM to cluster the surface sound speed in the study area, and the results are shown in the figure below ( Figure 5). We find that the position, width, and range of the KEF are different at different times.
In order to verify the validity of the results of FCM, we calculated the horizontal temperature gradient of the sea surface in the study area to obtain the detection results of the frontal zone under the traditional gradient method (the critical value is 0.05 • C/km) ( Figure 6). We found that the detection results of the frontal zone were similar to those of the FCM, especially in winter (Figure 6a,b). Compared with those of the gradient method, the detection results of FCM led to a more continuous frontal zone in the range, which is helpful for us to obtain the starting and ending positions of the frontal zone on the relevant sections and realize effective modeling. However, the continuity of the judgment results of the temperature gradient method is very poor, and this discontinuity may become worse due to different judgment standards.
vergence area in the KEF (Section 3.3).

Reconstruction Results of Ocean Front Sound Speed Field Model
We used FCM to cluster the surface sound speed in the study area, and the results are shown in the figure below ( Figure 5). We find that the position, width, and range of the KEF are different at different times. In order to verify the validity of the results of FCM, we calculated the horizontal temperature gradient of the sea surface in the study area to obtain the detection results of the frontal zone under the traditional gradient method (the critical value is 0.05 °C /km) ( Figure 6). We found that the detection results of the frontal zone were similar to those of the FCM, especially in winter (Figure 6a,b). Compared with those of the gradient method, the detection results of FCM led to a more continuous frontal zone in the range, which is helpful for us to obtain the starting and ending positions of the frontal zone on the relevant  Table 1). The melt parameters a of the model of the sound speed field of the ocean front in four periods were determined as 0.2, 0.3, 0, and −0.1, respectively.
Taking the experiment on 1 January as an example, we obtained the sound speed profiles of water masses on both sides of the frontal zone ((38 • N, 144 • E) and (36 • N, 144 • E)) and the sound speed profiles on the edge of the frontal zone (Figure 7). We found obvious differences between the sound speed profiles of water masses on both sides of the frontal zone. The depth of the SOFAR axis (the depth of the minimum sound speed) of the SSP on the north side was deeper, and the surface sound speed difference between the water masses on the south side and the north side was more than 18 m/s. Between the front edge, 1 and 2 are the sound speed profiles in the frontal zone.
Appl. Sci. 2021, 11, x FOR PEER REVIEW 9 of 23 sections and realize effective modeling. However, the continuity of the judgment results of the temperature gradient method is very poor, and this discontinuity may become worse due to different judgment standards.   Taking the experiment on 1 January as an example, we obtained the sound speed profiles of water masses on both sides of the frontal zone ((38° N, 144° E) and (36° N, 144° E)) and the sound speed profiles on the edge of the frontal zone (Figure 7). We found obvious differences between the sound speed profiles of water masses on both sides of the frontal zone. The depth of the SOFAR axis (the depth of the minimum sound speed) of the SSP on the north side was deeper, and the surface sound speed difference between the water masses on the south side and the north side was more than 18 m/s. Between the front edge, 1 and 2 are the sound speed profiles in the frontal zone. According to the melt parameters determined at different times and the sound speed profiles at points (38° N, 144° E) and (36° N, 144° E), the model of the sound speed field of the KEF was established quickly (Figure 8). We found that the distribution of sound speed obtained by the overall model was similar to that of the actual KEF ( Figure 2) in terms of the location of the front and the depth of the influence (especially in winter). However,

Comparison Results of TL Calculations
The BELLHOP ray model was used to calculate the TL of our model and the actual KEF environment. Compared with the normal mode model or parabolic equation model, the ray model can clearly describe the change of sound energy in the process of propagation, which is more suitable for the ocean front environment in this paper. The parameter settings in the sound field calculation process are as follows: the sound source frequencies were set to 50 Hz and 400 Hz, the glancing angle was −20°-20°, the sound source depth (SD) was 100 m, 300 m, and 500 m, and the receiving depth remained the same as the SD. The seabed sound speed is 1600 m/s, the density was 1.8 g/cm 3 , and the sound absorption coefficient of the seabed was 0.8 dB/λ.
We respectively calculated the TL under the model environment and the actual ocean front environment at different times (Figures 9-12). We found that when the SD was 100 m, the difference between the TL calculated by the model and the TL calculated by the actual ocean front was very small at low and high frequencies. The result of the model at 100 m can well replace the actual TL; at 300 m, only on February 1st (Figure 10b,e), the calculation results of low frequency and high frequency showed a small difference after passing the frontal zone, but in terms of the overall distribution trend, the results of the model were similar to the calculation results under the actual ocean front environment; at a deeper emission depth (SD = 500 m), since the SD is close to the SOFAR axis, there was

Comparison Results of TL Calculations
The BELLHOP ray model was used to calculate the TL of our model and the actual KEF environment. Compared with the normal mode model or parabolic equation model, the ray model can clearly describe the change of sound energy in the process of propagation, which is more suitable for the ocean front environment in this paper. The parameter settings in the sound field calculation process are as follows: the sound source frequencies were set to 50 Hz and 400 Hz, the glancing angle was −20 • -20 • , the sound source depth (SD) was 100 m, 300 m, and 500 m, and the receiving depth remained the same as the SD. The seabed sound speed is 1600 m/s, the density was 1.8 g/cm 3 , and the sound absorption coefficient of the seabed was 0.8 dB/λ.
We respectively calculated the TL under the model environment and the actual ocean front environment at different times (Figures 9-12). We found that when the SD was 100 m, the difference between the TL calculated by the model and the TL calculated by the actual ocean front was very small at low and high frequencies. The result of the model at 100 m can well replace the actual TL; at 300 m, only on February 1st (Figure 10b,e), the calculation results of low frequency and high frequency showed a small difference after passing the frontal zone, but in terms of the overall distribution trend, the results of the model were similar to the calculation results under the actual ocean front environment; at a deeper emission depth (SD = 500 m), since the SD is close to the SOFAR axis, there was no obvious convergence area or shadow area phenomenon in the distribution of TL, but in terms of the overall trend, the results of the model and the actual results also showed many similarities.
Appl. Sci. 2021, 11, x FOR PEER REVIEW 12 of 23 no obvious convergence area or shadow area phenomenon in the distribution of TL, but in terms of the overall trend, the results of the model and the actual results also showed many similarities.     In order to quantify the simulation results of the model, we smoothed the above results through the five-point smoothing method and calculated the root means square error (RMSE) and the percentage with the difference exceeding 10 dB compared with the actual In order to quantify the simulation results of the model, we smoothed the above results through the five-point smoothing method and calculated the root means square error (RMSE) and the percentage with the difference exceeding 10 dB compared with the actual TL (Table 2). It was found that before smoothing, the RMSE of the model was mostly controlled below 6 dB compared with the actual results. Compared with that in previous studies, the effect was further improved. Before smoothing, the percentage with the difference exceeding 10 dB was maintained near 10%, which indicated that the results of our model were less fluctuating than the actual results. In terms of RMSE, all smoothed data were reduced to less than 6 dB, and the percentage with the difference exceeding 10 dB was also further reduced. Except for individual data near 10%, most of the data had been reduced to less than 8%, which verified the effectiveness of the model. As to TL, the calculation results of the model were in good agreement with the actual situation. We calculated the RMSE under different sound source frequencies (Table 3) and different sound source depths (Table 4) to compare the differences between winter (the dates were 0101 and 0201) and summer (the dates were 0701 and 0801). We found that when frequencies were different, the RMSE of TL was close to the same in different seasons, and the seasonal difference of the RMSE caused by the change of frequency was not obvious. When sound source depths were different, the RMSE of TL was different between winter and summer. Specifically, when SD = 100 m, the RMSE in winter was less than that in summer (the difference of unsmoothed results was about 1 dB); When SD = 300 m, the unsmoothed results of the two seasons were equivalent; when SD = 500 m, the RMSE in winter was less than that in summer (the difference of unsmoothed results was still about 1 dB). Analyzing the reasons for the difference of the RMSE of TL under different sound source depths in the two seasons, we believed that in terms of the frontal zone position and influence depth, the sound speed distribution obtained by the model in winter was more similar to the actual KEF distribution than that in summer ( Figure 8). Therefore, when the sound source and reception depth were shallow, the RMSE of TL obtained in winter was smaller; at deeper depths, the phenomenon of convergence area is no longer obvious compared with that of shallow sound source depth (Figure 9c,f and Figure 12c,f), the RMSE may be greater in winter than in summer (this level was only less than 1 dB, which was acceptable, but we can not deny that the effectiveness of the model in winter was better). Table 3. Comparison between the TL calculated through the model and the actual situation at different seasons and different frequencies (the sound source depths were 100 m, 300 m, and 500 m, respectively).

Effectiveness of Forecasting the Changes of Depth in Convergence Areas
The convergence area is the area where the sound energy is concentrated near the sea surface during the sound propagation process [35]. The study of the convergence area is of practical significance for underwater acoustic detection [36,37]. Due to the existence of the ocean front, the depth change of the convergence area has become a concern. Taking the experiment on January 1st as an example, we set the sound source frequency to 50 Hz, and calculated the whole sea depth variation of the TL in the distance-independent environment (the environment without ocean fronts), model environment, and the actual ocean front environment at the SD of 100 m, 300 m, and 500 m ( Figure 13). It was found that due to the existence of the ocean front, there was a significant downward shift in the distribution of TL after its passing through the frontal zone in the convergence area, which was a very obvious phenomenon is in both the model and the actual ocean front environments. The change in the depth of the convergence area is related to not only the SD but also the location of the frontal zone and the range of its influence. Since the effectiveness of the model from the perspective of TL is verified in the previous section, in this section the effectiveness of melt function in the rapid forecast of the convergence area depth in the ocean front environment is verified by combining the melt function with the depth of the convergence area calculated by the model and the actual ocean front. surface of the outgoing layer respectively. According to Snell's law, when the sound speed changes with depth, the sound rays always bend to the place where the sound speed is lower. The convergence area is a manifestation of the concentrated area of sound ray turning, so the depth of the convergence area is not only related to the SD but also has a corresponding relationship with the distribution of sound speed.
Since only the sound speed profiles of the water masses on both sides of the frontal zone were used in our model as input values, we only obtained the data of these two The law of refraction followed by ray acoustics is Snell's law: where α 0 is the angle between the incident direction and the horizontal direction of the sound rays, called grazing angle, c 0 is the sound speed on the medium surface of the incident layer, α i and c i is the grazing angle and sound speed of the sound rays on the surface of the outgoing layer respectively. According to Snell's law, when the sound speed changes with depth, the sound rays always bend to the place where the sound speed is lower. The convergence area is a manifestation of the concentrated area of sound ray turning, so the depth of the convergence area is not only related to the SD but also has a corresponding relationship with the distribution of sound speed. Since only the sound speed profiles of the water masses on both sides of the frontal zone were used in our model as input values, we only obtained the data of these two sound speed profiles. Taking the experiment on January 1st as an example, the SOFAR axis depth of these two sound speed profiles was 700 m and 1250 m, respectively (Figure 7). When the SD is fixed, according to Snell's law, the depth of the convergence area formed by the sound rays emitted at a small glancing angle near the sound source should be basically the same as the SD (or slightly shallower than the SD) under a distance-independent environment. As shown in Figure 13, when the SD is 100 m, 300 m, and 500 m, the depth of the convergence area always stays near 0 m (Figure 13a), 200 m (Figure 13b), and 400 m (Figure 13c) under a distance-independent environment. In our model and actual ocean front environment, when the SD is located at 100 m, 300 m, and 500 m, the depth of the convergence area near the sound source is 0 m (Figure 13d,g), 200 m (Figure 13e,h), and 400 m (Figure 13f,i), respectively since the ocean environment not far away from the sound source is also a distance-independent environment. Based on the above analysis, we obtained a forecast curve of the depth change of the convergence area during the transmission process in the ocean front environment: where D 1 (z) is the depth of the convergence area near the sound source, D 2 (z) is the depth of the convergence area at the farthest distance (200 km). m(r; a) is the melt function, and the expression is consistent with Equation (4). In the ocean front environment, it is believed that the depth change of the convergence area is related to the change of the sound speed. Therefore, the depth of the SOFAR axis was extracted as the feedback depth and to estimate D 2 (z): where SOFAR 1 (z) and SOFAR 2 (z) are the depths of the SOFAR axis of two sound speed profiles. Take the experiment on January 1st as an example, the SOFAR 1 (z) and SOFAR 2 (z) are known to be 700 m and 1250 m. When the SD = 100 m, the depth of the convergence area near the sound source was about 0 m (D 1 (z) =0 m), so D 2 (z) should be 550 m. We completed the forecast of the depth of the convergence area at the SD of 100 m on 1 January and compared this result with the model calculated one (we selected two convergence areas passing through the frontal zone) and the actual results of the ocean front (corresponding to the two convergence areas of the model), to obtain the following results ( Figure 14). We found that the depth distribution of the convergence area in the model and the actual condition were in good agreement with the forecasted curve with only a little difference, which proves that the melt function was applied to effectively achieve the depth forecast of the convergence area in the ocean front environment.  By respectively analyzing the results of the depth of the convergence area during the four periods (Figures 14 and 15), it was concluded that since the two sound speed profiles were different in different periods when the SD was fixed, the depth of the convergence area is related to that of the SOFAR axis, the difference between them, as well as the SD. The forecasting results using the melt function based on the melt parameter a (which is related to the width and location of the frontal zone) were in good agreement with the model and the actual results, which proves the effectiveness of this method. Figure 14. The forecasting curve of the depth of the convergence area (the red dotted line) and the calculating results of the model (blue dots) and actual ocean front (black dots) on 1 January (SD = 100 m). We found that the depth distribution of the convergence area in the model and the actual condition were in good agreement with the forecasted curve with only a little difference. the melt parameter a (which is related to the width and location of the frontal zone) were in good agreement with the model and the actual results, which proves the effectiveness of this method. It proves that the melt function was applied to effectively achieve the depth forecast of the convergence area in the ocean front environment.
By calculating the forecasting results of other sound source depths together (Table 5), it was found that when the SD was shallow, the forecasting results were in good agreement with the actual results, but the degree of agreement was relatively poor at a deeper launch depth (SD = 500 m) because SD's closeness to the SOFAR axis makes most of the August. The forecasting results using the melt function based on the melt parameter a (which is related to the width and location of the frontal zone) were in good agreement with the model and the actual results, which proves the effectiveness of this method. It proves that the melt function was applied to effectively achieve the depth forecast of the convergence area in the ocean front environment.
By calculating the forecasting results of other sound source depths together (Table 5), it was found that when the SD was shallow, the forecasting results were in good agreement with the actual results, but the degree of agreement was relatively poor at a deeper launch depth (SD = 500 m) because SD's closeness to the SOFAR axis makes most of the sound energy converge near the SOFAR axis and propagates forward. It is impossible to form an obvious convergence area, which can not be detected during some periods. After calculation, the RMSE between the forecasting results and model calculation results in the first detection convergence area was 74.8 m, and the second detection convergence area was 45.0 m; The RMSE between the forecasting results and the actual calculation results in the first detection convergence area was 82.3 m, and the second detection convergence area was 43.3 m. It shows that the RMSE of the second detection convergence area was smaller and the effectiveness of forecasting was better.

Summary and Conclusions
In this paper, we first used the fuzzy C-means clustering method to perform fuzzy clustering on the surface sound speed of the Kuroshio extension study area and obtained the range and location of the frontal zone. The result of cluster analysis can be consistent with the result of surface temperature gradient judgment. Next, according to the position and range of the frontal zone at different times, we selected the corresponding melt parameters to construct a two-dimensional parametric model of ocean fronts based on the sound speed profile. The results showed that it had a good similarity with the actual two-dimensional sound speed distribution of the ocean front. In order to verify the effectiveness of the ocean front model in acoustic forecasting, we selected a characteristic section and used the BELLHOP ray model to obtain the comparison between the model and the actual transmission loss by setting different sound source depths and emission frequencies. It was found that the root means square error between the calculated results of the model and the actual results was mostly controlled below 6 dB. As far as the percentage with the difference exceeding 10 dB is concerned, this level was basically maintained at around 10%. In terms of the transmission loss, the calculation results of the model were in good agreement with the actual situation. Finally, we proposed a simple combination of the melt function and the depth of the convergence area to verify the effectiveness of applying the melt function to the rapid forecasting of the depth of the convergence area in the ocean front environment. After calculation, the root mean square error between the forecasting result and the actual calculation result in the first detection convergence area was 82.3 m, and the second detection convergence area was 43.3 m, the forecasting effect was better. It is hoped that this research may be to some degree of referential significance to acoustic detection in the ocean front environment.
As one of the most common mesoscale phenomena in the ocean, the impact of ocean fronts on underwater acoustic detection has become a hot spot for hydroacoustics. Because experimental data with high spatial resolution in the ocean front area can not be quickly and effectively obtained by oceanographic surveys, it is feasible to use theoretical models combined with typical data to reconstruct ocean fronts and analyze their impact on acoustic detection. Therefore, according to our method, future oceanographic surveys only require surface sound speed data (of course, this part of the data can also be obtained through the satellite remote sensing data of high spatial resolution) and sound speed profiles representing the typical properties of water masses (this part of the data can be obtained through the actual ocean surveys). Then we can construct the two-dimensional sound speed field distribution model in the ocean front environment. This method was proven effective in our research. This method can effectively predict the depth of the convergence area and can also provide us with relevant help in the acoustic detection of the relevant sea areas in the ocean front. For example, the decline of the depth of the convergence area is unfavorable to the acoustic detection under certain conditions, and this method can capture this disadvantage in advance.
It should be pointed out that although we have proved the applicability of this method to KEF in this paper, we have not proved whether this method is applicable in other ocean front environments. This is also one of the key points of future research. In addition, we also need to further study the interference of other ocean phenomena such as internal waves on the two-dimensional sound speed field modeling of ocean front.