New paleomagnetic results from Ceno-Mesozoic volcanic rocks along southern rim of the Tarim Basin,China*

Paleomagnetic samples were collected from four localities located in the southern rim of the Tarim basin. The samples were taken from volcanic rocks erupted between Jurassic and Quaternary. Detailed analysis of all samples has been carried out with progressive thermal demagnetization. A characteristic remanence (ChRM) with higher unblocking temperature has been isolated from all samples. The pole position from the middle Jurassic is at 52.5°N, 187.9°E(dp = 3.7°,dm =6.5°); the directions of the ChRM of Cretaceous correspond to a paleopole at 69.7°N, 211.6°E (dp = 9.8°,dm = 15.9°); the Quaternary pole from the Pulu site is at 79. 9°N, 183.1°E(dp = 1.6°.dm =2.4°). On the basis of these new paleomagnetic data, tectonic evolution of Tarim block is presented.

     

New Cambrian paleomagnetic pole for Yangtze Block

A total of 120 samples from 12 sites were collected from two flanks of a fold. Stepwise thermal demagnetization has successfully revealed characteristic magnetization components from the rocks in each case. A well-defined component determined from red fine-grained sandstone is clustered in the northeasterly direction with shallow upward inclination (D = 29.3°,I= -19.2°,k = 283.7, α₉₅ = 7.3°. tilt-corrected). The pole position (39.5°N, 247.3°E,dp = 4.0°,dm = 7.6°) derived from this component is close to the Permian pole for the Yangtze Block, indicating that the red fine-grained sandstone has been overprinted. The red mudstone reveals two characteristic components Component A with lower unblocking temperature, characterized by northerly declination and moderate to steep inclination corresponds to a pole position overlay with the present North Pole. Component B (D = 129.1°,I=-23.6°,k = 44.6, α₉₅ = 7.8°, tilt-corrected) with higher unblocking temperature, passes fold test, and yields a pole position (39.5°S, 185.l°E,dp = 4.4°,dm = 8.3°) different from the other poles for the Yangtze Block. It is therefore suggested that component B was probably a primary magnetization and the Yangtze Block was situated at low latitudes in the Southern Hemisphere in the Middle Cambrian.     

Paleomagnetic result of the Lower-Cretaceous from Luanping basin,Hebei Province and its tectonic implications*

A paleomagnetic study of about 95 samples from 16 sites sampled in the Early Cretaceous in Luanping basin in Hebei Province was reported. Stepwise thermal demagnetization was used to isolate magnetic components. Most samples have a characteristic direction with a high temperature component above 500°C. The tectonic-corrected data areD = 347.8°,I = 50.4°, α₉₅ = 7.l°, and the corresponding pole position is at 76.1°N, 346.3°E,with dp =6.4°,dm = 3.8°, paleolatitude λ = 31.1°N. This result indicates a counterclockwise post-Cretaceous rotation of 30.7° ±9.8° with respect to the stable Ordos basin in the west of North China Block, and a non-significant northward motion. This rotation could be related to local fault action or structural detachment, or regional NNW-NWWward motion and collision of Kula-Pacific plate with eastern China since the Early Cretaceous.     

Paleomagnetic results from Late Carboniferous to Early Permian rocks in the northern Qiangtang terrane,Tibet, China,and their tectonic implications

Results of a systematic paleomagnetic study are reported based on Late Carboniferous to Early Permian sedimentary rocks on the north slope of the Tanggula Mountains,in the northern Qiangtang terrane(NQT),Tibet,China.Data revealed that magnetic minerals in limestone samples from the Zarigen Formation(CP^z)are primarily composed of magnetite,while those in sandstone samples from the Nuoribagaribao Formation(Pnr)are dominated by hematite alone,or hematite and magnetite in combination.Progressive thermal,or alternating field,demagnetization allowed us to isolate a stable high temperature component(HTC)in 127 specimens from 16 sites which successfully passed the conglomerate test,consistent with primary remnance.The tilt-corrected mean direction for Late Carboniferous to Early Permian rocks in the northern Qiangtang terrane is D_s=30.2°,I_s=-40.9°,k_s=269.0,a_(95)=2.3°,N=16,which yields a corresponding paleomagnetic pole at 25.7°N,241.5°E(dp/dm=2.8°/1.7°),and a paleolatitude of 23.4°S.Our results,together with previously reported paleomagnetic data,indicate that:(1)the NQT in Tibet,China,was located at a low latitude in the southern hemisphere,and may have belonged to the northern margin of Gondwana during the Late Carboniferous to Early Permian;(2)the Paleo-Tethys Ocean was large during the Late Carboniferous to Early Permian,and(3)the NQT subsequently moved rapidly northwards,perhaps related to the fact that the Paleo-Tethys Ocean was rapidly contracting from the Late Permian to Late Triassic while the Bangong Lake-Nujiang Ocean,the northern branch of the Neo-Tethys Ocean,expanded rapidly during this time.     

