Apsidal motion and physical parameters of the eclipsing binary system AR Cas

Using the four-channel automatic photoelectric photometer of the Sternberg Astronomical Institute’s Tien Shan Mountain Observatory, we have acquired accurate (σ_obs≈0.004^m) W BV R brightness measurements for the eclipsing binary AR Cas during selected phases before eclipse ingress and after egress, as well as at the center of minima. A joint analysis of these measurements with other published data has enabled us to derive for the first time a self-consistent set of physical and geometrical parameters for the star and the evolutionary age of its components, t=(60±3)×10⁶ years. We have found the period of the apsidal motion (U_obs=1100±160 years, \(\dot \omega _{obs} = 0^\circ .327 \pm 0^\circ .049\) years^?1) and the apsidal parameter of the primary, logk _2,1 ^obs =?2.41±0.08, with the apsidal parameter being in good agreement with current models of stellar evolution. There is an ultraviolet excess in the primary’s radiation, Δ(U?B)=?0.12^m and Δ(B?V)=?0.06^m, possibly due to a metal deficiency in the star’s atmosphere.     

Cretaceous Rocks of the Penzhina Bay: Mineralogy,Petrography, and Geodynamic Sedimentation Conditions

Cretaceous sedimentary and volcanosedimentary rocks from northwestern Kamchatka are considered. The stadial analysis has revealed variable impacts of three major provenances upon the Cretaceous Penzhina sedimentary basin. The provenances were composed of volcanic and volcanosedimentary rocks (Uda–Murgal island arc and Okhotsk–Chukotka volcanic belt) and granitic–metamorphic rocks (the mature Asian continental margin). Sediments were largely accumulated owing to the erosion of island-arc volcanics during reactivation of the Uda–Murgal island arc (Hauterivian–Barremian) or the Okhotsk–Chukotka volcanic belt (middle Albian–Cenomanian). Eroded granitic–metamorphic rocks of the mature Asian continental margin (Berriasian–Valanginian) or Asian metamorphic–volcanic rocks (Santonian–Campanian) were supplied to the basin during tectonically quiet periods (Berriasian–Valanginian and late Cenomanian–Campanian). Compositional changes in provenances were related to active tectonic processes at the continental margin, including evolution and closure of the Uda–Murgal island-arc system and origin of the Okhotsk–Chukotka volcanic belt. The postsedimentary modification of Cretaceous rocks deposited in forearc trough beyond the tectonically active accretionary prism is characterized by a low degree of clastic component alteration.     

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.     

Using the Digital Models of Gravity Anomalies for Zoning of the Earth’s Gravity Field

This paper addresses the questions associated with using the models in regioning the Earth’s gravity field depending on its roughness and type of anomalies. A formalized technique for zoning based on the analysis of the intensity and variability of gravity anomalies is suggested. The suggested technique includes preprocessing the initial gravity anomalies aimed at eliminating the noise component; partitioning the studied region into the elementary segments; calculating the primary characteristics of the complexity and roughness (the intensity and variability of gravity anomalies) for each segment; and, based on these characteristics, classifying each segment into a particular category of complexity. The proposed system of classification of segments relies on the use of three classes of intensity and three classes of variability of gravity anomalies and four categories of complexity of the regions, which are derived from these classes. As a result of applying the technique, the mapped territory is subdivided into the regions that are uniform in terms of their geophysical properties. The developed technique is used for comparing the quality of the different digital models of gravity anomalies, including the Russian RGM model (A.P. Karpinsky Russian Geological Research Institute (VSEGEI)) and the WGM model (Bureau Gravimétrique International (BGI)) in the region of the Mongol–Okhotsk orogenic belt and Sikhote–Alin fold system. The results of the study can be used for zoning of the other geophysical fields and for planning the locations of the new survey networks for increasing the accuracy of the initial data used for compiling the maps.     

Low-frequency anisotropy in fractured and layered media

Computing effective medium properties is very important when upscaling data measured at small scale. In the presence of stratigraphic layering, seismic velocities and anisotropy parameters are scale and frequency dependent. For a porous layer permeated by aligned fractures, wave-induced fluid flow between pores and fractures can also cause significant dispersion in velocities and anisotropy parameters. In this study, we compare the dispersion of anisotropy parameters due to fracturing and layering at low frequencies. We consider a two-layer model consisting of an elastic shale layer and an anelastic sand layer. Using Chapman's theory, we introduce anisotropy parameters dispersion due to fractures (meso-scale) in the sand layer. This intrinsic dispersion is added to anisotropy parameters dispersion induced by layering (macro-scale) at low frequencies. We derive the series coefficients that control the behaviour of anisotropy parameters at low frequencies. We investigate the influences of fracture length and fracture density on fracturing effect, layering effect and combined effect versus frequency and volume fraction of sand layer. Numerical modelling results indicate that the frequency dependence due to layering is not always the dominant effect of the effective properties of the medium. The intrinsic dispersion is not negligible compared with the layering effect while evaluating the frequency-dependent properties of the layered medium.     

Frequency-dependent PP and PS reflection coefficients in fractured media

