Plane strain (2D) finite element models are used to examine factors contributing to basin effects observed for multiple seismic events at sites in the Thorndon basin of Wellington, New Zealand. The models consider linear elastic soil and rock response when subjected to vertically-propagating shear waves. Depth-dependent shear wave velocities are considered in the soil layers, and the effects of random variations of soil velocity within layers are modelled. Various rock shear wave velocity configurations are considered to evaluate their effect on the modelled surficial response. It is shown that these simple 2D models are able to capture basin reverberations and compare more favourably to observations from strong motion recordings than conventional 1D site response models. It is also shown that consideration of a horizontal impedance contrast across the Wellington Fault affects spectral response and amplification at longer periods, suggesting the importance of this feature in future ground motion modelling studies in the Wellington region.