Incipient fault detection in a hydraulic system is a challenge in the condition monitoring community. Existing research mainly monitors abnormal working conditions in hydraulic systems by separately detecting the key working parameter, which often causes a high miss warning rate for incipient faults due to the oversight of parameter dependence. A principal component analysis provides an effective method for incipient fault detection by taking the correlation of multiple parameters into consideration, but this technique assumes the systems are Gaussian-distributed, making it invalid for a dynamic non-Gaussian system. In this paper, we combine a canonical variable analysis (CVA) and adaptive kernel density estimation (AKDE) for the early fault detection of nonlinear dynamic hydraulic systems. The collected hydraulic system data set was used to construct the typical variable space, and the state space and residual space are divided to represent the characteristics of different correlations between the two variables, which are quantitatively described using Hotelling's T2 and Q. In order to investigate the proper upper control limits, AKDE was utilised to estimate the underlying probability density functions of T2 and Q by taking the nonlinearity of the hydraulic system variables into consideration. The advantages of the proposed approach for incipient fault detection are illustrated via a marine power plant lubrication system.