Comparison of SPT and <Emphasis Type="Italic">V</Emphasis><Subscript>s</Subscript>-based liquefaction analyses: a case study in Erciş (Van,Turkey)

Liquefaction which is one of the most destructive ground deformations occurs during an earthquake in saturated or partially saturated silty and sandy soils, which may cause serious damages such as settlement and tilting of structures due to shear strength loss of soils. Standard (SPT) and cone (CPT) penetration tests as well as the shear wave velocity (V _s)-based methods are commonly used for the determination of liquefaction potential. In this research, it was aimed to compare the SPT and V _s-based liquefaction analysis methods by generating different earthquake scenarios. Accordingly, the Erci? residential area, which was mostly affected by the 2011 Van earthquake (M _w = 7.1), was chosen as the model site. Erci? (Van, Turkey) and its surroundings settle on an alluvial plain which consists of silty and sandy layers with shallow groundwater level. Moreover, Çald?ran, Erci?–Kocap?nar and Van Fault Zones are the major seismic sources of the region which have a significant potential of producing large magnitude earthquakes. After liquefaction assessments, the liquefaction potential in the western part of the region and in the coastal regions nearby the Lake Van is found to be higher than the other locations. Thus, it can be stated that the soil tightness and groundwater level dominantly control the liquefaction potential. In addition, the lateral spreading and sand boiling spots observed after the 23rd October 2011 Van earthquake overlap the scenario boundaries predicted in this study. Eventually, the use of V _s-based liquefaction analysis in collaboration with the SPT results is quite advantageous to assess the rate of liquefaction in a specific area.     

Crustal Structure Along Sunda-Banda Arc Transition Zone from Teleseismic Receiver Functions

We analyzed receiver function of teleseismic events recorded at twelve Indonesian-GEOFON (IA-GE) broadband stations using nonlinear Neighbourhood Algorithm (NA) inversion and H-k stacking methods to estimate crustal thickness, V _p /V _s ratios and S-wave velocity structure along Sunda-Banda arc transition zone. We observed crustal thickness of 34–37 km in Timor Island, which is consistent with the previous works. The thick crust (> 30 km) is also found beneath Sumba and Flores Islands, which might be related to the arc-continent collision causing the thickened crust. In Timor and Sumba Islands, we observed high V _p /V _s ratio (> 1.84) with low velocity zone that might be associated with the presence of mafic and ultramafic materials and fluid filled fracture zone. The high V _p /V _s ratio observed at Sumbawa and Flores volcanic Islands might be an indication of partial melt related to the upwelling of hot asthenosphere material through the subducted slab.     

The investigation of soil–structure resonance of historical buildings using seismic refraction and ambient vibrations HVSR measurements: a case study from Trabzon in Turkey

In this study, two different historical structures built in Trabzon have been processed by ambient vibrations and seismic refraction measurements. One of the investigated historical structures is the Atatürk Pavilion built in the nineteenth century, and the other one is Hagia Sophia which was built in the thirteenth century. These two buildings are among the most important historical buildings in Trabzon and are very important for the tourism of the city. In order to determine peak/s frequency and amplitude from the horizontal-to-vertical spectral ratios (HVSRs), we have performed several measurements of ambient vibrations both inside (at different floors) and outside (on the ground) of structures. We have also conducted seismic prospecting to evaluate the vertical 1D and 2D profile of longitudinal and shear seismic waves, V_p and V_s, respectively. To this purpose, we have performed seismic refraction tomography and MASW. Ambient vibrations and seismic measurements were compared with each other. The results show that average predominant frequencies and HVSR amplitudes of inside and outside of Atatürk Pavilion are 4.0 Hz, 7.8 Hz and 2.6, 2.3, respectively. The V_p values vary from 300 to 2070 m/s, and the V_s for maximum effective depth is up to 790 m/s in Atatürk Pavilion. On the other hand, average predominant frequencies and HVSR amplitudes of inside and outside of Hagia Sophia and its tower are 4.7, 4.4 and 2.4 Hz and 1.6, 1.8 and 6.9, respectively. V_p values range from 450 to 2200 m/s, and V_s for maximum effective depth is also up to 1000 m/s in Hagia Sophia. The frequency values (F₀?=?V_s/4 h) calculated from the velocities up to the maximum effective depth for Atatürk Pavilion are in good agreement with the predominant frequency values determined from ambient vibrations. Atatürk Pavilion and Hagia Sophia soils have been classed according to Eurocode 8 by using V_S30 values. The class was defined as “B.” Moreover, the bedrock in studied area is basalt. The high V_p and V_s values are also compatible with the lithology. The HVSR curves measured at the Hagia Sophia show the presence of clear peaks when compared to the Atatürk Pavilion. At the same time, there are marked velocity changes in the V_s sections calculated in both areas. As a result, in both areas there are significant impedance contrasts in the subsoil. However, this impedance contrast is more evident in Hagia Sophia. This could be also compatible with a lithological transition. The possible soil–structure interaction was investigated by using all the results and evaluated in terms of resonance risk. It is thought that the probability of resonance risk at Atatürk Pavilion is low according to the ambient vibrations measurements. However, resonance risk should be taken into consideration at Hagia Sophia site since the predominant frequency values are very close to each other. Finally, this site should be investigated in detail and necessary precautions should be taken against the risk of resonance.     

