Triplications for the converted wave in transversely isotropic media with a tilted symmetry axis

In seismic data processing, serious problems could be caused by the existence of triplication and need to be treated properly for tomography and other inversion methods. The triplication in transversely isotropic medium with a vertical symmetry axis has been well studied and concluded that the triplicated traveltime only occurs for S wave and there is no triplication for P and converted PS waves since the P wave convexity slowness always compensates the S wave slowness concavity. Compared with the vertical symmetry axis model, the research of the triplication in transversely isotropic medium with a tilted symmetry axis is still keeping blank. In order to analyse the triplication for the converted wave in the tilted symmetry axis model, we examine the traveltime of the triplication from the curvature of averaged P and S wave slowness. Three models are defined and tested in the numerical examples to illustrate the behaviour of the tilted symmetry axis model for the triplicated traveltime with the change of the rotation angle. Since the orientation of an interface is related to the orientation of the symmetry axis, the triplicated traveltime is encountered for the converted wave in the tilted symmetry axis model assuming interfaces to be planar and horizontal. The triplicated region is influenced by the place and level of the concave curvature of the P and S wave slowness.     

Combined Gravimetric–Seismic Crustal Model for Antarctica

The latest seismic data and improved information about the subglacial bedrock relief are used in this study to estimate the sediment and crustal thickness under the Antarctic continent. Since large parts of Antarctica are not yet covered by seismic surveys, the gravity and crustal structure models are used to interpolate the Moho information where seismic data are missing. The gravity information is also extended offshore to detect the Moho under continental margins and neighboring oceanic crust. The processing strategy involves the solution to the Vening Meinesz-Moritz’s inverse problem of isostasy constrained on seismic data. A comparison of our new results with existing studies indicates a substantial improvement in the sediment and crustal models. The seismic data analysis shows significant sediment accumulations in Antarctica, with broad sedimentary basins. According to our result, the maximum sediment thickness in Antarctica is about 15 km under Filchner-Ronne Ice Shelf. The Moho relief closely resembles major geological and tectonic features. A rather thick continental crust of East Antarctic Craton is separated from a complex geological/tectonic structure of West Antarctica by the Transantarctic Mountains. The average Moho depth of 34.1 km under the Antarctic continent slightly differs from previous estimates. A maximum Moho deepening of 58.2 km under the Gamburtsev Subglacial Mountains in East Antarctica confirmed the presence of deep and compact orogenic roots. Another large Moho depth in East Antarctica is detected under Dronning Maud Land with two orogenic roots under Wohlthat Massif (48–50 km) and the Kottas Mountains (48–50 km) that are separated by a relatively thin crust along Jutulstraumen Rift. The Moho depth under central parts of the Transantarctic Mountains reaches 46 km. The maximum Moho deepening (34–38 km) in West Antarctica is under the Antarctic Peninsula. The Moho depth minima in East Antarctica are found under the Lambert Trench (24–28 km), while in West Antarctica the Moho depth minima are along the West Antarctic Rift System under the Bentley depression (20–22 km) and Ross Sea Ice Shelf (16–24 km). The gravimetric result confirmed a maximum extension of the Antarctic continental margins under the Ross Sea Embayment and the Weddell Sea Embayment with an extremely thin continental crust (10–20 km).     

Comparison of Bayesian estimates of peak ground acceleration (<Emphasis Type="Italic">A</Emphasis>max) with PSHA in Iran

Bayesian probability theory is an appropriate and useful method for estimating parameters in seismic hazard analysis. The analysis in Bayesian approaches is based on a posterior belief, also their special ability is to take into account the uncertainty of parameters in probabilistic relations and a priori knowledge. In this study, we benefited the Bayesian approach in order to estimate maximum values of peak ground acceleration (A_max) also quantiles of the relevant probabilistic distributions are figured out in a desired future interval time in Iran. The main assumptions are Poissonian character of the seismic events flow and properties of the Gutenberg-Richter distribution law. The map of maximum possible values of A_max and also map of 90% quantile of distribution of maximum values of A_max on a future interval time 100 years is presented. According to the results, the maximum value of the A_max is estimated for Bandar Abbas as 0.3g and the minimum one is attributed to Esfahan as 0.03g. Finally, the estimated values in Bayesian approach are compared with what was presented applying probabilistic seismic hazard (PSH) methods based on the conventional Cornel (1968) method. The distribution function of A_max for future time intervals of 100 and 475 years are calculated for confidence limit of probability level of 90%.     

New seismic array solution for earthquake observations and hydropower plant health monitoring

We present the novel fusion of seismic safety monitoring data of the hydropower plant in Chirkey (Caucasus Mountains, Russia). This includes new hardware solutions and observation methods, along with technical limitations for three types of applications: (a) seismic monitoring of the Chirkey reservoir area, (b) structure monitoring of the dam, and (c) monitoring of turbine vibrations. Previous observations and data processing for health monitoring do not include complex data analysis, while the new system is more rational and less expensive. The key new feature of the new system is remote monitoring of turbine vibration. A comparison of the data obtained at the test facilities and by hydropower plant inspection with remote sensors enables early detection of hazardous hydrodynamic phenomena.     

Anisotropy parameter inversion in vertical axis of symmetry media using diffractions

