Study on the earthquake fault of the Lulong <Emphasis Type="Italic">M</Emphasis><Subscript>L</Subscript>=6.2 earthquake by precise relocation of aftershock

We have selected 171 near-field records from 391 aftershock records of the Lulong, Hebei Province, earthquake in October 1982 and relocated the hypocenter of 45 aftershocks using the program Hypoinverse. The distribution of aftershocks reveals a set of earthquake faults: a WNW stretching fault truncates two NNE stretching faults. The two branches of faults show the conjugate structure which is often seen in brittle fracture. The NNE stretching faults are connected together. The Luanhe river valley near Lulong developed to a rudiment rift basin surrounded by a series of faults. The fault of Lulong earthquake is a strike-slip fault with tension component. This fault type matches with the activity of Zhangjiakou-Bohai seismic belt (Zhang-Bo belt) and also shows the action of Zhang-Bo belt as a boundary of two secondary active blocks that truncates the NNE fault.     

Estimation of velocity in the uppermost crust in a part of the western Gulf of Corinth,Greece, from the inversion of <Emphasis Type="Italic">P</Emphasis> and <Emphasis Type="Italic">S</Emphasis> arrival times using the neighbourhood algorithm

We determine the velocities in an upper crustal model, composed of three homogeneous layers, for one subregion of the western part of the Gulf of Corinth, NE of the town of Aigion, Greece. We have used local events that occurred there in the year 2001 and were recorded by the Corinth Rift Laboratory Network. Weighted P and S arrival time residuals are minimized using the Neighbourhood Algorithm of Sambridge (1999), combined with the grid search for source locations. The resolution of the inversion is tested by delete-one jackknifing. The model obtained is compared with some other models derived or applied to the subregion. A fast velocity increase between depths of 5 and 7 km is confirmed as the major structural element.     

Mechanism of formation of <Emphasis Type="Italic">Q</Emphasis>-disturbances in the <Emphasis Type="Italic">F</Emphasis><Subscript>2</Subscript> region of the equatorial ionosphere

Using model simulations, the morphological picture (revealed earlier) of the disturbances in the F ₂ region of the equatorial ionosphere under quiet geomagnetic conditions (Q-disturbances) is interpreted. It is shown that the observed variations in the velocity of the vertical E × B plasma drift, related to the zonal E _y component of the electric field, are responsible for the formation of Q-disturbances. The plasma recombination at altitudes of the lower part of the F ₂ region and the dependence of the rate of this process on heliogeophysical conditions compose the mechanism of Q-disturbance formation at night. The daytime positive Q-disturbances are caused exclusively by a decrease in the upward E × B drift, and this type of disturbances could be related to the known phenomenon of counter electrojet. Possible causes of formation of the daytime negative Q-disturbances are discussed.     

Improved method to determine the molar volume and compositions of the NaCl-H<Subscript>2</Subscript>O-CO<Subscript>2</Subscript> system inclusion

On the basis of Parry’s method (1986), an improved method was established to determine the molar volume (Vm) and compositions (X) of the NaCl-H2O-CO2 (NHC) system inclusion. To use this method, the determination of Vm-X only requires three microthermometric data of a NHC inclusion: partial homog-enization temperature (Th ,CO2), salinity (S) and total homogenization temperature (Th). Theoretically, four associated equations are needed containing four unknown parameters: X CO2, XNaCl, Vm and F (volume fraction of CO2 phase in total inclusion when occurring partial homogenization). When they are released, the Vm-X are determined. The former three equations, only correlated with Th ,CO2, S and F, have simplified expressions:XCO2=f1(Th,CO2,S,F),XNaCl=f2(Th,CO2,S,F),Vm=f3(Th,CO2,S,F). The last one is the thermodynamic relationship of X CO2, XNaCl, Vm and Th:f4(XCO2,XNaCl,Vm,Th)=0.Since the above four associated equations are complicated, it is necessary to adopt iterative technique to release them. The technique can be described by:(i) Freely input a F value (0≤F≤1),with Th ,CO2 and S, into the former three equations. As a result,X CO 2,XNaCl and the molar volume value recorded as Vm1 are derived. (ii) Input the X CO2 and XNaCl gotten in the step above into the last equation, and another molar volume value recorded as Vm2 is determined. (iii) If Vm1 is unequal to Vm2, the calculation will be restarted from “(i)”. The iteration is completed until Vm1 is equal to Vm2, which means that the four associated equations are released. Compared to Parry’s (1986) solution method, the improved method is more convenient to use, as well as more accurate to determine X CO 2. It is available for a NHC inlusion whose partial homogenization temperature is higher than clatherate melting temperature and there are no solid salt crystals in the inclusion at parital homogenization.     

