# Relationship Analysis Between Helicopter Gearbox Bearing Condition Indicators and Oil Temperature Through Dynamic ARDL and Wavelet Coherence Techniques

Lotfi Saidi, Eric Bechhofer, Mohamed Benbouzid

► To cite this version:

Lotfi Saidi, Eric Bechhofer, Mohamed Benbouzid. Relationship Analysis Between Helicopter Gearbox Bearing Condition Indicators and Oil Temperature Through Dynamic ARDL and Wavelet Coherence Techniques. Machines, 2025, 13 (8), pp.645. 10.3390/machines13080645. hal-05379583

HAL Id: hal-05379583

https://hal.univ-brest.fr/hal-05379583v1

Submitted on 24 Nov 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Relationship Analysis Between Helicopter Gearbox Bearing Condition Indicators and Oil Temperature Through Dynamic ARDL and Wavelet Coherence Techniques

Lotfi Saidi <fnref n="1" />,<fnref n="2" />,<fnref n="3" /> ID, Eric Bechhofer <fnref n="4" /> ID and Mohamed Benbouzid <fnref n="3" />,* ID

## Abstract

This study investigates the dynamic relationship between bearing gearbox condition indicators (BGCIs) and the lubrication oil temperature within the framework of health and usage monitoring system (HUMS) applications. Using the dynamic autoregressive distributed lag (DARDL) simulation model, we quantified both the short- and long-term responses of condition indicators to shocks in oil temperature, offering a robust framework for a counterfactual analysis. To complement the time-domain perspective, we applied a wavelet coherence analysis (WCA) to explore time-frequency co-movements and phase relationships between the condition indicators under varying operational regimes. The DARDL results revealed that the ball energy, cage energy, and inner and outer race indicators significantly increased in response to the oil temperature in the long run. The WCA results further confirmed the positive association between oil temperature and the condition indicators under examination, aligning with the DARDL estimations. The DARDL model revealed that the ball energy and the inner race energy have statistically significant long-term effects on the oil temperature, with *p*-values < 0.01. The adjusted *R*<fnref n="2" /> of 0.785 and the root mean square error (MSE) of 0.008 confirm the model's robustness. The wavelet coherence analysis showed strong time-frequency correlations, especially in the 8–16 scale range, while the frequency-domain causality (FDC) tests confirmed a bidirectional influence between the oil temperature and several condition indicators. The FDC analysis showed that the oil temperature significantly affected the BGCIs, with evidence of feedback effects, suggesting a mutual dependency. These findings contribute to the advancement of predictive maintenance frameworks in HUMSs by providing practical insights for enhancing system reliability and optimizing maintenance schedules. The integration of dynamic economic approaches demonstrates a robust methodology for monitoring critical mechanical components and encourages further research in broader aerospace and industrial contexts.

**Keywords:** bearing gearbox condition indicator; dynamic ARDL; frequency-domain causality; HUMS; lubrication oil temperature; wavelet coherence analysis

## 1. Introduction

Helicopter gearbox systems are critical mechanical components that operate under high stress and varying operational conditions, playing a central role in power transmission.

<figure><img src="image_2.png" /></figure>

check for
updates

Academic Editors: Qing Li, Yu Shao
and Xiang Chen

Received: 18 June 2025

Revised: 20 July 2025

Accepted: 21 July 2025

Published: 24 July 2025

**Citation:** Saidi, L.; Bechhofer, E.; Benbouzid, M. Relationship Analysis Between Helicopter Gearbox Bearing Condition Indicators and Oil Temperature Through Dynamic ARDL and Wavelet Coherence Techniques. *Machines* **2025**, *13*, 645. https://doi.org/10.3390/machines13080645

**Copyright:** © 2025 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

<footnote n="1">1 École Nationale Supérieure d'Ingénieurs de Tunis (ENSIT), University of Tunis, Labo SIME, Av Taha Hussein, Tunis 1008, Tunisia; lotfi.saidi@esstt.rnu.tn</footnote>

<footnote n="2">2 Ecole Supérieure des Sciences et de la Technologie de Hammam Sousse (ESSTHS), University of Sousse, Sousse 4011, Tunisia</footnote>

<footnote n="3">3 Institut de Recherche Dupuy de Lôme (UMR CNRS 6027 IRDL), University of Brest, 29238 Brest, France</footnote>

<footnote n="4">4 Green Power Monitoring Systems, LLC, 93 Pilgrim Park Rd, Waterbury, VT 05676, USA; eric@gpms-vt.com</footnote>

<footnote>* Correspondence: mohamed.benbouzid@univ-brest.fr</footnote>
Due to the demanding nature of their function, ensuring the health and reliability of these gearboxes is vital, especially during critical flight phases such as takeoff, landing, and adverse weather conditions [1–5]. The gearbox performance directly influences the flight safety, mission success, and operational costs. In this context, advanced predictive maintenance systems are increasingly replacing traditional maintenance strategies, which often rely on fixed schedules or failure-based models. These systems are integrated into HUMSs, which continuously monitor the status of critical components in real time, enabling more effective and timely maintenance decisions [4,6–8].

The analysis of the relationship between BGCIs and the oil temperature (OT) is of critical importance for enhancing the health monitoring, diagnostic precision, and predictive maintenance of helicopter gearbox systems [9,10]. While the OT serves as a thermal health indicator—signaling issues such as overheating, poor lubrication, or internal friction—BGCIs reflect mechanical degradation, particularly of components like the bearing inner race, outer race, cage, and rolling elements [11,12]. Studying these two condition indicators (CIs) in conjunction provides complementary insights into both thermal and mechanical health, allowing for early fault detection, degradation pattern recognition, and a more accurate estimation of the remaining useful life (RUL) [5–12]. Moreover, understanding the dynamic interaction between BGCIs and the OT is essential for developing robust prognostic models, especially in real-world environments where systems exhibit nonlinear and non-stationary behavior. This paper makes several original contributions to the field of gearbox diagnostics and health monitoring. While the individual analytical techniques employed—such as the DARDL model, the WCA, and FDC testing—are established in the literature, their combined application to real helicopter HUMS data for analyzing the interplay between the OT and BGCIs is, to our knowledge, novel. Specifically, this study is among the first to apply probabilistic DARDL simulations to model both positive and negative shocks in CIs, providing a more realistic and interpretable forecast of system behavior under stress. Moreover, the integration of both time-domain (DARDL) and frequency-domain frameworks enables a deeper understanding of how mechanical and thermal factors influence each other over multiple time scales. This dual perspective enhances the diagnostic sensitivity and system-level insights, extending the practical utility of condition-monitoring models in aerospace HUMS platforms.

Compared to traditional condition-monitoring approaches—which often rely on manual inspection, high-volume continuous data acquisition, or reactive diagnostics—HUMSs offer a more integrated and efficient alternative. Systems like GPMS enable automated, periodic in-flight monitoring that reduces data overhead while maintaining diagnostic sensitivity. By embedding sensors and analytics onboard, HUMSs facilitate real-time health tracking, early fault detection, and proactive maintenance planning, in contrast to the offline analysis or post-failure investigations typical of conventional systems.

One key challenge in gearbox health monitoring is the integration of multiple performance indicators that reflect different aspects of the system condition [9,12,13]. Among these, two critical parameters are the lubrication OT and BGCIs, both of which provide essential diagnostic information [14]. An elevated OT can indicate potential issues such as a lack of lubrication, overheating, or mechanical faults. BGCIs, on the other hand, directly reflect bearing wear and degradation, which are among the most failure-prone components in a gearbox. Understanding the relationship between the OT and BGCIs is therefore crucial for maintaining optimal gearbox performance and ensuring safe, efficient operations.

Although earlier studies have explored the individual effects of various CIs on the OT—such as the outer race energy (*OR*<sub>E</sub>), inner race energy (*IR*<sub>E</sub>), spectral kurtosis (*SK*), cage energy (*C*<sub>E</sub>), and ball energy (*B*<sub>E</sub>)—many treat these CIs in isolation. For instance, prior research has employed ARDL models to quantify the impact of condition indicators on the
OT trends in Bell 407 helicopters [9,13]. However, such approaches often assume linearity and stationarity and overlook the interaction effects across time scales. Recent work has proposed the use of wavelet-based methods, such as the autocorrelation of Morlet wavelet transforms [15] and continuous wavelet transforms [16], for early fault detection under non-stationary conditions. Additionally, a machine learning-assisted wavelet analysis has demonstrated strong potential for improving the health state classification in complex mechanical systems [17–19].

