Exploring timescale-specific functional brain networks and their associations with aging and cognitive performance in a healthy cohort without dementia.

INTRODUCTION: Functional brain networks (FBNs) coordinate brain functions and are studied in fMRI using blood-oxygen-level-dependent (BOLD) signal correlations. Previous research links FBN changes to aging and cognitive decline, but various physiological factors influnce BOLD signals. Few studies have investigated the intrinsic components of the BOLD signal in different timescales using signal decomposition. This study aimed to explore differences between intrinsic FBNs and traditional BOLD-FBN, examining their associations with age and cognitive performance in a healthy cohort without dementia. MATERIALS AND METHODS: A total of 396 healthy participants without dementia (men = 157; women = 239; age range = 20-85 years) were enrolled in this study. The BOLD signal was decomposed into several intrinsic signals with different timescales using ensemble empirical mode decomposition, and FBNs were constructed based on both the BOLD and intrinsic signals. Subsequently, network features-global efficiency and local efficiency values-were estimated to determine their relationship with age and cognitive performance. RESULTS: The findings revealed that the global efficiency of traditional BOLD-FBN correlated significantly with age, with specific intrinsic FBNs contributing to these correlations. Moreover, local efficiency analysis demonstrated that intrinsic FBNs were more meaningful than traditional BOLD-FBN in identifying brain regions related to age and cognitive performance. CONCLUSIONS: These results underscore the importance of exploring timescales of BOLD signals when constructing FBN and highlight the relevance of specific intrinsic FBNs to aging and cognitive performance. Consequently, this decomposition-based FBN-building approach may offer valuable insights for future fMRI studies.

The strengthened impact of water availability at interannual and decadal time scales on vegetation GPP.

Water availability (WA) is a key factor influencing the carbon cycle of terrestrial ecosystems under climate warming, but its effects on gross primary production (EWA-GPP ) at multiple time scales are poorly understood. We used ensemble empirical mode decomposition (EEMD) and partial correlation analysis to assess the WA-GPP relationship (RWA-GPP ) at different time scales, and geographically weighted regression (GWR) to analyze their temporal dynamics from 1982 to 2018 with multiple GPP datasets, including near-infrared radiance of vegetation GPP, FLUXCOM GPP, and eddy covariance-light-use efficiency GPP. We found that the 3- and 7-year time scales dominated global WA variability (61.18% and 11.95%), followed by the 17- and 40-year time scales (7.28% and 8.23%). The long-term trend also influenced 10.83% of the regions, mainly in humid areas. We found consistent spatiotemporal patterns of the EWA-GPP and RWA-GPP with different source products: In high-latitude regions, RWA-GPP changed from negative to positive as the time scale increased, while the opposite occurred in mid-low latitudes. Forests had weak RWA-GPP at all time scales, shrublands showed negative RWA-GPP at long time scales, and grassland (GL) showed a positive RWA-GPP at short time scales. Globally, the EWA-GPP , whether positive or negative, enhanced significantly at 3-, 7-, and 17-year time scales. For arid and humid zones, the semi-arid and sub-humid zones experienced a faster increase in the positive EWA-GPP , whereas the humid zones experienced a faster increase in the negative EWA-GPP . At the ecosystem types, the positive EWA-GPP at a 3-year time scale increased faster in GL, deciduous broadleaf forest, and savanna (SA), whereas the negative EWA-GPP at other time scales increased faster in evergreen needleleaf forest, woody savannas, and SA. Our study reveals the complex and dynamic EWA-GPP at multiple time scales, which provides a new perspective for understanding the responses of terrestrial ecosystems to climate change.

[Snoring noise removal method for bowel sound signal during sleep].

Monitoring of bowel sounds is an important method to assess bowel motility during sleep, but it is seriously affected by snoring noise. In this paper, the complete ensemble empirical mode decomposition with adaptive noise (CEEMDAN) method was applied to remove snoring noise from bowel sounds during sleep. Specifically, the noisy bowel sounds were first band-pass filtered, then decomposed by the CEEMDAN method, and finally the appropriate components were selected to reconstruct the pure bowel sounds. The results of semi-simulated and real data showed that the CEEMDAN method was better than empirical mode decomposition and wavelet denoising method. The CEEMDAN method is used to remove snoring noise from bowel sounds during sleep, which lays an important foundation for using bowel sounds to assess the intestinal motility during sleep.

Satellite observed reversal in trends of spring phenology in the middle-high latitudes of the Northern Hemisphere during the global warming hiatus.

