Calibration of Local Magnitude M<Subscript>L</Subscript> in the Azores Archipelago Based on Recent Digital Recordings

—The development of the digital seismic network in the Azores Archipelago during recent years made it possible to obtain the amplitudes (waveform) of recorded motion in a large set of stations. With this new data, maximum amplitudes of the Wood Anderson seismograph are computed, for each station/component, which, together with epicentral distances, allows for the estimation of local magnitude M _L . We used data recorded in 8 digital permanent three-component stations, with inter-stations distances up to 300 km, in the period June 1998 – June 2000, corresponding to a set of 1315 events with magnitude (M _L or M _D ) 2<M<5.8 and epicenters located in the Azores region, to estimate the coefficients of the equation to compute M _L , as well as to determine the corrections to be applied to each station. The new set of parameters, formed by attenuation coefficients and station corrections, were introduced in the calculations of the M _L , leading to smaller dispersions in the analyzed dataset. We also conclude that the attenuation in the first 150 km is similar to the California values, although higher for longer distances.     

Multi-scale Tomography for Crustal <Emphasis Type="Italic">P</Emphasis> and <Emphasis Type="Italic">S</Emphasis> Velocities in Southern California

Seismic tomography is a viable tool in building depth-velocity models in the presence of strong lateral velocity variations. In this study 3-D P- and S-velocity models for the crust of southern California are constrained using more than 1,000,000 P-wave first arrivals and 130,000 S-wave arrivals from local earthquakes. To cope with the uneven distribution of raypaths, a multi-scale tomography is applied with overlapping model cells of different sizes. Within the 300 × 480 × 39 km³ model volume, the smallest cell size is 10 × 10 × 3 km³. During the iterations of velocity updating, earthquake hypocenters are determined using both P and S arrivals, and full 3-D ray tracing is implemented. Except near the edges and in the lower crust, the resultant models are robust according to various tests on the effects of reference models, resolution and signal-to-noise ratio. The tomographic velocities at shallow depths correlate very well with the regional geology of southern California. In the upper crust the P-wave and S-wave models exhibit slow velocities in major sedimentary basins and fast velocities in areas of crystalline rocks. Mid-crustal low velocity zones are present under the Coso Range, San Gabriel Mountains, and a large portion of the Mojave Desert. P- and S-velocity patterns maintain their similarity in the lower crust though the models are less reliable there.     

Building and Testing an <Emphasis Type="Italic">a priori</Emphasis> Geophysical Model for Western Eurasia and North Africa

We construct and evaluate a new three-dimensional model of crust and upper mantle structure in Western Eurasia and North Africa (WENA) extending to 700 km depth and having 1° parameterization. The model is compiled in an a priori fashion entirely from existing geophysical literature, specifically, combining two regionalized crustal models with a high-resolution global sediment model and a global upper mantle model. The resulting WENA1.0 model consists of 24 layers: water, three sediment layers, upper, middle, and lower crust, uppermost mantle, and 16 additional upper mantle layers. Each of the layers is specified by its depth, compressional and shear velocity, density, and attenuation (quality factors, Q _P and Q _S ). The model is tested by comparing the model predictions with geophysical observations including: crustal thickness, surface wave group and phase velocities, upper mantle _n velocities, receiver functions, P-wave travel times, waveform characteristics, regional 1-D velocities, and Bouguer gravity. We find generally good agreement between WENA1.0 model predictions and empirical observations for a wide variety of independent data sets. We believe this model is representative of our current knowledge of crust and upper mantle structure in the WENA region and can successfully be used to model the propagation characteristics of regional seismic waveform data. The WENA1.0 model will continue to evolve as new data are incorporated into future validations and any new deficiencies in the model are identified. Eventually this a priori model will serve as the initial starting model for a multiple data set tomographic inversion for structure of the Eurasian continent.     

Stress Changes Modelled for the Sequence of Strong Earthquakes in the South Iceland Seismic Zone Since 1706

The South Iceland seismic zone is, roughly speaking, situated between two sections of the mid-Atlantic ridge, i.e., the Reykjanes Ridge southwest of Iceland and the Eastern Volcanic Zone on the island. It is a transform zone, where earthquakes are expected to occur on E-W-trending left-lateral shear faults, equivalent to conjugate, N-S-oriented right-lateral, rupture planes. In fact, earthquakes take place on en-échelon N-S-oriented faults, which is indicated by the distribution of main shock intensities, aftershocks as well as by surface fault traces. The stress field continuously generated in the fault zone by opening of the adjacent ridges is computed and superimposed on the stress field changes induced by a series of 13 earthquakes (M 6) between 1706 and 2000. The level of the pre-seismic stress field is analysed as well as the size of the area under high stress. Finally, the post-seismic stress field of June 2000 is analysed, to see where high stresses might have accumulated. The modelling indicates that the rupture planes located on separated parallel N-S-striking zones are dense enough to lead to an area-wide stress release by the series of events. The obtained pre-seismic stress level for most events is high and stable with the exception of situations when several strong shocks occur over a time span of several days, i.e., display typical main shock-aftershock patterns. The size of areas under high stress aside from of the rupture plane, i.e., where no event occurs at the specific time, is of medium to small size.     

New Methodology for Computation of the Earth’s Polarizability from Telluric Borehole Logging Measurements Taken at Different Frequencies

