## 1 Introduction  
Recently, multicomponent seismology has gained considerable attention due to providing more seismic wavefield information for imaging complex structures and predicting reservoirs [1]. Generally, three-component (3C) geophones simultaneously record two horizontal components and one vertical component of the incident seismic wave. Ground roll is coherent noise with characteristics of strong energy, low frequency, low velocity and elliptical polarization, which shield the seismic reflection wave from the middle-deep layer in near offset. Ground roll attenuation is a key step in seismic data processing. An adaptive matched filter is conventionally used to process each component separately for attenuating ground roll. This method consists of a prediction step and a robust subtraction step. Ground roll is predicted by data-driven interferometry method [2] or model-driven semi-analytic modeling [3] methods. Then, an optimization subtraction operation is used to adaptively separate the ground roll [2]. Filter methods in transform domain are also commonly used, such as the frequency-wavenumber (F-K) filter, Radon transform [4], Curvelet-domain filter [5] and time-frequency-domain filter [6], to attenuate the ground roll based on the different characteristics of frequency, velocity and polarization between the ground roll and seismic reflection waves. However, when the acquisition geometry is under sampled in space, the aliasing problem can decrease the effect of these methods. In land multicomponent seismic exploration, 3C seismic data are recorded and processed as a vector, which provides not only a robust analysis of each individual component but also valuable information about the coherency components [7, 8].  
Empirical mode decomposition (EMD) proposed by [PERSON] et al. [9] decomposes a signal into a series of intrinsic mode functions (IMFs). Each IMF has a relatively local constant frequency ranging from high frequency to low frequency. The EMD method does not require predetermined base functions and can adaptively separate nonlinear and non-stationary signals. EMD and its extensions have good performance on random and ground roll attenuation [10; 11; 12]. [PERSON] et al. [13] used complete ensemble empirical mode decomposition (CEEMD) to improve the standard EMD method and obtained better results when they were applied to attenuate the ground roll.  
The standard EMD method is conventionally conducted on each component separately, which results in misaligned IMFs for multicomponent data. The ground roll attenuation based on the standard EMD method damages the vector coherence of seismic data. Combining multiple components with a higher-dimensional signal could circumvent the mode alignment problem caused by the standard EMD method [14]. [PERSON] and [PERSON] [15] proposed a multivariate empirical mode decomposition (MEMD) method to process multicomponent signals and obtained the IMFs of aligned frequency range. These IMFs match well in the frequency scale and number properties for each component of the multicomponent data. In order to solve the mode mixing effect in standard EMD, extended EMD methods [13; 16] have been proposed. [PERSON] and [PERSON] [17] proposed noise-assisted MEMD (NA-MEMD) to improve MEMD by adding extra components containing independent white noise to the original multivariate signal. This method helps overcome the mode-mixing problem existing in the extracted IMFs.  
In this paper, we propose an alternative ground roll attenuation method for multicomponent seismic data in time domain. First, we start this article with a description of the MEMD and NA-MEMD methods. Then, we demonstrate the principle of the proposed method for ground roll attenuation of multicomponent seismic data. Finally, the synthetic and field seismic data tests demonstrate that the NA-MEMD method can effectively attenuate the ground roll. Compared to the F-K filter method, the proposed method preserves more low-frequency content of the seismic reflection wave and coherency information between components.