© 2015 World Scientific Publishing Company. Surrogate data methods have been widely applied to produce synthetic data, while maintaining the same statistical properties as the original. By using such methods, one can analyze certain properties of time series. In this context, Theiler's surrogate data methods are the most commonly considered approaches. These are based on the Fourier transform, limiting them to be applied only on stationary time series. Consequently, time series including nonstationary behavior, such as trend, produces spurious high frequencies with Theiler's methods, resulting in inconsistent surrogates. To solve this problem, we present two new methods that combine time series decomposition techniques and surrogate data methods. These new methods initially decompose time series into a set of monocomponents and the trend. Afterwards, traditional surrogate methods are applied on those individual monocomponents and a set of surrogates is obtained. Finally, all individual surrogates plus the trend signal are combined in order to create a single surrogate series. Using this method, one can investigate linear and nonlinear Gaussian processes in time series, irrespective of the presence of nonstationary behavior.