Site Specific Ground Motion Modeling and Seismic Response Analysis for Microzonation of Baku,Azerbaijan

We investigated ground response for Baku (Azerbaijan) from two earthquakes of magnitude M6.3 occurred in Caspian Sea (characterized as a near event) and M7.5 in Shamakhi (characterized as a remote extreme event). S-wave velocity with the average shear wave velocity over the topmost 30 m of soil is obtained by experimental method from the V _P values measured for the soils. The downtown part of Baku city is characterized by low V_S30 values (< 250 m/s), related to sand, water-saturated sand, gravel-pebble, and limestone with clay. High surface PGA of 240 gal for the M7.5 event and of about 190 gal for the M6.3 event, and hence a high ground motion amplification, is observed in the shoreline area, through downtown, in the north-west, and in the east parts of Baku city with soft clays, loamy sands, gravel, sediments.     

Optimal intensity measures for probabilistic seismic demand modeling of extended pile-shaft-supported bridges in liquefied and laterally spreading ground

Seismic intensity measures (IMs) perform a pivotal role in probabilistic seismic demand modeling. Many studies investigated appropriate IMs for structures without considering soil liquefaction potential. In particular, optimal IMs for probabilistic seismic demand modeling of bridges in liquefied and laterally spreading ground are not comprehensively studied. In this paper, a coupled-bridge-soil-foundation model is adopted to perform an in-depth investigation of optimal IMs among 26 IMs found in the literature. Uncertainties in structural and geotechnical material properties and geometric parameters of bridges are considered in the model to produce comprehensive scenarios. Metrics such as efficiency, practicality, proficiency, sufficiency and hazard computability are assessed for different demand parameters. Moreover, an information theory based approach is adopted to evaluate the relative sufficiency among the studied IMs. Results indicate the superiority of velocity-related IMs compared to acceleration, displacement and time-related ones. In particular, Housner spectrum intensity (HI), spectral acceleration at 2.0 s (S _a-20), peak ground velocity (PGV), cumulative absolute velocity (CAV) and its modified version (CAV ₅) are the optimal IMs. Conversely, Arias intensity (I _a ) and shaking intensity rate (SIR) which are measures often used in liquefaction evaluation or related structural demand assessment demonstrate very low correlations with the demand parameters. Besides, the geometric parameters do not evidently affect the choice of optimal IMs. In addition, the information theory based sufficiency ranking of IMs shows an identical result to that with the correlation measure based on coefficient of determination (R ²). This means that R ² can be used to preliminarily assess the relative sufficiency of IMs.     

Seismic behavior and zoning of the sabkha soils in Jubail industrial city,Saudi Arabia

This study aimed at the micro-level seismic behavior and zoning of the saline sabkha strata in Jubail industrial area in Eastern Saudi Arabia. It encompasses the evaluation of the site-specific seismic response parameters and the liquefaction potential for various possible subsurface conditions under the probable seismic event(s). The approach to achieve the objectives of this study included the following: analysis of geologic, hydrologic, and geotechnical data of the area; performance of field and laboratory dynamic testing; and dynamic modeling and analysis of the subsurface profiles. The results of the simulation have been used to develop liquefaction potential maps and site-specific spectra of the study area, consisting of ten seismic zones under a range of probable peak horizontal ground acceleration (PHA). Results do not show significant probability of liquefaction of the loose soil layers in the study area at the maximum possible design PHA of 0.035 g; however, liquefaction is anticipated at higher PHA values. Site-specific spectral response resulted in values of S _s and S ₁ spectral accelerations to be different as compared to those suggested by local standards. The resulting seismic micro-zonation maps and the corresponding parameters are very useful for the stability analysis of the existing and planned structures in the Jubail area.     

