Simulative Evaluation of the Underwater Geodetic Network Conﬁguration on Kinematic Positioning Performance

: The construction of underwater geodetic networks (UGN) is crucial in marine geodesy. To provide high-precision kinematic positioning for underwater submersibles, an underwater acoustic geodetic network conﬁguration of three seaﬂoor base stations, one subsurface buoy, and one sea surface buoy is proposed. The simulation results show that, for a 3 km-deep sea, based on the proposed UGN, the submersible positioning range and positioning accuracy are primarily affected by the size of the seaﬂoor base station array, while the height of the subsurface buoy has a greater impact on the submersible positioning accuracy than the positioning range. Considering current acoustic ranging technology, the kinematic positioning performance of the UGN is optimal when the seaﬂoor base stations are 9~13 km apart and the subsurface buoy is less than 2.5 km above the seaﬂoor, which can achieve a submersible positioning accuracy of less than 30 m within an underwater space of 25 km × 25 km × 3 km. The proposed cost-effective UGN conﬁguration can provide high-precision submersible kinematic positioning performance for seaﬂoor surveying and ocean precision engineering. The impact of the underwater environment on the acoustic transmission characteristics should be further investigated.


Introduction
Underwater geodetic networks (UGNs) are pivotal for ocean resources exploration, environmental monitoring, and underwater navigation [1]. Due to the attenuation of the electromagnetic field in water, the global navigation satellite system (GNSS) is unsuitable for underwater scenarios [2,3]. Based on the combination of GNSS and acoustic ranging technology (GNSS-A), acoustic underwater geodetic networks can accurately obtain seafloor displacement information [4][5][6][7][8][9]. Moreover, geodetic networks can provide location information for underwater targets, including autonomous underwater vehicles (AUV), remotely operated underwater vehicles (ROV), underwater gliders, underwater buoys, etc., which can be deployed in seafloor surveying and ocean precision engineering [2,3,10,11]. Therefore, the construction of UGNs is currently an important research topic in marine geodesy.
Acoustic positioning systems, such as Long Baseline (LBL) and Ultra-Short Baseline (USBL), are commonly used for underwater positioning [12]. The seafloor beacon array of LBL can be considered a simple underwater geodetic network. With the appropriate sonar and acoustic ranging technology, the positioning range of UGN can reach 10 km with the

Underwater Environmental Factors
The network configuration is the primary focus of our study. Generally, the seafloor topography, ocean currents, tides, internal waves, sea surface wind, sound speed variation, and other environmental factors all have an impact on the network kinematic positioning performance. The acoustic signal propagation bending always causes sound rays to sweep the seafloor, affecting acoustic signal propagation [15,33]. In the construction of an UGN, its height should be designed considering the seafloor topography and seabed sediment. Ocean currents and internal waves can affect the stability of the subsurface buoy [29], and the height of the subsurface buoy must be optimized. Similarly, waves and tides affect the stability of the sea surface buoy, and its position error also needs to be In this study, the basic network configuration, which consists of three seafloor base stations, one subsurface buoy, and one surface buoy, is expressed as the "3 seafloor base stations + 1 subsurface buoy + 1 sea surface buoy" ("3+1+1") configuration. The seafloor base stations are distributed evenly on a circle, and the distance between the seafloor base stations is set as D. The subsurface buoy is located above the horizontal center of the seafloor base station array with a vertical height of H (0 < H ≤ 3000 m), while the Remote Sens. 2022, 14, 1939 4 of 18 sea surface buoy is about 3000 m above the seafloor array center. Note that the basic configuration can be expanded or reduced to other configurations, such as the "4 seafloor base stations + 1 subsurface buoy + 1 sea surface buoy" configuration, "3 seafloor base stations + 0 subsurface buoy + 1 sea surface buoy" configuration, and so on. The network dimension, its positioning range, and kinematic positioning performance at different undersea points are closely related to the seafloor station separation and the subsurface buoy height. Therefore, the performance of the network configuration in the underwater space should be evaluated.

