[1] Forty-one wideband magnetotelluric (MT) soundings were collected in a 150-kmlong transect across the Southern Alps of the central South Island of New Zealand, an active compressional orogen. Decomposed MT impedance tensors, vertical magnetic field relations, and reconnaissance soundings at two locations off line imply an approximately two-dimensional geometry here with average regional geoelectric strike of $N40°E, similar to surface geologic trends. Two independent, two-dimensional inversion algorithms were applied to the MT data, and both imply a concave-upward (U-shaped), middle to lower crustal conductive zone beneath the west central portion of the island. The average conductivity of this zone in the strike direction appears to be much higher than that required across strike and may represent anisotropy or along-strike conductive strands narrower than the transverse magnetic (cross-strike) mode MT data can resolve. The deep crustal conductor under the Southern Alps is interpreted to represent mainly a volume of fluids arising from prograde metamorphism within a thickening crust. Fluid interconnection and electrical conduction are promoted by shear deformation. The conductor rises northwestward toward the trace of the Alpine Fault but attains a nearvertical configuration at a depth of $10 km and reaches close to the surface 5-10 km inland of the fault trace itself. The transition to vertical orientation at this depth is interpreted to occur as fluids ascend across the brittle-ductile transition in uplifting schist and approach the surface through induced hydrofractures. The high-grade schist becomes resistive after depletion of fluids and continues to extrude toward the Alpine Fault. Shallow extensions of the deep high conductivity are coincident with modern, hydrothermal veining and gold mineralization interpreted to be of deep crustal provenance. To the southeast, high conductivity also reaches the surface coincident with a major back thrust fault zone of the doubly vergent Southern Alps orogen, which also exhibits evidence for expulsion of high-temperature fluids. The higher conductivity inferred along strike (possible anisotropy) could reflect more efficient fluid interconnection in this higher-strain direction, as well as possible contributions by sheared, fluid-deposited graphite. Conductivity of the uppermost mantle of the South Island is low, consistent with advection of cold mantle lithosphere into the underlying asthenosphere as suggested by P wave delay studies.

No abstract

Convergent margin volcanism originates with partial melting, primarily of the upper mantle, into which the subducting slab descends. Melting of this material can occur in one of two ways. The flow induced in the mantle by the slab can result in upwelling and melting through adiabatic decompression. Alternatively, fluids released from the descending slab through dehydration reactions can migrate into the hot mantle wedge, inducing melting by lowering the solidus temperature. The two mechanisms are not mutually exclusive. In either case, the buoyant melts make their way towards the surface to reside in the crust or to be extruded as lava. Here we use magnetotelluric data collected across the central state of Washington, USA, to image the complete pathway for the fluid-melt phase. By incorporating constraints from a collocated seismic study into the magnetotelluric inversion process, we obtain superior constraints on the fluids and melt in a subduction setting. Specifically, we are able to identify and connect fluid release at or near the top of the slab, migration of fluids into the overlying mantle wedge, melting in the wedge, and transport of the melt/fluid phase to a reservoir in the crust beneath Mt Rainier.

Summary We report herein on a finite element algorithm for 2‐D magnetotelluric modelling which solves directly for secondary variations in the field parallel to strike, plus the subsequent vertical and transverse auxiliary fields, for both transverse electric and transverse magnetic modes. The governing Helmholtz equations for the secondary fields along strike are the same as those for total field algorithms with the addition of source terms involving the primary fields and the conductivity difference between the body and the host. Our approach has overcome a difficulty with numerical accuracy at low frequencies observed in total field solutions with 32‐bit arithmetic for the transverse magnetic mode especially, but also for the transverse electric mode. Matrix ill‐conditioning, which affects total field solutions, increases with the number of element rows with the square of the maximum element aspect ratio and with the inverse of frequency. In the secondary formulation, the field along strike and the auxiliary fields do not need to be extracted in the face of an approximately computed primary field which increasingly dominates the total field solution towards low frequencies. In addition to low‐frequency stability, the absolute accuracy of our algorithm is verified by comparison with the TM and the TE mode analytic responses of a segmented overburden model.

Newly forming subduction zones on Earth can provide insights into the evolution of major fault zone geometries from shallow levels to deep in the lithosphere and into the role of fluids in element transport and in promoting rock failure by several modes. The transpressional subduction regime of New Zealand, which is advancing laterally to the southwest below the Marlborough strike-slip fault system of the northern South Island, is an ideal setting in which to investigate these processes. Here we acquired a dense, high-quality transect of magnetotelluric soundings across the system, yielding an electrical resistivity cross-section to depths beyond 100 km. Our data imply three distinct processes connecting fluid generation along the upper mantle plate interface to rock deformation in the crust as the subduction zone develops. Massive fluid release just inland of the trench induces fault-fracture meshes through the crust above that undoubtedly weaken it as regional shear initiates. Narrow strike-slip faults in the shallow brittle regime of interior Marlborough diffuse in width upon entering the deeper ductile domain aided by fluids and do not project as narrow deformation zones. Deep subduction-generated fluids rise from 100 km or more and invade upper crustal seismogenic zones that have exhibited historic great earthquakes on high-angle thrusts that are poorly oriented for failure under dry conditions. The fluid-deformation connections described in our work emphasize the need to include metamorphic and fluid transport processes in geodynamic models.