Furthermore, studies in data-driven HUMS strategies have shown how integrating me-
chanical and thermal indicators—through dynamic modeling and spectral causality—enhances
the prognostic accuracy in rotating systems. These include the application of spectral
entropy, signal decomposition, and hybrid signal-processing techniques for gearbox diag-
nostics under real flight conditions [10,12,20].

Despite these advances, there remains a gap in the literature concerning the inte-
grated use of dynamic econometric models (e.g., DARDL), a wavelet coherence analysis
(WCA), and frequency-domain Granger causality to assess the coupled dynamics be-
tween the OT and BGCIs. Our work addresses this gap by combining these methods into
a unified framework for analyzing the temporal, spectral, and causal dependencies in
helicopter gearboxes.

Recent developments in fault detection frameworks further emphasize combin-
ing wavelet methods with advanced analytics such as adaptive filtering and deep
learning [21,22]. These studies highlight the growing integration of data-driven and hybrid
approaches to improve the early detection of complex mechanical failures.

Recent advancements in intelligent monitoring systems have also explored the use of federated learning for distributed fault diagnoses. For example, adaptive federated approaches have been proposed to address challenges like label heterogeneity and communication redundancy in collaborative industrial environments [23,24]. Additionally, client-level imbalances have been tackled to enhance model generalization across diverse operational conditions [25]. These techniques support the growing trend toward decentralized, privacy-preserving fault diagnosis frameworks, complementing the predictive analytics used in HUMS applications.

Building on these analytical advancements, a broader bibliometric perspective reveals how such methodologies align with the prevailing research trends in vertical flight HUMSs. Specifically, the integration of CIs like the OT and BGCIs within advanced diagnostic frameworks corresponds with the growing emphasis on data-driven techniques such as fault diagnoses, RUL estimation, and predictive maintenance, as illustrated in the keyword co-occurrence and thematic clusters shown in Figure 1.

Figure 1 provides a bibliometric overview supporting the relevance of this study’s methodology. Subfigure (a) groups research themes based on keyword co-occurrence, showing a dense cluster around “condition-based maintenance” and related terms. Subfigure (b) highlights the global distribution of research, with darker regions representing a higher level of publication activity. Subfigure (c) identifies the leading journals contributing to the HUMS literature, such as Reliability Engineering & System Safety and IEEE Access, while (d) maps the most cited authors, with the node size reflecting the citation count. The added legends clarify the visual encoding used in each network diagram.

The central theme, as evidenced by the prominence of the “condition-based mainte-
nance” node, underscores its foundational role in the field. Closely associated concepts
include “fault-diagnosis”, “deep learning”, “remaining useful life estimation”, and “prog-
nostics and health management” (PHM), indicating a strong research focus on advanced
diagnostic and prognostic methodologies. “Vibration monitoring” and “bearing” were also
significant, highlighting the importance of specific condition-monitoring techniques and
critical components in these systems. While "oil debris analysis" was not explicitly a large
node, its relevance is implied through the connections to "bearing" and "diagnosis", which
are integral to a comprehensive condition assessment. The field extensively leverages "artifi-
cial intelligence" and "neural networks" for developing robust "models" for "degradation"
and "system" analyses, ultimately contributing to "predictive maintenance" strategies.

<figure><img src="image_3.png" /><figcaption>**Figure 1.** HUMS bibliometric analysis results: (a) Keyword clusters, color-coded by thematic areas such as "predictive maintenance", "fault diagnosis", and "deep learning", with node size indicating keyword frequency and proximity representing co-occurrence strength. (b) Country-wise research contributions, where color intensity reflects publication volume. (c) Top publishing journals visualized by node size (citation impact) and color (disciplinary category). (d) Prominent authors, where larger nodes indicate higher citation frequency and edges denote co-authorship or thematic similarity. A legend and explanatory color map are included in each subplot to clarify data representation.</figcaption></figure>

In terms of publication outlets, Reliability Engineering & System Safety emerged as
a highly influential journal, serving as a central hub for research dissemination. Other
notable journals and proceedings included Sensors, IEEE Access, Mathematics, Journal of
Intelligent Manufacturing, Journal of Quality in Maintenance Engineering, Proceedings
of the Institution of Mechanical Engineers, IEEE Transactions on Reliability, and Ocean
Engineering, reflecting the interdisciplinary nature of the research, which encompasses en-
gineering, computer science, and applied sciences. Geographically, research contributions
were predominantly led by the "Peoples Republic of China" and the "USA", demonstrating
their significant impact and investment in this area. Other active countries include Italy,
France, England, Germany, Canada, and India, indicating a global collaborative effort.
Prominent authors such as A.K.S. Jardine, J. Lee, and Z. Li were identified as key contrib-
utors, whose work likely shapes the theoretical and practical advancements in HUMSs,
CBM, and the broader predictive maintenance landscape.
Building on this foundation, this study employs the DARDL model and a WCA to investigate the dynamic relationship between the OT and BGCIs in the Bell 407 helicopter model. The DARDL model, introduced by Jordan and Philips [19], is well-suited for helicopter gearbox data due to its robustness in capturing both short- and long-term dynamics, even with small sample sizes. Additionally, a WCA enables a time-frequency perspective, revealing how the relationship between BGCIs and the OT evolves over time and across various frequency bands. This is particularly valuable for detecting both short-term fluctuations and long-term trends that conventional time-domain models may overlook. Together, these methods provide a robust analytical framework for understanding the interactions that influence gearbox performance and contribute to system reliability and maintenance outcomes.

This study contributes to the literature in several key ways. First, while previous works—such as [9]—have utilized ARDL models to explore the relationships between CIs, this study employs the more advanced DARDL model, which accounts for actual variation in the independent CIs. This approach enables a more accurate assessment of how both positive and negative changes in the OT affect BGCIs and allows these effects to be simulated and visualized across different time horizons.

Second, most existing studies rely on time-domain Granger causality tests, which may miss causal interactions that vary across frequencies. Since the direction and strength of causality between the OT and BGCIs can change depending on the time scale—such as short-term operational shocks versus long-term wear patterns—a purely time-domain approach provides only a partial view of system behavior. To address this, we incorporated the FDC framework developed by Breitung and Candelon [20], which enables a more nuanced examination of causality across time scales, offering deeper insights into the interaction between lubrication performance and gearbox health.

The remaining structure of this paper is as follows: Section 2 provides a comprehensive review of the relevant literature. Section 3 outlines the materials and methodological approach employed in this study. Section 4 reports and analyzes the empirical findings. Lastly, Section 5 summarizes the main conclusions and discusses their implications for predictive maintenance practices, along with suggestions for future research.

## 2. Material and Methods

### 2.1. HUMS Data Collection

The HUMS provides a comprehensive aircraft health assessment by automating flight data monitoring, the rotor track and balance, the engine performance, and drivetrain diagnostics [26]. It uses ten smart accelerometers to gather vibration data, with an accessory drive sensor monitoring the duplex bearing at 4450 RPM. Data are collected at 46.875 samples per second over 2 s every four minutes during straight and level flight.

This 2 s sampling window at 46 Hz, taken every 4 min, reflects a common trade-off in HUMSs between a detailed resolution and a manageable data size for long-duration missions. It is designed to capture vibration under stable flight conditions, minimizing transient interference and ensuring consistent baseline monitoring. This approach is implemented in the GPMS Foresight MX system to optimize the diagnostic reliability while maintaining data efficiency.

Figures 1–3 show the GPMS Foresight MX Bell 407 kit, the variable reluctance sensors (VR sensors) and main rotor [26–28], and the OT and BGCIs for a Bell 407 helicopter, respectively.
<figure><img src="image_4.png" /><figcaption>**Figure 2.** The Foresight MX onboard HUMS system with labeled components: (1) smart accelerometers mounted on gearbox housings; (2) VR (variable reluctance) sensors monitoring the rotational speed; (3) the onboard processing unit that collects, filters, and compresses vibration signals; and (4) the data-transmission unit for remote health reporting. The system is STC-certified for various Bell 407 variants and ensures accurate in-flight condition monitoring [26].</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>**Figure 3.** VR sensors and the main rotor [26].</figcaption></figure>

