Abstract: A combination of the wave superposition method and the modified spectral representation method is applied to a wind field simulation for whole bridge including the deck and the towers. In wave superposition method, the variation of the auto-spectrum with respect to the elevation and the quadrature spectrum caused by asymmetric cross-correlation function as well as the frequency dependent phase lag are all taken into account. A practical decomposition method is presented for the case which involves complex target cross-spectral density matrix. In a spectrum comparison, it is found that the correlation functions of the simulated wind field generated by the above methods agree well with the target correlation functions. This indicates that the method presented in this paper is reliable. Some of the stochastic wind field samples of a long-span cable-stayed bridge are given.