Abstract:In recent decades, decomposition techniques have enabled increasingly more applications for dimension reduction, as well as extraction of additional information from geophysical time series. Traditionally, the principal component analysis (PCA)/empirical orthogonal function (EOF) method and more recently the independent component analysis (ICA) have been applied to extract, statistical orthogonal (uncorrelated), and independent modes that represent the maximum variance of time series, respectively. PCA and ICA can be classified as stationary signal decomposition techniques since they are based on decomposing the autocovariance matrix and diagonalizing higher (than two) order statistical tensors from centered time series, respectively. However, the stationarity assumption in these techniques is not justified for many geophysical and climate variables even after removing cyclic components, e.g., the commonly removed dominant seasonal cycles. In this paper, we present a novel decomposition method, the complex independent component analysis (CICA), which can be applied to extract non-stationary (changing in space and time) patterns from geophysical time series. Here, CICA is derived as an extension of real-valued ICA, where (a) we first define a new complex dataset that contains the observed time series in its real part, and their Hilbert transformed series as its imaginary part, (b) an ICA algorithm based on diagonalization of fourth-order cumulants is then applied to decompose the new complex dataset in (a), and finally, (c) the dominant independent complex modes are extracted and used to represent the dominant space and time amplitudes and associated phase propagation patterns. The performance of CICA is examined by analyzing synthetic data constructed from multiple physically meaningful modes in a simulation framework, with known truth. Next, global terrestrial water storage (TWS) data from the Gravity Recovery And Climate Experiment (GRACE) gravimetry mission (2003–2016), and satellite radiometric sea surface temperature (SST) data (1982–2016) over the Atlantic and Pacific Oceans are used with the aim of demonstrating signal separations of the North Atlantic Oscillation (NAO) from the Atlantic Multi-decadal Oscillation (AMO), and the El Niño Southern Oscillation (ENSO) from the Pacific Decadal Oscillation (PDO). CICA results indicate that ENSO-related patterns can be extracted from the Gravity Recovery And Climate Experiment Terrestrial Water Storage (GRACE TWS) with an accuracy of 0.5–1 cm in terms of equivalent water height (EWH). The magnitude of errors in extracting NAO or AMO from SST data using the complex EOF (CEOF) approach reaches up to ~50% of the signal itself, while it is reduced to ~16% when applying CICA. Larger errors with magnitudes of ~100% and ~30% of the signal itself are found while separating ENSO from PDO using CEOF and CICA, respectively. We thus conclude that the CICA is more effective than CEOF in separating non-stationary patterns.

Gravity signal extracting method based on independent component analysis with multiple reference signals

Independent Component Analysis(ica) for Denoising EMG Signal from ECG Signal

Study of Noise Removal Technique Based on Independent Component Analysis

Vibration sources identification with independent component analysis

A Fast Interpretation Method of Gravity Gradiometry Data Based on Magnetic Dipole Localization

Real-Time Joint Filtering of Gravity and Gravity Gradient Data Based on Improved Kalman Filter

Error analysis and filtering methods for absolute ocean gravity data

Geophysical Potential Field Anomaly Separation Method with Optimal Mother Wavelet and Spatial Locating Multiresolution Analysis (MRA)

Inversion of Gravity Anomaly Based on AlexNet Deep Neural Network

Noise Processing of Ocean Absolute Gravity Data Based on Savitzky-golay Filtering Algorithm

Developing a Complex Independent Component Analysis (CICA) Technique to Extract Non-stationary Patterns from Geophysical Time Series

Independent Component Analysis for Filtering Airwaves in Seabed Logging Application

Application research on adaptive Kalman filtering for airborne gravity anomaly determination

Methods for removal of unwanted signals from gravity time-series: Comparison using linear techniques complemented with analysis of system dynamics

An Atomic Gravimeter Dynamic Measurement Method Based on Kalman Filter

Unsupervised Noise Reductions for Gravitational Reference Sensors or Accelerometers Based on the Noise2Noise Method

A Cascade Gravity Matching Algorithm with Large Initial Position Error

A Fusion Algorithm of Underwater Dual-sensor Gravimeter

Unsupervised noise reductions for gravitational reference sensors or accelerometers based on Noise2Noise method

Scattering Suppression in Underwater LIDAR Signal Processing Based on Wavelet-ICA Method

Noise filtering of full-gravity gradient tensor data