  One of the most efficient ways to produce unconditional simulations is with
the spectral method using fast Fourier transform (FFT) [1]. But this approach
is not applicable to arbitrary surfaces because no regular grid exists.
However, points on the arbitrary surface can be generated randomly using
uniform distribution to replace a regular grid. This paper will describe a
nonstationary kernel convolution approach for data on arbitrary surfaces.
