P-waveform analysis for local subduction geometry south of Puget Sound,Washington

The local subduction geometry at a site south of Puget Sound in western Washington is investigated using teleseismicP-waveforms recorded on a three-component event triggered seismograph. The data are processed using source equalization deconvolution in order to isolate locally convertedP-to-S arrivals and stacked to improve the signal-to-noise ratio. Stable arrivals in the radial component indicate an oceanic Moho within the subducted slab at a depth of about 53 km beneath the station. Observed amplitude variations with azimuth in the radial data, as well as qualitative aspects of the tangential data, are used to establish a slab dip of 16° to the southeast. Our results are compatible with previous results from a site 60 km to the west, and further confirm a substantial warp in the regional geometry of the subducted Juan de Fuca plate.     

Fractal properties of simulated bed profiles in coarse-grained channels

Bed roughness characteristics in coarse-grained channels are fairly complex. A hierarchy of roughness elements can be observed, ranging from variable particle sizes and shapes and small-scale sedimentary structures, to large-scale bedforms such as riffle-pool sequences. The effects of these scales of roughness on the flow geometry still remain to be thoroughly investigated. The semivariogram has been suggested in the past as a means of quantifying bed roughness effects on streamflow, as well as for distinguishing between scales of roughness. However, field measurements are rather time-consuming. The low number of bed profiles measured in the field precludes the identification of generally applicable relationships between the statistical properties derived from the semivariograms (such as the Hausdorff dimensions and the scale of autocorrelation corresponding to each fractal band) and the bed configuration itself (geometrical and sedimentological properties). Simulation results of gravel-bed profiles are, therefore, presented in order to complement the original investigation of Robert (1988a). The simulation experiments, based on grain characteristics of sizes and shapes and on morphological properties of small-scale bedforms, yield very significant information on boundary roughness at the microscale and give insight into the interpretation of empirical semivariograms (derived from field measurements). Bed-material sorting, variable grain shapes, and height and spacing of cluster bedforms control the fractal dimensions obtained from the semivariograms, as well as the location of the break of slope and the range of the process.     

Ground variation of radon 222 for location of hidden structural features. Example of the south of France (Alpes Maritimes)

We report two examples from the south of the French Alps, showing that radon emanation monitored by alpha-sensitive film may be used to locate certain discrete structural features revealed in data collected by remote sensing from a satellite. The variations observed in our data, over a period of several months, are in accordance with atmospheric changes and might correlate with local seismic activity when the detectors are located directly above structural fractures and the magnitude of the seismic event is greater than 2.     

UPb zircon and rutile chronology of Archean greenstone formation and gold mineralization in the Val d'Or region, Quebec

UPb geochronology provides an absolute time framework for the evolution of the Sigma gold deposit and surrounding rocks at Val d'Or, southern Abitibi subprovince. The Bourlamaque batholith, the largest pluton in the area, gives a 2699.8 ± 1.0 Ma UPb zircon age. This pluton cuts the Val d'Or Formation which hosts the mineralization. A UPb zircon age of 2704.9 ± 1.1 Ma on a felsic volcanic rock, the Colombière “rhyolite”, 13 km east of the mine dates that formation. The gold-bearing quartz vein system at Sigma is hosted by andesites and two generations of porphyry intrusions, all metamorphosed to the greenschist facies. The oldest porphyry (“porphyritic diorite”) shows the same deformation as the volcanic rocks, and has a 2703.7 ± 2.5 Ma UPb zircon age. The porphyritic diorite and volcanic rocks are cut by feldspar-porphyry dykes which post-date regional folding and have a 2694.0 ± 2.2 Ma UPb zircon age.Regional greenschist metamorphism has been dated directly, with a UPb date of 2684 ± 7 Ma on rutile in the Colombière “rhyolite”. The mineralization and hydrothermal alteration in the mine are superimposed on the metamorphic minerals. Hydrothermal rutile, from an alteration halo around the veins in andesite, has a 2599 ± 9 Ma UPb age. Textural evidence clearly indicates that the wall-rock alteration and vein filling are contemporaneous, and hence the vein system and gold mineralization appear to have developed at least 80 m.y. after the formation and metamorphism of host greenstones.     

High-frequency sequences and their stacking patterns: sequence-stratigraphic evidence of high-frequency eustatic cycles

A hierarchy of interpreted eustatic cyclicity in siliciclastic sedimentary rocks has a pattern of superposed cycles with frequencies in the ranges of 9–10 m.y., 1–2 m.y., 0.1–0.2 m.y., and 0.01–0.02 m.y. (second- through fifth-order cyclicity, respectively). Stratigraphic units displaying this cyclicity include composite sequences, sequences, and parasequences. On the Exxon global cycle chart, fundamental third-order cycles (1–2 m.y. average duration) stack into related groups (second-order cycles: 9–10 m.y. duration). A much larger pattern (about 200 m.y.) is interpreted as tectonically controlled eustasy probably related to sea-floor spreading rates.

