Optimized Seismic Threshold Monitoring – Part 2: Teleseismic Processing

v--vThis second paper (Part 2) pertaining to optimized site-specific threshold monitoring addresses the application of the method to regions covered by a teleseismic or a combined regional-teleseismic network. In the first paper (Part 1) we developed the method for the general case, and demonstrated its application to an area well-covered by a regional network (the Novaya Zemlya nuclear test site). In the present paper, we apply the method to the Indian and Pakistani nuclear test sites, and show results during the periods of nuclear testing by these two countries in May 1998. Since the coverage by regional stations in these areas is poor, an optimized approach requires the use of selected, high-quality stations at teleseismic distances.¶To optimize the threshold monitoring of these test sites, we use as calibration events either one of the nuclear explosions or a nearby earthquake. From analysis of the calibration events we derive values for array beamforming steering delays, filter bands, short-term averages (STA) lengths, phase travel times (P waves), and amplitude-magnitude relationships for each station. By applying these parameters, we obtain a monitoring capability of both test sites ranging from m_b 2.8-3.0 using teleseismic stations only. When including the nearby Nilore station to monitor the Indian tests, we show that the threshold can be reduced by about 0.4 magnitude units. In particular, we demonstrate that the Indian tests on 13 May, 1998, which were not detected by any known seismic station, must have corresponded to a magnitude (m_b) of less than 2.4.¶We also discuss the effect of a nearby aftershock sequence on the monitoring capability for the Pakistani test sites. Such an aftershock sequence occurred in fact on the day of the last Pakistani test (30 May, 1998), following a large (m_b 5.5) earthquake in Afghanistan located about 1100 km from the test site. We show that the threshold monitoring technique has sufficient resolution to suppress the signals from these interfering aftershocks without significantly affecting the true peak of the nuclear explosion on the threshold trace.     

Least square approach to simulate wave propagation in irregular profiles using the indirect boundary element method

In earthquake engineering and seismology it is of interest to know the surface motion at a given site due to the incoming and scattered seismic waves by surface geology. This can be formulated in terms of diffraction of elastic waves and then the indirect boundary element method (IBEM) for dynamic elasticity is used. It is based on the explicit construction of diffracted waves at the boundaries from which they radiate. This provides the analyst with insight on the physics of diffraction. The IBEM has been applied to study the amplification of elastic waves in irregular soil profiles. From the strong or weak satisfaction of boundary conditions and a simple analytical discretization scheme a linear system of equations for the boundary sources is obtained. Here, we explore the use of a weak discretization strategy with more collocation points than force densities. The least squares enforcement of boundary conditions leads to a system with reduced number of unknowns. This approach naturally allows one to use both coarser and finer boundary discretizations for smooth and rapidly varying profiles, respectively. A well studied semicircular canyon under incident P or SV in-plane waves is used to calibrate this method. Several benefits are obtained using mixed meshing that leads to the least squares condensation of the IBEM.     

Cramer-Rao lower bounds for sonar broad-band modulation parameters

A role of passive sonar signal processing is the detection and estimation of the parameters associated with amplitude modulated broad-band signals. An example of such signals is propeller noise. Discrete frequency lines occur at the rotational frequency of the propulsion shaft and at the blade frequency. This correspondence provides expressions for the Cramer-Rao lower bounds for the estimates of broad-band signal power, modulation level, modulation frequency, and modulation phase. It is shown that for low broad-band-signal-to-broad-band-noise ratios, the estimates of power and modulation level are uncoupled from the estimates of modulation frequency and phase     

Edge waves on the New Zealand east coast

Tide gauge records from several New Zealand east coast ports show persistent oscillations in sea level with a period of 2.4–2.7 h. These are superimposed on local resonant oscillations, and appear to be associated with edge wave excitation on the continental shelf, possibly by meteorological disturbances. A shelf response is also excited by tsunamis.     

Exploratory seismic site response surveys in a complex geologic area: a case study from Mt. Etna volcano (southern Italy)