Underwater Environmental Factors
The network configuration is the primary focus of our study. Generally, the seafloor topography, ocean currents, tides, internal waves, sea surface wind, sound speed variation, and other environmental factors all have an impact on the network kinematic positioning performance. The acoustic signal propagation bending always causes sound rays to sweep the seafloor, affecting acoustic signal propagation [15,33]. In the construction of an UGN, its height should be designed considering the seafloor topography and seabed sediment. Ocean currents and internal waves can affect the stability of the subsurface buoy [29], and the height of the subsurface buoy must be optimized. Similarly, waves and tides affect the stability of the sea surface buoy, and its position error also needs to be considered [21,23]. In addition, the coordinate error and clock drift of network base stations also influence the kinematic positioning performance. Similar to GNSS, the impact of underwater environmental factors can be classified as User Equivalent Range Errors (UEREs), while the station coordinate and clock errors are considered as Signal-In-Underwater-Space Range Errors (SISREs) [34]. The impact of underwater environmental factors is evaluated in Section 3.4.

Kinematic Positioning Performance Assessment Methods
To evaluate the kinematic positioning performance of the UGN, an underwater space with a depth of 3000 m and a horizontal range of 25 km × 25 km was simulated in grids, and the grid spatial resolution was 500 m on the horizontal and 200 m on the vertical. For ocean near-bottom surveying and precision engineering applications, the user can distribute at any vertex of the grids 100 m above the seafloor, and the total number of grids was 39,015.
The underwater user's single-epoch kinematic positioning can be realized by receiving acoustic signals from UGN. The acoustic sonar properties and the underwater environment influence the kinematic positioning performance. Considering the current acoustic ranging technology [13,32], the omni-direction acoustic sonar's operating frequency is 8~16 kHz, and the maximum positioning range of one network station can reach 10 km with the same precision at a depth of 3000 m. Additionally, the time delay measurement error was set to 0.1 ms, which was equivalent to the ranging error of approximately 0.15 m. To simplify the impact of the underwater environment, only the constant sound speed error of 0.1 m/s was considered in the simulation [35].
In terms of the positioning range and precision, the Available Station Number (ASN), Effective Positioning Grid (EPG), and Position Dilution of Precision (PDOP) indicators are considered. For the user's random position, the Euclidean distance between the user and the station (seafloor base station, subsurface buoy, or sea surface buoy) is set to L. If 0 < L ≤ 10 km, then the user is within the effective positioning range of a station, and we consider this station as an available station for this user grid point; otherwise, it is unavailable. For different user grid points, the available station number (ASN) is variable.
Similar to GNSS, if the sonar signal delay is ignored and the clock of the synchronizationtype sonar is precisely synchronized, the user single-epoch kinematic positioning can be achieved by receiving acoustic signals from at least 3 stations (seafloor base stations and/or buoys). We defined the Effective Positioning Grid (EPG) as the user grid point with at least Remote Sens. 2022, 14, 1939 5 of 18 3 available station numbers (ASNs). Therefore, the EPG index is intended to evaluate the effective positioning range of the network configuration in the underwater space.

• Position Dilution of Precision (PDOP)
To realize user single-epoch kinematic positioning, the least-squares adjustment can be applied, and the basic positioning model is where l j i is the acoustic measurement in distance between user i and station j, L j i (x i ) is the Euclidean distance, v j i is the distance error, and x i and c j are the positions of user i and station j, respectively.
The linearization of the Equation (1) is where δx = x − x 0 , and x 0 is the initial value of x. Accumulating the acoustic measurements from all of the available stations, we obtain the error equations where δx is the parameter to be estimated, B is the design matrix, and n is the available station number. For user single-epoch kinematic positioning, we can evaluate the positioning accuracy σ X sound speed error (4) where N is the normal matrix, P is the weight matrix, which considers the ranging error and the constant sound speed error, and σ 0 2 is the variance of unit weight [26,36]. The a posteriori variance of the unit weight is assumed to equal to the a priori variance of the unit weight in the adjustment model.
The positioning accuracy σ X is related to the network geometric structure formed by the user and network stations, which can be evaluated by the position dilution of precision PDOP [26] We analyzed the PDOP in the underwater space using Equations (6)- (8), where m is the number of all grid points. The maximum PDOP determines the lower bound of the positioning accuracy.

