Analytical Model for Permeability Evolution in Microcracking Rock

— Analytical expressions to predict the enhancement of permeability due to stress-induced microcracking in initially low porosity rock are presented. A fracture mechanical microcrack model is employed to derive integrated effective hydraulic variables as a function of stress, which are then used to calculate the evolution of permeability using the statistically-based Dienes model. The model enables determination of permeability enhancement as a function of two loading parameters and three material parameters. Results are in reasonable agreement with experimental measurements and indicate that appreciable increases in permeability can be anticipated during brittle failure. The analytical nature of the model makes it easily incorporatable into numerical models that require quantification of the permeability evolution as a function of stress, for which there is currently no law.     

Evaluation of liquefaction potential using neural-networks and CPT results

In this research, a reliable Cone Penetration Test data set was gathered with a wide range of parameters. This data was incorporated in a Neural-Networks computer software called STATISTICA Neural-Networks. The back propagation algorithm with a multilayer perceptron network is utilized to analyze the liquefaction occurrence in different sites. In this study, different sets of effective parameters for the neural-network analyses are selected such that to reduce the noise and to obtain more accurate results.Considering the relative importance of effective parameters in liquefaction assessment, it is indicated that σ₀, σ′₀ together play a more important role than what previously was assumed and hence the relative importance of the q_c and seismic parameters are decreased compared with the previous works. The results presented here have more accuracy than previous works while at the same time, the range of the parameters used in this study is much wider than what was previously used. This range of parameters makes the proposed method applicable for practical purposes.     

Onshore and offshore seismic and lithostratigraphic analysis of a deeply incised Quaternary buried valley system in the Northern Netherlands

High-resolution seismic data (onshore and offshore), geophysical borehole data as well as detailed lithofacies from airlift boreholes were acquired in northern Netherlands on and around the island of Ameland. Marine and land seismic data combined with information from land boreholes have been explored with the objective of providing a sedimentary model. Qualitative seismic facies analysis of the valley fill commonly shows a thin unit with high amplitude reflectors at the base. Thick units of variable seismic facies (transparent to high amplitude) occur higher up in the sequence. Onlap is common at mid–upper levels within the sandy valley fill (with clay in mm layering), and a transparent seismic facies, corresponding to firm clays, is common at the top. Almost all lithological unit boundaries recognised within core parameters correspond with seismic unconformities within error margins. Subunits contain multiple cyclical trends in gamma ray and grain size. Cyclical trends show lower order fluctuations in gamma radiation on a scale of less than 1 m. Gamma-ray pattern variability between units, e.g. in general coarsening-up or fining-up units, suggests migration of subaqueous outwash fans or ice margin fluctuations. Seismic results could support a headward excavation and backfilling process suggested by Praeg [Morphology, stratigraphy and genesis of buried Elsterian tunnel valleys in the southern North Sea basin [PhD thesis]: University of Edinburgh, 207 pp.; Journal of Applied Geophysics, (this volume)] as being responsible for the formation of buried valleys. On a lithological scale, a more complicated, detailed and cyclical pattern arises. Catastrophic processes are considered unlikely as being responsible for the infill sequence because of the observed small-scale facies variability and because of the presence of diamicton layers. Diamicton layers at the base of basal unconformities as well as higher in sequence could suggest subglacial deformation by grounded ice before and during the valley-fill process.     

Estimating the orientation error of the Dahan Downhole Accelerometer using the maximum cross-correlation coefficient between the observed and synthetic waves

A new algorithm to correct the orientation error of the accelerometerat the Dahan Downhole Array, Hualien, Taiwan is developed. This algorithmconsists of three stages: (1) rotating two horizontal ground motions on thefree surface to the SH-SVdirection and SH axis offers a reference direction.(2) computing the synthetic downhole SH waves at a downhole station and (3)searching a rotation angle for downhole observation that yields a best waveformmatch between the synthetic and observed downhole seismograms. At this point, the rotated angle corresponding to the best waveform match can be considered as the orientation error. We selected five earthquakes with good data qualityfor analysis. Results show that this algorithm gives a more stable estimationthan a conventional method because it allows the selection of data from a wider time window for analyses. The estimated orientation error of the accelerometers at the Dahan Downhole Array after the 1999 reinstallation are40°, 114° and 285° at depths of 50, 100 and 200 m, respectively.     

An Mw based earthquake Catalogue for central, northern and northwestern Europe using a hierarchy of magnitude conversions

Data from 25 local catalogues and 30special studies of earthquakes in central,northern and northwestern Europe have beenincorporated into a Databank. The dataprocessing includes discriminating eventtypes, eliminating fake events and dupletsand converting different magnitudes andintensities to M_w if this is not givenby the original source. The magnitudeconversion is a key task of the study andimplies establishment of regressionequations where no local relations exist.The Catalogue contains tectonic events fromthe Databank within the area44°N–72°N,25°W–32°E and the time period1300–1993. The lower magnitude level forthe Catalogue entries is setat M_w == 3.50. The area covered by thedifferent catalogues are associated withpolygons. Within each polygon only datafrom one or a small number of the localcatalogues, supplemented by data fromspecial studies, enter the Catalogue. Ifthere are two or more such catalogues orstudies providing a solution for an event,a priority algorithm selects one entry forthe Catalogue. Then M_w is calculatedfrom one of the magnitude types, or frommacroseismic data, given by the selectedentry according to another priority scheme.The origin time, location, M_w magnitude and reference are specified for eachentry of the Catalogue. So is theepicentral intensity, I₀, if providedby the original source. Following thesecriteria, a total of about 5,000earthquakes constitute the Catalogue.Although originally derived for the purposeof seismic hazard calculation within GSHAP,the Catalogue provides a data base for manytypes of seismicity and seismic hazardstudies.     

