BRIEF RESEARCH REPORT article
Optimal Transponder Array and Survey Line Configurations for GNSS-A Observation Evaluated by Numerical Simulation
- 1Hydrographic and Oceanographic Department, Japan Coast Guard, Tokyo, Japan
- 2Institute of Industrial Science, University of Tokyo, Tokyo, Japan
The Global Navigation Satellite System-Acoustic ranging combination technique (GNSS-A) has enabled us to measure seafloor crustal deformation in the precision of centimeters, leading to numerous discoveries of subseafloor tectonic phenomena. The moving observation conducted by our research group allows us to measure both the horizontal and vertical absolute positions of a reference point on the seafloor. However, the observation frequency of our GNSS-A observation system is still insufficient to observe short-term phenomena. This paper focused on the possibility to reduce the observation time per a seafloor site by shrinking the seafloor transponder array size and the survey line radius, which were empirically defined to be equal to the seafloor site depth in the early research. We evaluated the effects of changing these sizes on the GNSS-A positioning accuracy by conducting a series of numerical experiments. The results of the numerical experiments indicated that for a seafloor site with a depth of 3,000 m, the positioning accuracy is rapidly degraded as the transponder array size and the survey line radius are reduced to less than 3,000 m. Additional experiments done for transponder array sizes and survey line radii around 2,000–4,000 m revealed that shrinking the survey line radius has a dominant effect on the decrease in positioning accuracy. Thus, shrinking the transponder array size and the survey line radius is not a suitable option for reducing observation time, and the empirically defined observation configurations are concluded to be quite optimal when regarding both the positioning accuracy and the observation time.
The Global Navigation Satellite System-Acoustic ranging combination technique (GNSS-A) has realized the precise measurement of a reference point on the seafloor, expanding the crustal observation network into the ocean (Spiess, 1985; Asada and Yabuki, 2001; Fujita et al., 2006). In the GNSS-A observation of the Japan Coast Guard (JCG), the precise position of a survey vessel sailing along a survey line is measured using kinematic GNSS positioning, and the distances between the onboard transducer and the mirror transponders installed on the seafloor are measured using acoustic ranging (Figure 1A). Combining these data with the vessel attitude data and XBT/XCTD observations, the global coordinates of the virtual reference point (the centroid of the transponders) is precisely determined in the centimeter level. GNSS-A observation has revealed numerous tectonic phenomena such as the coseismic and postseismic movements of the 2011 Tohoku-oki earthquake
FIGURE 1. Schematic images of the current GNSS-A observation system conducted by the Japan Coast Guard. (A) is the overview of the entire GNSS-A observation system, and (B) indicates the configurations of the transponder array and the survey line. As shown in (B), the seafloor site depth (D), transponder array size (L), and survey line radius (R) are set to be equal.
The observation frequency has been a crucial matter for GNSS-A observation over the past 20 years, and we have been making significant efforts to improve the observation frequency. Despite the increase in the number of our seafloor sites, the observation frequency of the JCG group has improved from approximately one time per year per site in the early 2000s to 4 times per year per site by 2018 (Ishikawa et al., 2020). However, this frequency is still inadequate for observing tectonic phenomena with time scales less than few months (e.g. shallow short-term slow slip events).
One of the key factors that affects the observation efficiency is the size of the transponder array and the survey line, as shown in Figure 1B. In the early research by Spiess (1985), the proposed stationary observation system required at least three transponders, and the transponder array size has been defined that the circumdiameter of the transponder array should be approximately equal to the depth of the seafloor site.
However, the stationary system using a triangular transponder array proposed by Spiess (1985) was unable to detect the vertical displacement. To overcome this problem of detecting the vertical displacement, our research group has taken a different approach. We have developed an observation system that uses a vessel to perform a moving survey (Figure 1A), which is capable to detect both the horizontal and vertical displacements. In our system, although it is possible to detect the movement of only a single transponder, we deploy four transponders to form a square array as shown in Figure 1A to increase the positioning accuracy and stability. The transponder array size in our observation system is configured to be equal to the site depth, as done in Spiess (1985).
Regarding the survey line, theoretical studies by a research group in China have indicated that the optimal survey line radius is
One possible approach to reduce the observation time for the deep sites is by reducing the transponder array size and the survey line radius which were empirically defined in the early studies (e.g. Spiess 1985). However, the effects of changing the transponder array size and the survey line radius on the accuracy of GNSS-A positioning have not been clearly investigated. In this study, we evaluated the effects of changing the transponder array size and the survey line radius on the accuracy of GNSS-A observation, by analyzing synthetically generated datasets.
We created a series of synthetic datasets with different transponder array sizes and survey line radii by using a numerical simulator (Yokota et al., 2016b) to evaluate the effects of changing these sizes on the positioning accuracy. In our numerical experiments, we assumed a simple square transponder array and assigned a symmetrical survey line. To focus only on the effects of changing the geometry of the observation system, we configured a simple sound speed structure with no horizontal and temporal variations. We analyzed these datasets like the actual GNSS-A datasets using the conventional linearized inversion method (Fujita et al., 2006).
The GNSS-A Simulator
We used a numerical simulator that has been developed to create synthetic GNSS-A observation datasets (Yokota et al., 2016b). In the GNSS-A simulator, the synthetic travel times are calculated between the synthetic vessel positions and the preconfigured transponder positions, and the outputs are formatted to be analyzed like an actual GNSS-A observation dataset. The GNSS-A simulator requires the following input files (also shown in Figure 2):
(I-1) Transponder position and error configuration file
(I-2) Survey line configuration file
(I-3) Sound speed profile configuration file (for synthetic acoustic travel time data)
(I-4) Sound speed profile configuration file (as pseudo-XBT data used in analysis)
FIGURE 2. The flow chart of our numerical experiment. The input files (I-1), (I-2), (I-3), and (I-4) are used for the configurations of the numerical simulation by the GNSS-A simulator. The outputs of the GNSS-A simulator (the synthesized GNSS-A data) are analyzed using the analysis software SGOBS version 4.0.2, in which the linearized inversion method by Fujita et al. (2006) is implemented (right side of the figure).
In our simulation, the three-dimensional position is expressed by local East-North-Up (ENU) coordinates, with the origin at the sea surface directly above the centroid of the transponder array. Configuration file (I-1) contains the three-dimensional positions of the transponders, which are written in the local ENU coordinates. The acoustic travel time error, pseudo-random walk of the GNSS, and Gaussian noise of the GNSS are also configured in this file. Survey line configuration file (I-2) contains the survey line geometry and size, number of acoustic shots and shot interval for each survey line, and the pseudo-random walk and Gaussian noises of the survey line (which express random distortions of the survey line geometry). In the sound speed profile configuration file for the theoretical ray path calculation (I-3), the spatiotemporal variation of the sound speed profile can be configured. There is also a sound speed profile configuration file (I-4) for the pseudo-XBT data used in the analysis. The sound speed profile in (I-3) and (I-4) are configured at several depth layers; the sound speed at depths between these configured depth layers are either linearly interpolated or extrapolated. Settings for the vessel attitude and the offset between the GNSS antenna and the transducer are not implemented in this simulator; thus, the calculated vessel positions are equivalent to the transducer positions.
The resulting outputs of the calculated data are the acoustic travel time file, vessel position file, and the sound speed observation file. These output files can be analyzed like an actual GNSS-A dataset using the method of Fujita et al. (2006), as described in Figure 2 and Analysis of the Synthetic Datasets.
Configurations of the Numerical Experiments
In our study, we mainly edited input files (I-1) and (I-2) to configure the transponder positions and the survey line radius for each experimental run, and kept the settings in files (I-3) and (I-4) constant throughout all of the numerical experiments.
Since the focus of our study is on the possibility to reduce the observation times for the deep sites, we have set the depth of the seafloor site in file (I-1) to 3,000 m, which is a representative depth of the deepest seafloor sites installed by our research group. In our numerical experiments, the depths of all transponders were set to 3,000 m. This value was kept constant for all of the experimental runs. The transponder positions were set so that the transponders form a square array, with the centroid at position (0 m, 0 m, −3,000 m) in local ENU coordinates. For example, when we set the transponder array size to L, we assigned the transponders at (L/2, 0 m, −3,000 m) (−L/2, 0 m, −3,000 m) (0 m, L/2, −3,000 m), and (0 m, −L/2, 3,000 m), in local ENU coordinates. The transponder array size was changed depending on the experimental run.
In the survey line configuration file (I-2), a symmetrical survey line configuration currently adopted in our GNSS-A observation, as shown in Figure 1B, is configured. Like the transponder positions, the survey line radius was also changed depending on the type of the experimental run. For each trial of data synthesis and analysis, a total of four sets with each having 784 acoustic shots were simulated to create a single GNSS-A dataset with a total of 3,136 shots. The acoustic shots were created at a constant time interval, to synthesize a geometrically well-balanced data, as we do in the actual observations. The number of acoustic shots per line and the shot intervals were tuned so that the total simulated acoustic shots is 3,136 for each trial, regardless of the survey line radius.
We configured a simple sound speed structure in files (I-3) and (I-4) that has no horizontal sound speed gradient nor temporal variation. Such a simplified sound speed structure was used instead of a more realistic model (e.g. Munk, 1974), since our focus is on how the GNSS-A positioning accuracy is affected by changing the transponder array size and the survey line radius. The configured values for the ray path calculation in (I-3) were 1,505.00 m/s at 0 m, 1,494.30 m/s at 400 m, and 1,490.00 m/s at depths below 600 m, as shown in the left of Figure 3A. The acoustic shots calculated using these values have round-trip travel times ranging roughly 4–10 s and takeoff angles ranging roughly from
FIGURE 3. The configurations of our numerical experiment. (A) shows the vertical sound speed profile (SSP) configured for the ray path calculation (left, blue line) and the vertical SSP used in the data analysis (right, red line). Note that the difference between these two profiles is so small that the difference cannot be distinguished in this figure (see Configurations of the Numerical Experiments for details). (B) indicates sample waveforms of the GNSS pseudo-random walk added to the synthesized data (from top to bottom, the east-west, north-south, and up-down components). (C) is the schematic diagram of the settings of EXP1, EXP2, and EXP3, from left to right (see Configurations of the Numerical Experiments for details).
We created high-rate pseudo-GNSS noises using the Box-Muller’s method (Box and Muller, 1958). For the pseudo-GNSS noise, we assigned a Gaussian noise and a random walk for the vertical and horizontal components. Gaussian noise and random walk were empirically known to be included in the GNSS noises (Bilich et al., 2008). The values of 1
We performed a series of experimental runs as shown in Figure 3C to evaluate the individual and the combined effects of changing the transponder array size and the survey line radius on the GNSS-A positioning accuracy. In our first experiment, only the transponder array size was changed from 1,000 m to 8,000 m at 1,000 m intervals, and the survey line radius was fixed to 3,000 m (EXP1 in Figure 3C). In our second experiment, only the survey line radius was changed from 1,000 m to 8,000 m at 1,000 m intervals, and the transponder array size was fixed to 3,000 m (EXP2 in Figure 3C). In the third experiment (EXP3 in Figure 3C), both the transponder array size and the survey line radius were simultaneously changed from 1,000 m to 8,000 m at 1,000 m intervals (e.g. when the transponder array size is configured to 8,000 m, the survey line radius is also configured to 8,000 m). Additionally, we have also performed an experiment with a transponder array size of 2,000 m and a survey line radius of 4,000 m, and another experiment with a transponder array size of 4,000 m and a survey line radius of 2,000 m, to closely examine the effects of transponder array sizes and survey line radii around 2,000–4,000 m. As stated above in this section, note that the site depth is fixed to 3,000 m for all experimental runs.
To evaluate the positioning accuracy, we conducted 300 trials of data synthesis and analysis for each experimental run. For each trial, we created pseudo-GNSS noises as assigned in file (I-1), and added the noises to the synthetic data, to create data variation among the trials.
Analysis of the Synthetic Datasets
For the analysis of the synthetic datasets, we used the analysis software SGOBS version 4.0.2, in which the method developed by Fujita et al. (2006) is implemented. SGOBS version 4.0.2 is regularly used in the observation and analysis routine of our group. In the method of Fujita et al. (2006), the sound speed structure and the three-dimensional transponder positions are alternately estimated using Bayesian inversion. As shown in the right of Figure 2, this algorithm first iterates to determine the positions of the seafloor transponders using the initial sound speed profile, and then the resulting residuals of the first loop are used in the second loop to determine the sound speed coefficient. Then, the transponder positions are determined again using the newly estimated sound speed coefficients, and so on. These iterations are repeated until both the sound speed and the transponder positions satisfy the convergence criteria. The centroid of the transponder positions is calculated from the final estimations of the transponder positions.
In the actual GNSS-A observations, acoustic shots with extremely small takeoff angles may be due to fake return signals, such as echoes from the seafloor. In order to remove such acoustic shots that may not be responses from the transponders, acoustic shots with takeoff angles less than
To evaluate the GNSS-A positioning accuracy, we defined three criteria as listed below:
The average horizontal error of the centroid (
The average horizontal error of the four transponders (
The average vertical error of the centroid (
The average horizontal error of the centroid (
Thus, the average horizontal error of the four transponders for
The average vertical error of the centroid
In Results and Discussion, the results will be evaluated by calculating and plotting the values of
Figures 4A–C indicate the values of
FIGURE 4. The results of the numerical experiments. Figure (A) is the plot of the average horizontal error of the centroid
Meanwhile, the values of
Figure 4C shows that the average error of the vertical position
Figures 4D–F shows the histograms of the horizontal error of the centroid, the average horizontal error of the four transponders, and the vertical error of the centroid, respectively. All histograms indicate the frequency of errors in the range of 0–4 cm. Figures 4D,E show that the distribution peaks of the horizontal errors become generally sharper, and the peak gradually shifts to 0 cm as the transponder array size and the survey line radius increase. This indicates that the horizontal positioning accuracy generally improves as these sizes increase. In Figure 4F, the histograms of the vertical error, the distribution peak shifts toward 0 cm as the transponder array size and the survey line radius increase. This result shows that the vertical positioning accuracy also improves by increasing these sizes.
Effects of Changing the Transponder Array Size
Figure 4B indicates that
The values of
However, the values of
Effects of Changing the Survey Line Radius
In general, larger survey line radius enhances the variation of the takeoff angles of the acoustic shots, which leads to improved positioning accuracy, as shown in Figure 4. This is analogous to a good constellation geometry in GNSS positioning; it is better that the GNSS satellites spread out evenly across the sky rather than being confined in a very small area.
Shrinking the survey line radius would eventually converge as a stationary observation at the sea surface directly above the centroid. The takeoff angle variation of the acoustic shots decreases since the survey range is confined to a smaller area. Unlike the transponder arrangement proposed by Tomita et al. (2019), our square transponder array cannot determine the vertical displacement by stationary observation. Our results shown in Figure 4C indicate that
When the survey line radius is increased, all of the errors tend to become very small. However, it must be noted that the availability of the acoustic data in the analyses for a very large survey line radius is very poor. For example, when the survey line radius is 8,000 m, the data availability is around 30–50%. This is due to the takeoff angle limitation in our GNSS-A observation to exclude acoustic signals transmitted from takeoff angles that are less than
The Optimal Transponder Array Size and Survey Line Radius
We can imply from Figure 4 that when regarding the horizontal and vertical positioning accuracies, the transponder array size and the survey line radius should be as large as possible. However, due to the takeoff angle limitation explained in Effects of Changing the Survey Line Radius, the optimal transponder array size and survey line radius for our observation system is limited to approximately 3,000–4,000 m. When regarding both the observation time and the positioning accuracy, the transponder array size and survey line radius of 3,000 m, which is equal to the empirically defined sizes for a depth of 3,000 m, seems to be quite optimal. When the sizes are less than 3,000 m, the positioning errors become too large for crustal monitoring. When the sizes are greater than 3,000 m, the improvement of the positioning accuracy is small compared to the significant increase in the observation time. It can be concluded that the trade-off between the positioning accuracy and the observation time is well balanced at our empirically defined sizes. Also, it must be noted that in the real world, increasing the transponder array size and the survey line radius would increase errors due to the spatial variability of the sound speed structure.
Figure 4A has indicated that the horizontal positioning error of the centroid of the transponders remains stable when the transponder array size and/or the survey line radius is changed. This result implies that it is possible to obtain the horizontal positions of a site with a transponder array size and a survey line radius that is smaller than our empirically defined sizes. One of the limitations when installing a site is the seafloor topography; there are cases where we cannot install a site with the empirically defined transponder array size, due to the lack of a smooth seafloor. When regarding only the horizontal positioning accuracy, we can shrink the transponder array size to some degree, which would give us more valuable options when installing a site. Also, the result that the horizontal positioning is not significantly affected by shrinking the survey line radius indicates a possibility that we can reduce the observation time. By reducing the observation time and increasing the observation frequency, we can improve the accuracy of the time series of the horizontal positions. An improved time series of the horizontal positions would provide us valuable information on the temporal variation of the seafloor crustal deformation.
Our simulation results may also be applied to the observations of other institutions, with different survey line geometry. As explained in Effects of Changing the Survey Line Radius, assigning a large survey line radius enables us to obtain acoustic data with larger variation of takeoff angles, which is necessary for achieving high positioning accuracies for both the horizontal and the vertical components. Thus, it may be possible that the optimal sizes indicated in our study can be applied to other types of survey line geometry, if the survey line geometry is configured to obtain acoustic data with large variation of the takeoff angles.
The results of the numerical experiments indicated that the positioning accuracy is degraded rapidly assuming transponder array size and survey line radius of less than 3,000 m for a site with 3,000 m of depth. Thus, shrinking the transponder array size and the survey line radius is not a suitable solution for improving observation frequency. In fact, the empirically defined sizes, which is equal to 3,000 m in the experiments of our study, are quite optimal values when regarding both the positioning accuracy and observation efficiency.
GNSS-A observation has numerous error sources other than those focused in this study. One of the main error sources is the spatiotemporal variation of the sound speed structure, which has been investigated in numerous studies (e .g. Yasuda et al., 2017; Honsho et al., 2019; Yokota et al., 2019). Assessing the effects of the temporal change of the sound speed structure and the horizontal sound speed gradient on the GNSS-A accuracy is one of the main topics for our future research using the GNSS-A simulator. Numerical experiments using different transponder arrangements such as those proposed by Tomita et al. (2019) is an important topic as well, since there is a possibility that the positioning accuracy may improve by changing the arrangement of the transponders. Also, a new analysis software that directly estimates the sound speed structure using a Bayesian approach is being developed (Watanabe et al., 2020). We can utilize the simulator to validate the capability of this new software to accurately determine the transponder positions and to properly estimate the sound speed structure. We plan to further investigate the GNSS-A observation errors by taking the approach using numerical simulation, as done in this study.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
YN designed the study and wrote this manuscript. YY developed and tuned the GNSS-A simulator. YN, YY, TI, SW discussed about the simulation configurations and the results, and commented to improving the manuscript.
The submission of this manuscript was funded by the Japan Coast Guard.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
We thank the many staff members from the Hydrographic and Oceanographic Department, Japan Coast Guard, including the crew of the survey vessels Takuyo, Shoyo, Meiyo, and Kaiyo for their support in our observations and technological developments. We especially thank the active senior staff members from the Geodesy and Geophysics Office, Hydrographic and Oceanographic Department, Japan Coast Guard, for their devoted maintenance and management of the equipment.
Fujita, M., Ishikawa, T., Mochizuki, M., Sato, M., Toyama, S., Katayama, M., et al. (2006). GPS/acoustic seafloor geodetic observation: method of data analysis and its application. Earth Planets Space 58, 265–275. doi:10.1186/BF03351923
Honsho, C., Kido, M., Tomita, F., and Uchida, N. (2019). Offshore postseismic deformation of the 2011 Tohoku earthquake revisited: application of an improved GPS-acoustic positioning method considering horizontal gradient of sound speed structure. J. Geophys. Res. Solid Earth 124. doi:10.1029/2018JB017135
Ishikawa, T., Yokota, Y., Watanabe, S., and Nakamura, Y. (2020). History of on-board equipment improvement for GNSS-A observation with focus on observation frequency. Front. Earth Sci. 8, 150. doi:10.3389/feart.2020.00150
Sato, M., Ishikawa, T., Ujihara, N., Yoshida, S., Fujita, M., Mochizuki, M., et al. (2011). Displacement above the hypocenter of the 2011 Tohoku-oki earthquake. Science 332, 1395. doi:10.1126/science.1207401
Tomita, F., Kido, M., Honsho, C., and Matsui, R. (2019). Development of a kinematic GNSS-Acoustic positioning method based on a state-space model. Earth Planets Space 71, 102. doi:10.1186/s40623-019-1082-y
Watanabe, S., Ishikawa, T., Yokota, Y., and Nakamura, Y. (2020). GARPOS: analysis software for the GNSS-A seafloor positioning with simultaneous estimation of sound speed structure. Front. Earth Sci. 8, 597532. doi:10.3389/feart.2020.597532
Watanabe, S., Sato, M., Fujita, M., Ishikawa, T., Yokota, Y., Ujihara, N., et al. (2014). Evidence of viscoelastic deformation following the 2011 Tohoku-oki earthquake revealed from seafloor geodetic observation. Geophys. Res. Lett. 41, 5789–5796. doi:10.1002/2014GL061134
Yasuda, K., Tadokoro, K., Taniguchi, S., Kimura, H., and Matsuhiro, K. (2017). Interplate locking condition derived from seafloor geodetic observation in the shallowest subduction segment at the Central Nankai Trough, Japan. Geophys. Res. Lett. 44, 3572–3579. doi:10.1002/2017GL072918
Yokota, Y., Ishikawa, T., and Watanabe, S. (2019). Gradient field of undersea sound speed structure extracted from the GNSS-A oceanography. Mar. Geophys. Res. 40 (4), 493–504. doi:10.1007/s11001-018-9362-7
Yokota, Y., Ishikawa, T., Watanabe, S., Tashiro, T., and Asada, A. (2016a). Seafloor geodetic constraints on interplate coupling of the Nankai Trough megathrust zone. Nature 534, 374–377. doi:10.1038/nature17632
Yokota, Y., Tashiro, T., Ishikawa, T., and Watanabe, S. (2016b). Development of numerical simulator for seafloor geodetic observation system [in Japanese]. Rep. Hydrogr. Oceanogr. Researches 53, 90–97.
Keywords: seafloor geodesy, simulator, transponder array, survey line, GNSS-A
Citation: Nakamura Y, Yokota Y, Ishikawa T and Watanabe S-i (2021) Optimal Transponder Array and Survey Line Configurations for GNSS-A Observation Evaluated by Numerical Simulation. Front. Earth Sci. 9:600993. doi: 10.3389/feart.2021.600993
Received: 31 August 2020; Accepted: 14 January 2021;
Published: 22 February 2021.
Edited by:Keiichi Tadokoro, Nagoya University, Japan
Reviewed by:Fumiaki Tomita, Japan Agency for Marine-Earth Science and Technology (JAMSTEC), Japan
Pierre Sakic, German Research Center for Geosciences, Helmholtz Center Potsdam, Germany
Copyright © 2021 Nakamura, Yokota, Ishikawa and Watanabe. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yuto Nakamura, email@example.com