Extension and Partitioning in an Oblique Subduction Zone, New Zealand: Constraints from Three-Dimensional Numerical Modeling

Contraction, strike slip, and extension displacements along the Hikurangi margin northeast of the North Island of New Zealand coincide with large lateral gradients in material properties. We use a finite- difference code utilizing elastic and elastic-plastic rheologies to build large- scale, three-dimensional numerical models which investigate the influence of material properties on velocity partitioning within oblique subduction zones. Rheological variation in the oblique models is constrained by seismic velocity and attenuation information available for the Hikurangi margin. We compare the effect of weakly versus strongly coupled subduction interfaces on the development of extension and the partitioning of velocity components for orthogonal and oblique convergence and include the effect of ponded sediments beneath the Raukumara Peninsula. Extension and velocity partitioning occur if the subduction interface is weak, but neither develops if the subduction interface is strong. The simple mechanical model incorporating rheological variation based on seismic observations produces kinematics that closely match those published from the Hikurangi margin. These include extension within the Taupo Volcanic Zone, uplift over ponded sediments, and dextral contraction to the south.

Along-Strike Differences in the Southern Alps of New Zealand: Consequences of Inherited Variation in Rheology

Long- and short-term strain variations along the Australian-Pacific plate boundary through the South Island of New Zealand, including a 300% increase in orogen width, coexistence of oblique thrusting on orthogonal structures, and variability in the locus of orogenic gold deposits, coincide with rheologically relevant geological variation. Our model investigates the consequences of thin, strong lower crust in the north and thick, weak lower crust in the south. Solution of the full 3-D mechanical equations reproduces the larger wavelength strain patterns of the orogen. A 3-D perturbation-based analytical solution leads to the identification of the sensitivity of displacement type to minor stress changes. Transition from boundary-normal thrusting to boundary-parallel thrusting occurs at the transition from strong to weak lower crust and is related to an increase in either tau(yz) (shear stress in the yz plane) or the ratio of the coordinate normal stresses, (sigma(yy)/sigma(xx)), where x and y are in the horizontal and z is vertical. Both mechanisms are compatible with the geologically dependent rheological variation employed in our model. Citation: Upton, P., P. O. Koons, D. Craw, C. M. Henderson, and R. Enlow (2009), Along-strike differences in the Southern Alps of New Zealand: Consequences of inherited variation in rheology, Tectonics, 28, TC2007, doi:10.1029/2008TC002353.