More information can be found in references [26–28].

Figure 4 illustrates the temporal evolution of key vibration-based indicators ($C_E$, $B_E$, $IR_E$, $OR_E$, and SK) alongside the oil temperature (T). Each subplot captures the sensor output over time during straight and level flight. Notably, $B_E$ and $IR_E$ exhibited multiple peaks that aligned with temperature spikes, suggesting strong thermal–mechanical coupling. $C_E$ and $OR_E$ also showed fluctuations that preceded or followed thermal events. These patterns support the hypothesis that gearbox heat buildup and mechanical wear are interrelated. SK, representing impulsiveness, fluctuated less predictably, but still showed minor co-variation. Understanding these trends is essential for HUMSs: elevated vibration levels paired with a rising OT can provide early warnings of bearing degradation or lubrication failure.
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Time-series plots of bearing gearbox condition indicators (BGCI s) and the oil temperature (OT) in the Bell 407 helicopter: (1) cage energy (CE), (2) ball energy (BE), (3) inner race energy (IRE), (4) outer race energy (ORE), (5) spectral kurtosis (SK), and (6) oil temperature (T). The x-axis represents the acquisition time, while the y-axis shows the normalized amplitude. A visual inspection revealed that spikes or upward shifts in CIs often correspond to increases in the OT, particularly for the BE and IRE. These time-aligned behaviors suggest correlated degradation and thermal stress, reinforcing their utility in early fault detection.</figcaption></figure>

## 2.2. Model Specification

The mathematical formulations used in this study originate from both the established literature and adaptations to suit our empirical framework. Specifically, Equations (1)–(3) are adapted from the ARDL modeling approach described by Pesaran et al. [29]. Equation (4) applies the dynamic ARDL (DARDL) simulation model introduced by Jordan and Philips [19]. Equations (5)–(9) represent standard continuous wavelet and coherence formulations, as developed by Goupillaud, Grossmann, and Morlet [29–33], and later operationalized by Torrence and Compo [34]. Where relevant, we have included direct citations with each formula to ensure clarity of origin.

This study analyzed the dynamic impacts of the $C_E$, $B_E$, $IR_E$, and $OR_E$, and SK on the OT, using the model presented in Equation (1):

$$OT_t = \alpha_0 + \alpha_1 C_{Et} + \alpha_2 B_{Et} + \alpha_3 IR_{Et} + \alpha_4 OR_{Et} + \alpha_5 SK_t + \varepsilon_t \quad (1)$$

where T is the OT and t denotes the acquisition time.

$\alpha_0$ is the constant term. $\alpha_1$, $\alpha_2$, $\alpha_3$, $\alpha_4$, and $\alpha_5$ represent the long-run elasticities of the explanatory CIs. $\varepsilon_t$ is the regression error term.

The econometric model in Equation (1) is not linear, which could lead to inconsistent results and complicate decision-making. To address this, we applied a natural logarithm (Ln) transformation to convert the CIs into comparable units, improving the results and interpretation. The log-linear form is shown in Equation (2):

$$LnOT_t = \alpha_0 + \alpha_1 LnC_{Et} + \alpha_2 LnB_{Et} + \alpha_3 LnIR_{Et} + \alpha_4 LnOR_{Et} + \alpha_5 LnSK_t + \varepsilon_t \quad (2)$$

Equations (1) and (2) define the foundational regression structure between the OT and selected CIs, using a log-linear form to ensure the comparability and interpretability of the coefficients. This formulation sets the base for assessing the long-term elasticity and informs subsequent modeling steps.

As mentioned above, this study addresses the challenges of using the ARDL model for analyzing the short- and long-term effects of independent CIs, including cointegration
testing, lag complexities, and software limitations. To overcome these, a dynamic simulation approach incorporating an error correction term is proposed. This innovative method enhances the traditional ARDL model, enabling a more effective analysis of exogenous CI impacts. Figure 5 outlines this methodology.

<figure><img src="image_7.png" /><figcaption>Figure 5. Empirical analysis roadmap.</figcaption></figure>

## 2.3. Algorithmic Procedure for Experimental Analysis

To improve the clarity and reproducibility of our methodological framework, we outline here the key steps followed in the experimental analysis. The overall procedure integrates DARDL simulations, a WCA, and FDC testing. These steps are illustrated in Figure 5 and can be summarized as follows:

- **Data Acquisition**: Vibration and OT data are collected from HUMS-equipped Bell 407 helicopters under steady-state flight conditions.

- **Preprocessing**: Raw signals are log-transformed and smoothed to ensure the stationarity and comparability of CIs.

- **Unit Root Testing**: The stationarity of the time series is assessed using the ADF and KPSS tests.

- **ARDL Bounds Testing**: Cointegration relationships between the OT and BGCIs are verified.

- **DARDL Simulation**: The DARDL model is applied to simulate the short- and long-term impacts of the BGCIs on the OT.

- **Wavelet coherence analysis**: A WCA is conducted to examine the time-frequency correlations between variables.

- **Frequency-domain causality**: Spectral Granger causality is tested across short-, medium-, and long-term bands.

- **Diagnostics and validation**: Statistical tests (LM, RESET, Jarque-Bera) validate the model robustness.

This structured approach enhances the transparency of our analysis and supports replication in future HUMS-related research.
All variable notations (e.g., frequency ω, lag terms, and log-transformed variables such as LnOT) have been standardized throughout the manuscript to ensure typographic and semantic consistency.

### 3. Empirical Analysis Procedure

To test for the presence of unit roots in the study, we applied the augmented Dickey-Fuller (ADF) test [30] and the Kwiatkowski-Phillips-Schmidt-Shin (KPSS) test [31]. Both tests assess the null hypothesis of a unit root against the alternative hypothesis of stationarity. The DARDL approach requires the dependent CI to be integrated at I(1).

To examine the long-term relationships between CIs, we employed the ARDL coin- tegration methodology proposed by Pesaran et al. [29]. This approach offers several advantages over the other traditional cointegration techniques.

Equation (3) shows the ARDL bounds testing method applied to analyze the long-run relationship between the CIs. This equation presents the ARDL bounds testing framework, specifying both lagged and differenced terms to estimate short- and long-term interactions between the CIs and the OT. This model forms the core of the empirical estimation strategy.

$$ \Delta LnOT_t = \alpha_0 + \sum_{k=1}^{n} \alpha_{1k} \Delta LnOT_{t-k} + \sum_{k=1}^{n} \alpha_{2k} \Delta LnC_{Et-k} + \sum_{k=1}^{n} \alpha_{3k} \Delta LnB_{Et-k} + \sum_{k=1}^{n} \alpha_{4k} \Delta LnIR_{Et-k} + \sum_{k=1}^{n} \alpha_{5k} \Delta LnOR_{Et-k} + \sum_{k=1}^{n} \alpha_{6k} \Delta LnSK_t + \delta_1 LnOT_{t-1} + \delta_2 LnC_{Et-1} + \delta_3 LnB_{Et-1} + \delta_4 LnIR_{Et-1} + \delta_5 LnOR_{Et-1} + \delta_6 LnSK_{t-1} + \varepsilon_t \quad (3) $$

These ARDL formulations follow the structure proposed by Pesaran et al. [29] and are adapted for our multi-indicator context.

Here, Δ represents the first difference operator, while t-1 refers to the optimal lag length determined using the Akaike information criterion (AIC). The parameters α and δ capture the short-run and long-run dynamics, respectively. When a long- run relationship among the CIs is established, the ARDL model is used to estimate both the short-term and long-term coefficients. The H₀ suggests no cointegration (H₀ = δ₁ = δ₂ = δ₃ = δ₄ = δ₅ = δ₆ = 0), while the alternative hypothesis indicates cointegra- tion (H₁ ≠ δ₁ ≠ δ₂ ≠ δ₃ ≠ δ₄ ≠ δ₅ ≠ δ₆ ≠ 0).

This hypothesis is tested using the F-statistic [32]. Cointegration is indicated if the F- value exceeds the upper bound, while no cointegration is present if it falls below the lower bound. If the F-value is between the bounds, the decision is inconclusive. A significant error correction term suggests long-term associations between the CIs.