Diffractions play a vital role in seismic processing as they can be utilized for high‐resolution imaging applications and analysis of subsurface medium properties like velocity. They are particularly valuable for anisotropic media as they inherently possess a wide range of dips necessary to resolve the angular dependence of velocity. However, until recently, the focus of diffraction imaging or inversion algorithms have been only on the isotropic approximation of the subsurface. Using diffracted waves, we develop a framework to invert for the effective η model. This effective model is obtained through scanning over possible effective η values and selecting the one that best fits the observed moveout curve for each diffractor location. The obtained effective η model is then converted to an interval η model using a Dix‐type inversion formula. The inversion methodology holds the potential to reconstruct the true η model with sufficiently high accuracy and resolution properties. However, it relies on an accurate estimation of diffractor locations, which in turn requires good knowledge of the background velocity model. We test the effectiveness and applicability of our method on the vertical transverse isotropic Marmousi model. The inversion results yield a reasonable match even for the complex Marmousi model.     

Kinematic parameters of pure‐ and converted‐mode waves in elastic orthorhombic media

I derive the kinematic properties of single‐mode P, S1, and S2 waves as well as converted PS1, PS2, and S1S2 waves in elastic orthorhombic media including vertical velocity, two normal moveout velocities defined in vertical symmetry planes, and three anelliptic parameters (two of them are defined in vertical symmetry plane and one parameter is the cross‐term one). I show that the azimuthal dependence of normal moveout velocity and anellipticity is different in phase and group domains. The effects on‐vertical‐axis singularity and on‐vertical‐axis triplication are considered for pure‐mode S1 and S2 waves and converted‐mode S1S2 waves. The conditions and properties of on‐vertical‐axis triplication are defined in terms of kinematic parameters. The results are illustrated in four homogeneous orthorhombic models and one multilayered orthorhombic model with no variation in azimuthal orientation for all the layers.     

The modified generalized moveout approximation: a new parameter selection

Non‐hyperbolic generalised moveout approximation is a powerful tool to approximate the travel‐time function by using information obtained from two rays. The standard approach for parameter selection is using three parameters defined from zero‐offset ray and two parameters obtained from a reference ray. These parameters include the travel time and travel‐time derivatives of different order. The original parameter selection implies more fit at zero offset compared with offset from a reference ray. We propose an alternative approach for parameter selection within the frame of generalised moveout approximation by transferring more fit from the zero offset to a reference ray by changing in parameter selection. The modified approximation is tested against the original one in few analytical model examples, including the multi‐layered model.     

Preserved travel‐time smoothing in orthorhombic media

Certain degree of smoothness of velocity models is required for most ray‐based migration and tomography. Applying conventional smoothing in model parameters results in offset‐dependent travel‐time errors for reflected events, which can be large even for small contrasts in model parameters between the layers. This causes the shift in both the depth and residual moveout of the migrated images. To overcome this problem in transversely isotropic medium with a vertical symmetry axis, the preserved travel‐time smoothing method was proposed earlier. We extend this method for orthorhombic media with and without azimuthal variation between the layers. We illustrate this method for a single interface between two orthorhombic layers and show that the smoothing‐driven errors in travel time are very small for practical application.     

Low-frequency layer-induced dispersion in a weak-contrast vertically heterogeneous orthorhombic medium

In this paper, we consider wave propagation in a periodically layered medium with orthorhombic symmetry. The weak-contrast approximation is utilized to derive the low-frequency dispersion in effective properties for P, S1 and S2 waves. We show that the dispersion term for all effective properties is controlled by the second-order contrasts in elastic properties from the layers. We also compute the sensitivity matrices for second- and fourth-order coefficients from eigenvalues of frequency-dependent system matrix associated with kinematic parameters for individual wave modes.     

Impacts as sources of the exosphere on Mercury

Chemical processes associated with meteoroid bombardment of Mercury are considered. Meteoroid impacts lead to production of metal atoms as well as metal oxides and hydroxides in the planetary exosphere. By using quenching theory, the abundances of the main Na-, K-, Ca-, Fe-, Al-, Mg-, Si-, and Ti-containing species delivered to the exosphere during meteoroid impacts were estimated. Based on a correlation between the solar photo rates and the molecular constants of atmospheric diatomic molecules, photolysis lifetimes of metal oxides and SiO are estimated. Meteoroid impacts lead to the formation of hot metal atoms (0.2-0.4 eV) produced directly during impacts and of very hot metal atoms (1-2 eV) produced by the subsequent photolysis of oxides and hydroxides in the exosphere of Mercury. The concentrations of impact-produced atoms of the main elements in the exosphere are estimated relative to the observed concentrations of Ca, assumed to be produced mostly by ion sputtering. Condensation of dust grains can significantly reduce the concentrations of impact-produced atoms in the exosphere. Na, K, and Fe atoms are delivered to the exosphere directly by impacts while Ca, Al, Mg, Si, and Ti atoms are produced by the photolysis of their oxides and hydroxides. The chemistry of volatile elements such as H, S, C, and N during meteoroid bombardment is also considered. Our conclusions about the temperature and the concentrations of impact-produced atoms in the exosphere of Mercury may be checked by the Messenger spacecraft in the near future and by BepiColombo spacecraft some years later.