Soil characterization of Tınaztepe region (İzmir/Turkey) using surface wave methods and nakamura (HVSR) technique

To determine the shear wave velocity structure and predominant period features of T?naztepe in ?zmir, Turkey, where new building sites have been planned, active–passive surface wave methods and single-station microtremor measurements are used, as well as surface acquisition techniques, including the multichannel analysis of surface waves (MASW), refraction microtremor (ReMi), and the spatial autocorrelation method (SPAC), to pinpoint shallow and deep shear wave velocity. For engineering bedrock (V _s > 760 m/s) conditions at a depth of 30 m, an average seismic shear wave velocity in the upper 30 m of soil (AV_s30) is not only accepted as an important parameter for defining ground behavior during earthquakes, but a primary parameter in the geotechnical analysis for areas to be classified by V _s30 according to the National Earthquake Hazards Reduction Program (NEHRP). It is also determined that Z1.0, which represents a depth to V _s = 1000 m/s, is used for ground motion prediction and changed from 0 to 54 m. The sediment–engineering bedrock structure for T?naztepe that was obtained shows engineering bedrock no deeper than 30 m. When compared, the depth of engineering bedrock and dominant period map and geology are generally compatible.     

Characterisation of ground motion recording stations in the Groningen gas field

The seismic hazard and risk analysis for the onshore Groningen gas field requires information about local soil properties, in particular shear-wave velocity (V_S). A fieldwork campaign was conducted at 18 surface accelerograph stations of the monitoring network. The subsurface in the region consists of unconsolidated sediments and is heterogeneous in composition and properties. A range of different methods was applied to acquire in situ V_S values to a target depth of at least 30 m. The techniques include seismic cone penetration tests (SCPT) with varying source offsets, multichannel analysis of surface waves (MASW) on Rayleigh waves with different processing approaches, microtremor array, cross-hole tomography and suspension P-S logging. The offset SCPT, cross-hole tomography and common midpoint cross-correlation (CMPcc) processing of MASW data all revealed lateral variations on length scales of several to tens of metres in this geological setting. SCPTs resulted in very detailed V_S profiles with depth, but represent point measurements in a heterogeneous environment. The MASW results represent V_S information on a larger spatial scale and smooth some of the heterogeneity encountered at the sites. The combination of MASW and SCPT proved to be a powerful and cost-effective approach in determining representative V_S profiles at the accelerograph station sites. The measured V_S profiles correspond well with the modelled profiles and they significantly enhance the ground motion model derivation. The similarity between the theoretical transfer function from the V_S profile and the observed amplification from vertical array stations is also excellent.     

Limits of transversely isotropic elastic parameters for the existence of classical Rayleigh waves

Solutions of P-SV equations of motion in a homogeneous transversely isotropic elastic layer contain a factor exp(±ν _j z), where z is the vertical coordinate and j?=?1, 2. For computing Rayleigh wave dispersion in a multi-layered half space, ν _j is computed at each layer. For a given phase velocity (c), ν _j becomes complex depending on the transversely isotropic parameters. When ν _j is complex, classical Rayleigh waves do not exist and generalised Rayleigh waves propagate along a path inclined to the interface. We use transversely isotropic parameters as α _H , β _V , ξ, ? and η and find their limits beyond which ν _j becomes complex. It is seen that ν _j depends on ? and η, but does not depend on ξ. The complex ν _j occurs when ? is small and η is large. For a given c/β _V , the region of complex ν _j in a ? -η plane increases with the increase of α _H /β _V . Further, for a given α _H /β _V , the complex region of ν _j increases significantly with the decrease of c/β _V . This study is useful to compute dispersion parameters of Rayleigh waves in a layered medium.     

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.     

Interaction of geometry and mechanical property of trapezoidal sedimentary basins with incident <Emphasis Type="Italic">SH</Emphasis> waves