The start of the growing season (SOS) is essential to track the responses of vegetation to climate change. However, recent findings on whether the SOS in the middle-high latitudes of the Northern Hemisphere (NH) continued to advance or reversed during the global warming hiatus were not consistent. It is necessary to investigate the causes of this controversy and to examine the relationship between the SOS and preseason temperature trends. To this end, we first applied four widely used phenology extraction methods to derive the SOS from the GIMMS NDVI3g dataset and then used the ensemble empirical modal decomposition (EEMD) method to extract the nonlinear trends of the SOS and preseason temperature. Our results clarify, for the first time, that the limitations of the linear assumption-based trend analysis methods are an important but overlooked cause of the discrepancies among existing studies on whether the SOS was advanced or delayed in the NH (>30° N) during the global warming hiatus. We further revealed the range of the mismatches between the SOS and preseason temperature trends at the latitude, altitude and biome levels. Specifically, we discovered that the SOS in the NH (>30° N) obtained by the four phenology extraction methods showed a significant reversal from advance to delay during the global warming hiatus, and the corresponding average rate of change was very small. The area showing increasing preseason temperatures decreased during the global warming hiatus, but it always occupied most of the NH (>30° N). However, delayed SOS trends were dominant in the NH from 50° N to 60° N, above 3000 m and in biomes other than TBMF and BF. Accordingly, using an EEMD-like approach to evaluate the changes in the SOS and preseason temperature is necessary for improving our understanding of the changes in the SOS and their association with climate.

Desert Abiotic Carbon Sequestration Weakening by Precipitation.

Desert carbon sequestration plays an active role in promoting carbon neutralization. However, the current understanding of the effect of hydrothermal interactions and soil properties on desert carbon sequestration after precipitation remains unclear. Based on the experiment in the hinterland of the Taklimakan Desert, we found that the heavy precipitation will accelerate the weakening of abiotic carbon sequestration in deserts under the background of global warming and intensified water cycle. The high soil moisture can significantly stimulate sand to release CO2 at an incredible speed by rapidly increasing microbial activity and organic matter diffusion. At this time, the CO2 flux in the shifting sand was synergistically affected by soil temperature and soil moisture. As far as soil properties are concerned, with less organic carbon substrate and stronger soil alkalinity, the carbon sequestration of shifting sand is gradually highlighted and strengthened at low temperature. On the contrary, the carbon sequestration of shifting sand is gradually weakened. Our study provides a new way to assess the contribution of desert to the global carbon cycle and improve the accuracy and scope of application.

Reliable Denoising Strategy to Enhance the Accuracy of Arrival Time Picking of Noisy Microseismic Recordings.

We propose a method to enhance the accuracy of arrival time picking of noisy microseismic recordings. A series of intrinsic mode functions (IMFs) of the microseismic signal are initially decomposed by employing the ensemble empirical mode decomposition. Subsequently, the sample entropy values of the obtained IMFs are calculated and applied to set an appropriate threshold for selecting IMFs. These are then reconstructed to distinguish between noise and useful signals. Ultimately, the Akaike information criterion picker is used to determine the arrival time of the denoised signal. Test results using synthetic noisy microseismic recordings demonstrate that the proposed approach can significantly reduce picking errors, with errors within the range of 1-3 sample intervals. The proposed method can also give a more stable picking result when applied to different microseismic recordings with different signal-to-noise ratios. Further application in real microseismic recordings confirms that the developed method can estimate an accurate arrival time of noisy microseismic recordings.

Research on Denoising Method for Hydroelectric Unit Vibration Signal Based on ICEEMDAN-PE-SVD.

Vibration monitoring and analysis play a crucial role in the fault diagnosis of hydroelectric units. However, accurate extraction and identification of fault features from vibration signals are challenging because of noise interference. To address this issue, this study proposes a novel denoising method for vibration signals based on improved complementary ensemble empirical mode decomposition with adaptive noise (ICEEMDAN), permutation entropy (PE), and singular value decomposition (SVD). The proposed method is applied for the analysis of hydroelectric unit sway monitoring. Firstly, the ICEEMDAN method is employed to process the signal and obtain several intrinsic mode functions (IMFs), and then the PE values of each IMF are calculated. Subsequently, based on a predefined threshold of PE, appropriate IMFs are selected for reconstruction, achieving the first denoising effect. Then, the SVD is applied to the signal after the first denoising effect, resulting in the SVD spectrum. Finally, according to the principle of the SVD spectrum and the variation in the singular value and its energy value, the signal is reconstructed by choosing the appropriate reconstruction order to achieve the secondary noise reduction effect. In the simulation and case analysis, the method is better than the commonly used wavelet threshold, SVD, CEEMDAN-PE, and ICEEMDAN-PE, with a signal-to-noise ratio (SNR) improvement of 6.9870 dB, 4.6789 dB, 8.9871 dB, and 4.3762 dB, respectively, and where the root-mean-square error (RMSE) is reduced by 0.1426, 0.0824, 0.2093 and 0.0756, respectively, meaning that our method has a better denoising effect and provides a new way for denoising the vibration signal of hydropower units.

