Reconnaissance paleomagnetism of Late Triassic blocks, Kuyul region, Northern Kamchatka Peninsula, Russia

The northernmost Kamchatka Peninsula is located along the northwestern margin of the Bering Sea and consists of complexly deformed accreted terranes. Progressing inland from the northwestern Bering Sea, the Olyutorskiy, Ukelayat and Koryak superterranes (OLY, UKL and KOR) are crossed. These terranes were accreted to the backstop Okhotsk-Chukotsk volcanic-plutonic belt (OChVB) in northernmost Kamchatka. A sedimentary sequence of Albian to Maastrichtian age overlaps the terranes and units of the Koryak superterrane, and constrains their accretion time. A paleomagnetic study of blocks within the Kuyul (KUY) terrane of the Koryak superterrane was completed at two localities (Camp 2: λ=61.83°N, φ=165.83°E and Camp 3: λ=61.67°N, φ=164.75°E). At both localities, paleomagnetic samples were collected from Late Triassic (225–208 Ma) limestone blocks (2–10 m in outcrop height) within a melange zone. Although weak in remanent magnetization, two components of remanent magnetization were observed during stepwise thermal demagnetization at 32 sites. The A component of magnetization was observed between room temperature and approximately 250 °C. This magnetic component is always of downward directed inclination and shows the best grouping at relatively low degrees of unfolding. Using McFadden–Reid inclination-only statistics and averaging all site means, the resulting A component mean is I_opt=60.3°, I₉₅=5.0° and n=36 (sites). The B magnetic component is observed up to 565 °C, at which temperature, most samples have no measurable remanent magnetization, or growth of magnetic minerals has disrupted the thermal demagnetization process. Combining sites with Fisher estimates of kappa (k-value)≥13 and n (sites)≥3, where bedding orientation differs within a block, most of these sites show the best grouping of B component directions at 100% unfolding, and two of the blocks display remanent magnetizations of both upward and downward directed magnetic inclination. Combining sites with Fisher estimates of kappa (k-value)≥13 and n (sites)≥3, the resulting overall B component paleolatitude and associated uncertainty are λ_obs=30.4°N or S, λ₉₅=8.9° and n=19 (sites). When compared with the expected North America paleolatitude of λ_{APWP expected}=57.9°N, our data support a model in which blocks within the Koryak superterrane are allochthonous and far travelled.     

Evaluation of models for Fourier amplitude spectra for the Taiwan region

In the Taiwan region, the empirical spectral models for estimating ground-motion parameters were obtained recently on the basis of recordings of small to moderate (5.0≤M_L≤6.5) earthquakes. A large collection of acceleration records from the M_L=7.3 Chi-Chi earthquake (21 September, 1999) makes it possible to test the applicability of the established relationships in the case of larger events. The comparison of ground-motion parameters (Fourier amplitude spectra, peak accelerations and response spectra), which were calculated using the models, and the observed data demonstrates that the models could provide an accurate prediction for the case of the Chi-Chi earthquake and the largest aftershocks. However, there are some peculiarities in the ground-motion frequency content and attenuation that, most probably, are caused by the features of the rupture process of the large shallow earthquake source.     

Chironomid responses to long-term metal contamination: a paleolimnological study in two bays of Lake Imandra, Kola Peninsula, northern Russia

Short sediment cores from two gulfs, Monche Bay and Kunchast Bay, of Lake Imandra (Kola Peninsula, northern Russia) were analysed for sediment chemistry and chironomid head capsule remains. Monche Bay has been receiving metals from the Severonikel copper-nickel smelter since the late 1930's. Kunchast Bay was selected in the remotest lake basin as an internal reference site. There were no pronounced changes in the chironomid assemblages with the beginning of slight metal contamination of Kunchast Bay. Based on the reconstructed environmental variables and chironomid assemblages, three developmental stages were distinguished from the chironomid fauna history of Monche Bay: (1) A natural development stage; (2) the early warning stage; and (3) the developing crises stage. During the first period, the changes in the chironomid fauna reflect an anthropogenically undisturbed assemblage, with Micropsectra insignilobus dominating (17–23%). The changes during the second period reflect the initial phase of anthropogenic succession associated with the beginning of metal pollution. The main species showed opposite distributional patterns in this period: the abundance of M. insignilobus decreased, whereas the abundance of Chironomus, Procladius and Sergentia coracina increased. At the same time, maximal numbers were attained for species richness (45) and Shannon-Weaver diversity (4.85) of chironomid assemblages, and the highest head capsule concentration (75 head capsules · g^–1 of dry sediment). The third period was characterized by a major shift in the faunal assemblages, from M. insignilobus to other dominant species, including Chironomus (22–44%), Procladius (10–30 %) and S. coracina (15–18%). Besides fauna changes, assemblages of the third period are distinguished by the occurrence of mouthpart deformities in Chironomus head capsules.     

