Analysis of the resolution function in seismic prestack depth imaging

We consider the problem of estimating subsurface quantities such as velocity or reflectivity from seismic measurements. Because of a limited aperture and band-limited signals, the output from a seismic prestack reconstruction method is a distorted or blurred image. This distortion can be computed using the concept of resolution function, which is a quantity readily accessible in the Fourier space of the model. The key parameter is the scattering wavenumber, which at a particular image point is defined by the incident and scattered ray directions in a given background model. Any location in any background model can be considered. In general, the resolution function will depend on the following four quantities: the background velocity model, the frequency bandwidth, the wavefield type and the acquisition geometry.
We first establish the resolution function for a general scattering model assuming local reaction. We then adapt this result for two well-known scattering models: Born and Kirchhoff. For each of these approximations the corresponding resolution function is derived and discussed. Finally, by employing a simple synthetic data example we demonstrate the ability of the resolution function to predict the image distortions.     

A methodology for the evaluation of long-term volcanic risk from pyroclastic flows in Campi Flegrei (Italy)

The volcanological history of Campi Flegrei suggests that the most frequent eruptions are characterized by the emplacement of pyroclastic flow and surge deposits erupted from different vents scattered over a 150-km² caldera. The evaluation of volcanic risk in volcanic fields is complex because of the lack of a central vent. To approach this problem, we subdivided the entire area of Campi Flegrei into a regular grid and evaluated the relative spatial probability of opening of vents based on geological, geophysical and geochemical data. We evaluated the volcanic risk caused by pyroclastic flows based on the formula proposed by UNESCO (1972), R=H×V×Va, where H is the hazard, V is the vulnerability and Va is the value of the elements at risk. The product H×V was obtained by performing simulations of type eruptions centered in each cell of the grid. The simulation is based on the energy cone scheme proposed by Sheridan and Malin [J. Volcanol. Geotherm. Res. 17 (1983) 187–202], hypothesizing a column collapse height of 100 m for eruptions of VEI=3 and 300 m for eruptions of VEI=4 with a slope angle of 6°. Each simulation has been given the relative probability value associated with the corresponding cell. We made use of the GIS software ArcView 3.2 to evaluate the intersection between the energy cone and the topography. The superposition of the areas invaded by pyroclastic flows (124 simulations for VEI=3 and 37 for VEI=4) was used to obtain the relative hazard map of the area. The relative volcanic risk map is obtained by superimposing the urbanization maps.     

A rifted inside corner massif on the Mid-Atlantic Ridge at 5°S

The structure of the Mid-Atlantic Ridge at 5°S was investigated during a recent cruise with the FS Meteor. A major dextral transform fault (hereafter the 5°S FZ) offsets the ridge left-laterally by 80 km. Just south of the transform and to the west of the median valley, the inside corner (IC – the region bounded by the ridge and the active transform) is marked by a major massif, characterized by a corrugated upper surface. Fossil IC massifs can also be identified further to the west. Unusually, a massif almost as high as the IC massif also characterizes the outside corner (OC) south of the inactive fracture zone and to the east of the median valley. This OC massif has axis-parallel dimensions identical to the IC massif and both are bounded on their sides closest to the spreading axis by abrupt, steep slopes. An axial volcanic ridge is well developed in the median valley both south of the IC/OC massifs and in an abandoned rift valley to the east of the OC massif, but is absent along the new ridge-axis segment between the IC and OC massifs. Wide-angle seismic data show that between the massifs, the crust of the median valley thins markedly towards the FZ. These observations are consistent with the formation of the OC massif by the rifting of an IC core complex and the development of a new spreading centre between the IC and OC massifs. The split IC massif presents an opportunity to study the internal structure of the footwall of a detachment fault, from the corrugated fault surface to deeper beneath the fault, without recourse to drilling. Preliminary dredging recovered gabbros from the scarp slope of the rifted IC massif, and serpentinites and gabbros from the intersection of this scarp with the corrugated surface. This is compatible with a concentration of serpentinites along the detachment surface, even where the massif internally is largely plutonic in nature.     

A method to estimate hydraulic conductivity while drilling

