For further improvements of gravity field mod- els based on Gravity Recovery and Climate Experiment (GRACE) observations, it is necessary to identify the error sources within the recovery process. Observation residuals obtained during the gravity field recovery contain most of the measurement and modeling errors and thus can be con- sidered a realization of actual errors.In this work, we investigate the ability of wavelets to help in identifying specific error sources in GRACE range-rate residuals. The multiresolution analysis (MRA) using discrete wavelet transform (DWT) is applied to decompose the resid- ual signal into different scales with corresponding frequency bands. Temporal, spatial, and orbit-related features of each scale are then extracted for further investigations.The wavelet analysis has proven to be a practical tool to find the main error contributors. Besides the previously known sources such as K-band ranging (KBR) system noise and systematic attitude variations, this method clearly shows effects which the classic spectral analysis is hardly able or unable to represent. These effects include long-term signa- tures due to satellite eclipse crossings and dominant ocean tide errors.