Jordan and Philips [19] introduced a DARDL simulation approach to simplify the analysis of short- and long-run impacts in ARDL models. This method simulates and plots probabilistic shifts in the dependent CI due to changes in a single independent CI, with other regressors held constant. It requires the data to be integrated of order one and cointegrated. Using 5000 parameter vector simulations from a multivariate normal distribution, the DARDL error correction algorithm ensures a robust analysis. The ARDL bounds test for error correction follows the methodology outlined in [19] and is given in Equation (4). Equation (4) introduces the DARDL simulation, where shocks in individual CIs are simulated to project the probabilistic responses of the OT over time. This tech- nique enhances the ARDL estimation by allowing Monte Carlo simulations and graphical scenario analyses.

Equation (4) implements the DARDL simulation model as outlined by Jordan and Philips [19].

$$ \Delta LnT_t = \alpha_0 + \theta_0 LnT_{t-1} + \beta_1 \Delta LnC_{Et} + \theta_1 LnC_{Et-1} + \beta_2 \Delta LnB_{Et} + \theta_2 LnB_{Et-1} + \beta_3 \Delta LnIR_{Et} + \theta_3 LnIR_{Et-1} + \beta_4 \Delta LnOR_{Et} + \theta_4 LnOR_{Et-1} + \beta_5 \Delta LnSK_t + \theta_5 LnSK_{t-1} + \varepsilon_t \quad (4) $$
This equation illustrates the DARDL simulation model, where $\beta_i$ ($i$ = 1 to 5) represents the long-run coefficients and $\theta_i$ ($i$ = 1 to 5) corresponds to the short-run coefficients. The error correction term (ECT) captures the speed at which the system returns to equilibrium following a disturbance.

To ensure the robustness of the findings, several diagnostic tests were performed. These included the Breusch-Godfrey Lagrange multiplier (LM) test to detect serial correlations, the Cameron and Trivedi decomposition of the Im-Pesaran-Shin (IM) test for heteroscedasticity, skewness and kurtosis tests for normality, and the cumulative sum of recursive residuals (CUSUM) and cumulative sum of recursive residual squares (CUSUMSQ) tests to evaluate the model stability [34].

This investigation employed the wavelet approach to explore the time-frequency dependency among the selected CIs. Specifically, this research identified the temporal and frequency-based interactions between inclusive growth and its explanatory factors. To achieve this, we applied a WCA that included phase differences developed by [33,35]. This method reveals the identification of spectral properties in the time series, especially in a manner in which various periodic characteristics of the time-series change over time. Equation (5) presents the mathematical representation of the Morlet wavelet function:

$$ \omega(t) = \pi^{-\frac{1}{4}}e^{-j\omega t}e^{-\frac{1}{2}t^2} \qquad (5) $$

where $\omega$ represents the frequency used for a limited time series and $j$ is the imaginary unit.

Wavelets are characterized by two fundamental factors: time (denoted by $k$) and frequency (denoted by $f$). The CWT, using time-series data $p(t)$, is given by Equation (6):

$$ \omega_p(k, f) = \int_{-\infty}^{\infty} p(t) \frac{1}{\sqrt{f}} \omega\left(\frac{\overline{t-k}}{f}\right) dt \qquad (6) $$

The wavelet power spectrum (WPS) is used to illustrate the sensitivity of both the dependent and explanatory CIs, as demonstrated in the following expression:

$$ WPS_p(k, f) |W_p(k, f)|^2 \qquad (7) $$

The CWT method was used to convert the time-series CI in Equation (7) into the form presented in Equation (7). In order to explore any relationship between two time series ($p(t)$ and $q(t)$) in the context of time-frequency-based causalities, we utilized the WCA. Equation (8) describes the application of the CWT to the time series:

$$ W_{pq}(k, f) = W_p(k, f) \overline{W_q(k, f)} \qquad (8) $$

The CWT of the OT and the BGCIs are represented by $W_p(k, f)$ and $W_q(k, f)$, respectively. The squared wavelet coherence is presented in Equation (9):

$$ R^2(k, f) = \frac{|S(f^{-1}W_{pq}(k, f))|^2}{S(f^{-1}|W_p(k, f)|^2)S(f^{-1}|W_q(k, f)|^2)} \qquad (9) $$

A value of $R^2(k, f)$ approaching one suggests a strong connection between the time-series components, potentially indicating a causal link between the CIs at a specific frequency. This relationship is typically represented by a black contour outlined in red.

Torrence and Compo [34] introduced a technique to detect discrepancies in wavelet coherence by employing deferral indicators to analyze the fluctuations between two time
series. The corresponding formulation for the wavelet coherence difference is presented in Equation (10) below:

$$ \phi_{pq}(k, f) = \tan^{-1} \left( \frac{L\{S(f^{-1}W_{pj}(k, f))\}}{O\{S(f^{-1}W_{pj}(k, f))\}} \right) \quad (10) $$

Equations (5)–(9) are standard in the wavelet coherence literature and are based on the formulations of Goupillaud et al. [33] and Torrence and Compo [34].

These equations introduce the wavelet coherence methodology, which complements the time-domain analysis by identifying co-movement and directional dependence between the OT and CIs across various frequency bands. This dual-domain approach strengthens the robustness of the experimental framework shown in Figure 5.

We applied the FDC method by Breitung and Candelon [20] to examine the spectral Granger causality between the OT and BGCI values. This approach predicts CI values at specific time frequencies, identifying historical patterns to guide maintenance and operations. Unlike conventional tests, it decomposes causality into short-term (high-frequency) and long-term (low-frequency) components, providing insights into how factors like a fluctuating OT impact the gearbox condition over varying time horizons. For example, it reveals whether short-term temperature changes affect BGCIs differently from long-term thermal trends. The equation from [20] is expressed as follows:

$$ x_t = \alpha_1 X_{t-1} + \dots + \alpha_p X_{t-p} + \beta_1 Y_{t-1} + \dots + \beta_p Y_{t-p} + \varepsilon_t \quad (11) $$

The null hypothesis of Equation (8) is $M_{y→x}(\omega) = 0$. However, $\alpha$ and $\beta$ represent the assessed parameters at time $t$, with a lag $p$, and $\varepsilon_t$ denotes the error term.

# 4. Empirical Results and Discussion

## 4.1. Unit Root Test

Prior to analyzing the interrelationships among the CIs, it is crucial to assess their stationarity. Table 1 presents the outcomes of the unit root tests employed in this study. The KPSS test, which assumes stationarity as the null hypothesis, and the ADF test, which assumes the presence of a unit root, were both applied. The findings revealed that most CIs—specifically, LnC<sub>E</sub>, LnB<sub>E</sub>, LnOR<sub>E</sub>, and LnSK—are stationary at level, while LnIR<sub>E</sub> and LnOT become stationary only after first differencing. Consequently, the CIs exhibited a mixed integration order, with some being I(0) and others I(1), as consistently indicated by the results of both tests.

<table><caption>**Table 1.** Results of ADF and KPSS unit root tests.</caption><thead><tr><th rowspan="2">CIs</th><th colspan="2">ADF</th><th colspan="2">KPSS</th></tr><tr><th>Level</th><th>First Difference</th><th>Level</th><th>First Difference</th></tr></thead><tbody><tr><td>LnOT</td><td>-0.460</td><td>-32.875 ***</td><td>-0.013</td><td>-23.324 ***</td></tr><tr><td>LnC<sub>E</sub></td><td>-3.619 ***</td><td>-32.923 ***</td><td>-9.183 ***</td><td>-32.202 ***</td></tr><tr><td>LnB<sub>E</sub></td><td>-5.725 ***</td><td>-30.881 ***</td><td>-8.252 ***</td><td>-35.617 ***</td></tr><tr><td>LnIR<sub>E</sub></td><td>-0.703</td><td>-35.604 ***</td><td>-1.356</td><td>-31.267 ***</td></tr><tr><td>LnOR<sub>E</sub></td><td>-4.548 ***</td><td>-30.650 ***</td><td>-15.142 ***</td><td>-30.262 ***</td></tr><tr><td>LnSK</td><td>-3.817 ***</td><td>-34.203 ***</td><td>-10.461 ***</td><td>-30.166 ***</td></tr></tbody></table>