This paper investigated the effects of basin geometry and material property on the response of 2D trapezoidal sediment-filled basin to incident plane SH waves. Ten basin configurations with different geometries were developed, and then their seismic responses to both Ricker wavelets and seismic records were simulated by using an explicit finite difference scheme. The definition of deep/shallow basin, the precondition for the observation of prominent surface waves and the influential area of edge effects of the shallow basin were discussed quantitatively in this study. The followings were concluded: in the common velocity contrast range (v _s1/v _s2 < 10), the fundamental frequency a basin with W/H > 3.0 can be estimated approximately by 1D theory. The complexity of peak ground acceleration distribution pattern, the width of the most affected section as well as the amplitude of ground motion in the Edge Region increase with incident frequency. Prominent surface waves can only be observed when the incident wavelength is shorter than the critical wavelength λ _c . The interaction between incident wave and basin dynamic property plays a dominant role on the peak ground acceleration amplitude while the interaction between incident wave and geometry plays a more significant role on the peak ground acceleration distribution. For very shallow basin, different areas along the basin width are affected to different extents. Only a limited area close to the basin edge is influenced significantly. It is more feasible to propose spectral aggravation factor for different surface zones respectively than a uniform constant as a tool to calibrate the 1D-based design spectrum so as to take the basin effects into account.     

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.     

Evapotranspiration from citrus trees growing in sandy soil under drip irrigation with saline water

An experiment on evapotranspiration from citrus trees under irrigation with saline water was carried out for 4 months. Two lysimeters planted with a citrus tree in the green house were used. One lysimeter was irrigated with saline water (NaCl and CaCl₂ of 2000 mg/L equivalence,EC = 3.8 dS/m, SAR = 5.9) and the other was irrigated with freshwater using drip irrigation. The applied irrigation water was 1.2 times that of the evapotranspiration on the previous day. Evapotranspiration was calculated as the change in lysimeter weight recorded every 30 minutes. The lysimeters were filled with soil with 95.8% sand. The results of the experiment were as follows. (i) The evapotranspiration from citrus tree was reduced after irrigation with saline water. The evapotranspiration returns to normal after leaching. However it takes months to exhaust the salt from the tree. (ii) To estimate the impact of irrigation with saline water on the evapotranspiration from citrus trees, the reduction coefficient due to salt stress (K_s) was used in this experiment. Evapotranspiration under irrigation with saline water (ET _s ) can be calculated from evapotranspiration under irrigation with freshwater (ET) by the equationET _s =K _s × ET. K_s can be expressed as a function ofEC _sw . (iii) The critical soil-water electrical conductivity (EC _sw ) is 9.5 dS/m, beyond which adverse effects on evapotranspiration begin to appear. IfEC _sw can be controlled at below 9.5 dS/m, saline water can be safely used for irrigation.     

Three-dimensional admittance analysis of lithospheric elastic thickness over the Louisville Ridge

Using bathymetry and altimetric gravity anomalies, a 1° 9 1° lithospheric effective elastic thickness(Te) model over the Louisville Ridge and its adjacent regions is calculated using the moving window admittance technique. For comparison, three bathymetry models are used: general bathymetric charts of the oceans, SIO V15.1,and BAT_VGG. The results show that BAT_VGG is more suitable for calculating T e than the other two models. T e along the Louisville Ridge was re-evaluated. The southeast of the ridge has a medium Te of 10–20 km, while Te increases dramatically seaward of the Tonga-Kermadec trench as a result of the collision of the Pacific and IndoAustralian plates.     

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.     

Experimental study on the shear behavior of the interface between cushion materials and the concrete raft

Cushion is a layer of granular materials between the raft and the ground. The shear behavior of the interface between the cushion and the raft may influence the seismic performance of the superstructure. In order to quantify such influences, horizontal shear tests on the interfaces between different cushion materials and concrete raft under monotonic and cyclic loading were carried out. The vertical pressure P_v, material type and cushion thickness h_c were taken as variables. Conclusions include: 1) under monotonic loading, P_v is the most significant factor; the shear resistance P_(hmax) increases as P_v increases, but the normalized factor of resistance μ_n has an opposite tendency; 2) for the materials used in this study, μ_n varies from 0.40 to 0.70, the interface friction angle δ_s varies from 20° to 35°, while u_(max) varies from 3 mm to 15 mm; 3) under cyclic loading, the interface behavior can be abstracted as a "three-segment" back-bone curve, the main parameters include μ_n, the displacement u_1 and stiffness K_1 of the elastic stage, the displacement u_2 and stiffness K_2 of the plastic stage; 4) by observation and statistical analysis, the significance of different factors, together with values of K_1, K_2 and μ_n have been obtained.     

