Thin-shell finite element methods that incorporate faults, realistic rheology, laterally varying heat flow and topography, and plate velocity boundary conditions have been used to model the neotectonics of New Zealand. We find that New Zealand's faults have effective friction of ~0.17, comparable to that found in other Pacific Rim regions. The long-term average slip rate of the Alpine fault varies along strike, generally increasing northeastward until slip is partitioned among the strands of the Marlborough system. The average slip rate, ~30 mm/yr, when combined with published geodetic results and historical seismicity, strongly suggests a high probability of future large earthquakes. Tectonic deformation of North Island is controlled by a balance between differential topographic pressure and traction from the Hikurangi subduction thrust. The Hikurangi forearc is an independent plate sliver moving relative to the Pacific and Australian plates. There is a complicated zone of slip partitioning in the transition from the Alpine fault to the Puysegur trench. An offshore thrust fault, the southern segment of which may correspond to the Waipounamou fault system, parallels to the SE coast of South Island and needs to be included in seismic hazard estimates.