Note: *** p < 0.01.

## 4.2. ARDL Cointegration Test

Following the confirmation of stationarity among the CIs, this study applied the ARDL bounds testing approach to assess potential cointegration relationships among the
CIs. The results of the bounds test are summarized in Table 2. Specifically, when LnOT was treated as the dependent variable in the ARDL framework, the analysis revealed evidence of cointegration. This is supported by the fact that the computed F-statistic (5.534) exceeded both the lower and upper critical bounds, indicating a long-run relationship among the variables.

<table><caption>**Table 2.** Results of estimated ARDL model for cointegration.</caption><thead><tr><th>Estimated Model</th><th colspan="2">F-Statistics</th></tr></thead><tbody><tr><td>LnOT<sub>t</sub> = f(LnC<sub>Et</sub>, LnB<sub>Et</sub>,<br/>LnIR<sub>Et</sub>, LnOR<sub>Et</sub>, LnSK<sub>t</sub>)</td><td colspan="2">5.534 *</td></tr><tr><td>Significance level</td><td>Lower bound</td><td>Upper bound</td></tr><tr><td>1%</td><td>3.09</td><td>3.86</td></tr><tr><td>5%</td><td>2.93</td><td>3.83</td></tr><tr><td>10%</td><td>2.101</td><td>3.869</td></tr></tbody></table>

Note: * *p* < 0.05.

### 4.3. *Dynamic Autoregressive Distributed Lag Simulations*

The results of the DARDL simulations are presented in Table 3.

<table><caption>**Table 3.** DARDL simulations results.</caption><thead><tr><th>CIs</th><th>Coefficient</th><th>St. Error</th><th>t-Value</th></tr></thead><tbody><tr><td>ECT</td><td>-0.778</td><td>0.163</td><td>-2.96 ***</td></tr><tr><td>LnOR<sub>E</sub></td><td>0.205</td><td>-5.287</td><td>0.000 ***</td></tr><tr><td>LnOR<sub>Et-1</sub></td><td>0.122</td><td>2.366</td><td>0.000 ***</td></tr><tr><td>LnSK</td><td>0.025</td><td>0.636</td><td>0.250</td></tr><tr><td>LnSK<sub>t-1</sub></td><td>0.020</td><td>0.211</td><td>0.704</td></tr><tr><td>LnIR<sub>E</sub></td><td>0.262</td><td>3.139</td><td>0.000 ***</td></tr><tr><td>LnIR<sub>Et-1</sub></td><td>0.292</td><td>3.766</td><td>0.000 ***</td></tr><tr><td>LnC<sub>E</sub></td><td>0.374</td><td>0.712</td><td>0.010</td></tr><tr><td>LnC<sub>Et-1</sub></td><td>-0.020</td><td>-0.110</td><td>0.028</td></tr><tr><td>LnB<sub>E</sub></td><td>0.120</td><td>4.448</td><td>0.000 ***</td></tr><tr><td>LnB<sub>Et-1</sub></td><td>0.266</td><td>4.883</td><td>0.000 ***</td></tr><tr><td>Cons</td><td>3.281</td><td>0.069</td><td>0.000 ***</td></tr><tr><td>Adj R-squared</td><td>0.7850</td><td>Root MSE</td><td>0.008</td></tr><tr><td>R<sup>2</sup></td><td>0.772</td><td></td><td></td></tr><tr><td>Simulation</td><td>5000</td><td></td><td></td></tr></tbody></table>

Note: *** denote statistical significance at the 10% level.

The DARDL method can simulate, predict, and plot projections of counterfactual shifts in one regressor and their effects on the dependent CI while maintaining the stability of the model. A key advantage of the DARDL simulated model is its ability to provide instant forecasts, robust estimations, and highly accurate predictions. This approach also enables the visualization of actual changes in explanatory and dependent CIs, offering valuable insights through dynamic estimations.

In terms of computational complexity, the DARDL simulation involves 5000 iterations of a Monte Carlo process using a multivariate normal distribution. The estimated time complexity is approximately $O(n \times p \times s)$, where n is the number of regressors, p is the lag order, and s is the number of simulation draws (typically 5000). The WCA has a computational complexity of $O(N \log N)$ for each pair of time series, where N is the length of the signal. These operations remain tractable and suitable for embedded HUMS platforms with appropriate data preprocessing and hardware acceleration.
Graphical illustrations of the DARDL were employed to depict the precise relationships between the explanatory CIs and their impacts on the dependent CI. The DARDL simulations illustrate the forecasted effects of changes in the actual regressors on the LnOT.

Specifically, the impacts of LnC<sub>Et</sub>, LnB<sub>Et</sub>, LnIR<sub>Et</sub>, LnOR<sub>Et</sub>, and LnSK<sub>t</sub> on LnT were predicted under scenarios of a 10% increase and a 10% decrease in these CIs.

It is also important to consider the role of potential sensor defects or measurement errors. Although the GPMS smart accelerometers and VR sensors are STC-certified and undergo regular calibration, they are not immune to noise, drift, or minor malfunctioning. These imperfections may propagate through the signal-processing pipeline and affect the accuracy of BGCI calculations. Since our model uses log-transformed and smoothed data to reduce the noise impact, some robustness is retained. However, we acknowledge that such sensor-induced uncertainty could influence DARDL outputs, especially over long prediction horizons. Future work may include incorporating uncertainty quantification techniques, such as bootstrapping or Bayesian inference, to formally assess how sensor variability affects model confidence.

The graphical results of these simulations are presented in Figures 6–10.

<figure><img src="image_8.png" /><figcaption>**Figure 6.** The impulse response of LnOR<sub>E</sub> to a 10% increase and decrease over time. The x-axis represents time in minutes, and the y-axis shows the percentage deviation in oil temperature (LnT). Blue shaded bands represent 75%, 90%, and 95% confidence intervals. (a) Impact of a 10% positive change in the ORE on the OT in both the short and long term. (b) Effect of a 10% negative change in the ORE on the OT in the short and long term.</figcaption></figure>

Figure 6 illustrates the effect of both positive and negative OR<sub>E</sub> variations on the OT. Figure 6a shows the impact of a 10% positive change in the OR<sub>E</sub> on the OT in both the short and long term. Figure 6b depicts the effect of a 10% negative change in the OR<sub>E</sub> on the OT in the short and long term. Figure 7 presents the impulse response plot illustrating the nexus between the IR<sub>E</sub> and the OT. Graph (a) demonstrates that a 10% rise in the IR<sub>E</sub> positively affects the OT in the short and long run. Figure 8a indicates that a 10% increase in the C<sub>E</sub> positively impacts the OT in the short and long term. Conversely, a 10% decrease in the C<sub>E</sub> is associated with an increased OT in both the short and long run (Figure 8b). Figure 9b indicates that a 10% decrease in the SK negatively influences the OT in both the short and long run.

This result may appear counterintuitive, but suggests that a reduction in the B<sub>E</sub>, typically indicative of less vibration or frictional excitation, may signal sensor underperformance or hidden mechanical anomalies, leading to a compensatory rise in the OT.
<figure><img src="image_9.png" /><figcaption>**Figure 7.** The impulse response of LnIR<sub>E</sub> to a 10% increase and decrease over time. The time (minutes) is plotted on the horizontal axis, showing how the oil temperature (OT, LnT) reacts dynamically. Confidence intervals are displayed using shaded regions. (a) A 10% rise in the IRE positively affects the OT in both the short and long run. (b) A 10% decrease in the IRE shows a corresponding negative impact on the OT in the short and long term.</figcaption></figure>

<figure><img src="image_10.png" /><figcaption>**Figure 8.** The impulse response of LnCE over time following a 10% change. The x-axis represents time in minutes. The graphs illustrate the short- and long-run effects on oil temperature (LnT). (a) A 10% increase in CE positively impacts the OT in both the short and long term. (b) A 10% decrease in CE is also associated with an increase in OT in both the short and long run.</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>**Figure 9.** The dynamic response of the oil temperature (OT) to changes in the $B_E$. The x-axis shows time in minutes. The graphs illustrate the system's behavior following 10% upward and downward shifts in the $B_E$. (a) A 10% increase in the BE and its effect on the OT in both the short and long term. (b) A 10% decrease in the BE, which may positively or negatively influence the OT in both the short and long run.</figcaption></figure>
<figure><img src="image_12.png" /><figcaption>**Figure 10.** The impulse response plot for LnSK, with time labeled in minutes on the x-axis. The graphs illustrate changes in oil temperature (OT) in response to SK perturbations, providing insight into the system's dynamic behavior. (a) A 10% increase in SK leads to a marginal but delayed rise in OT. (b) A 10% decrease in SK results in a faster and more pronounced increase in OT.</figcaption></figure>

