A Gaussian random process with a given power spectral density (PSD) function can be modeled as a sum of sinusoids (SOS), and has been widely used to simulate Rayleigh-fading communication channels. In practice, a finite number of sinusoids can be used to approximate a Gaussian process, which reduces complexity. For this reason, the method of sum of sinusoids (SOS) to simulate fading channel has been extensively investigated [7], [8], [10]-[12], [14], [17]. In [11], the method of SOS is studied in more detail and several approaches to determining the sinusoidal frequencies are developed. The very popular Jakes' channel simulation model calculates the frequencies of the sinusoids based on a physical channel model with uniformly distributed scatterrers around a circle. The conventional one-dimensional (1-D) channel model cannot capture the spatial correlation of shadowing processes. Here we develop a two-dimensional (2-D) SOS-based channel model to simulate the shadowing process.