Displacement Monitoring of a Bridge Based on BDS Measurement by CEEMDAN-Adaptive Threshold Wavelet Method.

From the viewpoint of BDS bridge displacement monitoring, which is easily affected by background noise and the calculation of a fixed threshold value in the wavelet filtering algorithm, which is often related to the data length. In this paper, a data processing method of Complete Ensemble Empirical Mode Decomposition with Adaptive Noise (CEEMDAN), combined with adaptive threshold wavelet de-noising is proposed. The adaptive threshold wavelet filtering method composed of the mean and variance of wavelet coefficients of each layer is used to de-noise the BDS displacement monitoring data. CEEMDAN was used to decompose the displacement response data of the bridge to obtain the intrinsic mode function (IMF). Correlation coefficients were used to distinguish the noisy component from the effective component, and the adaptive threshold wavelet de-noising occurred on the noisy component. Finally, all IMF were restructured. The simulation experiment and the BDS displacement monitoring data of Nanmao Bridge were verified. The results demonstrated that the proposed method could effectively suppress random noise and multipath noise, and effectively obtain the real response of bridge displacement.

A Multi-Target Localization and Vital Sign Detection Method Using Ultra-Wide Band Radar.

Life detection technology using ultra-wideband (UWB) radar is a non-contact, active detection technology, which can be used to search for survivors in disaster rescues. The existing multi-target detection method based on UWB radar echo signals has low accuracy and has difficulty extracting breathing and heartbeat information at the same time. Therefore, this paper proposes a new multi-target localization and vital sign detection method using ultra-wide band radar. A target recognition and localization method based on permutation entropy (PE) and K means++ clustering is proposed to determine the number and position of targets in the environment. An adaptive denoising method for vital sign extraction based on ensemble empirical mode decomposition (EEMD) and wavelet analysis (WA) is proposed to reconstruct the breathing and heartbeat signals of human targets. A heartbeat frequency extraction method based on particle swarm optimization (PSO) and stochastic resonance (SR) is proposed to detect the heartbeat frequency of human targets. Experimental results show that the PE-K means++ method can successfully recognize and locate multiple human targets in the environment, and its average relative error is 1.83%. Using the EEMD-WA method can effectively filter the clutter signal, and the average relative error of the reconstructed respiratory signal frequency is 4.27%. The average relative error of heartbeat frequency detected by the PSO-SR method was 6.23%. The multi-target localization and vital sign detection method proposed in this paper can effectively recognize all human targets in the multi-target scene and provide their accurate location and vital signs information. This provides a theoretical basis for the technical system of emergency rescue and technical support for post-disaster rescue.

A Frequency Estimation Scheme Based on Gaussian Average Filtering Decomposition and Hilbert Transform: With Estimation of Respiratory Rate as an Example.

Frequency estimation plays a critical role in vital sign monitoring. Methods based on Fourier transform and eigen-analysis are commonly adopted techniques for frequency estimation. Because of the nonstationary and time-varying characteristics of physiological processes, time-frequency analysis (TFA) is a feasible way to perform biomedical signal analysis. Among miscellaneous approaches, Hilbert-Huang transform (HHT) has been demonstrated to be a potential tool in biomedical applications. However, the problems of mode mixing, unnecessary redundant decomposition and boundary effect are the common deficits that occur during the procedure of empirical mode decomposition (EMD) or ensemble empirical mode decomposition (EEMD). The Gaussian average filtering decomposition (GAFD) technique has been shown to be appropriate in several biomedical scenarios and can be an alternative to EMD and EEMD. This research proposes the combination of GAFD and Hilbert transform that is termed the Hilbert-Gauss transform (HGT) to overcome the conventional drawbacks of HHT in TFA and frequency estimation. This new method is verified to be effective for the estimation of respiratory rate (RR) in finger photoplethysmography (PPG), wrist PPG and seismocardiogram (SCG). Compared with the ground truth values, the estimated RRs are evaluated to be of excellent reliability by intraclass correlation coefficient (ICC) and to be of high agreement by Bland-Altman analysis.