When a porous layer is permeated by mesoscale fractures, wave-induced fluid flow between pores and fractures can cause significant attenuation and dispersion of velocities and anisotropy parameters in the seismic frequency band. This intrinsic dispersion due to fracturing can create frequency-dependent reflection coefficients in the layered medium. In this study, we derive the frequency-dependent PP and PS reflection coefficients versus incidence angle in the fractured medium. We consider a two-layer vertical transverse isotropy model constituted by an elastic shale layer and an anelastic sand layer. Using Chapman's theory, we introduce the intrinsic dispersion due to fracturing in the sand layer. Based on the series coefficients that control the behaviour of velocity and anisotropy parameters in the fractured medium at low frequencies, we extend the conventional amplitude-versus-offset equations into frequency domain and derive frequency-dependent amplitude-versus-offset equations at the elastic–anelastic surface. Increase in fracture length or fracture density can enlarge the frequency dependence of amplitude-versus-offset attributes of PP and PS waves. Also, the frequency dependence of magnitude and phase angle of PP and PS reflection coefficients increases as fracture length or fracture density increases. Amplitude-versus-offset type of PP and PS reflection varies with fracture parameters and frequency. What is more, fracture length shows little impact on the frequency-dependent critical phase angle, while the frequency dependence of the critical phase angle increases with fracture density.     

The kinematics of S waves in acoustic orthorhombic media

Orthorhombic models are often used in the seismic industry nowadays to describe azimuthal and polar anisotropy and reasonably realistic in capturing the features of the earth interior. It is challenging to handle so many model parameters in the seismic data processing. In order to reduce the number of the parameters for P wave, the acoustic orthorhombic medium is proposed by setting all on-axis S wave velocities to zero. However, due to the coupled behaviour for P and S waves in the orthorhombic model, the ‘S wave artefacts’ are still remained in the acoustic orthorhombic model, which kinematics needs to be defined and analysed. In this paper, we analyse the behaviour of S wave in acoustic orthorhombic media. By analysis of the slowness surface in acoustic orthorhombic media, we define the S waves (or S wave artefacts) that are more complicated in shape comparing to the one propagating in an acoustic transversely isotropic medium with a vertical symmetry axis. The kinematic properties of these waves are defined and analysed in both phase and group domain. The caustics, amplitude and the multi-layered case for S wave in acoustic orthorhombic model are also discussed. It is shown that there are two waves propagating in this acoustic orthorhombic medium. One of these waves is similar to the one propagating in acoustic vertical symmetry axis media, whereas another one has a very complicated shape consisting of two crossing surfaces.     

Water retention effects on elastic properties of Opalinus shale

Shales play an important role in many engineering applications such as nuclear waste, CO₂ storage and oil or gas production. Shales are often utilized as an impermeable seal or an unconventional reservoir. For both situations, shales are often studied using seismic waves. Elastic properties of shales strongly depend on their hydration, which can lead to substantial structural changes. Thus, in order to explore shaly formations with seismic methods, it is necessary to understand the dependency of shale elastic properties on variations in hydration. In this work, we investigate structural changes in Opalinus shale at different hydration states using laboratory measurements and X-ray micro-computed tomography. We show that the shale swells with hydration and shrinks with drying with no visible damage. The pore space of the shale deforms, exhibiting a reduction in the total porosity with drying and an increase in the total porosity with hydration. We study the elastic properties of the shale at different hydration states using ultrasonic velocities measurements. The elastic moduli of the shale show substantial changes with variations in hydration, which cannot be explained with a single driving mechanism. We suggest that changes of the elastic moduli with variations in hydration are driven by multiple competing factors: (1) variations in total porosity, (2) substitution of pore-filling fluid, (3) change in stiffness of contacts between clay particles and (4) chemical hardening/softening of clay particles. We qualitatively and quantitatively analyse and discuss the influence of each of these factors on the elastic moduli. We conclude that depending on the microstructure and composition of a particular shale, some of the factors dominate over the others, resulting in different dependencies of the elastic moduli on hydration.     

A new parameterization for generalized moveout approximation,based on three rays

A simple and accurate traveltime approximation is important in many applications in seismic data processing, inversion and modelling stages. Generalized moveout approximation is an explicit equation that approximates reflection traveltimes in general two-dimensional models. Definition of its five parameters can be done from properties of finite offset rays, for general models, or by explicit calculation from model properties, for specific models. Two versions of classical finite-offset parameterization for this approximation use traveltime and traveltime derivatives of two rays to define five parameters, which makes them asymmetrical. Using a third ray, we propose a balance between the number of rays and the order of traveltime derivatives. Our tests using different models also show the higher accuracy of the proposed method. For acoustic transversely isotropic media with a vertical symmetry axis, we calculate a new moveout approximation in the generalized moveout approximation functional form, which is explicitly defined by three independent parameters of zero-offset two-way time, normal moveout velocity and anellipticity parameter. Our test shows that the maximum error of the proposed transversely isotropic moveout approximation is about 1/6 to 1/8 of that of the moveout approximation that had been reported as the most accurate approximation in these media. The higher accuracy is the result of a novel parameterization that do not add any computational complexity. We show a simple example of its application on synthetic seismic data.