Distribution and significance of C~_(40) alkanes in the extracts of Cretaceous source rocks from the Songliao Basin

Source rock extracts and crude oils from the Songliao Basin were analyzed by high-temperature gas chromatography (HTGC), gas chromatography-mass spectrometry (HTGC-MS) and gas chromatography-isotope ratio-mass spectrometry (GC-IRMS), for high molecular-weight alkanes. The distributions of n-alkanes in the Nenjiang Formation extracts are in the C14―C63 range; a bimodal distribution occurs in the C-21 and C21―40 regions. The C30―C37 n-alkanes are accompanied by C29―C35 hopanes, whereas the high molecular-weight C45―C47 n-alkanes co-occur with abundant isoalkanes, alkylcyclohexanes and alkylcyclopentanes. The high δ 13C values of the n-alkanes and the microscopic maceral compositions indicate a highly diversified organic source input for the Nenjiang Formation source rocks, ranging from aquatic plants, blue alge-bacteria, to land plant material. In contrast, n-alkanes in the rock extracts of the Qingshankou Formation are characterized by a single modal distribution, with relatively low abundances of C29―C35 hopanes, but high molecular-weight isoalkanes, alkylcyclohexanes and alkylcyclopentanes. The relatively low δ 13C values of C22―C44 n-alkanes and organic material compositions indicate that the source rocks in the Qingshankou Formation contain dominantly type I algal organic matter. The relative abundance of C 40 compounds in source rocks changes little at low maturity stage, but decreases drastically at higher maturity levels, with a concurrent reduction in the odd/even carbon predominance. In crude oils, in contrast, the relative abundance of C 40 compounds appears to relate closely with the oil source and oil viscosity.     

Numerical simulation of the structure of the high-latitude ionospheric <Emphasis Type="Italic">F</Emphasis> region during meridional HF propagation

A previously developed model of the high-latitude ionosphere is used to calculate the distribution of the ionospheric parameters in the polar region. A specific method for specifying input parameters of the mathematical model, using the experimental data obtained by the method of satellite radio tomography, is used in this case. The spatial distributions of the ionospheric parameters characterized by a complex inhomogeneous structure in the high-latitude region, calculated with the help of the mathematical model, are used to simulate the HF propagation along the meridionally oriented radio paths extending from middle to high latitudes. The method for improving the HF communication between a midlatitude transmitter and a polar-cap receiver is proposed.     

Mechanism of the M<Subscript><Emphasis Type="Italic">L</Emphasis></Subscript>4.0 25 April 2016 earthquake in southwest of France in the vicinity of the Lacq gas field

The source mechanism of the M_L 4.0 25 April 2016 Lacq earthquake (Aquitaine Basin, South-West France) is analyzed from the available public data and discussed with respect to the geometry of the nearby Lacq gas field. It is one of the biggest earthquakes in the area in the past few decades of gas extraction and the biggest after the end of gas exploitation in 2013. The routinely obtained location shows its hypocenter position inside the gas reservoir. We first analyze its focal mechanism through regional broad-band seismograms recorded in a radius of about 50 km epicentral distances and obtain EW running normal faulting above the reservoir. While the solution is stable using regional data only, we observe a large discrepancy between the recorded data on nearby station URDF and the forward modeling up to 1 Hz. We then look for the best epicenter position through performing wave propagation simulations and constraining the potential source area by the peak ground velocity (PGV). The resulting epicentral position is a few to several km away to the north or south direction with respect to station URDF such that the simulated particle motions are consistent with the observation. The initial motion of the seismograms shows that the epicenter position in the north from URDF is preferable, indicating the north-east of the Lacq reservoir. This study is an application of full waveform simulations and characterization of near-field ground motion in terms of an engineering factor such as PGV. The finally obtained solution gives a moment magnitude of M_w 3.9 and the best focal depth of 4 km, which corresponds to the crust above the reservoir rather than its interior. This position is consistent with the tendency of Coulomb stress change due to a compaction at 5 km depth in the crust. Therefore, this earthquake can be interpreted as a relaxation of the shallow crust due to a deeper gas reservoir compaction so that the occurrence of similar events cannot be excluded in the near future. It would be necessary to continue monitoring such local induced seismicity in order to better understand the reservoir/overburden behavior and better assess the local seismic hazard even after the end of gas exploitation.     