Evaluation of Network Kinematic Positioning Performance
To evaluate the kinematic positioning performance of the basic network configuration "3+1+1", various distances between the seafloor base stations and heights of the subsurface buoy were simulated. As the user can obtain single-epoch kinematic positioning by receiving acoustic signals from at least three stations, we analyzed the number of effective positioning grids (EPGs) and the position dilution of precision (PDOP) with at least three available station numbers (ASN). Figure 2a shows the percentage of EPGs in the total grid number of 39,015. For the same subsurface buoy height, the EPG varied significantly with the seafloor base station distance. When the distance was less than 17 km, the EPG decreased slowly as the distance increased. However, when the distance was greater than 17 km, the circumscribed circle radius of the seafloor base stations array was greater than 10 km, and the EPG decreased rapidly. For the same seafloor station distance, the EPG changed gently with the subsurface buoy height. The EPG values were relatively larger when the subsurface buoy height was above 1.5 km. Overall, the EPG reached its maximal value of 18,486 when the seafloor base station distance was 1 km and the subsurface buoy height was 1.5 km.

Impact of Network Dimension on Positioning Performance
For the "3+1+1" configuration, the dimensions relate to the distance between the se floor base stations and the height of the subsurface buoy. The EPGs of different dime  Figure 2b shows the PDOP value variation with the seafloor base station distance and the subsurface buoy height. According to Equation (5), due to the inversion of N, outliers of PDOP value may occur. We detected and removed the PDOP outlier when its value exceeded ±1.5 × IQR (interquartile range), which amounted to approximately 7% of all of the PODPs [37]. Figure 2b indicates that, when the subsurface buoy height is lower than 1.5 km, the PDOP value changes slowly with the base station distance. However, when the height is higher than 1.5 km and the distance changes from 7 km to 9 km, the PDOP values change notably. Overall, the results show that, when the subsurface buoy height is less than 1.5 km, the PDOP mainly changes with the seafloor base station distance. Otherwise, the PDOP fluctuates with the subsurface buoy height and the station distance.