Effect of zonal E × B plasma drift on electron density in the low-latitude ionospheric <Emphasis Type="Italic">F</Emphasis> region at a solar activity maximum near vernal equinox

The variations in the density of the ionospheric F2 layer maximum (NmF2) under the action of the zonal plasma drift perpendicularly to the magnetic (B) and electric (E) fields in the direction geomagnetic west-geomagnetic east have been studied using the three-dimensional nonstationary theoretical model of electron and ion densities (N _e and N _i ) and temperatures (T _e and T _i ) in the low-latitude and midlatitude ionospheric F region and plasmasphere. The method of numerical calculations of N _e , N _i , T _e , and T _i , including the advantages of the Lagrangian and Eulerian methods, is used in the model. A dipole approximation of the geomagnetic field (B), taking into account the non-coincidence of the geographic and geomagnetic poles and differences between the positions of the Earth’s and geomagnetic dipole centers, is accepted in the calculations. The calculated NmF2 and altitudes of the F2 layer maximum (hmF2) have been compared with these quantities measured at 16 low-latitude ionospheric sounding stations during the geomagnetically quiet period October 11–12, 1958. This comparison made it possible to correct the input model parameters: the NRLMSISE-00 model [O], the meridional component of the neutral wind velocity according to the HWW90 model, and the meridional component of the equatorial plasma drift due to the electric field specified by the empirical model. It has been indicated that the effect of the zonal E × B plasma drift on NmF2 can be neglected under daytime conditions and changes in NmF2 and hmF2 under the action of this drift are insignificant under nighttime conditions north of 25° and south of ?26° geomagnetic latitude. The effect of the zonal E × B plasma drift on NmF2 and hmF2 is most substantial in the nightside ionosphere approximately from ?20° to 20° geomagnetic latitude, and the neglect of this drift results in an up to 2.4-fold underestimation of NmF2. The found dependence of the effect of the zonal E × B plasma drift on NmF2 and hmF2 on geomagnetic latitude is related to the longitudinal asymmetry of B, asymmetry of the neutral wind about the geomagnetic equator, and changes in the meridional E × B plasma drift at a change in geomagnetic longitude.     

Spatial inhomogeneity of variations in the <Emphasis Type="Italic">F</Emphasis><Subscript>2</Subscript> layer critical frequencies under the conditions of winter 1995

The daily samples of the hourly measurements of the foF2 critical frequency, obtained on January 5–21, 1995, at the midlatitude and high-latitude automated ionospheric stations (geographic latitude higher than 60°), are considered. The {fo} sets are transformed into the {δfoF2} sets of relative variations, for which asymmetry (A) and excess (E) are calculated. The selected stations are grouped into 20 pairs of automated ionospheric stations (AISs) located at distances of 200–10 000 km from one another. Sign estimates of the cross-correlation coefficients between the sets of 16 A and E values for different pairs of stations are used. Two types of structures of the statistical invariant spatial distribution are established: the structures with a scale of about 300 km, invariant with respect to latitude, and with a scale of about 6000 km (for only high latitudes).     

Spherical cap harmonic analysis of regional magnetic anomalies based on CHAMP satellite data

We used CHAMP satellite vector data and the latest IGRF12 model to investigate the regional magnetic anomalies over mainland China. We assumed satellite points on the same surface (307.69 km) and constructed a spherical cap harmonic model of the satellite magnetic anomalies for elements X, Y, Z, and F over Chinese mainland for 2010.0 (SCH2010) based on selected 498 points. We removed the external field by using the CM4 model. The pole of the spherical cap is 36N° and 104°E, and its half-angle is 30°. After checking and comparing the root mean square (RMS) error of ΔX, ΔY, and ΔZ and X, Y, and Z, we established the truncation level at K ^max = 9. The results suggest that the created China Geomagnetic Referenced Field at the satellite level (CGRF2010) is consistent with the CM4 model. We compared the SCH2010 with other models and found that the intensities and distributions are consistent. In view of the variation of F at different altitudes, the SCH2010 model results obey the basics of the geomagnetic field. Moreover, the change rate of X, Y, and Z for SCH2010 and CM4 are consistent. The proposed model can successfully reproduce the geomagnetic data, as other data-fitting models, but the inherent sources of error have to be considered as well.     