<Superscript>40</Superscript>Ar/<Superscript>39</Superscript>Ar geochronology of the Faku tectonites: Implications for the tectonothermal evolution of the Faku block,northern Liaoning

For lack of reliable isotopic chronological data, the metamorphic rock series in the Faku region of northern Liaoning has long been regarded as the platform basement. Recent studies reveal that these deformed and metamorphosed rocks, with a variety of protoliths of plutonic intrusions and supracrustal volcanic and sedimentary rocks, were genetically related to later ductile shearing events, and they, together with the syntectonic intrusions, constituted the large-scale Faku tectonites. In this paper, we report new ⁴⁰Ar/³⁹Ar data on hornblende, biotite, and K-feldspar from typical granitic mylonites in this suite of tectonites. The plateau age 256 Ma of FK53 hornblende and the high-temperature plateau age 262 Ma of Fk51-1 biotite should represent the cooling ages when the granites, formed as a result of Paleozoic oceanic crustal subduction beneath the continental crust or collision of multiple micro-continental blocks, were exhumed into shallow crustal levels. The plateau age 231 Ma of FK51-1 boitite and the apparent age 227 Ma of Fk51-2 K-feldspar are interpreted to record the time of ductile deformation occurring under greenschist facies conditions, i.e. the formation age of the Faku tectonites, while the age gradient from 197 Ma to 220 Ma of Fk51-2 K-feldspar probably record the subsequent stable uplift-cooling process. The tectonic exhumation event indicated by the plateau age 180 Ma of Fk51-2 K-feldspar may be associated with the onset of paleo-Pacific subduction beneath the North China plate. In addition, the U-Pb dating of FK54 zircon from later-intruded granite yields the age of crystallization of this super-unit intrusion at 159 Ma, thus establishing an upper limit for the formation age of the Faku tectonites, while the plateau age 125 Ma of Fk54 K-feldspar most likely corresponds to the rapid cooling and tectonic denudation event associated with the final collision between the Siberian plate and the North China plate. These isotopic ages provide important geochronological constraints for re-evaluating the tectonic essence of the Faku Faulted Convex and ascertaining the suturing boundary between the North China Platform and the Xingmeng Fold System.     

Numerical simulation of the <Emphasis Type="Italic">F</Emphasis>2-layer stratification and appearance of the <Emphasis Type="Italic">F</Emphasis>3 and <Emphasis Type="Italic">G</Emphasis> layers in the equatorial ionosphere: The morphology of the phenomena

This work is devoted to a numerical simulation of the equatorial ionosphere, performed using the GSM TIP model completed with a new block for calculating the electric field. It has been indicated that the usage of the wind system calculated according to the MSIS-90 model makes it possible to reproduce the electromagnetic drift velocities at the equator, the effect of the F2-layer stratification, and the appearance of the F3 layer in the equatorial ionosphere. The calculations performed using the modified GSM TIP model made it possible to detect a maximum in the electron density vertical profile at an altitude of ∼1000 km, formed by H⁺ ions, which we called the G layer. If this layer actually exists, it can be observed during sounding the low-latitude ionosphere from satellites during dark time of day.     

Variations in <Emphasis Type="Italic">foF</Emphasis>2 and <Emphasis Type="Italic">hmF</Emphasis>2 at the end of the 1990s and the beginning of the 2000s

This article considers sparse available data on variations in the main parameters of the ionospheric F2 layer foF2(ss + 2) and hmF2(ss + 2) at the end of the 1990s and the beginning of 2000s. It is shown that the monotonous behavior of hmF2(ss + 2) obtained in earlier publications for the period after 1980 is violated. The hmF2(ss + 2) behavior obtains a more complicated nature by time with a tendency towards a decrease in hmF2(ss + 2) at the beginning of a new century. A statistically significant relationship between foF2(ss + 2) and hmF2(ss + 2) is discovered confirming the Rishbeth statement that during the first hours after sunset, the critical frequency foF2 is governed by dynamical processes via changes in the F2-layer height. It is found that at the end of the interval in question, there is a tendency towards deviations from the above-mentioned dependence. The latter could be a manifestation of the fact that changes in the aeronomical parameters caused by the general cooling and contraction of the thermosphere begin influencing the foF2 value. It is found that in the summer months, the foF2(ss + 2) value demonstrated a systematic decline tendency from the “boundary date” towards the beginning of the 2000s.     