Seismic scattering at the top of the mantle Transition Zone

We document strong seismic scattering from around the top of the mantle Transition Zone in all available high resolution explosion seismic profiles from Siberia and North America. This seismic reflectivity from around the 410 km discontinuity indicates the presence of pronounced heterogeneity in the depth interval between 320 and 450 km in the Earth’s mantle. We model the seismic observations by heterogeneity in the form of random seismic scatterers with typical scale lengths of kilometre size (10-40 km by 2-10 km) in a 100-140 km thick depth interval. The observed heterogeneity may be explained by changes in the depths to the α-β-γ spinel transformations caused by an unexpectedly high iron content at the top of the mantle Transition Zone. The phase transformation of pyroxenes into the garnet mineral majorite probably also contributes to the reflectivity, mainly below a depth of 400 km, whereas we find it unlikely that the presence of water or partial melt is the main cause of the observed strong seismic reflectivity. Subducted oceanic slabs that equilibrated at the top of the Transition Zone may also contribute to the observed reflectivity. If this is the main cause of the reflectivity, a substantial amount of young oceanic lithosphere has been subducted under Siberia and North America during their geologic evolution. Subducted slabs may have initiated metamorphic reactions in the original mantle rocks.     

Explosion seismic reflections from the Earth’s core

Recent seismological studies have presented evidence for the existence of a layer with ultra-low seismic velocities at the core-mantle boundary at ca. 2900 km depth. We report high-amplitude, high-frequency, and laterally coherent seismic arrivals from three nuclear explosions in Siberia. With recording station intervals of 15 km, the seismic phases are readily correlated and show the presence of a thin, ultra-low velocity zone in a region where it was not previously reported. The duration and complexity of the arrivals are inconsistent with a simple core-mantle boundary and require a hitherto unidentified, kilometre-scale, fine structure in the ultra-low velocity zone. The observations may be explained by a ca. 7 km thick, two-layer, ultra-low velocity zone with exceptional low velocities, which indicate the presence of high percentages of melt (>15%), in particular in the lower part of the zone. Waveform variation implies lateral change in the thickness and physical properties of the ultra-low velocity zone with a wavelength of less than 100 km.     

Evaluating foundation mass,damping and stiffness by the least‐squares method

This paper discusses how to use the three‐dimensional (3D) time‐domain finite‐element method incorporating the least‐squares method to calculate the equivalent foundation mass, damping and stiffness matrices. Numerical simulations indicate that the accuracy of these equivalent matrices is acceptable when the applied harmonic force of 1+sine is used. Moreover, the accuracy of the least‐squares method using the 1+sine force is not sensitive to the first time step for inclusion of data. Since the finite‐element method can model problems flexibly, the equivalent mass, damping and stiffness matrices of very complicated soil profiles and foundations can be established without difficulty using this least‐squares method. Copyright © 2003 John Wiley & Sons, Ltd.     

An Efficient Numerical Method for Computing Synthetic Seismograms for a Layered Half-space with Sources and Receivers at Close or Same Depths

— It is difficult to compute synthetic seismograms for a layered half-space with sources and receivers at close to or the same depths using the generalized R/T coefficient method (Kennett, 1983; Luco and Apsel, 1983; Yao and Harkrider, 1983; Chen, 1993), because the wavenumber integration converges very slowly. A semi-analytic method for accelerating the convergence, in which part of the integration is implemented analytically, was adopted by some authors (Apsel and Luco, 1983; Hisada, 1994, 1995). In this study, based on the principle of the Repeated Averaging Method (Dahlquist and Björck, 1974; Chang, 1988), we propose an alternative, efficient, numerical method, the peak-trough averaging method (PTAM), to overcome the difficulty mentioned above. Compared with the semi-analytic method, PTAM is not only much simpler mathematically and easier to implement in practice, but also more efficient. Using numerical examples, we illustrate the validity, accuracy and efficiency of the new method.     

Environmental Emission of Mercury During Gold Mining by Amalgamation Process and its Impact on Soils of Gympie, Australia

— The aims of this study were to estimate the total amount of mercury released to the environment during 60 years of gold mining (1867–1926) at Gympie, Queensland, Australia and to measure the mercury levels in soil samples surrounding the mining activity. We estimated that 1902 tonnes of mercury was released to the environment and about 1236 tonnes of which was released to the air. The mean mercury in the soil samples in the vicinity of the Scottish battery varied from 1.07 to 99.26 μg g^?1 as compared to 0.075 μg g^?1 as background mercury concentrations. The maximum mercury concentration measured in sediments of the Langton Gully was 6.12 μg g^?1. These results show that large amount of mercury was used in this area during gold mining. Since mining is active in the area and Langton Gully flows into Mary River, we therefore, recommend that mercury concentration in air and fish should be monitored.