Impact of Network Dimension on Positioning Performance
For the "3+1+1" configuration, the dimensions relate to the distance between the seafloor base stations and the height of the subsurface buoy. The EPGs of different dimensions can be scattered or concentrated in the underwater space. We need to analyze the impact of the seafloor base station distance and subsurface buoy height on the user positioning range, available station number, and positioning accuracy. Figure 3a shows the EPG values with various available stations in the "3+1+1" network configuration. We defined EPGs with 3, 4, and 5 available station numbers as EPG ASN=3 , EPG ASN=4 , and EPG ASN=5 , respectively. For the convenience of comparison, we also expressed EPGs with 0, 1, and 2 available station numbers as EPG ASN=0 , EPG ASN=1 , and EPG ASN=2 , respectively.
From Figure 3a, in general, the seafloor base station distance has a much greater impact on the EPG with different available station numbers than the subsurface buoy height. When the height was constant, EPG ASN=0 decreased with the increase in distance; when the distance was 13~19 km, EPG ASN=0 tended to be stable. However, EPG ASN=0 gradually increased when the distance was greater than 19 km. Due to the limited acoustic positioning range (10 km), there were always grid points with unavailable base stations. EPG ASN=1 increased with the distance at first; when the distance was close to 11 km, EPG ASN=1 was stable; when the distance was greater than 19 km, EPG ASN=1 showed a decreasing trend. EPG ASN=2 was significantly lower than the others, and EPG ASN=2 began to increase when the distance was greater than 19 km. EPG ASN=3 increased with the distance at first; when the distance approached 19 km, EPG 3 reached the maximum value and then decreased. EPG ASN=4 was similar to EPG ASN=3 , while the maximum value of EPG ASN=4 occurred when the distance was close to 11 km, and EPG ASN=4 disappeared with distances greater than 19 km. EPG ASN=5 showed a downward trend with the distance. When the distance was greater than 17 km, the distance between the seafloor base stations exceeded the maximum acoustic range (10 km), and user positioning with five stations was no longer possible. Figure 3b shows the average, standard deviation, and maximum values of PDOP. PDOP was inversely proportional to the geometric structure formed by the user and stations, and the user positioning accuracy based on the network configuration was affected by the seafloor base station distance and the subsurface buoy height. After the removal of outliers with the IQR criterion, most PDOP values were less than 30, except for when the height was above 2.5 km. We found that PDOP outliers appeared when users were located near the seafloor base stations. Moreover, when the height was lower than 1.5 km, the PDOP value changed slowly with the seafloor base station distance. When the distance increased from 7 km to 9 km, the variation of PDOP was obvious.
Since the subsurface buoy height had little impact on EPG, we investigated the spatial distribution of the EPG with the fixed subsurface buoy height and various seafloor base station distances (Figure 3c,d). The height was fixed at 2.9 km, and the distances were 7 km and 9 km, respectively. In both scenarios, grids with five available stations covered the station array, and grids with three and four available stations were outside the station array. It is obvious that they were all within a circle in the plane, and the radius was Remote Sens. 2022, 14,1939 8 of 18 approximately 10 km. Because the subsurface buoy and surface buoy were above the horizontal center of the seafloor base station array, grids with at least three available base stations were widely distributed. While grids with two available stations were only sporadically distributed, grids with one or unavailable base stations were distributed at the edge of the underwater space, where the acoustic ranging signal of the subsurface buoy and the surface buoy could not reach. When the station distance increased from 7 km to 9 km, EPG ASN=0 and EPG ASN=2 decreased by 27.63% and 3.82%; EPG ASN=1 , EPG ASN=3 , and EPG ASN=4 increased by 19.40%, 45.83%, and 10.85%, respectively; while EPG ASN=5 decreased by 86.63%, which was primarily linked to the obvious variation in the PDOP. From Figure 3a, in general, the seafloor base station distance has a much greater impact on the EPG with different available station numbers than the subsurface buoy height. When the height was constant,  Based on the above analysis, for the "3+1+1" network configuration, the seafloor base station distance of 9~13 km and the subsurface buoy height of lower than 2.5 km were optimal in terms of positioning range and positioning accuracy. Considering Equation (4), the submersible positioning accuracy was less than 30 m within an underwater space of 25 km × 25 km × 3 km.
We analyzed EPGs with at least three available stations (EPG ASN≥3 ). Figure 4a shows that the EPG ASN≥3 of the "3+0+0" network configuration gradually decreased as the distance increased. When the distance was greater than 17 km, the available station number was not enough for positioning. For the "3+0+1" network configuration with the sea surface buoy, the EPG ≥3 values were much larger than those of the "3+0+0" network configuration. However, when the distance was greater than 19 km, the available station number was also not enough for positioning. Figure 4b shows that, for the "3+1+0" network configuration, by replacing the sea surface buoy with the subsurface buoy, the EPG ≥3 values had a similar situation. Compared with Figure 2a, for the "3+1+1" network configuration, it can be found that the EPG ASN≥3 values increased significantly when the sea surface buoy and subsurface buoy were configured simultaneously. Moreover, when the distance was greater than 19 km, the proportion of EPG ASN≥3 was still greater than >21%.   PDOP outliers were removed if the PDOP value exceeded ±1.5 × IQR. For the "3+0+0" network configuration, the PDOP statistics first decreased and then increased as the distance increased. For the "3+0+1" network configuration, the PDOP values also had similar trends, but decreased significantly with the surface buoy. When we replaced the surface buoy with the subsurface buoy, and when the height of the buoy was greater than 1 km, the PDOP values decreased first and then increased as the distance increased; when the height was less than 1 km, the PDOP values also decreased with increasing distance. When the buoy height was less than 2.5 km, the PDOP average and standard deviation of the "3+1+0" network configuration were greater than those of the "3+1+1" network configuration. Moreover, when the subsurface buoy of the "3+1+0" network configuration was close to the sea surface, the positioning indicators were close to the "3+0+1" network configuration, and there was no significant difference between the functions of the subsurface buoy and the sea surface buoy.
In general, configurations without subsurface or surface buoys were unable to achieve positioning when the seafloor station distance exceeded 19 km, which can be avoided in the "3+1+1" network configuration with a subsurface buoy and sea surface buoy. In addition, the height of the subsurface buoy can be adjusted to improve the performance of the geodetic network. Therefore, the subsurface buoy is a necessary and important component to make UGNs more resilient and precise.
Although the sea surface buoy and the subsurface buoy had similar performances in kinematic positioning, the surface buoy could be equipped with a satellite communication module for transmitting the UGN data, receiving remote-control commands, etc. Moreover, the GNSS equipment on the surface buoy can be used for the clock synchronization of acoustic equipment, and the acoustic signal from the buoy can send time synchronization information to other stations in the UGN, which is crucial to maintaining accurate system time synchronization.

