Distribution of seismic anisotropy in the subduction zone beneath the Wellington region, New Zealand

Shear wave splitting measurements from S arrivals of local earthquakes recorded at the Incorporated Research Institutions for Seismology (IRIS) broadband sensor SNZO are used to determine a basic anisotropic structure for the subduction zone in the Wellington region. With the use of high-frequency filters, fast anisotropic polarization ( φ ) and splitting time ( δt ) measurements typical of crustal anisotropy are evident, but the larger splitting expected from the mantle is often not resolved. The small splitting seen agrees well with the results of previous studies concerning shallow crustal anisotropy. With the use of lower-frequency filters, measurements more consistent with mantle anisotropy are made. Anisotropy of 4.4 ± 0.9 per cent with a fast polarization of 29° ± 38° is calculated for the subducting slab, from 20 to 70  km depth. Using this result in addition to the results of previous studies, a model is proposed. The model requires a frequency-dependent anisotropy of less than 1.4 per cent when measured with a period of ~2  s to be present in the sub-slab mantle.
Separate from this population, a band of events in northern Cook Strait with an 86° ± 10° fast polarization is seen. This is at about 40° from the strike of the Hikurangi margin, and suggests a source of shear strain 40° removed from that found in the majority of the region. The cause of this is probably a deformation in the subducting slab in this region, as it moves towards a greater incline to the south.     

Effects of active crustal movements on thermal structure in subduction zones

In young suduction zones we observe steady uplift of island arcs. The steady uplift of island arcs is always accompanied by surface erosion. The long duration of uplift and erosion effectively transports heat at depth to shallower parts by advection. If the rates of uplift and erosion are sufficiently large, such a process of heat transportation will strongly affect thermal structure in subduction zones. First, we quantitatively examine the effects of uplift and erosion on thermal structure by using a simple 1-D heat conduction model, based on the assumption that the initial thermal state is in equilibrium. The results show that temperature increase, Δ T  , due to uplift and erosion can be approximately evaluated by Δ T  = ν _e tβ at depth, where ν _e is the rate of uplift (erosion), t is the duration of uplift (erosion), and β is the gradient of the geotherm in the initial state. Next, considering the effects of vertical crustal movements such as uplift and erosion in island arcs and subsidence and sedimentation in ocean trenches, in addition to the effects of radioactive heat generation in the crust, frictional heating at plate boundaries and accretion of oceanic sediments to overriding continental plates, we numerically simulate the evolution process of the thermal structure in subduction zones. The result shows that the temperature beneath the island arc gradually increases as a result of uplift and erosion as plate subduction progresses. Near the ocean trench, on the other hand, the low-temperature region gradually expands as a result of sedimentation and accretion in addition to direct cooling by the cold descending slab. The surface heat flow expected from this model is low in fore-arc basins, high in island arcs and moderately high in back-arc regions.