Two-dimensional modelling of subduction zone anisotropy with application to southwestern Japan

We present a series of 2-D numerical models of viscous flow in the mantle wedge induced by a subducting lithospheric plate. We use a kinematically defined slab geometry approximating the subduction of the Philippine Sea plate beneath Eurasia. Through finite element modelling we explore the effects of different rheological and thermal constraints (e.g. a low-viscosity region in the wedge corner, power law versus Newtonian rheology, the inclusion of thermal buoyancy forces and a temperature-dependent viscosity law) on the velocity and finite strain field in the mantle wedge. From the numerical flow models we construct models of anisotropy in the wedge by calculating the evolution of the finite strain ellipse and combining its geometry with appropriate elastic constants for effective transversely isotropic mantle material. We then predict shear wave splitting for stations located above the model domain using expressions derived from anisotropic perturbation theory, and compare the predictions to ∼500 previously published shear wave splitting measurements from seventeen stations of the broad-band F-net array located in southwestern Japan. Although the use of different model parameters can have a substantial effect on the character of the finite strain field, the effect on the average predicted splitting parameters is small. However, the variations with backazimuth and ray parameter of individual splitting intensity measurements at a given station for different models are often different, and rigorous analysis of details in the splitting patterns allows us to discriminate among different rheological models for flow in the mantle wedge. The splitting observed in southwestern Japan agrees well with the predictions of trench-perpendicular flow in the mantle wedge along with B-type olivine fabric dominating in a region from the wedge corner to about 125 km from the trench.     

Evolution of optical properties of vitrinite, sporinite and semifusinite in response to heating under inert conditions

The objective of the study was to characterize changes of reflectance, reflectance anisotropy and reflectance indicating surface (RIS) shape of vitrinite, sporinite and semifusinite subjected to thermal treatment under inert conditions. Examination was performed on vitrinite, liptinite and inertinite concentrates prepared from channel samples of steam coal (R_r = 0.70%) and coking coal (R_r = 1.25%), collected from seam 405 of the Upper Silesian Coal Basin. The concentrates were heated at temperatures of 400–1200 °C for 1 h time in an argon atmosphere.All components examined in this study: vitrinite, sporinite and semifusinite as well as matrix of vitrinite and liptinite cokes, despite of rank of their parent coal, show, in general, the most important changes of reflectance value and optical anisotropy when heated at 500 °C, 800 °C (with the exception of bireflectance value of sporinite) and 1200 °C.After heating the steam coal at 1200 °C, the vitrinite and the semifusinite reveal similar reflectances, whereas the latter a slightly stronger anisotropy. Sporinite and matrix of liptinite coke have lower reflectances but anisotropy (R_bi and R_am values) similar to those observed for vitrinite and semifusinite. However, at 1000 °C sporinite and matrix of liptinite coke have the highest reflectivity of the studied components. The RIS at 1200 °C is the same for all components.The optical properties of the three macerals in the coking coal become similar after heating at 1000 °C. Coke obtained at 1200 °C did not contain distinguishable vitrinite grains. At 1200 °C semifusinite and vitrinite coke matrix have highest R_r values among the examined components. Maximum reflectance (R_max) reach similar values for vitrinite and sporinite, slightly lower for semifusinite. Matrix of liptinite coke and matrix of vitrinite coke have considerably stronger anisotropy (R_bi and R_am values) than other components. RIS at 1200 °C is also similar for all components.     

Continuation of the San Andreas fault system into the upper mantle: Evidence from spinel peridotite xenoliths in the Coyote Lake basalt, central California

The Coyote Lake basalt, located near the intersection of the Hayward and Calaveras faults in central California, contains spinel peridotite xenoliths from the mantle beneath the San Andreas fault system. Six upper mantle xenoliths were studied in detail by a combination of petrologic techniques. Temperature estimates, obtained from three two-pyroxene geothermometers and the Al-in-orthopyroxene geothermometer, indicate that the xenoliths equilibrated at 970–1100 °C. A thermal model was used to estimate the corresponding depth of equilibration for these xenoliths, resulting in depths between 38 and 43 km. The lattice preferred orientation of olivine measured in five of the xenolith samples show strong point distributions of olivine crystallographic axes suggesting that fabrics formed under high-temperature conditions. Calculated seismic anisotropy values indicate an average shear wave anisotropy of 6%, higher than the anisotropy calculated from xenoliths from other tectonic environments. Using this value, the anisotropic layer responsible for fault-parallel shear wave splitting in central California is less than 100 km thick. The strong fabric preserved in the xenoliths suggests that a mantle shear zone exists below the Calaveras fault to a depth of at least 40 km, and combining xenolith petrofabrics with shear wave splitting studies helps distinguish between different models for deformation at depth beneath the San Andrea fault system.     

粘弹性裂隙型单斜介质多方位地面三分量记录模拟

针对粘弹性裂隙型单斜介质,本文应用高阶交错网格差分法,对其传播快照以及多方位地面三分量记录进行模拟,结果证明了裂隙填充物的性质引起的各向异性与粘弹性性质对单斜介质波场有明显的影响,总结了不同观测方位地震记录呈现出的规律性,这为进一步了解波场在介质中的传播情况,提出更合理的理论模型,提供一种尝试。     

Shear wave anisotropy in D＂ region beneath the western Pacific

Using seismic shear phases from 47 Tonga-Fiji and its adjacent region events recorded by the CENC and IRIS, and from 26 northeast Asia and north Pacific events recorded by IRIS, we studied the shear wave anisotropy in D" region beneath the western Pacific utilizing the ScS-S differential travel time method and obtained the splitting time values between the radial and transverse components of each ScS wave corresponding to each core-mantle boundary (CMB) reflection point. We found that most shear waves involved horizontally polarized shear wave components traveling faster than vertically polarized shear wave components through the D" region. The splitting time values of ScS wave range from ?0.91 s to 3.21 s with an average value of 1.1 s. The strength of anisotropy varies from ?0.45% to 1.56% with an average value of 0.52%. The observations and analyses show that in the D" region beneath the western Pacific the lateral flow is expected to be dominant and the vertical transverse isotropy may be the main anisotropic structure. This structure feature may be explained by the shape preferred orientation of the CMB chemical reaction products or partial melt and the lattice preferred orientation of the lower mantle materials caused by the lateral flow at lowermost mantle.     

ulti-azimuth three-component surface seismic modeling for viscoelastic cracked monoclinic media

针对粘弹性裂隙型单斜介质,本文应用高阶交错网格差分法,对其传播快照以及多方位地面三分量记录进行模拟,结果证明了裂隙填充物的性质引起的各向异性与粘弹性性质对单斜介质波场有明显的影响,总结了不同观测方位地震记录呈现出的规律性,这为进一步了解波场在介质中的传播情况,提出更合理的理论模型,提供一种尝试.