Configuration with More Seafloor Base Stations
Generally, four seafloor base stations are used in traditional networks [3,28], and the basic "3+1+1" network configuration can be expanded to the "4 seafloor base stations + 1 subsurface buoy + 1 sea surface buoy" ("4+1+1") network configuration. Figure 5 illustrates this expanded configuration. We analyzed two kinds of seafloor base station array geometries, the triangular-center-form array geometry and the square-form array geometry. Figure 5a,b show that, for the "4+1+1" network configurations, the EPG ASN≥3 values fluctuated slightly when the base station distance was less than 19 km. Compared with the "3+1+1" network configuration (Figure 2a), when the distance was greater than 19 km, the EPG ASN≥3 of the "4+1+1" network configuration (triangular-center-form array geometry) showed an obvious increase, and the percentage of EPG ASN≥3 changed from 21.62% to 46.51%. In contrast, the EPG ASN≥3 values of the "4+1+1" network configuration (squareform array geometry) decreased significantly, and the minimum percentage of EPG ASN≥3 was 6.72%. Figure 5c,d show that, when the distance was less than 13 km, the PDOP values of the triangular-center-form array geometry were slightly greater than those of the squareform array geometry. When the distance was greater than 13 km, the PDOP values of the square-form array geometry increased significantly and were greater than those of the triangular-center-form array. Moreover, compared with the "3+1+1" network configuration (Figure 3b), the PDOP values of the expanded network configuration (square-form array geometry) showed similar trends, but decreased slightly.
In general, the performance of positioning range and positioning accuracy can be improved by increasing the number of seafloor base stations. The improvement of kinematic positioning performance is also closely related to the array geometry. However, more seafloor base stations always mean more costs in equipment, construction, and maintenance.
To optimize the cost and efficiency of the UGN, the basic "3+1+1" network configuration is more appropriate.
form array geometry) decreased significantly, and the minimum percentage of ASN 3 EPG  was 6.72%. Figure 5c,d show that, when the distance was less than 13 km, the PDOP values of the triangular-center-form array geometry were slightly greater than those of the squareform array geometry. When the distance was greater than 13 km, the PDOP values of the square-form array geometry increased significantly and were greater than those of the triangular-center-form array. Moreover, compared with the "3+1+1" network configuration (Figure 3b), the PDOP values of the expanded network configuration (square-form array geometry) showed similar trends, but decreased slightly.

Impact of Acoustic Transmission Characteristics on Acoustic Signal Propagation
In the simulation section, the Euclidean distance between the base station and the user is adopted. However, the acoustic ray bending in the acoustic signal propagation path should be considered. Additionally, the acoustic transmission characteristics, such as the multipath effect caused by the acoustic signal reflected from the seabed or the sea surface, the Doppler effect caused by the motion of the acoustic source, and the transmission loss of acoustic energy, also affect acoustic signal propagation in the underwater environment. Herein, we simulated the acoustic transmission characteristics by BELLHOP [38] and studied the impact on acoustic signal propagation. The frequency of the omni-direction acoustic sonar was set to 8 kHz, and the users were randomly distributed in an underwater space of 3000 m in depth. To illustrate the impact on acoustic signal propagation, the signal service range was limited to 30 km in the horizontal, which was different from the signal positioning range in Section 2.3. Note that it is complex to simulate the real underwater environment, such as the temporal-spatial-varying ocean sound speed, sea surface wind and waves, and the seafloor topography. In the simulation, the typical deep-sea Munk model was used as the reference sound velocity profile [39], and a calm sea surface and flat seafloor were considered, as shown in Figure 6. tion, the signal service range was limited to 30 km in the horizontal, which was different from the signal positioning range in Section 2.3. Note that it is complex to simulate the real underwater environment, such as the temporal-spatial-varying ocean sound speed, sea surface wind and waves, and the seafloor topography. In the simulation, the typical deep-sea Munk model was used as the reference sound velocity profile [39], and a calm sea surface and flat seafloor were considered, as shown in Figure 6.

Acoustic Ray Bending
We first studied the direct acoustic signal from the base station to the user, excluding the reflected signal from the seafloor or the sea surface. Due to acoustic ray bending, there was an acoustic shadow where the sound signal could not reach. By finding the critical reversal point of the bending ray, where the ray grazed the horizontal, we referred to the horizontal distance between the reversal point and the base station as the maximum service range. For the fixed base station, when the user was within the maximum range, the acoustic signal could reach the user directly by adjusting the initial grazing angle; when the user was outside the maximum range, there existed the acoustic shadow. Therefore, we could investigate the impact of acoustic ray bending on the maximum service range.

