This paper proposes the multiresolution mode decomposition (MMD) as a novel model for adaptive time series analysis. The main conceptual innovation is the introduction of the multiresolution intrinsic mode function (MIMF) of the formto model nonlinear and non-stationary data with time-dependent amplitudes, frequencies, and waveforms. The multiresolution expansion coefficients {a n }, {b n }, and the shape function series {s cn (t)} and {s sn (t)} provide innovative features for adaptive time series analysis. For complex signals that are a superposition of several MIMFs with well-differentiated phase functions φ(t), a new recursive scheme based on Gauss-Seidel iteration and diffeomorphisms is proposed to identify these MIMFs, their multiresolution expansion coefficients, and shape function series. Numerical examples from synthetic data and natural phenomena are given to demonstrate the power of this new method.