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

A multiphase seismic investigation of the shallow subduction zone, southern North Island, New Zealand

The shallow structure of the Hikurangi margin, in particular the interface between the Australian Plate and the subducting Pacific Plate, is investigated using the traveltimes of direct and converted seismic phases from local earthquakes. Mode conversions take place as upgoing energy from earthquakes in the subducted slab crosses the plate interface. These PS and SP converted arrivals are observed as intermediate phases between the direct P and S waves. They place an additional constraint on the depth of the interface and enable the topography of the subducted plate to be mapped across the region. 301 suitable earthquakes were recorded by the Leeds (Tararua) broad-band seismic array, a temporary line of three-component short-period stations, and the permanent stations of the New Zealand national network. This provided coverage across the land area of southern North Island, New Zealand, at a total of 17 stations. Rays are traced through a structure parametrized using layered B-splines and the traveltime residuals inverted, simultaneously, for hypocentre relocation, interface depth and seismic velocity. The results are consistent with sediment in the northeast of the study region and gentle topography on the subducting plate. This study and recent tectonic reconstructions of the southwest Pacific suggest that the subducting plate consists of captured, oceanic crust. The anomalous nature of this crust partly accounts for the unusual features of the Hikurangi margin, e.g. the shallow trench, in comparison with the subducting margin further north

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

A three-dimensional crustal seismic velocity model for southern California from a composite event method

We present a new crustal seismic velocity model for southern California derived from P and S arrival times from local earthquakes and explosions. To reduce the volume of data and ensure a more uniform source distribution, we compute "composite event" picks for 2597 distributed master events that include pick information for other events within spheres of 2 km radius. The approach reduces random picking error and maximizes the number of S wave picks. To constrain absolute event locations and shallow velocity structure, we also use times from controlled sources, including both refraction shots and quarries. We implement the SIMULPS tomography algorithm to obtain three-dimensional (3-D) V_p and V_p/V_s structure and hypocenter locations of the composite events. Our new velocity model in general agrees with previous studies, resolving low-velocity features at shallow depths in the basins and some high-velocity features in the midcrust. Using our velocity model and 3-D ray tracing, we relocate about 450,000 earthquakes from 1981 to 2005. We observe a weak correlation between seismic velocities and earthquake occurrence, with shallow earthquakes mostly occurring in high P velocity regions and midcrustal earthquakes occurring in low P velocity regions. In addition, most seismicity occurs in regions with relatively low V_p/V_s ratios, although aftershock sequences following large earthquakes are often an exception to this pattern

High fidelity: extra-pair fertilisations in eight Charadrius plover species are not associated with parental relatedness or social mating system

Extra-pair paternity is a common reproductive strategy in many bird species. However, it remains unclear why extra-pair paternity occurs and why it varies among species and populations. Plovers (Charadrius spp.) exhibit considerable variation in reproductive behaviour and ecology, making them excellent models to investigate the evolution of social and genetic mating systems. We investigated inter- and intra-specific patterns of extra-pair parentage and evaluated three major hypotheses explaining extra-pair paternity using a comparative approach based on the microsatellite genotypes of 2049 individuals from 510 plover families sampled from twelve populations that constituted eight species. Extra-pair paternity rates were very low (0 to 4.1% of chicks per population). No evidence was found in support of the sexual conflict or genetic compatibility hypotheses, and there was no seasonal pattern of extra-pair paternity (EPP). The low prevalence of EPP is consistent with a number of alternative hypotheses, including the parental investment hypothesis, which suggests that high contribution to care by males restricts female plovers from engaging in extra-pair copulations. Further studies are needed to critically test the importance of this hypothesis for mate choice in plovers

Polarization filtering for automatic picking of seismic data and improved converted phase detection

Data-adaptive polarization filtering is used to improve the detection of converted seismic phases. Both direct waves and mode-converted PS and SP arrivals may be more easily picked on the filtered records. An autopicking routine is applied that cuts the polarization filtered traces according to the modelled traveltime of each phase through an initial structure. Use of forward-modelled, source–receiver times reduces the likelihood of an automatic pick being incorrectly made on spurious spikes in the polarization filtered trace. It is therefore a realistic way of automatically picking multiphase data sets or, more generally, linearly polarized phases where low signal-to-noise ratios may be encountered. The method is suitable for any three-component seismic data and is here applied to local earthquakes recorded in North Island, New Zealand. Intermediate energy is observed between the direct P and S arrivals due to phase conversion at the interface between the Indo-Australian and subducting Pacific plates. The amplitudes of these converted arrivals are often too low for them to be identified above the P-wave coda but polarization filtering of the records enables the yield of converted phase picks to be greatly increased

