Common procedures used to account for spatial variability of shear wave velocity (Vs) in one-dimensional (1D) ground response analyses (GRAs), such as stochastic randomization of Vs or increasing small-strain damping, have been shown to improve seismic site response predictions relative to 1D GRAs where no attempts are made to account for spatial variability. However, even after attempting to account for spatial variability using common procedures, 1D GRAs often still yield results that are different than ground motions recorded at many downhole array sites. When 1D predictions differ from observations, the site is typically considered to be too spatially variable to effectively use 1D GRAs. While there is no doubt that some sites are indeed too variable for 1D GRAs, it is also possible that simple 1D analyses could still be effectively used at many sites if spatial variability is accounted for via a more rational, site-specific approach. In this study, an H/V geostatistical approach for building pseudo-3D Vs models is implemented to account for spatial variability in 1D GRAs. The geostatistical approach is used to generate a uniform grid of Vs profiles that have been scaled to match fundamental site frequency estimates from horizontal-to-vertical spectral ratio (H/V) noise measurements. In this paper, 1D GRAs are performed for each grid-point and the results are statistically combined to reflect the average site response and its variability. This 1D application is demonstrated at the Treasure Island and Delaney Park Downhole Array sites, where it is shown to produce superior fits to the small-strain recorded site response relative to existing approaches used to account for spatial variability in 1D GRAs. Using the proposed approach, we also investigate the lateral area that is likely influencing site response at each site and show that it could extend to significant distances (as much as 1 km) from the boreholes.