The use of uniaxial recordings in moment tensor inversions for induced seismic sources

This study shows that the use of uniaxial recordings maintains the linear dependence of the low frequency displacements, with first polarities included, on the moment tensor components, and thus can be easily incorporated in the evaluation of failure mechanism using a moment tensor approach. Synthetic data analysis demonstrates that the errors in the orientation of the principal strain axes of the mechanism solution can be reliably evaluated using a linear error propagation approximation for up to 30% variance in the amplitude data. Case studies at the Kidd Mine indicate that the general characteristics of the mechanism solution are well retrieved, regardless of the type of data used, as relatively small disorientation angles of 15–20° are found between the solutions derived from rotated triaxial data and those from unrotated triaxial data used as uniaxial data, and of 25–30° between the solutions based on rotated triaxial data and those on independent uniaxial data.     

阿穆尔板块相对欧亚板块的欧拉运动

利用中国东北地区1998～2005年的GPS数据,结合俄罗斯、韩国和蒙古的部分GPS观测资料,用加权最小二乘方法以迭代的方式解算了阿穆尔板块相对欧亚板块的欧拉运动矢量。最终的33个测站的位移速率加权均方根残差(WRMS)为0.98mm/a。计算结果表明阿穆尔相对欧亚板块运动的欧拉极位于阿穆尔板块的北部边界外兴安岭东端(54.06±1.13°N,135.87±7.16°E),欧拉旋转量为0.099±0.00735°/Ma。     

Earthquake hazard and risk assessment based on Unified Scaling Law for Earthquakes: Greater Caucasus and Crimea

We continue applying the general concept of seismic risk analysis in a number of seismic regions worldwide by constructing regional seismic hazard maps based on morphostructural analysis, pattern recognition, and the Unified Scaling Law for Earthquakes (USLE), which generalizes the Gutenberg-Richter relationship making use of naturally fractal distribution of earthquake sources of different size in a seismic region. The USLE stands for an empirical relationship log₁₀N(M, L)?=?A?+?B·(5 – M)?+?C·log₁₀L, where N(M, L) is the expected annual number of earthquakes of a certain magnitude M within a seismically prone area of linear dimension L. We use parameters A, B, and C of USLE to estimate, first, the expected maximum magnitude in a time interval at seismically prone nodes of the morphostructural scheme of the region under study, then map the corresponding expected ground shaking parameters (e.g., peak ground acceleration, PGA, or macro-seismic intensity). After a rigorous verification against the available seismic evidences in the past (usually, the observed instrumental PGA or the historically reported macro-seismic intensity), such a seismic hazard map is used to generate maps of specific earthquake risks for population, cities, and infrastructures (e.g., those based on census of population, buildings inventory). The methodology of seismic hazard and risk assessment is illustrated by application to the territory of Greater Caucasus and Crimea.     

Generalized Tikhonov's algorithm for accurate calculation of one-dimensional transient responses directly in time domain