Acoustic Ray Bending
We first studied the direct acoustic signal from the base station to the user, excluding the reflected signal from the seafloor or the sea surface. Due to acoustic ray bending, there was an acoustic shadow where the sound signal could not reach. By finding the critical reversal point of the bending ray, where the ray grazed the horizontal, we referred to the horizontal distance between the reversal point and the base station as the maximum service range. For the fixed base station, when the user was within the maximum range, the acoustic signal could reach the user directly by adjusting the initial grazing angle; when the user was outside the maximum range, there existed the acoustic shadow. Therefore, we could investigate the impact of acoustic ray bending on the maximum service range.
The relationship between the initial grazing angle α source of the acoustic signal at the base station and the grazing angle α reversal at the reversal point is as follows [40]: where p is the ray constant, c is the sound speed, and α is the angle of the ray with respect to the horizontal, for which the upward angle is negative and the downward angle is positive. Using the initial grazing angle α source , the horizontal distance x between the reversal point and the base station is calculated by the 2-D planar ray-trace method, as follows [40]: where g is the gradient of the sound speed, n is the number of layers from the source to the reversal point. Figure 7a shows the acoustic propagation path of the base station at different depths, and acoustic rays graze the sea surface and the seafloor, respectively. For the rays that grazed the sea surface, fixing the base station and decreasing the initial grazing angle of the upward acoustic signal, the ray could reach the area within the maximum service range on the sea surface. In contrast, when the initial grazing angle of the upward acoustic signal increased, the reversal point was away from the sea surface, and the ray could not reach the area between the reversal point and the sea surface. For the ray that grazed the seafloor, the situation was similar to that of the ray that grazed the sea surface. When the base station was fixed, the ray could reach the area within the maximum service range on the seafloor by increasing the initial grazing angle of the downward acoustic signal. In contrast, when the initial grazing angle of the upward acoustic signal decreased, the reversal point was away from the seafloor, and the ray could not reach the area between the reversal point and the seafloor. Figure 7b shows the maximum service range on the sea surface and the seafloor, and the maximum range varied with the depth of the base station. For the maximum service range on the sea surface, it gradually increased with increasing depth. Decreasing the height of the subsurface buoy could increase the maximum service range on the sea surface, and the maximum service range of the seafloor base station was close to 15 km. For the maximum service range on the seafloor, the ray grazed on the seafloor when the sound speed satisfied c source ≤ c reversal according to Equation (9). For the Munk sound velocity profile, the sound speed in the water shallower than 400 m could not satisfy the condition, and the rays from the depth of 0~400 m were reflected or absorbed at the seafloor. When the depth was greater than 400 m, the sound speed satisfied the condition, and the ray could graze on the seafloor. The maximum service range on the seafloor gradually decreased as the depth increased, and the ray was refracted upward when the base station was on the seafloor, for which the maximum service range on the seafloor was 0 m and not plotted.
In general, the effect of ray bending on the service range can be reduced by lifting the station height. For the seafloor base station, increasing the erection height is costly and laborious, and weakens the stability of the station. Therefore, subsurface buoys can be an alternative station to the seafloor base station to increase the service range. In addition, this study only considered the case of flat seafloor, and the effects of seafloor topographic features, such as steep slopes and sea hills, will be considered in future studies.