Figure 10 provides insight into the dynamic behavior of the OT in response to SK variation. In Figure 10a, a 10% increase in the SK led to a marginal, but delayed, rise in the OT, while Figure 10b reveals that a 10% decrease in the SK resulted in a faster and more pronounced increase in the OT. This implies that a lower SK, which typically represents less impulsiveness in the signal, could mask early-stage faults, leading to delayed detection and increased thermal buildup over time.

These findings emphasize the importance of jointly monitoring the SK and $B_E$ in addition to traditional energy-based indicators.

Table 4 summarizes the outcomes of various diagnostic tests performed in this study. To assess the serial correlation among the selected CIs, the LM test was applied, and the findings confirmed the absence of a serial correlation in the model. The Breusch-Pagan-Godfrey test was used to evaluate the heteroscedasticity, revealing that the model does not suffer from any heteroscedasticity issues. Furthermore, the Ramsey RESET test was conducted to verify the correctness of the DARDL model specification, and the results suggest that the model is appropriately specified. Lastly, the Jarque-Bera test confirmed that the residuals follow a normal distribution.

<table><caption>**Table 4.** Diagnostic tests.</caption><thead><tr><td>Diagnostic Test</td><td>X<sup>2</sup> (p-Value)</td><td>Result</td></tr></thead><tbody><tr><td>Breusch-Godfrey LM</td><td>0.16 (0.863)</td><td>No evidence of serial correlations</td></tr><tr><td>Breusch-Pagan-Godfrey</td><td>0.388 (0.960)</td><td>No evidence of heteroscedasticity</td></tr><tr><td>Ramsey RESET test</td><td>2.839 (0.213)</td><td>Model specified correctly</td></tr><tr><td>Jarque-Bera test</td><td>0.144 (0.930)</td><td>Residuals are normally estimated</td></tr></tbody></table>

### 4.4. Wavelet Coherence Analysis

The WCA simultaneously captures the correlation and potential causality between two CIs across both the time and frequency domains [33–39]. In this study, we applied a WCA to examine the time-varying relationships and directional causality between the BGCIs and the lubrication OT.

Figure 11 illustrates the analysis over short (0–4 scale), medium (4–8 scale), and long (8–16 scale) periods. The strength of the association between the two time series is represented by color: blue (cool) indicates a weak or no correlation (a correlation close to zero), while red (warm) indicates a strong and significant relationship within the cone of influence. Areas outside the cone are not considered reliable for interpretation. In
the wavelet coherence diagrams, arrows within the thick black contour lines indicate the direction of causality and the phase relationship between the CIs.

Figure 11a presents the WCA of the relationship between the *OT* and the *OR*<sub>E</sub>. The presence of upward arrows indicates a positive relationship at low frequencies (short-term), which becomes stronger as the frequency increases (longer-term). This suggests that the *OR*<sub>E</sub> influences the *OT* in the short, medium, and long run. Figure 11b provides evidence of a positive relationship between the *OT* and the *SK*, with most arrows pointing upward or upward-right, especially at high frequencies in the long run. The relationship also appears more stable in the long run compared to the short and medium term.

Figure 11c,d illustrate the positive relationship between the *OT* and the *IR*<sub>E</sub>. In Figure 11c, the rightward-pointing arrows indicate that the *IR*<sub>E</sub> influenced the *OT*. Similarly, Figure 11d,e show that a positive relationship exists between the *OT* and the *C*<sub>E</sub> and *B*<sub>E</sub>, respectively, particularly at higher frequencies.

4.5. *Frequency-Domain Causality*

This study employed the FDC test introduced by [20] to examine potential links be-
tween the OT and the BGCI values in the Bell 407 helicopter. The results of the FDC test
provide valuable insights into the dynamics of the causal relationships between the CIs
across different temporal frequencies. By analyzing causality across long-term, medium-
term, and short-term frequencies, this method offers a more detailed understanding of
interrelationships that may not be captured by traditional causality tests. The null hypoth-
esis posits that there is no Granger causality between the independent and dependent
CIs. This hypothesis can be rejected at specific frequencies: ωᵢ = 0→1 for the long term,
ωᵢ = 1→2 for the medium term, and ωᵢ = 2→3 for the short term. The results of the FDC
test are presented in Figures 12–16.

In Figures 12–16, the x-axis represents the frequency (ω) in Hertz [Hz] and the y-axis displays the corresponding F-statistic. The horizontal dashed line indicates the 5% critical value for significance. F-statistic values above this threshold indicate the rejection of the null hypothesis of no Granger causality from the CIs to the OT at that frequency.

Figure 12 indicates a statistically significant causal relationship between the *OT* and the *OR*<sub>E</sub> in the long term. However, as the frequency increases, the test statistic decreases, falling below the critical values. This implies that the causal relationship from the *OT* to the *OR*<sub>E</sub> is not significant in the short term. Figure 13 demonstrates the presence of causality in the medium and short run, confirming the interdependence between the *B*<sub>E</sub> and the *OT*. Figure 14 illustrates that there is causality from the *C*<sub>E</sub> to the *OT* in the short and medium terms at the 5% significance level. However, in the long term, the causal relationship is significant at the frequency range ω ∈ [0.8, 1]. Figure 15 indicates a significant causal relationship from the *OT* to the *IR*<sub>E</sub> in the long term. Figure 16 demonstrates the presence of causality in both the medium and short run, confirming the causality relationship from the *SK* to the *OT*.
<figure><img src="image_13.png" /><figcaption>**Figure 11.** WCT analysis between the *OT* and the explanatory CIs: (a) WCT between the *OT* and the *OR*<sub>*E*</sub>, (b) WCT between the *OT* and the *SK*, (c) WCT between the *OT* and the *IR*<sub>*E*</sub>, (d) WCT between the *OT* and the *C*<sub>*E*</sub>, and (e) WCT between the *OT* and the *B*<sub>*E*</sub>.</figcaption></figure>
<figure><img src="image_14.png" /><figcaption>Figure 12. Breitung–Candelon Spectral Granger Causality Test between *OT* and *OR*<sub>*E*</sub>.</figcaption></figure>

<figure><img src="image_15.png" /><figcaption>Figure 13. Breitung–Candelon Spectral Granger Causality Test between *OT* and *B*<sub>*E*</sub>.</figcaption></figure>

<figure><img src="image_16.png" /><figcaption>Figure 14. Breitung–Candelon Spectral Granger Causality Test between *OT* and *C*<sub>*E*</sub>.</figcaption></figure>
<figure><img src="image_17.png" /><figcaption>Figure 15. Breitung–Candelon Spectral Granger Causality Test between OT and *IR*<sub>E</sub>.</figcaption></figure>

<figure><img src="image_18.png" /><figcaption>Figure 16. Breitung–Candelon Spectral Granger Causality Test between *OT* and *SK*.</figcaption></figure>

**5. Conclusions and Policy Implications**

This study offers valuable insights into the relationship between the *OT* and the BGCIs in the Bell 407 helicopter. It effectively applies the DARDL model to examine both the short- and long-term relationships between the *OT* and the BGCIs, improving the understanding of how changes in condition indicators influence the *OT*. This addresses gaps in earlier research, which mainly focused on the effects of individual CIs. In addition, the wavelet coherence method provides a time-frequency analysis that captures how the *OT–BGCI* relationship evolves over time and across different frequencies. This is especially useful for detecting short-term fluctuations and long-term trends that traditional time-domain models may overlook. The WCA also proved effective in identifying lead-lag relationships and structural shifts, offering deeper insight into the dynamic interactions between the *OT* and the BGCIs. By revealing localized time-frequency dependencies, this method is particularly well-suited for assessing how variations in the *OT* impact the mechanical performance under different operational conditions. Moreover, this study also explored the causality between the *OT* and the BGCIs using a frequency-domain approach. This method reveals that the direction of causality may vary across different frequencies, providing a more comprehensive understanding of the interactions between these CIs compared to traditional time-domain analyses.