In the last two decades, forward modelling for the time domain (transient) electromagnetic method has concentrated almost entirely on multi-dimensional models and algorithms. At the same time, the interpretation of real field data is still mainly one dimensional. This is caused by the lack of an efficient multi-dimensional acquisition procedure supported by sufficiently fast and reliable inversion software, on the one hand, and by the great efficiency of one-dimensional field set up and interpretation of the data on the other hand. The latter is particularly true for the short offset transient electromagnetic method, which is much less sensitive to multi-dimensional effects, compared to long offset methods. The most commonly used one-dimensional forward modelling algorithms are based on the spectral method, which requires calculating rapidly oscillating Fourier–Bessel (Hankel) integrals. Due to the very fast decay of short offset responses, the integrals become computationally unstable at late times of the transient process. Although this problem has been successfully solved for practically feasible measurement times of conventional short offset systems using transverse electric and mixed transverse electric and transverse magnetic fields, it turned out crucial for novel methods based on the use of unimodal transverse magnetic fields. These methods are much more sensitive to geoelectric parameters of the Earth in general and those of resistive targets, in particular, but they generate responses, which drop at late times significantly faster than those of conventional methods. Such behaviour of transverse magnetic fields represents severe computational problem for the spectral method, but is successfully solved by direct time domain algorithms. This article describes a generalization of the well-known Tikhonov's solution to a boundary value problem directly in time domain, which is applied to an arbitrary one-dimensional earth model excited by an arbitrary source. Contrary to existing spectral algorithms, the described method allows accurate calculations of both transverse electric and transverse magnetic transient responses at arbitrarily late times. On the other hand, it is more time efficient than finite-difference/finite element direct time domain algorithms and provides analytical late-stage asymptotic solutions.     

Image-domain wavefield tomography for tilted transversely isotropic media

Transversely isotropic models with a tilted symmetry axis have become standard for imaging beneath dipping shale formations and in active tectonic areas. Here, we develop a methodology of wave-equation-based image-domain tomography for acoustic tilted transversely isotropic media. We obtain the gradients of the objective function using an integral wave-equation operator based on a separable dispersion relation that takes the symmetry-axis tilt into account. In contrast to the more conventional differential solutions, the integral operator produces only the P-wavefield without shear-wave artefacts, which facilitates both imaging and velocity analysis. The model is parameterized by the P-wave zero-dip normal-moveout velocity, the Thomsen parameter δ, anellipticity coefficient η and the symmetry-axis tilt θ. Assuming that the symmetry axis is orthogonal to reflectors, we study the influence of parameter errors on energy focusing in extended (space-lag) common-image gathers. Distortions in the anellipticity coefficient η introduce weak linear defocusing regardless of reflector dip, whereas δ influences both the energy focusing and depth scale of the migrated section. These results, which are consistent with the properties of the P-wave time-domain reflection moveout in tilted transversely isotropic media, provide important insights for implementation of velocity model-building in the image-domain. Then the algorithm is tested on a modified anticline section of the BP 2007 benchmark model.     

Seismicityof the Novaya Zemlya archipelago: relocated event catalog from 1974 to 2014

We have relocated seismic events registered in the area of the Novaya Zemlya Archipelago, which are not identified as nuclear explosions but are probably of tectonic nature. For the relocation, we collected all available seismic bulletins and waveform data. The location was performed by applying a modified method of generalized beamforming. Verification of the modified method and selection of the travel time model were performed using data on two Novaya Zemlya nuclear explosions that occurred on 02.11.1974 and 24.10.1990. The modified method and the BARENTS travel time model provide sufficient accuracy for event location in the region. The relocation procedure was applied for 9 seismic events registered in the area of the Novaya Zemlya Archipelago. As a result, the new coordinates of the five events turned out to be significantly different from those that were defined previously.     

Role of the Rossiter-McLaughlin effect in the study of close binaries

This paper is devoted to binary stars belonging to the class of eclipsing-variable systems.Photometric and spectroscopic analysis of eclipses allows us to determine geometric parameters of the orbit and physical characteristics of stellar components as well as inclinations of stellar equators to the orbital plane. Estimations of inclinations can be obtained from measurement of the Rossiter-McLaughlin effect, which is discussed using examples of some eccentric binaries with an anomalous apsidal effect. Our task is to find the complete spectrum of solutions of the equation of apsidal motion, depending on the inclinations of the polar axes of the components to the orbital one for these systems, based on their individual spectroscopic and photometric observational data. The matrix of solutions allows us to select those pairs of polar inclinations that provide agreement with the observational apsidal period.