An important approach for computation of the earths polarizability is achieved by developing a novel method with new parameters called Natural Percent Frequency Effect (NPFE) and Natural Metal Factor (NMF) resulting from telluric logging data taken at different frequencies. The essential advantage of the developed approach relies on using natural electrical currents instead of an artificial electrical source usually needed in the standard frequency domain applied in the induced polarization method. Furthermore, a good qualitative correlation has been found between the new parameters and those of the traditional induced polarization method obtained in time and frequency domains. The new method has been tested in well 32 in the Rouez mine in France, where sulphide mineralization is well known to occur.     

Ground Surface Temperature History in Poland in the 16th–20th Centuries Derived from the Inversion of Geothermal Profiles

Ground Surface Temperature (GST) history in Poland was derived from the inversion of temperature-depth profiles in over 20 wells. Temperature histories for the period 1500 A.D. through 1977 A.D. agree well with the instrumental record of the surface-air temperature available for the last two centuries. A statistical correlation of the reconstructed histories (from the well temperature data) with the instrumental record (air temperature) from the homogeneous Warsaw series is high (>0.8). Functional space inversion (FSI) of the temperature data with depth shows that beginning in the early 19th century, temperatures warmed by 0.9 ± 0.1°C following a long period of colder climate before. The last number could be a minimal as higher warming was calculated using a simple model based on surface temperature for the observational period (homogenized Warsaw surface temperature series, Lorenc, 2000) and POM (pre-observational mean; Harris and Chapman, 1998) of –1.53^oC below the 1951–1980 mean temperature level.     

Nonlinear Time Series Analysis of Volcanic Tremor Events Recorded at Sangay Volcano,Ecuador

The assumption that volcanic tremor may be generated by deterministic nonlinear source processes is now supported by a number of studies at different volcanoes worldwide that clearly demonstrate the low-dimensional nature of the phenomenon. We applied methods based on the theory of nonlinear dynamics to volcanic tremor events recorded at Sangay volcano, Ecuador in order to obtain more information regarding the physics of their source mechanism. The data were acquired during 21–26 April 1998 and were recorded using a sampling interval of 125 samples s^–1 by two broadband seismometers installed near the active vent of the volcano. In a previous study Johnson and Lees (2000) classified the signals into three groups: (1) short duration (<1 min) impulses generated by degassing explosions at the vent; (2) extended degassing chugging events with a duration 2–5 min containing well-defined integer overtones (1–5 Hz) and variable higher frequency content; (3) extended degassing events that contain significant energy above 5 Hz. We selected 12 events from groups 2 and 3 for our analysis that had a duration of at least 90 s and high signal-to-noise ratios. The phase space, which describes the evolution of the behavior of a nonlinear system, was reconstructed using the delay embedding theorem suggested by Takens. The delay time used for the reconstruction was chosen after examining the first zero crossing of the autocorrelation function and the first minimum of the Average Mutual Information (AMI) of the data. In most cases it was found that both methods yielded a delay time of 14–18 samples (0.112–0.144 s) for group 2 and 5 samples (0.04 s) for group 3 events. The sufficient embedding dimension was estimated using the false nearest neighbors method which had a value of 4 for events in group 2 and was in the range 5–7 for events in group 3. Based on these embedding parameters it was possible to calculate the correlation dimension of the resulting attractor, as well as the average divergence rate of nearby orbits given by the largest Lyapunov exponent. Events in group 2 exhibited lower values of both the correlation dimension (1.8–2.6) and largest Lyapunov exponent (0.013–0.022) in comparison with the events in group 3 where the values of these quantities were in the range 2.4–3.5 and 0.029–0.043, respectively. Theoretically, a nonlinear oscillation described by the equation ++g(x)=fcost can generate deterministic signals with characteristics similar to those observed in groups 2 and 3 as the values of the parameters ,,f, are drifting, causing instability of orbits in the phase space.     

An investigation of analysis of variance as a tool for exploring regional differences in strong ground motions

The statistical technique known as analysis of variance is applied to a large set of European strong-motion data to investigate whether strong ground motions show a regional dependence. This question is important when selecting strong-motion records for the derivation of ground motion prediction equations and also when choosing strong-motion records from one geographical region for design purposes in another. Five regions with much strong-motion data (the Caucasus region, central Italy, Friuli, Greece and south Iceland) are investigated here. For the magnitude and distance range where there are overlapping data from the five areas (2.50 M_s 5.50, 0 d 35 km) and consequently analysis of variance can be performed, there is little evidence for a regional dependence of ground motions. There is a lack of data from moderate and large magnitude earthquakes (M_s > 5.5) so analysis of variance cannot be performed there. Since there is uncertainty regarding scaling ground motions from small to large magnitudes whether ground motions from large earthquakes are significantly different in different parts of Europe is not known. Analysis of variance has the ability to complement other techniques for the assessment of regional dependence of ground motions.     

Criteria for defining and recognizing the various orders of sequences in outcrop sequence stratigraphy

Sequence stratigraphy which began in the late seventies of the previous century has not only enriched and widened the scope of stratigraphy, but also has been widely used in sedimentary geology, basin analy-sis, exploration of oil and gas, and other fields of ge-ology. In the last two decades and more, sequence stratigraphy, as a new branch of stratigraphical disci-pline, has become a cornerstone of modern stratigra-phy, and has produced more profound influences than all others. The currently …