The DARDL model results indicate that the condition indicators *B*<sub>E</sub>, *C*<sub>E</sub>, *IR*<sub>E</sub>, and *ORE* have a significant positive impact on the *OT* in the long term. Additionally, the WCA highlights the dynamic interactions between these CIs and the *OT* across various time scales and frequencies. The overall WCA findings support the conclusion that the *B*<sub>E</sub>, *C*<sub>E</sub>, *IR*<sub>E</sub>,
and *OR*<sub>E</sub> drive changes in the *OT* in both the short and long run, aligning with the DARDL estimates and confirming a positive relationship. Furthermore, the FDC analysis revealed that the *OT* has a significant influence on the BGCIs, with clear evidence of bidirectional causality, indicating mutual interdependence.

The study reinforces how HUMS frameworks differ from and improve upon conven-
tional condition monitoring. Unlike threshold-based systems that generate alerts post-
deviation, HUMSs continuously track mechanical and thermal health parameters like BG-
CIs and the OT during normal operations. By applying advanced models such as DARDL
and wavelet coherence, we can move beyond descriptive monitoring to predictive diagnos-
tics, enabling the simulation of failure scenarios and deeper causality analyses—features
not achievable with traditional methods. This evolution is crucial for minimizing downtime
and optimizing maintenance in aerospace applications.

This study contributes valuable insights for predictive maintenance strategies. For
maintenance engineers and operational managers, this shows the importance of closely
monitoring these CIs as part of predictive maintenance strategies. Elevated values in
these parameters can serve as early warning signs of a rising *OT*, which may compromise
lubrication effectiveness and lead to increased wear or the potential failure of gearbox
components. Implementing real-time monitoring systems that track these CIs alongside
the *OT* can improve the diagnostic accuracy and maintenance scheduling. Additionally,
integrating these insights into maintenance protocols can help reduce unplanned downtime,
extend equipment life, and enhance the overall operational safety and efficiency in aviation
and similar high-performance mechanical systems [40–45].

However, the snapshot-based data acquisition method—where 2 s samples are taken
every 4 min—introduces a known tradeoff. While it offers efficiency in terms of the data
volume and transmission bandwidth, it may miss short-lived transient fault signatures that
occur outside the sampling window. This limitation is particularly relevant in dynamic
systems such as helicopters, where rapid load changes or short-duration events may pro-
vide early signs of faults. To partially mitigate this, the HUMS employs long-term repeated
snapshots under steady-state conditions, and our use of a wavelet-based time-frequency
analysis enhances the sensitivity to transient patterns within those windows. Nonetheless,
we acknowledge that continuous monitoring or adaptive, higher-frequency snapshot in-
tervals could further improve the system’s diagnostic resolution and responsiveness to
intermittent faults.

**Author Contributions:** Conceptualization, L.S., E.B., and M.B.; methodology, L.S. and E.B.; software, L.S.; validation, M.B. and E.B.; formal analysis, L.S.; investigation, E.B.; resources, E.B.; data curation, E.B.; writing—original draft preparation, L.S.; writing—review and editing, M.B.; visualization, E.B. and M.B.; supervision, L.S., E.B., and M.B.; project administration, L.S. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research received no external funding.

**Data Availability Statement:** The dataset is available from the authors upon request.

**Conflicts of Interest:** Eric Bechhofer is employed by Green Power Monitoring Systems, LLC—USA. The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as potential conflicts of interest.
Abbreviations

<table><tr><td>BGCI</td><td>Bearing Gearbox Condition Indicator</td></tr><tr><td>CWT</td><td>Continuous Wavelet Transformation</td></tr><tr><td>OT</td><td>Oil Temperature</td></tr><tr><td>HUMS</td><td>Health and Usage Monitoring System</td></tr><tr><td>DARDL</td><td>Dynamic Autoregressive Distributed Lag</td></tr><tr><td>CIs</td><td>Condition Indicators</td></tr><tr><td>WCA</td><td>Wavelet Coherence Analysis</td></tr><tr><td>ORE</td><td>Outer Race Energy</td></tr><tr><td>IRE</td><td>Inner Race Energy</td></tr><tr><td>CE</td><td>Cage Energy</td></tr><tr><td>BE</td><td>Ball Energy</td></tr><tr><td>SK</td><td>Spectral Kurtosis</td></tr><tr><td>Ln</td><td>Natural Logarithm</td></tr><tr><td>KPSS</td><td>Kwiatkowski-Phillips-Schmidt-Shin</td></tr><tr><td>AIC</td><td>Akaike Information Criterion</td></tr><tr><td>ECT</td><td>Error Correction Terms</td></tr><tr><td>CUSUM</td><td>Cumulative Sum of Recursive Residuals</td></tr><tr><td>CUSUMSQ</td><td>Cumulative Sum of Recursive Residual Squares</td></tr><tr><td>WPS</td><td>Wavelet Power Spectrum</td></tr><tr><td>FDC</td><td>Frequency-Domain Causality</td></tr><tr><td>LM</td><td>Lagrange Multiplier</td></tr></table>

References

1. Ranasinghe, K.; Kapoor, R.; Gardi, A.; Sabatini, R.; Wickramanayake, V.; Ludovici, D. Vehicular Sensor Network and Data Analytics for a Health and Usage Management System. *Sensors* **2020**, *20*, 5892. [CrossRef] [PubMed]

2. Zhang, Q.; Wu, Z.; An, B.; Sun, R.; Cui, Y. Digital Twin-Based Technical Research on Comprehensive Gear Fault Diagnosis and Structural Performance Evaluation. *Sensors* **2025**, *25*, 2775. [CrossRef] [PubMed]

3. Kannan, V.; Zhang, T.; Li, H. A Review of the intelligent condition monitoring of rolling element bearings. *Machines* **2024**, *12*, 484. [CrossRef]

4. Qiu, S.; Cui, X.; Ping, Z.; Shan, N.; Li, Z.; Bao, X.; Xu, X. Deep Learning Techniques in Intelligent Fault Diagnosis and Prognosis for Industrial Systems: A Review. *Sensors* **2023**, *23*, 1305. [CrossRef] [PubMed]

5. Tiboni, M.; Remino, C.; Bussola, R.; Amici, C. A Review on Vibration-Based Condition Monitoring of Rotating Machinery. *Appl. Sci.* **2022**, *12*, 972. [CrossRef]

6. Xie, F.; Liu, H.; Dong, J.; Wang, G.; Wang, L.; Li, G. Research on the Gearbox Fault Diagnosis Method Based on Multi-Model Feature Fusion. *Machines* **2022**, *10*, 1186. [CrossRef]

7. Wei, L.; Cheng, Z.; Cheng, J.; Hu, N.; Yang, Y. A Fault Detection Method Based on an Oil Temperature Forecasting Model Using an Improved Deep Deterministic Policy Gradient Algorithm in the Helicopter Gearbox. *Entropy* **2022**, *24*, 1394. [CrossRef] [PubMed]

8. Saidi, L.; Benbouzid, M. Prognostics and Health Management of Renewable Energy Systems: State of the Art Review, Challenges, and Trends. *Electronics* **2021**, *10*, 2732. [CrossRef]

9. Babay, O.; Saidi, L.; Bechhofer, E.; Benbouzid, M. Investigation of the dynamic relationship between oil temperature and bearing gearbox condition indicator values for the Bell 407 helicopter based on cointegration analysis. In Proceedings of the Annual Conference of the PHM Society, Salt Lake City, UT, USA, 4–7 November 2024.

10. Wu, J.; Wang, Z.; Cao, H.; Huang, Y.; Liu, Y. Deep Clustering Variational Network for Helicopter Regime Recognition in HUMS. *Aerosp. Sci. Technol.* **2022**, *124*, 107553. [CrossRef]

11. Modaresahmadi, S.; Chittenden, R.; Mueller, M. Convective Heat Transfer Analysis of a Laminated Flux Focusing Magnetic Gearbox. *Therm. Sci. Eng. Prog.* **2020**, *18*, 100552. [CrossRef]