One and probably two higher orders of cyclicity (fourth-order: 0.1–0.2 m.y.; and fifth-order: 0.01–0.02 m.y.) are now observed in work with well logs, cores, and outcrops in areas of very rapid deposition. These frequencies are in the range of Milankovitch cycles, and may represent part of the Milankovitch hierarchy which has been widely interpreted for cyclical units in carbonate rocks.

High-frequency (fourth-order) sequences, which form at a 0.1–0.2 m.y. cyclicity, have all the stratal attributes of conventional sequences, including constituent parasequences and systems tracts, and play a dominant role controling reservoir, source, and sealing rock distribution. A consistent hierarchy of stratigraphy is observed. Parasequences (probable fifth-order cyclicity) stack into sets to form systems tracts in fourth-order sequences. Groups (sets) of fourth-order sequences are deposited between major third-order boundaries within third-order composite sequences. Sequences in these sets stack in prograding and backstepping patterns to form third-order lowstand, transgressive, and highstand sequence sets.

Third-order sequence boundaries are marked by greater basinward shifts in facies, by larger more widespread incised valleys, and by more extensive onlap than are fourth-order sequence boundaries. Third-order condensed sections commonly are widespread, faunally rich, and widely correlated biozone and mapping markers. Fourth-order sequence analysis helps to understand reservoir, source, and seal distribution at the play and prospect scale. An example from the Gulf of Mexico is discussed.     

Glacial meltwater cooling of the Gulf of Mexico: GCM implications for Holocene and present-day climates

An atmospheric general circulation model, the NCAR CCM, has been used to investigate the possible effects that reduced Gulf of Mexico sea surface temperatures (SST) could have on regional and hemispheric climates. ¹⁸O records and terrestrial evidence indicate at least two major glacial meltwater discharges into the Gulf of Mexico subsequent to the last glacial maximum. It is probable that these discharges reduced Gulf of Mexico SST. We have conducted three numerical experiments, with imposed gulf-wide SST coolings of 3°C, 6°C, and 12°C, and find in all three experiments significant reductions in the North Atlantic storm-track intensity, along with a strong decrease in transient eddy water vapor transport out of the Gulf of Mexico. Surface pressures are higher over the North Atlantic, indicating a reduction of the climatological Icelandic low. The region is generally cooler and drier, with a reduction in precipitation that agrees well with evidence from Greenland ice cores. Other statistically significant changes occur across the Northern Hemisphere, but vary between the three experiments. In particular, warmer, wetter conditions are found over Europe for both the 6°C and 12°C SST reductions, but cooler conditions are found for the 3°C reduction. This indicates a dependence, in both the sign and magnitude of the model response, on the magnitude of the imposed SST anomaly. The results suggest that the present-day North Atlantic storm track is dependent on warm Gulf of Mexico SST for much of its intensity. They also suggest that meltwater-induced coolings may help account, in part, for some of the climatic oscillations that occurred during the last glacial/interglacial transition.     

Nonlinear least-squares method via an isomorphic geometrical setup

The resolution of a nonlinear parametric adjustment model is addressed through an isomorphic geometrical setup with tensor structure and notation, represented by a u-dimensional “model surface” embedded in a flat n-dimensional “observational space”. Then observations correspond to the observational-space coordinates of the pointQ, theu initial parameters correspond to the model-surface coordinates of the “initial” pointP, and theu adjusted parameters correspond to the model-surface coordinates of the “least-squares” point . The least-squares criterion results in a minimum-distance property implying that the vector Q must be orthogonal to the model surface. The geometrical setup leads to the solution of modified normal equations, characterized by a positive-definite matrix. The latter contains second-order and, optionally, thirdorder partial derivatives of the observables with respect to the parameters. This approach significantly shortens the convergence process as compared to the standard (linearized) method.     

Inversion of arrival-times in a region of dilatancy anisotropy

Shear-wave splitting has been identified in many three-component seismograms from two separate field experiments on a section of the North Anatolian Fault in North-West Turkey. These observations are consistent with shear-wave propagation through a zone of extensive-dilitancy anisotropy. A preliminary attempt has been made to confirm this interpretation by simultaneously inverting suites of arrival-times for hypocentral locations and for parameters describing an anisotropic halfspace. Although the inversion procedure is not globally convergent, it is possible to recognize the true solution by systematically varying the initial conditions. Applied to selected data sets, the inversion defines several anisotropic models that fit the data significantly better than a simple isotropic model, and display the anisotropy required by the shear-wave splitting. However, most of these anisotropic models are not superior when they are used to individually locate events in a much larger data set. However, for each experiment, there is a single model that produces clearly superior locations for the larger data sets than those of other anisotropic or simple isotropic models. Both models display similar velocity variations which are characteristic of propagation through distributions of biplanar cracks displaying orthorhombic symmetry. The principal axes of the two models are oriented in similar directions and are within 20° of the principal axis of regional stress derived from fault-plane solutions. The solutions indicate low velocities close to the tensional axis, as would be expected in extensive-dilatancy anisotropy.