Chemical fingerprints encode mother-offspring similarity, colony membership, relatedness and genetic quality in fur seals

Chemical communication underpins virtually all aspects of vertebrate social life, yet remains poorly understood because of its highly complex mechanistic basis. We therefore used chemical fingerprinting of skin swabs and genetic analysis to explore the chemical cues that may underlie mother–offspring recognition in colonially breeding Antarctic fur seals. By sampling mother–offspring pairs from two different colonies, using a variety of statistical approaches and genotyping a large panel of microsatellite loci, we show that colony membership, mother–offspring similarity, heterozygosity, and genetic relatedness are all chemically encoded. Moreover, chemical similarity between mothers and offspring reflects a combination of genetic and environmental influences, the former partly encoded by substances resembling known pheromones. Our findings reveal the diversity of information contained within chemical fingerprints and have implications for understanding mother–offspring communication, kin recognition, and mate choice

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

Source parameters of the great Sumatran megathrust earthquakes of 1797 and 1833 inferred from coral microatolls

Large uplifts and tilts occurred on the Sumatran outer arc islands between 0.5° and 3.3°S during great historical earthquakes in 1797 and 1833, as judged from relative sea level changes recorded by annually banded coral heads. Coral data for these two earthquakes are most complete along a 160-km length of the Mentawai islands between 3.2° and 2°S. Uplift there was as great as 0.8 m in 1797 and 2.8 m in 1833. Uplift in 1797 extended 370 km, between 3.2° and 0.5°S. The pattern and magnitude of uplift imply megathrust ruptures corresponding to moment magnitudes (M_w) in the range 8.5 to 8.7. The region of uplift in 1833 ranges from 2° to at least 3.2°S and, judging from historical reports of shaking and tsunamis, perhaps as far as 5°S. The patterns and magnitude of uplift and tilt in 1833 are similar to those experienced farther north, between 0.5° and 3°N, during the giant Nias-Simeulue megathrust earthquake of 2005; the outer arc islands rose as much as 3 m and tilted toward the mainland. Elastic dislocation forward modeling of the coral data yields megathrust ruptures with moment magnitudes ranging from 8.6 to 8.9. Sparse accounts at Padang, along the mainland west coast at latitude 1°S, imply tsunami runups of at least 5 m in 1797 and 3–4 m in 1833. Tsunamis simulated from the pattern of coral uplift are roughly consistent with these reports. The tsunami modeling further indicates that the Indian Ocean tsunamis of both 1797 and 1833, unlike that of 2004, were directed mainly south of the Indian subcontinent. Between about 0.7° and 2.1°S, the lack of vintage 1797 and 1833 coral heads in the intertidal zone demonstrates that interseismic submergence has now nearly equals coseismic emergence that accompanied those earthquakes. The interseismic strains accumulated along this reach of the megathrust have thus approached or exceeded the levels relieved in 1797 and 1833

Non-monotonicity of the frictional bimaterial effect

Sliding along frictional interfaces separating dissimilar elastic materials is qualitatively different from sliding along interfaces separating identical materials due to the existence of an elastodynamic coupling between interfacial slip and normal stress perturbations in the former case. This bimaterial coupling has important implications for the dynamics of frictional interfaces, including their stability and rupture propagation along them. We show that while this bimaterial coupling is a monotonically increasing function of the bimaterial contrast, when it is coupled to interfacial shear stress perturbations through a friction law, various physical quantities exhibit a non-monotonic dependence on the bimaterial contrast. In particular, we show that for a regularized Coulomb friction, the maximal growth rate of unstable interfacial perturbations of homogeneous sliding is a non-monotonic function of the bimaterial contrast, and provide analytic insight into the origin of this non-monotonicity. We further show that for velocity-strengthening rate-and-state friction, the maximal growth rate of unstable interfacial perturbations of homogeneous sliding is also a non-monotonic function of the bimaterial contrast. Results from simulations of dynamic rupture along a bimaterial interface with slip-weakening friction provide evidence that the theoretically predicted non-monotonicity persists in non-steady, transient frictional dynamics.Comment: 14 pages, 5 figure