Three-dimensional lithospheric structure below the New Zealand Southern Alps

Uppermost mantle seismic structure below the Southern Alps in South Island, New Zealand, is investigated by teleseismic P wave travel time residual inversion. The three-dimensional tomographic images show a near-vertical, high-velocity (2–4%) structure in the uppermost mantle that directly underlies thickened crust along the NNESSW axis of the Southern Alps. The center of the high-velocity anomaly lies to the east of the Alpine fault which bounds Pacific and Australian plate rocks. The oblique collision of these plates resulted in the uplift of the Southern Alps during the past 5–7 m.y. Also, a high-velocity anomaly (3–5%) corresponding to the Hikurangi subduction zone lies to the northeast of the Southern Alps anomaly, and low-velocity anomalies (-3%) underlying parts of northwestern and southern South Island may be signatures of late Tertiary extension and volcanism. The data consist of teleseismic arrival times from the New Zealand National Seismograph Network and arrival times recorded during the 1995–1996 Southern Alps Passive Seismic Experiment. Crustal heterogeneity was accounted for by back projecting the rays through an independently obtained three-dimensional crustal velocity and Moho depth model. The Southern Alps uppermost mantle velocity anomalies are most simply explained by lithospheric thickening below the center of convergence accompanied by thinning and asthenospheric upwelling adjacent to the region of convergence

Creating realistic models based on combined forward modeling and tomographic inversion of seismic profiling data

Amplitudes and shapes of seismic patterns derived from tomographic images often are strongly biased with respect to real structures in the earth. In particular, tomography usually provides continuous velocity distributions, whereas major velocity changes in the earth often occur on first-order interfaces. We propose an approach that constructs a realistic structure of the earth that combines forward modeling and tomographic inversion (FM&TI). Using available a priori information, we first construct a synthetic model with realistic patterns. Then we compute synthetic times and invert them using the same tomographic code and the same parameters as in the case of observed data processing. We compare the reconstruction result with the tomographicimage of observed data inversion. If a discrepancy is observed, we correct the synthetic model and repeat the FM&TI process. After several trials, we obtain similar results of synthetic and observed data inversion. In this case, the derived synthetic model adequately represents the real structure of the earth. In a working scheme of this approach, we three authors used two different synthetic models with a realistic setup. One of us created models, but the other two performed the reconstruction with no knowledge of the models. We discovered that the synthetic models derived by FM&TI were closer to the true model than the tomographic inversion result. Our reconstruction results from modeling marine data acquired in the Musicians Seamount Province in the Pacific Ocean indicate the capacity and limitations of FM&TI

Intermediate-Depth Earthquakes in a Region of Continental Convergence: South Island, New Zealand

It is rare to find earthquakes with depths greater than 30 km in continent–continent collision zones because the mantle lithosphere is usually too hot to enable brittle failure. However, a handful of small, intermediate-depth earthquakes (30–97 km) have been recorded in the continental collision region in central South Island, New Zealand. The earthquakes are not associated with subduction but all lie within or on the margins of thickened crust or uppermost mantle seismic high-velocity anomalies. The largest of the earthquakes has M_L 4.0 corresponding to a rupture radius of between 100 and 800 m, providing bounds on the upper limit to the rupture length over which brittle failure is taking place in the deep brittle–plastic transition zone. The earthquake sources may be controlled by large shear strain gradients associated with viscous deformation processes in addition to depressed geotherms

Internal deformation of the subducted Nazca slab inferred from seismic anisotropy

Within oceanic lithosphere a fossilized fabric is often preserved originating from the time of plate formation. Such fabric is thought to form at the mid-ocean ridge when olivine crystals align with the direction of plate spreading1, 2. It is unclear, however, whether this fossil fabric is preserved within slabs during subduction or overprinted by subduction-induced deformation. The alignment of olivine crystals, such as within fossil fabrics, can generate anisotropy that is sensed by passing seismic waves. Seismic anisotropy is therefore a useful tool for investigating the dynamics of subduction zones, but it has so far proved difficult to observe the anisotropic properties of the subducted slab itself. Here we analyse seismic anisotropy in the subducted Nazca slab beneath Peru and find that the fast direction of seismic wave propagation aligns with the contours of the slab. We use numerical modelling to simulate the olivine fabric created at the mid-ocean ridge, but find it is inconsistent with our observations of seismic anisotropy in the subducted Nazca slab. Instead we find that an orientation of the olivine crystal fast axes aligned parallel to the strike of the slab provides the best fit, consistent with along-strike extension induced by flattening of the slab during subduction (A. Kumar et al., manuscript in preparation). We conclude that the fossil fabric has been overprinted during subduction and that the Nazca slab must therefore be sufficiently weak to undergo internal deformation

Evaluating the importance of metamorphism in the foundering of continental crust