Anomalous variations in the ionospheric <Emphasis Type="Italic">F</Emphasis><Subscript>2</Subscript>-layer structure at geomagnetic midlatitudes of the Southern and Northern hemispheres at the transition from summer to winter conditions under low solar activity

The structure and dynamics of the ionosphere and plasmasphere at low solar activity under quiet geomagnetic conditions on January 15–17, 1985, and July 10–13, 1986, over Millstone Hill station and Argentine Islands ionosonde, the locations of which are approximately magnetically conjugate, have been theoretically calculated. The detected correction of the model input parameters makes it possible to coordinate the measured and calculated anomalous variations in the electron density NmF2 at the height hmF2 of the ionospheric F2 layer over Argentine Islands ionosonde as well as the calculated and measured values of NmF2 and electron temperature at the hmF2 height over Millstone Hill station. It has been shown that vibrationally excited N₂ and O₂ molecules almost do not influence the formation of the winter anomaly under the conditions of low solar activity. A difference between the influence of electronically excited O⁺ on N _e ions under winter and summer conditions forms not more than 11% of the N _e winter anomaly event in the F ₂ layer and topside ionosphere. The model without electronically excited O⁺ ions reduces the duration of the N _e winter anomaly event. It has been shown that the seasonal variations in the composition of the neutral atmosphere form mainly the NmF2 winter anomaly event over the Millstone Hill radar at low solar activity.     

Mantle upwelling beneath Madagascar: evidence from receiver function analysis and shear wave splitting

Crustal receiver functions have been calculated from 128 events for two three-component broadband seismomenters located on the south coast (FOMA) and in the central High Plateaux (ABPO) of Madagascar. For each station, crustal thickness and V _p /V _s ratio were estimated from H- κ plots. Self-consistent receiver functions from a smaller back-azimuthal range were then selected, stacked and inverted to determine shear wave velocity structure as a function of depth. These results were corroborated by guided forward modeling and by Monte Carlo error analysis. The crust is found to be thinner (39 ± 0.7 km) beneath the highland center of Madagascar compared to the coast (44 ± 1.6 km), which is the opposite of what would be expected for crustal isostasy, suggesting that present-day long wavelength topography is maintained, at least in part, dynamically. This inference of dynamic support is corroborated by shear wave splitting analyses at the same stations, which produce an overwhelming majority of null results (>96 %), as expected for vertical mantle flow or asthenospheric upwelling beneath the island. These findings suggest a sub-plate origin for dynamic support.     

Ionospheric <Emphasis Type="Italic">N</Emphasis><Subscript>e</Subscript> disturbances before 2007 Pu’er,Yunnan, China,earthquake

Selecting three half orbits near the epicenter of Pu’er earthquake, we analyzed the N_e data recorded in their revisited orbits during a year before this earthquake, and extracted N_e precursors. The results show that: ① There are significant seasonal variations of ionospheric N_e in night time, which exhibit different shapes respectively in four seasons; ② There are three main shapes of N_e: single-peak, saddle-shaped and even-shaped, all of which may occur in four seasons, but each season with its typical shape relatively; ③ Spatial images of N_e showed high values near the epicenter in 30 days before the earthquake, and there is a good correlation between anomaly and distribution of earthquake in space and time, which reflects that these spatial anomalies were indeed concerned with the earthquake; ④ There shows a certain similarity of the N_e curves among revisited orbits, which can provide background information for distinguishing and identification of seismic anomaly.     

Shear Wave Velocity Estimates through Combined Use of Passive Techniques in a Tectonically Active Area

We made an attempt to assess the shear wave velocity values V _S and, to a lesser extent, the V _P values from ambient noise recordings in an array configuration. Five array sites were situated in the close proximity to borehole sites. Shear wave velocity profiles were modeled at these five array sites with the aid of two computational techniques, viz. spatial autocorrelation (SPAC) and H/V ellipticity. Out of these five array sites, velocity estimates could be reliably inferred at three locations. The shear wave velocities estimated by these methods are found to be quite consistent with each other. The computed V _S values up to 30 m depth are in the range from 275 to 375 m/s in most of the sites, which implies prevalence of a low velocity zone at some pocket areas. The results were corroborated by evidence of site geology as well as geotechnical information.