Upper mantle structure beneath central Western Europe from data on phase velocities of surface waves

The deep structure of the upper mantle is determined from data on phase velocities of Love and Rayleigh waves measured by a differential method on traces between two stations in central Western Europe. One-dimensional velocity structures are first constructed from data of each pair of stations, after which two-dimensional distributions of SH and SV velocities are calculated by the method of two-dimensional tomography from S wave velocities at fixed depths. The results are presented in the form of 2-D vertical structures of the average S wave velocity (S = (SV + SH)/2) constructed along profiles crossing the region in directions of the best resolution. The main structural features are a higher velocity zone at depths of 60–80 km in the area (48°–50°N, 9°–11°E) and a lower velocity zone in the western part of the region at depths of 100–150 km, probably extending farther beyond the studied area.     

A study of Guptkashi,Uttarakhand earthquake of 6 February 2017 (<Emphasis Type="Italic">M</Emphasis><Subscript><Emphasis Type="Italic">w</Emphasis></Subscript>5.3) in the Himalayan arc and implications for ground motion estimation

The 2017 Guptkashi earthquake occurred in a segment of the Himalayan arc with high potential for a strong earthquake in the near future. In this context, a careful analysis of the earthquake is important as it may shed light on source and ground motion characteristics during future earthquakes. Using the earthquake recording on a single broadband strong-motion seismograph installed at the epicenter, we estimate the earthquake’s location (30.546° N, 79.063° E), depth (H?=?19 km), the seismic moment (M₀?=?1.12×10¹⁷ Nm, M _w 5.3), the focal mechanism (φ?=?280°, δ?=?14°, λ?=?84°), the source radius (a?=?1.3 km), and the static stress drop (Δσ _s ~22 MPa). The event occurred just above the Main Himalayan Thrust. S-wave spectra of the earthquake at hard sites in the arc are well approximated (assuming ω^?2 source model) by attenuation parameters Q(f)?=?500f^0.9, κ?=?0.04 s, and f_max?=?infinite, and a stress drop of Δσ?=?70 MPa. Observed and computed peak ground motions, using stochastic method along with parameters inferred from spectral analysis, agree well with each other. These attenuation parameters are also reasonable for the observed spectra and/or peak ground motion parameters in the arc at distances ≤?200 km during five other earthquakes in the region (4.6?≤?M _w ?≤?6.9). The estimated stress drop of the six events ranges from 20 to 120 MPa. Our analysis suggests that attenuation parameters given above may be used for ground motion estimation at hard sites in the Himalayan arc via the stochastic method.     

Petromagnetic features of sediments at the Mesozoic-Cenozoic boundary: Results from the Gams section

The paper continues a cycle of petromagnetic investigations of epicontinental deposits at the Mesozoic-Cenozoic (K/T) boundary and is devoted to the study of the Gams section (Austria). Using thermomagnetic analysis, the following magnetic phases are identified: goethite (T _C = 90–150°C), hemoilmenite (T _C = 200?300°C), metallic nickel (T _C = 350–360°C), magnetite and titanomagnetite (T _C = 550–610°C), Fe-Ni alloy (T _C = 640–660°C), and metallic iron (T _C = 740–770°C). Their concentrations are determined from M(T). In all samples, ensembles of magnetic grains have similar coercivity spectra and are characterized by a high coercivity. An exception is the lower coercivity of the boundary clay layer due to grains of metallic nickel and iron. With rare exceptions, the studied sediments are anisotropic and generally possess a magnetic foliation, which indicates a terrigenous accumulation of magnetic minerals. Many samples of sandy-clayey rocks have an inverse magnetic fabric associated with the presence of acicular goethite. The values of paramagnetic and diamagnetic components in the deposits are calculated. According to the results obtained, the K/T boundary is marked by a sharp increase in the concentration of Fe hydroxides. The distribution of titanomagnetite reflects its dispersal during eruptive activity, which is better expressed in the Maastrichtian and at the base of the layer J. The along-section distribution of metallic iron, most likely of cosmic origin, is rather uniformly chaotic. The presence of nickel, most probably of impact origin, is a particularly local phenomenon as yet. The K/T boundary is not directly related to an impact event.     

