Closed-form solutions are analytically derived for stochastic properties of earthquake ground motion fields, which are conditioned by an observed time series at certain observation sites and are characterized by spectra with uncertainties. The theoretical framework presented here can estimate not only the expectations of such simulated earthquake ground motions, but also the prediction errors which offer important information for the field of engineering. Before these derivations are made, the theory of conditional random fields is summarized for convenience in this study. Furthermore, a method for stochastic interpolation of power spectra is explained.