  Most processes in nature are coupled; however, extensive null models for
generating such processes still lacks. We present a new method to generate two
coupled Gaussian stochastic processes with arbitrary correlation functions.
This method is developed by modifying the Fourier filtering method. The
robustness of this method is proved by generating two coupled fractional
Brownian motions and extending its range of application to Gaussian random
fields.