Behavior of Parameters of Enhancements in the Nighttime Electron Concentration in the Ionospheric <Emphasis Type="Italic">F</Emphasis>2 Layer

We have analyzed the behavior of the F2 layer parameters during nighttime periods of enhanced electron concentration by the results of vertical sounding of the ionosphere carried out with five-minute periodicity in Almaty (76°55′ E, 43°15′ N) in 2001–2012. The results are obtained within the frameworks of the unified concept of different types of ionospheric plasma disturbances manifested as variations in the height and half-thickness of the layer accompanied by an increase and decrease of N _m F2 at the moments of maximum compression and expansion of the layer. A good correlation is found between height h _Am , which corresponds to the maximum increase, and layer peak height h _m F, while h _Am is always less than h _m F. The difference between h _Am and h _m F linearly increases with increasing h _m F. Whereas the difference is ~38 km for h _m F = 280 km, it is ~54 km for h _m F = 380 km. Additionally, the correlation is good between the increase in the electron concentration in the layer maximum ΔN _m and the maximum enhancement at the fixed height ΔN; the electron concentration enhancement in the layer maximum is about two to three times lower than its maximum enhancement at the fixed height.     

Ionospheric Time-delay over Akure Using Global Positioning System Observations

Ionospheric time delay (VΔt) variability using Global Positioning System (GPS) data over Akure (7.15°N, 5.12°E), Nigeria, has been studied. The observed variability of VΔt in comparison to older results of vertical total electron content (TEC) across similar regions has shown equivalent signatures. Higher monthly mean values of VΔt (MVΔt) were observed during daytime as compared to nighttime (pre- and post-midnight) hours in all months. The highest MVΔt observed in September during daytime hours range between ~6 and ~21 ns (~1.80 and ~6.30 m) and at post-midnight, they are in the range of ~1 to ~6 ns (~0.3 to ~1.80 m). The possible mechanisms responsible for this variability were discussed. Seasonal VΔt were investigated as well.     

Distribution of Lg coda <Emphasis Type="Italic">Q</Emphasis> in Xinjiang and its adjacent regions

In this study, we collected 1 156 broadband vertical components records at 22 digital seismic stations in Xinjiang region, Ürümqi station, and 7 stations in the adjacent regions during the period of 1999–2003. The records were firstly processed by the stacked spectral ratio method to obtain Q ₀ (Q at 1 Hz) and the frequency correlation factor η corresponding to each path. Based on the results, the distribution images of Q ₀ and η in 1°×1° grids for Xinjiang region were gained by the back-projection technique. The results indicate that Q ₀ is high (300–450) in the Tarim platform and marginal Siberian platform, while Q ₀ is low (150–250) in the southern regions as west Kunlun fold system and Songpan-Ganzi fold system. In the northern regions as Junggar fold system and Tianshan fold system, Q ₀ is also low (250–300) and η varies between 0.5 and 0.9.     

The geomagnetic field at the Proterozoic-Paleozoic boundary and lower-mantle plumes

Data on the amplitude of variations in the direction of the geomagnetic field and the frequency of reversals in the Vendian-Cambrian are presented. It has been established from these data that (a) distributions of variations in the direction of the geomagnetic field S _p are bimodal (modes 9° and 11°); (b) the maximum of the average amplitude S _p takes place by 5–10 Myr later than the Vendian-Cambrian boundary; (c) S _p tends to increase as plume epicenters are approached; and (d) the plume formation is more often confined to intervals with different frequencies of geomagnetic reversals than to the interval of a stable state of the geomagnetic field without reversals (Vendian hyperchron). The listed features of the geomagnetic field behavior are repeated near all boundaries of geological eras of the Phanerozoic.     

Palaeomagnetism of Middle Ordovician Carbonate Sequence,Vaivara Sinimäed Area,Northeast Estonia,Baltica

The hill range of Vaivara Sinimäed in northeast Estonia consists of several narrow east- to northeast-trending glaciotectonic fold structures. The folds include tilted (dips 4–75°) Middle Ordovician (early Darriwilian) layered carbonate strata that were studied by mineralogical, palaeomagnetic, and rock magnetic methods in order to specify the postsedimentational history of the area and to obtain a better control over the palaeogeographic position of Baltica during the Ordovician. Mineralogical studies revealed that (titano)magnetite, hematite, and goethite are carriers of magnetization. Based on data from 5 sites that positively passed a DC tilt test, a south-easterly downward directed component A (D^ref = 154.6°± 15.3°, I^ref = 60.9°± 9.7°) was identified. The component is carried by (titano)magnetite, dates to the Middle Ordovician (Plat = 17.9°, Plon = 47.3°, K = 46.7, A95 = 11.3°), and places Baltica at mid-southerly latitudes. Observations suggest that in sites that do not pass the tilt test, the glaciotectonic event has caused some rotation of blocks around their vertical axis.     