12. Rashid, H.S.J.; Yang, Z.; Wang, J. Reliability Model for Helicopter Main Gearbox Lubrication System Using Influence Diagrams. *Reliab. Eng. Syst. Saf.* **2015**, *139*, 50–57. [CrossRef]

13. Saidi, L.; Bechhofer, E.; Benbouzid, M. Dynamic Relationship Between Oil Temperature and BGCI in Bell 407 Helicopter: A Cointegration and Autoregressive Distributed Lag Approach. *Int. J. Progn. Health Manag.* **2025**, *16*, 1. [CrossRef]

14. Bechhoefer, E.; Hajimohammadali, M. Process for Turboshaft Engine Performance Trending. *Annu. Conf. PHM Soc.* **2023**, *15*, 1. [CrossRef]
15. Ayad, M.; Saoudi, K.; Rezki, M.; Benziane, M.; Arabi, A.; Chikouche, D. Early detection of defects in gear systems using autocorrelation of Morlet wavelet transforms. *Stroj. Vestn. J. Mech. Eng.* **2022**, *68*, 56–65. [CrossRef]

16. Praveen, H.M.; Sabareesh, G.R. Wavelet based real-time planetary gearbox health monitoring under non-stationary operation. *Exp. Tech.* **2022**, *46*, 861–875. [CrossRef]

17. Lupea, I.; Lupea, M. Continuous Wavelet Transform and CNN for Fault Detection in a Helical Gearbox. *Appl. Sci.* **2025**, *15*, 950. [CrossRef]

18. Majeed, R.; Haddar, M.; Chaari, F.; Haddar, M. Exponentially Weighted Moving Average Control Chart for Fault Detection of the Spur Gear Transmission System. In *Proceedings of the International Conference on Acoustics and Vibration*; Springer: Cham, Switzerland, 2022; pp. 123–132. [CrossRef]

19. Jordan, S.; Philips, A.Q. Cointegration testing and dynamic simulations of autoregressive distributed lag models. *Stata J.* **2018**, *18*, 902–923. [CrossRef]

20. Breitung, J.; Candelon, B. Testing for short-and long-run causality: A frequency-domain approach. *J. Econom.* **2006**, *132*, 363–378. [CrossRef]

21. Łuczak, D. Machine Fault Diagnosis through Vibration Analysis: Continuous Wavelet Transform with Complex Morlet Wavelet and Time–Frequency RGB Image Recognition via Convolutional Neural Network. *Electronics* **2024**, *13*, 452. [CrossRef]

22. Habbouche, H.; Benkedjouh, T.; Amirat, Y.; Benbouzid, M. A Wavelet Transform-Based Transfer Learning Approach for Enhanced Shaft Misalignment Diagnosis in Rotating Machinery. *Electronics* **2025**, *14*, 341. [CrossRef]

23. Lee, J.-H.; Okwuosa, C.N.; Hur, J.-W. Extruder Machine Gear Fault Detection Using Autoencoder LSTM via Sensor Fusion Approach. *Inventions* **2023**, *8*, 140. [CrossRef]

24. Gao, Z.-W.; Zhao, R.; Zhang, J.; Zhang, X.; Zhang, Y. An optimized updating adaptive federated learning for pumping units collaborative diagnosis with label heterogeneity and communication redundancy. *Eng. Appl. Artif. Intell.* **2025**, *152*, 110724. [CrossRef]

25. Zhou, F.; Yang, Y.; Wang, C.; Hu, X. Federated Learning Based Fault Diagnosis Driven by Intra-Client Imbalance Degree. *Entropy* **2023**, *25*, 606. [CrossRef] [PubMed]

26. GPMS. Available online: https://www.gpms-vt.com/ (accessed on 15 June 2025).

27. GPMS International. GPMS Foresight MX HUMS for MD530F. Heli-Expo 2023. 2023. Available online: https://www.helis.com/database/news/heliexpo23-gpms-md530f-hums (accessed on 17 June 2025).

28. GPMS International. GPMS Foresight MX HUMS for Airbus Models. Heli-Expo 2023. 2023. Available online: https://www.helis.com/database/news/heliexpo23-gpms-as350b2-ec145 (accessed on 17 June 2025).

29. Pesaran, M.H.; Shin, Y.; Smith, R.J. Bounds testing approaches to the analysis of level relationships. *J. Appl. Econom.* **2001**, *16*, 289–326. [CrossRef]

30. Dickey, D.A.; Fuller, W.A. Distribution of the estimators for autoregressive time series with a unit root. *J. Am. Stat. Assoc.* **1979**, *74*, 427–431. [CrossRef] [PubMed]

31. Kwiatkowski, D.; Phillips, P.C.; Schmidt, P.; Shin, Y. Testing the null hypothesis of stationarity against the alternative of a unit root: How sure are we that economic time series have a unit root? *J. Econom.* **1992**, *54*, 159–178. [CrossRef]

32. Narayan, P.K. The saving and investment nexus for China: Evidence from cointegration tests. *Appl. Econ.* **2005**, *37*, 1979–1990. [CrossRef]

33. Goupillaud, P.; Grossmann, A.; Morlet, J. Cycle-octave and related transforms in seismic signal analysis. *Geoexploration* **1984**, *23*, 85–102. [CrossRef]

34. Hainmueller, J.; Hazlett, C. Kernel regularized least squares: Reducing misspecification bias with a flexible and interpretable machine learning approach. *Polit. Anal.* **2014**, *22*, 143–168. [CrossRef]

35. Torrence, C.; Compo, G.P. A practical guide to wavelet analysis. *Bull. Am. Meteorol. Soc.* **1998**, *79*, 61–78. [CrossRef]

36. Huang, W.; Wang, H.; Li, C.; Zi, Y.; He, Z. Feature Extraction for Rolling Element Bearing Faults Using Resonance Sparse Signal Decomposition. *Exp. Tech.* **2017**, *41*, 251–265. [CrossRef]

37. Bechhoefer, E.; He, D. Case Study of a Faulted Planet Bearing. *Proc. Annu. Conf. PHM Soc.* **2016**, *8*, 1. [CrossRef]

38. Saidi, L.; Benbouzid, M.; Diallo, D.; Amirat, Y.; Elbouchikhi, E.; Wang, T. Higher-order spectra analysis-based diagnosis method of blades biofouling in a PMSG driven tidal stream turbine. *Energies* **2020**, *13*, 2888. [CrossRef]

39. Saidi, L.; Bechhoefer, E.; Ben Ali, J.; Benbouzid, M. Wind turbine high-speed shaft bearing degradation analysis for run-to-failure testing using spectral kurtosis. In *Proceedings of the 16th International Conference on Sciences and Techniques of Automatic Control and Computer Engineering (STA), Monastir, Tunisia*, 16–18 December 2015.

40. Vasić, M.; Stojanović, B.; Blagojević, M. Failure analysis of idler roller bearings in belt conveyors. *Eng. Fail. Anal.* **2020**, *117*, 104898. [CrossRef]

41. Pastukhov, A.; Timashov, E. Procedure for simulation of stable thermal conductivity of bearing assemblies. *Adv. Eng. Lett.* **2023**, *2*, 58–63. [CrossRef]
42. Vertical Aviation International (VAI); GPMS International. VAI Partners with GPMS to Expand Access to Advanced HUMS. 2025. Available online: https://verticalavi.org/vai-daily/vai-partners-with-gpms-to-expand-access-to-advanced-hums/ (accessed on 17 June 2025).

43. Airbus Helicopters. Connected Services. 2023. Available online: https://www.airbus.com/en/products-services/helicopters/hcare-services/connected-services (accessed on 17 June 2025).

44. SciELO Brazil. Addressing Gearbox Health Monitoring Challenges for Helicopters: A Machine Learning Approach. 2023. Available online: https://www.scielo.br/j/aabc/a/ZBmSydzGg6GdMWycGn8RDjL/ (accessed on 17 June 2025).

45. Helicopted. The Role of Predictive Maintenance in Helicopter Fleet Management. 2023. Available online: https://helicopted.com/the-role-of-predictive-maintenance-in-helicopter-fleet-management/ (accessed on 17 June 2025).

**Disclaimer/Publisher's Note:** The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.