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.     

Optimal placement of dampers and actuators based on stochastic approach

A general method is developed for optimal application of dampers and actuators by installing them at optimal location on seismic-resistant structures. The study includes development of a statistical criterion, formulation of a general optimization problem and establishment of a solution procedure. Numerical analysis of the seismic response in time-history of controlled structures is used to verify the proposed method for optimal device application and to demonstrate the effectiveness of seismic response control with optimal device location. This study shows that the proposed method for the optimal device application is simple and general, and that the optimally applied dampers and actuators are very efficient for seismic response reduction.     

Models of passive margin inversion: implications for the Rhenohercynian fold-and-thrust belt, Belgium and Germany

Positive tectonic inversion is related to the transmission of compressional stresses along a décollement into the foreland of an orogenic zone. This stress and strain concentration in regions remote from the main orogenic front is commonly related to the presence of pre-existing rheological heterogeneities such as normal syn-depositional faults. During inversion, these pre-existing normal faults are reactivated as reverse faults. Tectonic inversion in the Rhenohercynian fold-and-thrust belt during the Variscan Orogeny shows that inversion is likely synchronous with the onset of collision in the hinterland. Here, we present the results of a simplified thermo-mechanical model (STM) which allows one to study strain partitioning between two orogenic zones. We show that, if the two orogenic zones have the same mechanical properties, the viscosity of the décollement, which links them, controls the initial strain partitioning. During subsequent finite shortening, erosional processes determine the partitioning of strain rate. The presence of a weak structure in the inverted zone and of a low-viscosity décollement leads to initial strain concentration in the inverted track rather than in the collision zone and a progressive decrease in strain partitioning between the two orogenic zones. The STM results are in good agreement with results of a 2D finite-element model. We conclude that, in the western part of the Rhenohercynian Massif, simultaneous uplift and deformation within the Mid-German Crystalline Rise (the main collision zone) and the Ardenne Anticlinorium (the inverted zone) lead to interpreting this orogenic event as a case of vice tectonic rather than the propagation of a ‘wave of folding’ towards the Variscan front, as suggested by previous authors.     

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.     

Model-independent Travel-time Attributes for 2-D, Finite-offset Multicoverage Reflections

— In this paper, we provide a 5-parameter stacking formula to transform 2-D prestack data into a particular common-offset section. This requires the knowledge of the near-surface velocity only and it is expected that ray theory holds to describe primary reflections. The earth model can be arbitrarily inhomogeneous. The new stacking approach can be viewed as a generalization of the 3-parameter common-reflection-surface (CRS) stack, by which 2-D multicoverage data are stacked into a simulated zero-offset section. The new 5-parameter formula can handle P-P, P-S and S-S reflections.     

Microseismicity Induced by Heavy Rainfall Around Flooded Vertical Ore Veins

—?Microseismicity (M?M?M?M??6 were also monitored at a crustal movement monitoring station located several hundred meters from the veins. It was found that the opening of the vertical ore veins primary led to significant strain and tilt, but not to seismicity, because the delay and the longer duration of the seismicity were significant. Most seismic events involve thrusting mechanisms that are consistent with the present stress state of E-W-oriented tectonic compression, but are not consistent with the opening of the deepest ore vein. Interstingly, all the events within a few months of the heavy rainfall occurred near the faults that offset the deepest ore veins, wheareas all those events located away from the deepest ore veins occurred many months after the heavy rainfall. Consequently, the delayed diffusion of water appears to have played a dominant role in reducing rock strength, which led to seismicity in the Ikuno mine.     

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.