Multipath and Doppler Effects
Multipath signals, such as reflected signals from the seabed, the sea surface, the equipment, etc., affect acoustic signal propagation [19], and misrecognizing reflected signals as direct signals also leads to acoustic time delay errors [41]. The time deviation of the reflected signal from the direct signal was simulated by BELLHOP for base station depths of 500 m and 2998 m, as shown in Figure 8a,b. The results show that the multipath signals caused time deviations of milliseconds to seconds. In general, methods such as the signal cross-correlation function (CCF) accompanied by the Energy Ratio (ER) algorithm and Phase-Correlation-Only (PCO) method can be utilized to identify the direct acoustic signal from the multipath signal [19][20][21]41,42]. Figure 8a,b also show that the propagation time deviation decreased with increasing distance between the base station and the user. For a fixed user, the propagation time deviations of the shallow base station are significantly less than those of the deep base station. Therefore, the subsurface buoy helps to reduce the propagation time deviation. However, for the impact of decreasing the time deviation, on the one hand, the time deviation can be neglected when its value is close to a microsecond; on the other hand, it's hard to distinguish the multipath signal when a sub-millisecond time deviation exists [19,21,43]. Nevertheless, these cases were not present in the simulation.
In addition, the acoustic signal cannot reach the user directly due to ray bending; e.g., with the seafloor base station at 2998 m depth and the user at 2900 m depth, there was no Figure 7. Simulation of acoustic ray bending. (a) Acoustic propagation path of the base station at different depths, with the rays grazing the sea surface or the seafloor. Blue and orange colors denote the rays grazing the sea surface and the seafloor, respectively; (b) base station maximum service range on the sea surface and the seafloor. Blue and orange colors denote the rays grazing the sea surface and the seafloor, respectively.

Multipath and Doppler Effects
Multipath signals, such as reflected signals from the seabed, the sea surface, the equipment, etc., affect acoustic signal propagation [19], and misrecognizing reflected signals as direct signals also leads to acoustic time delay errors [41]. The time deviation of the reflected signal from the direct signal was simulated by BELLHOP for base station depths of 500 m and 2998 m, as shown in Figure 8a,b. The results show that the multipath signals caused time deviations of milliseconds to seconds. In general, methods such as the signal cross-correlation function (CCF) accompanied by the Energy Ratio (ER) algorithm and Phase-Correlation-Only (PCO) method can be utilized to identify the direct acoustic signal from the multipath signal [19][20][21]41,42]. Figure 8a,b also show that the propagation time deviation decreased with increasing distance between the base station and the user. For a fixed user, the propagation time deviations of the shallow base station are significantly less than those of the deep base station. Therefore, the subsurface buoy helps to reduce the propagation time deviation. However, for the impact of decreasing the time deviation, on the one hand, the time devia-tion can be neglected when its value is close to a microsecond; on the other hand, it's hard to distinguish the multipath signal when a sub-millisecond time deviation exists [19,21,43]. Nevertheless, these cases were not present in the simulation.
In addition, the acoustic signal cannot reach the user directly due to ray bending; e.g., with the seafloor base station at 2998 m depth and the user at 2900 m depth, there was no direct acoustic signal when the horizontal distance was more than 5 km. In contrast, the acoustic signal from the subsurface buoy of 500 m depth could reach the user directly under the same conditions, which also shows that the deployment of a subsurface buoy could increase the service range.
Due to the relative motion between the target and the base station, the acoustic signal generated a Doppler shift. During acoustic signal detection, the optimal frequency shift rate of the signal was generally estimated by searching the maximum cross-correlation coefficient to correct for the effect of the Doppler shift [20].
In general, both multipath and Doppler effects can be weakened or corrected during acoustic signal detection, which is beyond the scope of this paper and will be considered in further studies.

Acoustic Transmission Loss
Multipath signals also result in acoustic energy transmission loss ( TL ), and the transmission loss of the acoustic signal can be assessed as follows [38]:  (11) where d is the distance between the station and the user, f is the signal frequency, and coefficient k represents the propagation geometry space. Figure 9a,b show the transmission loss when the base station was at a depth of 500 m. Figure 9a shows that the transmission loss increased from 40 dB to 80 dB within the horizontal distance of 5 km. The transmission loss then increased slowly and approached 100 dB when the horizontal distance was 10 km. However, there was an obvious caustic curve above the depth of 2500 m-the transmission loss changed significantly from the left side of the curve to the right side. On the right side of the caustic curve, the rays reflected by the seafloor or the sea surface with large reflection loss, and the zone of large transmission loss decreased with increasing station depth. Therefore, for the station in shallow water, such as a subsurface buoy, the acoustic signals had a large transmission loss in the shallow water. Figure 9b shows the rapid increase in the transmission loss. When the user was at the depth of 100 m, there was a drastic change in the transmission loss at a horizontal distance of 6 km, while the drastic change occurred at a horizontal distance of 10.5 km when the user was at the depth of 1000 m. Figure 9c,d show the transmission loss when the base station was at a depth of 2998 m. Figure 9c shows that, when the horizontal distance was within 2.5 km, the transmission