Seasonal to decadal spatiotemporal variations of the global ocean carbon sink.

The global ocean has absorbed approximately 30% of anthropogenic CO2  since the beginning of the industrial revolution. However, the spatiotemporal evolution of this important global carbon sink varies substantially on all timescales and has not yet been well evaluated. Here, based on a reconstructed observation-based product of surface ocean pCO2 and air-sea CO2  flux (the MPI-SOMFFN method), we investigated seasonal to decadal spatiotemporal variations of the ocean CO2  sink during the past three decades using an adaptive data analysis method. Two predominant variations are modulated annual cycles and decadal fluctuations, which account for approximately 46% and 25% of all extracted components, respectively. Although the whole summer to non-summer seasonal difference pattern is determined by the Southern Ocean, the non-summer CO2  sink at mid-latitudes in both hemispheres shows an increasing trend (a total increase of approximately 1.0 PgC during the period 1982-2019), while it is relatively stable in summer. On decadal timescales for the global ocean carbon sink, unlike the weakening decade (1990-1999) and the reinvigoration decade (2000-2009) in which the Southern Ocean plays the dominant role, the reinforcement decade (2010-2019) is mainly the result from the weakening source effect in the equatorial Pacific Ocean. Our results suggest that except for the Southern Ocean's role in the global ocean carbon sink, the strengthening non-summer's sink at mid-latitudes in both hemispheres and the decadal or longer timescales of equatorial Pacific Ocean dynamics should be fully considered in understanding the oceanic carbon cycle on a global scale.

Multiscale Temporal Variations of Atmospheric Mercury Distinguished by the Hilbert-Huang Transform Analysis Reveals Multiple El Niño-Southern Oscillation Links.

Atmospheric mercury (Hg) cycling is sensitive to climate-driven changes, but links with various teleconnections remain unestablished. Here, we revealed the El Niño-Southern Oscillation (ENSO) influence on gaseous elemental mercury (GEM) concentrations recorded at a background station in East Asia using the Hilbert-Huang transform (HHT). The timing and magnitude of GEM intrinsic variations were clearly distinguished by ensemble empirical mode decomposition (EEMD), revealing the amplitude of the GEM concentration interannual variability (IAV) is greater than that for diurnal and seasonal variability. We show that changes in the annual cycle of GEM were modulated by significant IAVs at time scales of 2-7 years, highlighted by a robust GEM IAV-ENSO relationship of the associated intrinsic mode functions. With confirmation that ENSO modulates the GEM annual cycle, we then found that weaker GEM annual cycles may have resulted from El Niño-accelerated Hg evasion from the ocean. Furthermore, the relationship between ENSO and GEM is sensitive to extreme events (i.e., 2015-2016 El Niño), resulting in perturbation of the long-term trend and atmospheric Hg cycling. Future climate change will likely increase the number of extreme El Niño events and, hence, could alter atmospheric Hg cycling and influence the effectiveness evaluation of the Minamata Convention on Mercury.

Data Processing Approaches to Measure Velocity of Electromagnetic Gun on Laser Screen in Complex Environment.

The exit velocity of the armature is an important indicator in measuring the launching performance of the electromagnetic gun. The non-contact photoelectric detection technology with the use of a laser screen was applied to the measurement of the armature velocity of the electromagnetic gun. By means of taking the signals that pass through the laser screen obtained by the velocity measurement system as the research object, we solved problems such as the harsh test environment of the launch armature velocity of the electromagnetic gun, the interferences on the armature signal passing through the laser screen unavoidably caused by various factors such as vibration, electromagnetic interference, shock wave, flare, smoke and fragments, and even the non-recognition of the signal passing through the laser screen in severe cases. A data-processing algorithm that combines the Ensemble Empirical Mode Decomposition (EEMD) with Correlation Algorithm (CA) was proposed, with the aim of processing the signals passing through the laser screen, while using the maximum slope point as the time passing through the laser screen so as to calculate the velocity of the armature passing the laser screen. This method can effectively reduce the influence of interference on the test results, and the test results from two sets of velocity measuring systems show that the velocity obtained by the proposed approach is highly consistent.

Using Complementary Ensemble Empirical Mode Decomposition and Gated Recurrent Unit to Predict Landslide Displacements in Dam Reservoir.