Activation of seismicity in Central and South Asia after the Makran earthquakes: Possible acceleration of preparation of large seismic events in the Tien Shan region

Two zones of seismicity (ten events with M _w = 7.0–7.7) stretching from Makran and the Eastern Himalaya to the Central and EasternTien Shan, respectively, formed over 11 years after the great Makran earthquake of 1945 (M _w = 8.1). Two large earthquakes (M _w = 7.7) hit theMakran area in 2013. In addition, two zones of seismicity (M ≥ 5.0) occurred 1–2 years after theMakran earthquake in September 24, 2013, stretching in the north-northeastern and north-northwestern directions. Two large Nepal earthquakes struck the southern extremity of the “eastern” zone (April 25, 2015, M _w = 7.8 and May 12, 2015, M _w = 7.3), and the Pamir earthquake (December 7, 2015, M _w = 7.2) occurred near Sarez Lake eastw of the “western” zone. The available data indicate an increase in subhorizontal stresses in the region under study, which should accelerate the possible preparation of a series of large earthquakes, primarily in the area of the Central Tien Shan, between 70° and 79° E, where no large earthquakes (M _w ≥ 7.0) have occurred since 1992.     

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.     

Anomalous Variation in GPS TEC,Land and Ocean Parameters Prior to 3 Earthquakes

The present study reports the analysis of GPS TEC prior to 3 earthquakes (M > 6.0). The earthquakes are: (1) Loyalty Island (22°36′S, 170°54′E) on 19 January 2009 (M = 6.6), (2) Samoa Island (15°29′S, 172°5′W) on 30 August 2009 (M = 6.6), and (3) Tohoku (38°19′N, 142°22′E) on 11 March 2011 (M = 9.0). In an effort to search for a precursory signature we analysed the land and ocean parameters prior to the earthquakes, namely SLHF (Land) and SST (Ocean). The GPS TEC data indicate an anomalous behaviour from 1–13 days prior to earthquakes. The main purpose of this study was to explore and demonstrate the possibility of any changes in TEC, SST, and SLHF before, during and after the earthquakes which occurred near or beneath an ocean. This study may lead to better understanding of response of land, ocean, and ionosphere parameters prior to seismic activities.     

Dynamics of the midlatitude ionospheric <Emphasis Type="Italic">F</Emphasis> region at sunrise

The behavior of the F2 layer at sunrise has been studied based on vertical-incidence ionospheric sounding data in Almaty (76°55′E, 43°15′N). Records with small amplitudes of electron density background fluctuations were selected in order to exactly estimate the onsets of a pronounced increase in the electron density at different altitudes. It has been indicated that the electron density growth rate is a function of altitude; in this case, the growth rate at the F2 layer maximum is much lower than such values at fixed altitudes of ～30–55 km below the layer maximum. The solar zenith angle (χ) and the blanketing layer thickness (h ₀) at the beginning of a pronounced increase in the electron density at altitude h are linearly related to the h value, and these quantities vary within ～90° < χ < 100° and 180 km < h ₀ < 260 km, respectively.     

A way of estimating the characteristic slip displacement

During the ruptures of an earthquake,the strain energy.△E,.will be transferred into,at least,three parts,i.e..the seismic radiation energy(E_s),fracture energy(E_g),and frictional energy(E_f),that is,△E = E_s + E_g + E_f.Friction,which is represented by a velocity- and state-dependent friction law by some researchers,controls the three parts.One of the main parameters of the law is the characteristic slip displacement.D_c.It is significant and necessary to evaluate the reliable value of D_c from observed and inverted seismic data.Since D_c controls the radiation efficiency.η_R = E_s/(E_s+ E_g),the value of η_r is a good constraint of estimating D_c.Integrating observed data and inverted results of source parameters from recorded seismograms.the values of E_s and E_g of an earthquake can be measured,thus leading to the value of η_R.The constraint used to estimate the reliable value of D_c will be described in this work.An example of estimates of D_c.based on the observed and inverted values of source parameters of the September 20,1999 M_S 7.6 Chi-Chi(Ji-Ji).Taiwan region,earthquake will be presented.