Acoustic Transmission Loss
Multipath signals also result in acoustic energy transmission loss (TL), and the transmission loss of the acoustic signal can be assessed as follows [38]: where d is the distance between the station and the user, f is the signal frequency, and coefficient k represents the propagation geometry space. Figure 9a,b show the transmission loss when the base station was at a depth of 500 m. Figure 9a shows that the transmission loss increased from 40 dB to 80 dB within the horizontal distance of 5 km. The transmission loss then increased slowly and approached 100 dB when the horizontal distance was 10 km. However, there was an obvious caustic curve above the depth of 2500 m-the transmission loss changed significantly from the left side of the curve to the right side. On the right side of the caustic curve, the rays reflected by the seafloor or the sea surface with large reflection loss, and the zone of large transmission loss decreased with increasing station depth. Therefore, for the station in shallow water, such as a subsurface buoy, the acoustic signals had a large transmission loss in the shallow water. Figure 9b shows the rapid increase in the transmission loss. When the user was at the depth of 100 m, there was a drastic change in the transmission loss at a horizontal distance of 6 km, while the drastic change occurred at a horizontal distance of 10.5 km when the user was at the depth of 1000 m. Figure 9c,d show the transmission loss when the base station was at a depth of 2998 m. Figure 9c shows that, when the horizontal distance was within 2.5 km, the transmission loss gradually increased from 40 dB to 70 dB, and then the increase in the transmission loss slowed down. When the horizontal distance was 15 km, the transmission loss was close to 100 dB. Two obvious caustic curves could be found, and the depth of the turning point was about 500 m. For the caustic curve above 500 m, the zone of large transmission loss decreased with increasing station depth. In contrast, the zone of large transmission loss increased with increasing station depth for the other caustic curve. Figure 9d shows the rapid increase in the transmission loss when the base station was at a depth of 2998 m. When the user was at a depth of 100 m, the horizontal distance of drastic change increased to 16.5 km, and when the user was at a depth of 1000 m, the horizontal distance increased to 21.5 km. Therefore, for the base station in deep water, the acoustic signals had a large transmission loss in both shallow and deep water.
From the above analysis, when the base stations were located at 500 m and 2998 m depths simultaneously, they could compensate for the signal transmission loss in shallow and deep water, which further indicates the necessity of deploying the subsurface buoy and seafloor base stations.

Conclusions
Our study proposes an underwater geodetic network configuration consisting of three seafloor base stations, one underwater subsurface buoy, and one sea surface buoy. The main purpose of the UGN is to provide kinematic positioning for regional underwater submersibles for ocean near-bottom surveying and precision engineering. To evaluate the performance of the network configuration, the positioning range and accuracy were analyzed by changing the seafloor base station separation and subsurface buoy height. The

Conclusions
Our study proposes an underwater geodetic network configuration consisting of three seafloor base stations, one underwater subsurface buoy, and one sea surface buoy. The main purpose of the UGN is to provide kinematic positioning for regional underwater submersibles for ocean near-bottom surveying and precision engineering. To evaluate the performance of the network configuration, the positioning range and accuracy were analyzed by changing the seafloor base station separation and subsurface buoy height. The simulation results indicated that, in a deep-sea environment with a water depth of 3000 m, the positioning accuracy and range were mainly affected by the distance between the seafloor base stations, and the height of the subsurface buoy had a greater impact on the positioning accuracy than on the positioning range. In general, when the distance between seafloor base stations was 9~13 km and the height of the subsurface buoy above the seafloor was lower than 2.5 km, the configuration had the optimal kinematic positioning performance, which could achieve a submersible positioning accuracy of less than 30 m within an underwater space of 25 km × 25 km × 3 km. Compared with the network configurations with fewer or more stations, the proposed configuration with the underwater subsurface buoy and the sea surface buoy made the kinematic positioning performance of the UGN more precise and cost-effective.
The purpose of our study was to evaluate the configuration of UGNs on kinematic positioning performance. Based on the proposed evaluation method in this study, an optimal network configuration with the best positioning performance can be obtained under various conditions, such as the effective sonar ranging distance and accuracy, the water depth, and the underwater base station configuration. However, a variety of other factors, including the underwater acoustic environment, the base station coordinate accuracy, and stability, affect the user kinematic positioning performance, which should be thoroughly investigated in future studies.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.