A field test and analysis method has been developed to estimate the vertical distribution of hydraulic conductivity in shallow unconsolidated aquifers. The field method uses fluid injection ports and pressure transducers in a hollow auger that measure the hydraulic head outside the auger at several distances from the injection point. A constant injection rate is maintained for a duration time sufficient for the system to become steady state. Exploiting the analogy between electrical resistivity in geophysics and hydraulic flow two methods are used to estimate conductivity with depth: a half-space model based on spherical flow from a point injection at each measurement site, and a one-dimensional inversion of an entire dataset.

The injection methodology, conducted in three separate drilling operations, was investigated for repeatability, reproducibility, linearity, and for different injection sources. Repeatability tests, conducted at 10 levels, demonstrated standard deviations of generally less than 10%. Reproducibility tests conducted in three, closely spaced drilling operations generally showed a standard deviation of less than 20%, which is probably due to lateral variations in hydraulic conductivity. Linearity tests, made to determine dependency on flow rates, showed no indication of a flow rate bias. In order to obtain estimates of the hydraulic conductivity by an independent means, a series of measurements were made by injecting water through screens installed at two separate depths in a monitoring pipe near the measurement site. These estimates differed from the corresponding estimates obtained by injection in the hollow auger by a factor of less than 3.5, which can be attributed to variations in geology and the inaccurate estimates of the distance between the measurement and the injection sites at depth.     

Geoecological Characteristics of Lake Seliger

Data of hydrological and hydrochemical observations conducted in August 2000 at Lake Seliger and its major tributaries are presented. The analyzed characteristics included: O₂, H₂S, CH₄, microelements; mineral and organic compounds of P and N, ionic composition components (HCO₃ ^–, Cl^–, Na⁺, K⁺, dry residue, pH; color index, electrical conductivity. Bottom sediment samples were analyzed for microelement concentrations, CH₄, oil hydrocarbons. The influence of the basin on the ecological state of the lake was assessed, and tendencies in the evolution of its ecosystem are analyzed. It is shown that the ecological conditions of the lake has not changed significantly since 1960–1991, however, nutrient concentrations in the lake water was found to increase.     

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.     

3-D Viscoelastic FEM Modeling of Crustal Deformation in Northeast Japan

-- As a first step toward establishing a standard earthquake cycle model in Japan, we simulate the crustal deformation during the past 100 years in northeast Japan, using a 3-D FEM based on the kinematic model. Then, we compare the computed results with the observed long-term leveling data and the recent GPS data. On the whole, although the effect of the subducting PAC is dominant, coseismic deformation of the interplate earthquakes can be clearly seen in the inland. Moreover, the postseismic deformation of the earthquakes due to the viscoelastic upper mantle seriously affects the inland movements, and continues for a few decades. Our modeling, including the effects of the interplate earthquakes and the three-dimensional viscoelastic inhomogeneity, reasonably explains the observed movement. Finally, we stress that the viscoelastic effect should be taken into consideration in the analyses, even if no earthquakes occur in the analyzed period.     

The null distribution of sample serial correlation coefficient

 The null distribution of the lag-k sample serial correlation coefficient (r _k , k=1,2,3) was investigated by Monte Carlo simulation. For a time series with normal, exponential, Pearson 3, EV1 (Gumbel), or generalized Pareto (GP) distribution type, the null distribution of its r _k can be approximated by the normal distribution with mean −1/(n−k) and variance 1/(n−1). But for a time series with the lognormal, EV2 or EV3 (Weibull) distribution type, the null distribution of r _k is skewed distributed. In such cases, a simulation technique is suggested to construct percentile confidence intervals at a given significance level.     

Analysis of Time Series of the Ufa River Water Pollution by Benz[a]pyrene

Mathematical analysis of time series is shown to be applicable to assessing the factors governing benz[a]pyrene concentration in the Ufa River water. The contributions of the deterministic and random components to the total benz[a]pyrene concentration. Seasonal increase (October–January, March–April) and decrease (February, May–September) in benz[a]pyrene concentrations in the Ufa River are established.     

Reduction of the Trophic Status of a Water Body with the Use of Large-Size Zooplankton

It is shown that large-size zooplankton can be used as a natural filter in flow-through plants for reducing the trophic status of water bodies and improve water transparency. Experiments in Chistye Prudy pond in the period of blue-green algae blooming show large-size zooplankton to reduce phytoplankton biomass in the pond throughout the period of phytoplankton growth. Daphnia magna specimens (with a biomass of 92 to 3450 mg/l) in lothic ecosystems are found to reduce the biomass of dominating Anabaena spirodes blue-green algae when this biomass attains its maximum.