Compressed delayed-information projection for six-degree-of-freedom underwater vehicle navigation under delayed acoustic positioning
Authors: Shuyue Li, Miguel López-Benítez, Eng Gee Lim, Fei Ma, Qian Dong, Mengze Cao, Limin Yu, Xiaohui Qin
Organizations: Academy of Artificial Intelligence and Advanced Technology, Xi’an Jiaotong-Liverpool University, No. 111, Ren’ai Road, Suzhou Industrial Park, Suzhou, Jiangsu Province 215123, China · Department of Electrical Engineering and Electronics, School of Engineering, University of Liverpool, Brownlow Hill, Liverpool L69 3GJ, United Kingdom · Jiangsu JITRI Tsingunited Intelligent Control Technology Co., Ltd., D1 Building, No. 999 Gaolang East Road, Wuxi, Jiangsu Province, China · School of Computer Science and Informatics, University of Liverpool, Ashton Building, Liverpool L69 3DR, United Kingdom · ARIES Research Centre, Polytechnic School, Universidad Antonio de Nebrija, C/ Pirineos, 55, 28040 Madrid, Spain · School of Mathematics and Physics, Xi’an Jiaotong-Liverpool University, No. 111, Ren’ai Road, Suzhou Industrial Park, Suzhou, Jiangsu Province 215123, China · College of Mechanical and Vehicle Engineering, Hunan University, Lushan Road (S), Yuelu District, Changsha, Hunan Province 410082, China
Delayed acoustic positioning packets constrain historical navigation states, but a current-time update evaluates them against a mismatched state, whereas exact rewind/replay re-executes the intervening estimator history. This paper introduces compressed delayed-information projection (CDIP), a causal 15-state error-state Kalman filter (ESKF) treatment for delayed-acoustic unmanned underwater vehicle (UUV) navigation. CDIP retains a source-epoch snapshot and the historical-to-current cross-covariance, then projects the delayed source-epoch acoustic correction directly to the current state without full rewind/replay. Exact fixed-lag rewind/replay out-of-sequence-measurement (OOSM) processing serves as a high-fidelity accuracy reference. In 154 usable paired recordings at a fixed 1.5-s acoustic delay without an outage, CDIP reduced mean trajectory-position root-mean-square error (RMSE) from 1.062 m for the baseline to 0.456 m (57.1%). Its 0.456-m mean was 1.03% higher than the 0.451-m replay mean, while its measured mean per-update runtime was 99.2% lower (approximately 127-fold). A separate predeclared sweep across six fixed delays, with 30 paired recordings per delay, and a truth-supported 9-D consistency analysis bound the interpretation. Additional targeted experiments showed near-replay trajectory accuracy across 50-300-s acoustic outages while preserving sub-millisecond update cost. CDIP therefore provides a compact delayed-information treatment with an empirical accuracy-computation trade-off under the evaluated configuration; the evidence does not establish statistical equivalence or non-inferiority relative to replay.
Figures & tables
Figure 1: Causal sensing and delayed-acoustic context. IMU, DVL, and depth observations are locally available to the causal ESKF. An acoustic position packet generated at source epoch j becomes available at reception epoch k and therefore constrains the historical state rather than the current state.
Figure 2: Conceptual contrast between the two historical treatments. (a) M2 (CDIP) retains a source snapshot and the historical-to-current cross-covariance needed to project the source-epoch correction directly to the current state. (b) M3 (Exact rewind/replay) restores a checkpoint before the source epoch, applies the historical acoustic update, and replays the intervening causal event sequence. The replay path is an accuracy-reference contrast, not an equivalence claim.
Method
Traj. RMSE (m)
Final err. (m)
Peak err. (m)
Vel. RMSE (m/s)
Runtime (ms/update)
M1 (Baseline)
1.062±0.421
1.469±0.625
1.825±0.682
0.095±0.035
0.0273
M2 (CDIP)
0.456±0.171
0.346±0.220
1.116±0.401
0.057±0.020
0.0743
M3 (Exact rewind/replay)
0.451±0.171
0.347±0.218
1.088±0.371
0.057±0.019
9.469
Table 1: Final-main descriptive performance at a fixed 1.5-s acoustic delay without an outage. Accuracy and error metrics are mean ± standard deviation; runtime/update is reported as a mean only. N=154 paired recordings per method; all failure counts were zero.
Comparison
Metric
Mean difference (95% CI)
p ( t /Wilcoxon)
Cohen’s dz
M1–M2
Traj. RMSE (m)
0.606[0.549,0.664]
<0.001/<0.001
1.68
M1–M2
Final err. (m)
1.123[1.025,1.222]
<0.001/<0.001
1.81
M1–M2
Peak err. (m)
0.709[0.616,0.802]
<0.001/<0.001
1.21
M1–M2
Vel. RMSE (m/s)
0.038[0.033,0.042]
<0.001/<0.001
1.34
M2–M3
Traj. RMSE (m)
0.005[−0.001,0.010]
0.0912/0.0187
0.137
M2–M3
Final err. (m)
−0.001[−0.005,0.003]
0.615/0.906
−0.041
Table 2: Paired final-main comparisons ( N=154 ). For M1 minus M2, positive differences favor CDIP. For M2 minus M3, negative differences favor CDIP. Trajectory RMSE is the primary performance metric in the manuscript presentation. The paired t -test and Wilcoxon signed-rank p values are the saved unadjusted values from the frozen final-main analysis; no multiplicity-adjusted p values are claimed, and secondary endpoint tests are interpreted descriptively. No equivalence or non-inferiority inference is made.
Figure 3: Final-main trajectory-position RMSE at the fixed 1.5-s acoustic delay without an outage. M1 denotes the baseline, M2 CDIP, and M3 exact rewind/replay. Open markers show all N=154 paired recordings per method; boxes show the interquartile range with median bars, whiskers extend to the non-outlier range, and filled diamonds denote means.
Figure 4: Final-main accuracy–runtime display at the fixed 1.5-s acoustic delay without an outage. Points are saved means and horizontal and vertical error bars are the saved 95% confidence intervals of those means. The horizontal axis is logarithmic. This is an implementation-specific accuracy–computation comparison, not a significance or equivalence display.
Method
Mean 9-D ANEES
Below / inside / above
Mean DVL NIS
Mean acoustic NIS
M1 (Baseline)
18.110
0.0%/11.7%/88.3%
2.995
2.084
M2 (CDIP)
8.728
50.0%/40.0%/10.0%
2.990
1.993
M3 (Exact rewind/replay)
8.691
50.0%/40.0%/10.0%
2.990
1.992
Table 3: Truth-supported 9-D consistency and innovation diagnostics for the final main experiment. Fractions denote the 60-point aggregate time grid below, inside, and above the 95% ANEES bounds [8.342,9.682] . The expected NIS dimensions are three for DVL and two for acoustic position.
Figure 5: Truth-supported 9-D consistency summary for the final-main experiment. M1 denotes the baseline, M2 CDIP, and M3 exact rewind/replay. (a) Saved mean 9-D ANEES values relative to the saved 95% expected range. (b) Fractions of the 60-point aggregate time grid below, inside, and above that range. This diagnostic complements the trajectory-error results and is not a 15-D ANEES analysis or a performance endpoint.
Delay (s)
M1 RMSE
M2 RMSE
M3 RMSE
M2–M3 difference (95% CI)
M2 runtime
M3 runtime
0
0.455
0.455
0.455
—
0.0307
9.565
0.5
0.535
0.457
0.462
−0.005[−0.015,0.005]
0.0416
9.207
1.5
1.053
0.469
0.477
−0.008[−0.023,0.006]
0.0877
10.774
3
1.934
0.511
0.501
0.010[−0.007,0.027]
0.1141
7.980
5
2.948
0.531
0.537
−0.006[−0.035,0.023]
0.2409
7.664
10
4.354
0.674
0.631
0.043[−0.008,0.094]
0.6525
6.834
Table 4: Separate fixed-delay sweep ( N=30 paired recordings per delay). M1 denotes the baseline, M2 CDIP, and M3 exact rewind/replay. Entries are mean trajectory RMSE in metres and mean runtime in ms/update. The M2 minus M3 RMSE difference is shown with its 95% confidence interval for nonzero delays; negative values favor CDIP. All such intervals include zero.
Figure 6: Separate fixed-delay sweep ( N=30 paired recordings at each predeclared delay). M1 denotes the baseline, M2 CDIP, and M3 exact rewind/replay. (a) Mean trajectory-position RMSE with saved 95% confidence intervals. (b) Mean runtime/update with saved 95% confidence intervals; the vertical axis is logarithmic.
Outage (s)
M1 Baseline ESKF
M2 Proposed CDIP
M3 Exact rewind/replay
CDIP reduction
50
1.816
0.412
0.421
77.3%
100
1.825
0.412
0.421
77.4%
200
1.826
0.435
0.443
76.2%
300
1.920
0.470
0.478
75.5%
Table 5: Targeted prolonged-outage results at a fixed 1.5-s acoustic delay. Entries are mean trajectory-position RMSE in metres over five paired scenarios per duration. CDIP reduction is relative to M1 Baseline ESKF.
Figure 7: Targeted prolonged-outage study at a fixed 1.5-s acoustic delay ( n=5 paired scenarios per duration). Points denote mean trajectory-position RMSE, and bars show the saved 95% confidence intervals. Packets generated during each outage were permanently lost, with no batch recovery.
In Global Navigation Satellite System (GNSS)-denied underwater environments, individual unmanned underwater vehicles (UUVs) suffer from unbounded dead-reckoning drift, making collaborative navigation (CN) crucial for accurate state estimation. However, the severe communication delay inherent in underwater acoustic channels poses serious challenges to real-time state estimation. Traditional filters, such as Extended Kalman Filters (EKFs) or Unscented Kalman Filters (UKFs), usually block the main control loop while waiting for delayed data, or effectively discard Out-of-Sequence Measurements (OOSMs), resulting in serious drift. To address this, we propose an Asynchronous Two-Speed Kalman Filter (TSKF) enhanced by a novel projection mechanism, which we term Variational History Distillation (VHD). The proposed architecture decouples the estimation process into two parallel threads: a fast-rate thread that utilizes Gaussian Process (GP) compensated dead reckoning to guarantee high-frequency real-time control, and a slow-rate thread dedicated to processing asynchronously delayed collaborative information. By introducing a Finite-Length Circular State Buffer (FLCSB), the algorithm applies delayed measurements to their corresponding historical states, and utilizes a VHD-based projection to fast-forward the correction to the current time without computationally heavy recalculations. Simulation results demonstrate that the proposed TSKF maintains a trajectory error comparable to computationally intensive batch-optimization methods under severe delays (up to 30,s). Executing in sub-millisecond time, it significantly outperforms standard EKF/UKF. The results demonstrate an effective control, communication, and computing (3C) co-design that significantly enhances the resilience of autonomous marine automation systems.
Shuyue Li, Miguel López-Benítez, Eng Gee Lim +5
School of Advanced Technology, Xi’an Jiaotong-Liverpool University, Suzhou, China · School of Computer Science and Informatics, University of Liverpool, Liverpool, UK · ARIES Research Centre, Universidad Antonio de Nebrija, Madrid, Spain +3
We present a safety-oriented framework for autonomous underwater vehicles (AUVs) that improves localization accuracy, enhances trajectory prediction, and supports efficient search operations during communication loss. Acoustic signals emitted by the AUV are detected by a network of fixed buoys, which compute Time-Difference-of-Arrival (TDOA) range-difference measurements serving as position observations. These observations are subsequently fused with a Kalman-based prediction model to obtain continuous, noise-robust state estimates. The combined method achieves significantly better localization precision and trajectory stability than TDOA-only baselines. Beyond real-time tracking, our framework offers targeted search-and-recovery capability by predicting post-disconnection motion and explicitly modeling uncertainty growth. The search module differentiates between continued navigation and propulsion failure, allowing search resources to be deployed toward the most probable recovery region. Our framework fuses multi-buoy acoustic data with Kalman filtering and uncertainty propagation to maintain navigation accuracy and yield robust search-region definitions during communication loss.
Zizhan Tang, Yao Liu, Jessica Liu
Department of Electrical and Computer Engineering, University of Southern California, Los Angeles, CA, USA · Department of Computer and Information Science, University of Pennsylvania, Philadelphia, PA, USA
The reliable operation of Unmanned Underwater Vehicle (UUV) clusters is highly dependent on continuous acoustic communication. However, this communication method is highly susceptible to intermittent interruptions. When communication outages occur, standard state estimators such as the Unscented Kalman Filter (UKF) will be forced to make open-loop predictions. If the environment contains unmodeled dynamic factors, such as unknown ocean currents, this estimation error will grow rapidly, which may eventually lead to mission failure. To address this critical issue, this paper proposes a Variational History Distillation (VHD) approach. VHD regards trajectory prediction as an approximate Bayesian reasoning process, which links a standard motion model based on physics with a pattern extracted directly from the past trajectory of the UUV. This is achieved by synthesizing ``virtual measurements'' distilled from historical trajectories. Recognizing that the reliability of extrapolated historical trends degrades over extended prediction horizons, an adaptive confidence mechanism is introduced. This mechanism allows the filter to gradually reduce the trust of virtual measurements as the communication outage time is extended. Extensive Monte Carlo simulations in a high-fidelity environment demonstrate that the proposed method achieves a 91% reduction in prediction Root Mean Square Error (RMSE), reducing the error from approximately 170 m to 15 m during a 40-second communication outage. These results demonstrate that VHD can maintain robust state estimation performance even under complete communication loss.
Shuyue Li, Miguel López-Benítez, Eng Gee Lim +5
School of Advanced Technology, Xi’an Jiaotong-Liverpool University, Suzhou, China · School of Engineering, University of Liverpool, Liverpool, UK · Jiangsu JITRI Tsingunited Intelligent Control Technology Co., Ltd., Wuxi, China +4