Time-Aware Validation of Machine Learning Fuel Consumption Models: Evidence from 1,Hz Operational Data, CCGS \textit{Sir Wilfrid Laurier}
Organizations: National Research Council Canada, Ocean, Coastal and River Engineering, P.O. Box 12093, St. John’s, A1B 3T5, NL, Canada · Department of Computer Science, Memorial University of Newfoundland, 40 Arctic Avenue, St. John’s, A1B 3X5, NL, Canada
Abstract
Ship fuel consumption (SFC) prediction supports vessel operation optimisation, emissions estimation, and decision support systems (DSS) for sustainable maritime transportation. Numerous data-driven fuel models have been developed over the past two decades, but a critical and often overlooked limitation lies in their validation practices: most studies evaluate performance using random train--test splits, which, applied to high-frequency records, admit temporal leakage and yield optimistic results that do not reflect deployment conditions. This paper examines that gap using time-aware evaluation, specifically Time Series Cross-Validation (TSCV) and Blocked TSCV (BTSCV). Using the Canadian Coast Guard Ship (CCGS) \textit{Sir Wilfrid Laurier} as a case study, six regression models and a physics baseline are tuned under three time-aware schemes and three feature configurations, then evaluated on a common chronological hold-out set drawn from approximately 3.88 million steady-state 1,Hz records.
Figures & tables
| Literature | Ship Type | Data Source (Size) | Method(s) | Validation | Observations |
| Coraddu et al. (2017) | Handymax chemical/product tanker | Onboard 15 s data aggregated to 15-min intervals; 100,000 data points | RLS, Lasso, RF | Random – 30 random train–test splits across varying data sizes; Bayesian optimization for hyperparameter tuning | Compared WBM, BBM, and GBM for trim optimization. Used 30 random splits with Bayesian optimization for tuning. |
| Bui-Duy et al. (2020) | Container ship | Operational + secondary (wind, TEU); 10,853 records over 2 years | Deep learning neural network (DLNN) | Random – 80:20 train–test split; 7:1 internal validation | Model used in ATSP framework for route optimization; MAPE = 5.89% |
| Kim et al. (2021) | Container ship (13,000 TEU) | 1-min AMS data (6 months); size not specified | ANN, MLR, SVM, Ensemble | Random + voyage-level – 70:30 random split + voyage-level evaluation | Evaluated both randomly and voyage-wise; emphasis on real-world deployment performance |
| Agand et al. (2022) | Passenger ferry | 1-min sensor data (2019–2021); 1M records | MLR, DT, RF, ANN, XGBoost | Random – 70:30 random split + 10-fold CV | XGBoost performed best; consistent route; 34 features used |
| Xie et al. (2023) | Oil tanker (trial) | 1 Hz sensor data; 378,468 records (4.38 days), 496 features | XGBoost, RF, BP-ANN, MLR, polynomial regression | Random – 10 random 70:30 splits, 10-fold CV | Applied Kwon method for acceleration/deceleration filtering; reached 0.9977 (XGBoost); robust to sensor noise and data quality issues |
| Hu et al. (2023) | Bulk carrier | Flowmeter, navigation, meteorological sensors (300 s); 24,300 points | Stacking (hybrid ML), trim optimization | Random – random 80:20 split (5 seeds) | Two-voyage real-world case study; hybrid stacking model outperformed all single models; trim optimization reduced fuel/emissions by up to 1.82% |
| Category | Column Name | Description | Unit | Source |
| Fuel Flow | ME1_Flow | Fuel flow rate, Main Engine 1 | L/h | OpDAQ |
| ME2_Flow | Fuel flow rate, Main Engine 2 | L/h | OpDAQ | |
| ME3_Flow | Fuel flow rate, Main Engine 3 | L/h | OpDAQ | |
| Aux_Gen_Flow | Fuel flow rate, auxiliary generator | L/h | OpDAQ | |
| Speed | OpDAQ GPS SOG | Speed over ground (SOG) | kn | VTG |
| VHW Speed Relative to Water Kn | Speed through water (STW) | kn | VHW |
| Parameter | Variable Code | Description | Unit |
| Significant wave height | swh | Average height of the highest one-third of ocean waves | m |
| Mean wave direction | mwd | Mean direction from which the combined wind-sea and swell wave field is propagating | degrees |
| Mean wave period | mwp | Mean period of the wave spectrum | s |
| Peak wave period | pp1d | Period corresponding to the maximum of the wave spectrum | s |
| Significant swell height | shts | Average height of the swell component of waves | m |
| Mean swell direction | mdts | Direction from which the swell component is propagating | degrees |
| Feature | Derivation | Unit |
| Total_ME_Flow | ME1_Flow + ME2_Flow + ME3_Flow | L/h |
| Relative_wave_direction | ( mwd Heading + 360) mod 360 | degrees |
| Relative_swell_direction | ( mdts Heading + 360) mod 360 | degrees |
| Wind_Direction_Class | mod 8, MWV Wind Angle | – |
| Wave_Direction_Class | mod 8, Relative_wave_direction | – |
| Swell_Direction_Class | mod 8, Relative_swell_direction | – |
| Variant | Speed Channel | Purpose |
| A | Speed over Ground (SOG) | Available at the voyage-planning stage; the configuration evaluated for deployment. |
| B | Speed through Water (STW) | Relates fuel consumption directly to hull resistance, independent of surface current and drift. |
| C | SOG + STW | Tests whether combining both signals improves accuracy beyond either alone, using their difference as an implicit proxy for current and drift. |
| Model | Hyperparameter | Search Grid | Search Iterations |
| Ridge | Regularization strength ( ) | 30-point logarithmic grid, to | 30 a |
| Lasso | Regularization strength ( ) | 20-point logarithmic grid, to | 20 a |
| ElasticNet | Regularization strength ( ) | 15-point logarithmic grid, to | 30 |
| Mixing parameter ( ratio) | 9-point linear grid, 0.1 to 0.9 | ||
| Random Forest | Number of trees | 50, 100, 150, 200 | 20 |
| Maximum tree depth | 10, 20, 30 |
| Model Pair | Feature Variant | Cross-Validation Scheme | Model with Lower Error | Mean Squared-Error Reduction | Blocks with Lower Error (of 100) | Wilcoxon Signed-Rank -value |
| Ridge vs XGBoost | A | Sequential K-Fold | Ridge | < 0.0001 | ||
| A | Custom TSCV | Ridge | ||||
| A | Blocked TSCV | Ridge | < 0.0001 | |||
| B | Sequential K-Fold | Ridge | ||||
| B | Custom TSCV | Ridge | ||||
| B | Blocked TSCV | Ridge |
Appendix figures & tables5 assets
Supplementary material from the paper’s appendix.
Appendix
| Model | Train | Train RMSE | Train MAE | Train Accuracy (%) | Train PA15 (%) | Test | Test RMSE | Test MAE | Test Accuracy (%) | Test PA15 (%) |
| MLR | ||||||||||
| Ridge | ||||||||||
| Lasso | ||||||||||
| ElasticNet | ||||||||||
| Random Forest | ||||||||||
| XGBoost |
| Model | Train | Train RMSE | Train MAE | Train Accuracy (%) | Train PA15 (%) | Test | Test RMSE | Test MAE | Test Accuracy (%) | Test PA15 (%) |
| MLR | ||||||||||
| Ridge | ||||||||||
| Lasso | ||||||||||
| ElasticNet | ||||||||||
| Random Forest | ||||||||||
| XGBoost |
| Model | Train | Train RMSE | Train MAE | Train Accuracy (%) | Train PA15 (%) | Test | Test RMSE | Test MAE | Test Accuracy (%) | Test PA15 (%) |
| Physics | ||||||||||
| MLR | ||||||||||
| Ridge | ||||||||||
| Lasso | ||||||||||
| ElasticNet | ||||||||||
| Random Forest |
| Variant | Model | CV Scheme | Train | Train RMSE | Train MAE | Val | Val RMSE | Val MAE |
| A | Physics | Sequential K-Fold | ||||||
| Custom TSCV | ||||||||
| Blocked TSCV | ||||||||
| Ridge | Sequential K-Fold | |||||||
| Custom TSCV | ||||||||
| Blocked TSCV |
| Variant | Model | CV Scheme | Train | Train RMSE | Train MAE | Train Accuracy (%) | Train PA15 (%) | Test | Test RMSE | Test MAE | Test Accuracy (%) | Test PA15 (%) |
| A | Physics | – | ||||||||||
| MLR | – | |||||||||||
| Ridge | Sequential K-Fold | |||||||||||
| Custom TSCV | ||||||||||||
| Blocked TSCV | ||||||||||||
| Lasso | Sequential K-Fold |