Seasonal variation of the <Emphasis Type="Italic">M</Emphasis><Subscript>2</Subscript> tide

The seasonal cycle of the main lunar tidal constituent M ₂ is studied globally by an analysis of a high-resolution ocean circulation and tide model (STORMTIDE) simulation, of 19 years of satellite altimeter data, and of multiyear tide-gauge records. The barotropic seasonal tidal variability is dominant in coastal and polar regions with relative changes of the tidal amplitude of 5–10 %. A comparison with the observations shows that the ocean circulation and tide model captures the seasonal pattern of the M ₂ tide reasonably well. There are two main processes leading to the seasonal variability in the barotropic tide: First, seasonal changes in stratification on the continental shelf affect the vertical profile of eddy viscosity and, in turn, the vertical current profile. Second, the frictional effect between sea-ice and the surface ocean layer leads to seasonally varying tidal transport. We estimate from the model simulation that the M ₂ tidal energy dissipation at the sea surface varies seasonally in the Arctic (ocean regions north of 60°N) between 2 and 34 GW, whereas in the Southern Ocean, it varies between 0.5 and 2 GW. The M ₂ internal tide is mainly affected by stratification, and the induced modified phase speed of the internal waves leads to amplitude differences in the surface tide signal of 0.005–0.0150 m. The seasonal signals of the M ₂ surface tide are large compared to the accuracy demands of satellite altimetry and gravity observations and emphasize the importance to consider seasonal tidal variability in the correction processes of satellite data.     

Mechanisms of stratification of the <Emphasis Type="Italic">F</Emphasis>2 layer and formation of the <Emphasis Type="Italic">F</Emphasis>3 and <Emphasis Type="Italic">G</Emphasis> layers in the equatorial ionosphere

The paper is dedicated to the studies of formation mechanisms of additional layers in the equatorial ionosphere carried out using numerical simulations with use of the Global Self-Consistent Model of the Thermosphere, Ionosphere, and Protonosphere (GSM TIP) modified in the part of the solution of the electric field equation in the Earth’s ionosphere. Calculations were preformed for quiet geomagnetic conditions using the MSIS-90 model for the calculation of thermospheric parameters. The obtained spatio-temporal pattern of thermospheric circulation and the variations in the dynamo electric field obtained on its basis make it possible to reproduce the stratification effect of the F2 layer and the appearance of the F3 layer in the equatorial ionosphere due to the action of the nonuniform in height zonal electric field at the geomagnetic equator. On the basis of the earlier presented results of calculations using the modified GSM TIP model, the appearance of a maximum in the vertical profile of the electron density at a height of ∼1000 km formed by H+ ions, which we called the G layer, has been predicted. Numerical simulations showed that this layer is formed by the meridional component of the thermospheric wind and is related to the formation of the nighttime midlatitude maximum at heights of the ionospheric F region.     

Seismotectonic studies of the pleistoseist area of the <Emphasis Type="Italic">M</Emphasis> <Subscript> <Emphasis Type="Italic">s</Emphasis> </Subscript> = 6.9 Ilin-Tass earthquake in northeast Yakutia

The complex seismotectonic studies of the pleistoseist area of the Ilin-Tas earthquake (M_s = 6.9), one of the strongest seismic events ever recorded by the regional seismic network in northeastern Russia, are carried out. The structural tectonic position, morphotectonic features of present-day topography, active faults, and types of Cenozoic deformations of the epicentral zone are analyzed. The data of the instrumental observations are summarized, and the manifestations of the strong seismic events in the Yana–Indigirka segment of the Cherskii seismotectonic zone are considered. The explanation is suggested for the dynamical tectonic setting responsible for the Andrei-Tas seismic maximum. This setting is created by the influence of the Kolyma–Omolon indenter, which intrudes into the Cherskii seismotectonic zone from the region of the North American lithospheric plate and forms the main seismogenic structures of the Yana–Indigirka segment in the frontal zone (the Ilin-Tas anticlinorium). The highest seismic potential is noted in the Andrei- Tas block—the focus of the main tectonic impacts from the Kolyma–Omolon superterrane. The general trend of this block coincides with the orientation of the major axis of isoseismal ellipses (azimuth 50°–85°), which were determined from the observations of macroseismic effects on the ground after the Uyandina (M_s = 5.6), Andrei-Tas (M_s = 6.1), and Ilin-Tas (M_s = 6.9) earthquakes.