A new method is presented that simulates nonstationary multipoint earthquake accelerograms at local field points, when the observed ground motions are specified at several locations. This process effectively avoids the stochastic decomposition of the spectral matrix. In addition, it does not use the target spectrum and spatial coherence function to define the information of each point. First, the observed ground motions are reconstructed using an autoregressive model. Second, a new cross-regressive model is proposed based on the formulation of the auto-regressive model. The cross-regressive model can reconstruct different observed ground motions by using recordings from any nearby area. To consider the nonstationary characteristics and spatial coherence, the ground motions of unobserved points between two observed points can be generated using an interpolation algorithm. This new method simulates the ground motions of target points in the Chiba array using partially observed earthquake accelerograms. The simulated earthquake accelerograms satisfy the nonstationary characteristics and spatial coherence of the Chiba array. This study is promising for engineering applications and hazard analyses.