The metamorphic conditions and mechanisms required to induce foundering in deep arc crust are assessed using an example of representative lower crust in SW New Zealand. Composite plutons of Cretaceous monzodiorite and gabbro were emplaced at ~1.2 and 1.8 GPa are parts of the Western Fiordland Orthogneiss (WFO); examples of the plutons are tectonically juxtaposed along a structure that excised ~25 km of crust. The 1.8 GPa Breaksea Orthogneiss includes suitably dense minor components (e.g. eclogite) capable of foundering at peak conditions. As the eclogite facies boundary has a positive dP/dT, cooling from supra-solidus conditions (T > 950 ºC) at high-P should be accompanied by omphacite and garnet growth. However, a high monzodioritic proportion and inefficient metamorphism in the Breaksea Orthogneiss resulted in its positive buoyancy and preservation. Metamorphic inefficiency and compositional relationships in the 1.2 GPa Malaspina Pluton meant it was never likely to have developed densities sufficiently high to founder. These relationships suggest that the deep arc crust must have primarily involved significant igneous accumulation of garnet–clinopyroxene (in proportions >75%). Crustal dismemberment with or without the development of extensional shear zones is proposed to have induced foundering of excised cumulate material at P > 1.2 GPa

Seismically invisible fault zones: Laboratory insights into imaging faults in anisotropic rocks

Phyllosilicate‐rich rocks which commonly occur within fault zones cause seismic velocity anisotropy. However, anisotropy is not always taken into account in seismic imaging and the extent of the anisotropy is often unknown. Laboratory measurements of the velocity anisotropy of fault zone rocks and gouge from the Carboneras fault zone in SE Spain indicate 10–15% velocity anisotropy in the gouge and 35–50% anisotropy in the mica‐schist protolith. Greater differences in velocity are observed between the fast and slow directions in the mica‐schist rock than between the gouge and the slow direction of the rock. This implies that the orientation of the anisotropy with respect to the fault is key in imaging the fault seismically. For example, for fault‐parallel anisotropy, a significantly greater velocity contrast between fault gouge and rock will occur along the fault than across it, highlighting the importance of considering the foliation orientation in design of seismic experiments

Campylobacteriosis, Eastern Townships, Québec

Independent risk factors for campylobacteriosis (eating raw, rare, or undercooked poultry; consuming raw milk or raw milk products; and eating chicken or turkey in a commercial establishment) account for <50% of cases in Québec. Substantial regional and seasonal variations in campylobacteriosis were not correlated with campylobacter in chickens and suggested environmental sources of infection, such as drinking water

Fore-arc deformation and underplating at the northern Hikurangi margin, New Zealand

Geophysical investigations of the northern Hikurangi subduction zone northeast of New Zealand, image fore‐arc and surrounding upper lithospheric structures. A seismic velocity (Vp) field is determined from seismic wide‐angle data, and our structural interpretation is supported by multichannel seismic reflection stratigraphy and gravity and magnetic modeling. We found that the subducting Hikurangi Plateau carries about 2 km of sediments above a 2 km mixed layer of volcaniclastics, limestone, and chert. The upper plateau crust is characterized by Vp = 4.9–6.7 km/s overlying the lower crust with Vp > 7.1 km/s. Gravity modeling yields a plateau thickness around 10 km. The reactivated Raukumara fore‐arc basin is >10 km deep, deposited on 5–10 km thick Australian crust. The fore‐arc mantle of Vp > 8 km/s appears unaffected by subduction hydration processes. The East Cape Ridge fore‐arc high is underlain by a 3.5 km deep strongly magnetic (3.3 A/m) high‐velocity zone, interpreted as part of the onshore Matakaoa volcanic allochthon and/or uplifted Raukumara Basin basement of probable oceanic crustal origin. Beneath the trench slope, we interpret low‐seismic‐velocity, high‐attenuation, low‐density fore‐arc material as accreted and recycled, suggesting that underplating and uplift destabilizes East Cape Ridge, triggering two‐sided mass wasting. Mass balance calculations indicate that the proposed accreted and recycled material represents 25–100% of all incoming sediment, and any remainder could be accounted for through erosion of older accreted material into surrounding basins. We suggest that continental mass flux into the mantle at subduction zones may be significantly overestimated because crustal underplating beneath fore‐arc highs have not properly been accounted for

A meta-analysis of long-term effects of conservation agriculture on maize grain yield under rain-fed conditions

Conservation agriculture involves reduced tillage, permanent soil cover and crop rotations to enhance soil fertility and to supply food from a dwindling land resource. Recently, conservation agriculture has been promoted in Southern Africa, mainly for maize-based farming systems. However, maize yields under rain-fed conditions are often variable. There is therefore a need to identify factors that influence crop yield under conservation agriculture and rain-fed conditions. Here, we studied maize grain yield data from experiments lasting 5 years and more under rain-fed conditions. We assessed the effect of long-term tillage and residue retention on maize grain yield under contrasting soil textures, nitrogen input and climate. Yield variability was measured by stability analysis. Our results show an increase in maize yield over time with conservation agriculture practices that include rotation and high input use in low rainfall areas. But we observed no difference in system stability under those conditions. We observed a strong relationship between maize grain yield and annual rainfall. Our meta-analysis gave the following findings: (1) 92% of the data show that mulch cover in high rainfall areas leads to lower yields due to waterlogging; (2) 85% of data show that soil texture is important in the temporal development of conservation agriculture effects, improved yields are likely on well-drained soils; (3) 73% of the data show that conservation agriculture practices require high inputs especially N for improved yield; (4) 63% of data show that increased yields are obtained with rotation but calculations often do not include the variations in rainfall within and between seasons; (5) 56% of the data show that reduced tillage with no mulch cover leads to lower yields in semi-arid areas; and (6) when adequate fertiliser is available, rainfall is the most important determinant of yield in southern Africa. It is clear from our results that conservation agriculture needs to be targeted and adapted to specific biophysical conditions for improved impact