It is crucial to predict landslide displacement accurately for establishing a reliable early warning system. Such a requirement is more urgent for landslides in the reservoir area. The main reason is that an inaccurate prediction can lead to riverine disasters and secondary surge disasters. Machine learning (ML) methods have been developed and commonly applied in landslide displacement prediction because of their powerful nonlinear processing ability. Recently, deep ML methods have become popular, as they can deal with more complicated problems than conventional ML methods. However, it is usually not easy to obtain a well-trained deep ML model, as many hyperparameters need to be trained. In this paper, a deep ML method-the gated recurrent unit (GRU)-with the advantages of a powerful prediction ability and fewer hyperparameters, was applied to forecast landslide displacement in the dam reservoir. The accumulated displacement was firstly decomposed into a trend term, a periodic term, and a stochastic term by complementary ensemble empirical mode decomposition (CEEMD). A univariate GRU model and a multivariable GRU model were employed to forecast trend and stochastic displacements, respectively. A multivariable GRU model was applied to predict periodic displacement, and another two popular ML methods-long short-term memory neural networks (LSTM) and random forest (RF)-were used for comparison. Precipitation, reservoir level, and previous displacement were considered to be candidate-triggering factors for inputs of the models. The Baijiabao landslide, located in the Three Gorges Reservoir Area (TGRA), was taken as a case study to test the prediction ability of the model. The results demonstrated that the GRU algorithm provided the most encouraging results. Such a satisfactory prediction accuracy of the GRU algorithm depends on its ability to fully use the historical information while having fewer hyperparameters to train. It is concluded that the proposed model can be a valuable tool for predicting the displacements of landslides in the TGRA and other dam reservoirs.

Temperature Prediction of Seasonal Frozen Subgrades Based on CEEMDAN-LSTM Hybrid Model.

Improving the temperature prediction accuracy for subgrades in seasonally frozen regions will greatly help improve the understanding of subgrades' thermal states. Due to the nonlinearity and non-stationarity of the temperature time series of subgrades, it is difficult for a single general neural network to accurately capture these two characteristics. Many hybrid models have been proposed to more accurately forecast the temperature time series. Among these hybrid models, the CEEMDAN-LSTM model is promising, thanks to the advantages of the long short-term memory (LSTM) artificial neural network, which is good at handling complex time series data, and its combination with the broad applicability of the complete ensemble empirical mode decomposition with adaptive noise (CEEMDAN) in the field of signal decomposition. In this study, by performing empirical mode decomposition (EMD), ensemble empirical mode decomposition (EEMD), and CEEMDAN on temperature time series, respectively, a hybrid dataset is formed with the corresponding time series of volumetric water content and frost heave, and finally, the CEEMDAN-LSTM model is created for prediction purposes. The results of the performance comparisons between multiple models show that the CEEMDAN-LSTM model has the best prediction performance compared to other decomposed LSTM models because the composition of the hybrid dataset improves predictive ability, and thus, it can better handle the nonlinearity and non-stationarity of the temperature time series data.

Non-Contact Detection of Vital Signs Based on Improved Adaptive EEMD Algorithm (July 2022).

Non-contact vital sign detection technology has brought a more comfortable experience to the detection process of human respiratory and heartbeat signals. Ensemble empirical mode decomposition (EEMD) is a noise-assisted adaptive data analysis method which can be used to decompose the echo data of frequency modulated continuous wave (FMCW) radar and extract the heartbeat and respiratory signals. The key of EEMD is to add Gaussian white noise into the signal to overcome the mode aliasing problem caused by original empirical mode decomposition (EMD). Based on the characteristics of clutter and noise distribution in public places, this paper proposed a static clutter filtering method for eliminating ambient clutter and an improved EEMD method based on stable alpha noise distribution. The symmetrical alpha stable distribution is used to replace Gaussian distribution, and the improved EEMD is used for the separation of respiratory and heartbeat signals. The experimental results show that the static clutter filtering technology can effectively filter the surrounding static clutter and highlight the periodic moving targets. Within the detection range of 0.5 m~2.5 m, the improved EEMD method can better distinguish the heartbeat, respiration, and their harmonics, and accurately estimate the heart rate.

Multi-scale causality analysis between COVID-19 cases and mobility level using ensemble empirical mode decomposition and causal decomposition.