A preliminary study targeting to evaluate the local seismic response was performed in the eastern flank of Mt. Etna (southern Italy) using ambient noise measurements. The obtained spectral ratios were subdivided through cluster analysis into different classes of fundamental frequency permitting to draw an iso-frequency contour map. The analysis set into evidence the extreme heterogeneity of lava sequences, which makes difficult to identify a single seismic bedrock formation. Another important outcome, concerning the local seismic effects in terms of frequency and azimuth, is the important role played by the fracture fields associated with the main structural systems of the area. The existence of two zones with strong directional effects striking WNW–ESE and NW–SE, nearly orthogonal to the orientation of the main fracture fields, corroborate such hypothesis.     

The Y-3 tephra: A Last Glacial stratigraphic marker for the central Mediterranean basin

The paper reviews the existing data on the Y-3 tephra layer, first recognised in the Ionian Sea (Mediterranean basin). The collection and collation of old and new data on distal tephra occurrences in terrestrial, marine and lacustrine successions indicate that the Y-3 layer is dispersed over a wide area of the central Mediterranean basin. The peculiar homogeneous chemical composition of this layer makes its recognition rather straightforward and permits it being distinguished from other stratigraphically adjacent tephras. The best age estimate for the Y-3 layer of ca 30–31 cal ka BP, its peculiar stratigraphic position close to the Marine Isotope Stage 3/2 transition or Heinrich Event 3 onset, as well as its wide dispersion makes this layer an important marker to link and date late Pleistocene terrestrial and marine archives of the central Mediterranean basin.     

Estimation of Site Effects in Beijing City

— For the realistic modeling of the seismic ground motion in lateral heterogeneous anelastic media, the database of 3-D geophysical structures for Beijing City has been built up to model the seismic ground motion in the City, caused by the 1976 Tangshan and the 1998 Zhangbei earthquakes. The hybrid method, which combines the modal summation and the finite-difference algorithms, is used in the simulation. The modeling of the seismic ground motion, for both the Tangshan and the Zhangbei earthquakes, shows that the thick Quaternary sedimentary cover amplifies the peak values and increases the duration of the seismic ground motion in the northwestern part of the City. Therefore the thickness of the Quaternary sediments in Beijing City is the key factor controling the local ground effects. Four zones are defined on the base of the different thickness of the Quaternary sediments. The response spectra for each zone are computed, indicating that peak spectral values as high as 0.1 g are compatible with past seismicity and can be well exceeded if an event similar to the 1697 Sanhe-Pinggu occurs.     

Nonlinear transfer function modelling of beach morphology at Duck,North Carolina

This paper presents a simple nonlinear data-based modelling approach for predicting the beach profile volume at Duck, North Carolina, USA. The state-dependent parameter form of the general transfer function (SDP TF) model is used to describe nonlinearity influencing these morphological data in two case examples. Case 1 investigates the nonlinearity associated with the dependency of wave forcing on the preceding beach volume. Case 2 investigates the ability to model the variables within the well-known diffusion equation for beach volume using this data-based approach. The results of this study show that the SDP TF approach can be used successfully to develop statistically robust models for describing nonlinearity in beach morphological systems. Furthermore, these models are shown to predict the beach volumes over both short (1 month ahead) and long (2 years ahead) time periods, and thus show great potential for practical applications in coastal zone management and engineering.     

Nitrate contamination in a shallow urban aquifer in East Ukraine: evidence from hydrochemical,stable isotopes of nitrate and land use analysis

A combined hydrochemical and stable isotope approach was used to investigate the origin of nitrate in the shallow unconfined groundwater of Kharkiv city, Eastern Ukraine. The contamination was investigated in the context of land use within the catchment area. The observed enrichment of sulfate, chloride and nitrate suggests significant groundwater contamination in the shallow urban aquifer, which is widely used as drinking water source for the urban population. Characteristic nitrate/chloride ratios as well as stable isotope ratios (N and O) of nitrate in the most contaminated springs confirmed that septic waste from leaky sewer systems was the main source of nitrate contamination in the groundwater. Nitrate contamination is linked to the type of land use and sewage treatment regime in the catchment area. It is also modulated by the regional hydrogeology, which determines the susceptibility of a given aquifer toward groundwater pollution. A more quantitative assessment of nitrate sources based on the nitrate isotope analysis alone is rather difficult. However, our study confirms that the combination of hydrochemical tracers, robust land-use analysis and nitrate stable isotope measurements represents a valuable approach to identify the origin of the nitrate contamination.