The global spread of the coronavirus disease 2019 (COVID-19) pandemic has affected the world in many ways. Due to the communicable nature of the disease, it is difficult to investigate the causal reason for the epidemic's spread sufficiently. This study comprehensively investigates the causal relationship between the spread of COVID-19 and mobility level on a multi time-scale and its influencing factors, by using ensemble empirical mode decomposition (EEMD) and the causal decomposition approach. Linear regression analysis investigates the significance and importance of the influential factors on the intrastate and interstate causal strength. The results of an EEMD analysis indicate that the mid-term and long-term domain portrays the macroscopic component of the states' mobility level and COVID-19 cases, which represents overall intrinsic characteristics. In particular, the mobility level is highly associated with the long-term variations of COVID-19 cases rather than short-term variations. Intrastate causality analysis identifies the significant effects of median age and political orientation on the causal strength at a specific time-scale, and some of them cannot be identified from the existing method. Interstate causality results show a negative association with the interstate distance and the positive one with the airline traffic in the long-term domain. Clustering analysis confirms that the states with the higher the gross domestic product and the more politically democratic tend to more adhere to social distancing. The findings of this study can provide practical implications to the policymakers that whether the social distancing policies are effectively working or not should be monitored by long-term trends of COVID-19 cases rather than short-term.

CEGH: A Hybrid Model Using CEEMD, Entropy, GRU, and History Attention for Intraday Stock Market Forecasting.

Intraday stock time series are noisier and more complex than other financial time series with longer time horizons, which makes it challenging to predict. We propose a hybrid CEGH model for intraday stock market forecasting. The CEGH model contains four stages. First, we use complete ensemble empirical mode decomposition (CEEMD) to decompose the original intraday stock market data into different intrinsic mode functions (IMFs). Then, we calculate the approximate entropy (ApEn) values and sample entropy (SampEn) values of each IMF to eliminate noise. After that, we group the retained IMFs into four groups and predict the comprehensive signals of those groups using a feedforward neural network (FNN) or gate recurrent unit with history attention (GRU-HA). Finally, we obtain the final prediction results by integrating the prediction results of each group. The experiments were conducted on the U.S. and China stock markets to evaluate the proposed model. The results demonstrate that the CEGH model improved forecasting performance considerably. The creation of a collaboration between CEEMD, entropy-based denoising, and GRU-HA is our major contribution. This hybrid model could improve the signal-to-noise ratio of stock data and extract global dependence more comprehensively in intraday stock market forecasting.

Multiscale Analysis of Solar Loading Thermographic Signals for Wall Structure Inspection.

Infrared thermography has been widely adopted in many applications for material structure inspection, where data analysis methods are often implemented to elaborate raw thermal data and to characterize material structural properties. Herein, a multiscale thermographic data analysis framework is proposed and applied to building structure inspection. In detail, thermograms are first collected by conducting solar loading thermography, which are then decomposed into several intrinsic mode functions under different spatial scales by multidimensional ensemble empirical mode decomposition. At each scale, principal component analysis (PCA) is implemented for feature extraction. By visualizing the loading vectors of PCA, the important building structures are highlighted. Compared with principal component thermography that applies PCA directly to raw thermal data, the proposed multiscale analysis method is able to zoom in on different types of structural features.

A Hydraulic Pump Fault Diagnosis Method Based on the Modified Ensemble Empirical Mode Decomposition and Wavelet Kernel Extreme Learning Machine Methods.

To address the problem that the faults in axial piston pumps are complex and difficult to effectively diagnose, an integrated hydraulic pump fault diagnosis method based on the modified ensemble empirical mode decomposition (MEEMD), autoregressive (AR) spectrum energy, and wavelet kernel extreme learning machine (WKELM) methods is presented in this paper. First, the non-linear and non-stationary hydraulic pump vibration signals are decomposed into several intrinsic mode function (IMF) components by the MEEMD method. Next, AR spectrum analysis is performed for each IMF component, in order to extract the AR spectrum energy of each component as fault characteristics. Then, a hydraulic pump fault diagnosis model based on WKELM is built, in order to extract the features and diagnose faults of hydraulic pump vibration signals, for which the recognition accuracy reached 100%. Finally, the fault diagnosis effect of the hydraulic pump fault diagnosis method proposed in this paper is compared with BP neural network, support vector machine (SVM), and extreme learning machine (ELM) methods. The hydraulic pump fault diagnosis method presented in this paper can diagnose faults of single slipper wear, single slipper loosing and center spring wear type with 100% accuracy, and the fault diagnosis time is only 0.002 s. The results demonstrate that the integrated hydraulic pump fault diagnosis method based on MEEMD, AR spectrum, and WKELM methods has higher fault recognition accuracy and faster speed than existing alternatives.