USGS Science⌕ Search

SEARCH · USGS Science

Results for “Geophysical Journal”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2Linked to original sources
USGS2009

Offshore double-planed shallow seismic zone in the NE Japan forearc region revealed by sP depth phases recorded by regional networks

We detected the sP depth phase at small epicentral distances of about 150 km or more in the seismograms of shallow earthquakes in the NE Japan forearc region. The focal depths of 1078 M > 3 earthquakes that occurred from 2000 to 2006 were precisely determined using the time delay of the sP phase from the initial P-wave arrival. The distribution of relocated hypocentres clearly shows the configuration of a double-planed shallow seismic zone beneath the Pacific Ocean. The upper plane has a low dip angle near the Japan Trench, increasing gradually to ???30?? at approximately 100 km landward of the Japan Trench. The lower plane is approximately parallel to the upper plane, and appears to be the near-trench counterpart of the lower plane of the double-planed deep seismic zone beneath the land area. The distance between the upper and lower planes is 28-32 km, which is approximately the same as or slightly smaller than that of the double-planed deep seismic zone beneath the land area. Focal mechanism solutions of the relocated earthquakes are determined from P-wave initial motion data. Although P-wave initial motion data for these offshore events are not ideally distributed on the focal sphere, we found that the upper-plane events that occur near the Japan Trench are characterized by normal faulting, whereas lower-plane events are characterized by thrust faulting. This focal mechanism distribution is the opposite to that of the double-planed deep seismic zone beneath the land area. The characteristics of these focal mechanisms for the shallow and deep doubled-planed seismic zones can be explained by a bending-unbending model of the subducting Pacific plate. Some of relocated earthquakes took place in the source area of the 1933 Mw8.4 Sanriku earthquake at depths of 10-23 km. The available focal mechanisms for these events are characterized by normal faulting. Given that the 1933 event was a large normal-fault event that occurred along a fault plane dipping landward, the earthquakes that currently occur just beneath or oceanwards of the Japan Trench are probably its aftershocks, suggesting that aftershock activity continues to the present day, 70 years after the main shock. ?? 2009 The Authors, Journal compilation ?? 2009 RAS.

Geophysical Journal International↗
USGS2008

Effect of 3-D viscoelastic structure on post-seismic relaxation from the 2004 M = 9.2 Sumatra earthquake

The 2004 M=9.2 Sumatra-Andaman earthquake profoundly altered the state of stress in a large volume surrounding the ???1400 km long rupture. Induced mantle flow fields and coupled surface deformation are sensitive to the 3-D rheology structure. To predict the post-seismic motions from this earthquake, relaxation of a 3-D spherical viscoelastic earth model is simulated using the theory of coupled normal modes. The quasi-static deformation basis set and solution on the 3-D model is constructed using: a spherically stratified viscoelastic earth model with a linear stress-strain relation; an aspherical perturbation in viscoelastic structure; a 'static'mode basis set consisting of Earth's spheroidal and toroidal free oscillations; a "viscoelastic" mode basis set; and interaction kernels that describe the coupling among viscoelastic and static modes. Application to the 2004 Sumatra-Andaman earthquake illustrates the profound modification of the post-seismic flow field at depth by a slab structure and similarly large effects on the near-field post-seismic deformation field at Earth's surface. Comparison with post-seismic GPS observations illustrates the extent to which viscoelastic relaxation contributes to the regional post-seismic deformation. ?? Journal compilation ?? 2008 RAS.

Geophysical Journal International↗
USGS2007

Lithology-derived structure classification from the joint interpretation of magnetotelluric and seismic models

Magnetotelluric and seismic methods provide complementary information about the resistivity and velocity structure of the subsurface on similar scales and resolutions. No global relation, however, exists between these parameters, and correlations are often valid for only a limited target area. Independently derived inverse models from these methods can be combined using a classification approach to map geologic structure. The method employed is based solely on the statistical correlation of physical properties in a joint parameter space and is independent of theoretical or empirical relations linking electrical and seismic parameters. Regions of high correlation (classes) between resistivity and velocity can in turn be mapped back and re-examined in depth section. The spatial distribution of these classes, and the boundaries between them, provide structural information not evident in the individual models. This method is applied to a 10 km long profile crossing the Dead Sea Transform in Jordan. Several prominent classes are identified with specific lithologies in accordance with local geology. An abrupt change in lithology across the fault, together with vertical uplift of the basement suggest the fault is sub-vertical within the upper crust. ?? 2007 The Authors Journal compilation ?? 2007 RAS.

Geophysical Journal International↗
USGS1970

Annotated bibliography on artificial recharge of ground water, 1955-67

Artificial ground-water recharge has become more important as water use by agriculture, industry, and municipalities increases. Water management agencies are increasingly interested in potential use of recharge for pollution abatement, waste-water disposal, and re-use and reclamation of locally available supplies. Research projects and theoretical analyses of operational recharge systems show increased scientific emphasis on the practice. Overall ground-water basin management systems generally now contain considerations of artificial recharge, whether by direct or indirect methods. Artificial ground-water recharge is a means of conserving surface runoff for future use in places where it would otherwise be lost, of protecting ground-water basins from salt-water encroachment along coastal areas, and of storing and distributing imported water. The biblio-graphy emphasizes technology; however, annotations of articles on waste-water reclamation, ground-water management and ground-water basin management are included. Subjects closely related to artificial recharge, including colloidal flow through porous media, field or laboratory instrumentation, and waste disposal by deep well injection are included where they specifically relate to potential recharge problems. Where almost the same material has been published in several journals, all references are included on the assumption that some publications may be more readily available to interested persons than others. Other publications, especially those of foreign literature, provided abstracts that were used freely as time limitations precluded obtaining and annotating all materials. Abstracts taken from published sources are noted. These are: "Abstracts of North American Geology," U.S. Department of the Interior, Geological Survey; "Abstracts of Recent Published Material on Foil and Water Conservation," ARS-41 series, Agricultural F.esearch Service, U.S. Department of Agriculture; "Water and1 Water Engineering," published by Fuel and Metallurgical Journals, Ltd., London, England; "Journal of Geophysical Research," American Geophysical Union, Washington, D.C.; "American Society of Civil Engineers Transactions," New York; "Selected Bibliography of Hydrology, United Kingdom, for the Years 1955-59," International Association of Scientific Hydrology; "Water Wells, an Annotated Bibliography," California University Water Resources Center Archives Report 13; "Re-use of Effluent in the Future With an Annotated Bibliography," by G. A. Whetstone, Texas Water Development Board Report 8, Austin, Tex.; "Journal of Water Pollution Control Federation," Washington, D.C.; and "A List of Selected Technical References on Artificial Recharge of Ground-Water Reservoirs," compiled by Roy W. Graves, Tulsa University, Information Services Department, Tulsa, Okla. Other notations are self-explanatory, and initials are those of the authors (DCS, DJG, WK). An unpublished compilation of recharge references by Arnon Arad sponsored by the United Nations Educational, Scientific, and Cultural Organization during a training period with the U.S. Geological Survey was also used. The bibliography is arranged alphabetically by author. Where an author has more than one publication, the arrangement is chronological; where an author has more than one publication in a given year, a, b, c, . . . are added. The indexing is by subject and geographic location. Each article was assigned the key words or phrases to best characterize its contents. Units of measure are as they were in the original article; abbreviations retained are generally those in common use such as mg/1 (milligrams per liter), ppm (parts per million), gpm (gallons per minute), km (kilometers), m (meters), cu m per hr (cubic meters p^r hour), cfs (cubic feet per second), me/1 (milliequivalents per liter), psi (pounds per square inch), BOD (biochemical oxygen demand), sq m (square meters), gpd (gallons per day), and mgd (million gallons per day). The bibliography was prepared because of the worldwide interest in the field of artificial recharge and the need for a single source of references to the literature published since 1954. The work is a sequel to the "Annotated Bibliography on Artificial Recharge of Ground Water Through 1954," by D. K. Todd, U.S. Geological Survey Water-Supply Paper 1477, published in 1959.

Water Supply Paper↗
USGS1980

Comment on 'A reinterpretation of the linear heat flow and heat production relationship for the exponential model of the heat production in the crust' by R.N. Singh & J.G. Negi.

In their recent paper, Singh & Negi, (This journal, 57, 741-744) contend that if thd slope of the empirical linear relation between heat flow and heat production is interpreted as the decay-length of an exponential depth-distribution of sources, a discrepancy rises, whereas if it is interpreted as the depth of a step distribution, it does not. I should like to point out that their discrepancy follows from their arbitrary assumption of one of a range of physical possibilities unconstrained by the observations; with an equally valid alternate assumption (Lachenbruch 1970) the discrepancy disappears. In any case such discrepancies are probably minor compared to physical difficulties that arise from the step model, and to uncertainties introduced by other assumptions in any simple model.-Author

Geophysical Journal of the Royal Astronomical Soci↗
USGS2008

Summit CO 2 emission rates by the CO 2 /SO 2 ratio method at Kīlauea Volcano, Hawaiʻi, during a period of sustained inflation

The emission rate of carbon dioxide escaping from the summit of Kīlauea Volcano, Hawaiʻi, proved highly variable, averaging 4900 ± 2000 metric tons per day (t/d) in June–July 2003 during a period of summit inflation. These results were obtained by combining over 90 measurements of COSPEC-derived SO 2 emission rates with synchronous CO 2 /SO 2 ratios of the volcanic gas plume along the summit COSPEC traverse. The results are lower than the CO 2 emission rate of 8500 ± 300 t/d measured by the same method in 1995–1999 during a period of long-term summit deflation [Gerlach, T.M., McGee, K.A., Elias, T., Sutton, A.J. and Doukas, M.P., 2002. Carbon dioxide emission rate of Kīlauea Volcano: Implications for primary magma and the summit reservoir. Journal of Geophysical Research-Solid Earth, 107(B9): art. no.-2189.]. Analysis of the data indicates that the emission rates of the present study likely reflect changes in the magma supply rate and residence time in the summit reservoir. It is also likely that emission rates during the inflation period were heavily influenced by SO 2 pulses emitted adjacent to the COSPEC traverse, which biased CO 2 /SO 2 ratios towards low values that may be unrepresentative of the global summit gas plume. We conclude that the SO 2 pulses are consequences of summit re-inflation under way since 2003 and that CO 2 emission rates remain comparable to, but more variable than, those measured prior to re-inflation.

Hawaii↗
USGS1986

Field observations of bed shear stress and sediment resuspension on continental shelves, Alaska and California

Bed shear stress was estimated using wave and current measurements obtained with the GEOPROBE bottom-tripod system during resuspension events in Norton Sound, Alaska, and on the northern California shelf. The boundary-layer model of Grant and Madsen (1979, Journal of Geophysical Research, 84, 1797-1808) was used to compute the bed shear stress under combined wave-generated and quasi-steady currents. Resuspension events were identified by sudden, large increases in light scattering at 1.9 m above the sea floor. The shear-stress values were used to compute the Shields parameter (??). The results for Norton Sound are in excellent agreement with the Shields threshold criterion; the data for the California shelf plot somewhat above the Shields threshold curve, though generally within the scatter envelope. Although the surface sediments in each area contain substantial fine-grained fractions (mean diameters were 0.007 cm in Norton Sound and 0.002 cm on the California shelf), the results do not indicate significant cohesion, because the sediment was entrained at bed shear-stress values close to those predicted by the modified Shields curve for cohesionless fine-grained particles. We suspect that frequent wave stirring and observed plowing of the surface sediment by benthonic animals maintain a high water content and contribute to the ease with which these materials are resuspended.

Alaska, California↗
USGS1990

Sediment movement along the U.S. east coast continental shelf-II. Modelling suspended sediment concentration and transport rate during storms

Long-term near-bottom wave and current observations and a one-dimensional sediment transport model are used to calculate the concentration and transport of sediment during winter storms at 60-80 m water depth along the southern flank of Georges Bank and in the Mid-Atlantic Bight. Calculations are presented for five stations, separated by more than 600 km alongshelf, that have different bottom sediment texture, bedforms and current conditions. A modified version of the sediment transport model presented by Grant and Glenn (1983, Technical Report to the American Gas Association), Glenn (1983, D.Sc. Thesis, M.I.T.), and Glenn and Grant (1987, Journal of Geophysical Research, 92, 8244-8264) is used to examine the influence of wave-current interaction, sediment stratification, and limitations on the erodibility of the bottom sediments on the concentration of sediment in the water column and on transport. Predicted suspended sediment concentrations are higher than observed, based on beam transmissometer measurements, unless an erosion limit of order a few millimeters for sediments finer than 94 ??m is imposed. The agreement between predicted and measured beam attenuation is better at stations that have significant amounts of silt plus clay in the surficial sediments than for stations with sandy sediments. Sediment concentrations during storms estimated by Moody et al. (1987, Continental Shelf Research, 7, 609-628) are within 50% of the model predictions. Sediment transport rates for sediments 94 ??m and finer are determined largely by the concentrations in the surficial sediment and the erosion depth limit. Large alongshelf transports in the direction of storm-driven currents are inferred for stations in the Mid-Atlantic Bight. During a 115-day period in winter 1979-1980, the net transport of sediment along the shelf was westward; benthic storms (defined as periods when the bottom wave stress exceeded the current stress by 2 dyn cm-2) occurred between 23 and 73% of the time, and greater than 91% of the net alongshelf transport was during storms. ?? 1990.

Georges Bank, Mid-Atlantic Bight↗
USGS2008

Time-averaged near-bed suspended sediment concentrations under waves and currents: Comparison of measured and model estimates

Profiles of suspended sediment concentration and velocity were measured over a 15-day period at a near-shore site off Santa Cruz, CA in Monterey Bay. The concentration and velocity data were collected with an Acoustic Backscattering System (ABS) and Acoustic Current Profiler (ACP) that were mounted on a bottom tripod. High-resolution bottom scanning sonar was also attached to the tripod to provide images of bed features during the experiment. Hourly time-averaged near-bed concentrations of suspended sediment were calculated from three models and compared with the measurements. Surface waves and currents that were generated by a storm of moderate intensity caused bed stresses that exceeded threshold stress for D 50 =0.02 cm, the median size of the moderately well-sorted bottom sediment, over a period of about 7 days. Estimates of the concentration at 1 cm above the bottom, C a1 , were obtained using the ABS measurements. These observations have been compared with predictions for the concentration at 1 cm above the bottom, C 1 . Nielsen's models for reference concentration C o [Nielsen, P., 1986. Suspended sediment concentrations under waves. Coastal Engineering 10, 32-31; Nielsen, P., 1992. Coastal Bottom Boundary Layers and Sediment Transport, Advanced Series on Ocean Engineering. World Scientific, Hackensack, NJ.] are purely wave-based and do not include effects of bottom currents on bed stress and bedform scales. C 1 calculated from this model compared well with measured C a1 when currents were weak and small oscillatory ripples were observed in the sonar images. However, during the 3-day period of highest bottom stresses modeled C 1 did not compare well to C a1 . The other two models for C 1 , Glenn and Grant [Glenn, S.M., Grant, W.D., 1987. A suspended sediment stratification correction for combined wave and current flows. Journal of Geophysical Research 92(C8), 8244-8264.] and van Rijn and Walstra [Van Rijn, L.C., Walstra, D.J.R., 2004. Description of TRANSPOR2004 and implementation in Delft3D-ONLINE. Interim Report prepared for DG Rijkswaterstaat, Rijksinstituut voor Kust en Zee. Delft Hydraulics Institute, The Netherlands.], accounted for combined wave-current stresses and included different formulations for predicting bedform scales. C 1 calculated from these models compared well with measurements throughout the early and middle period of the storm, but underpredicted measured values during the latter part of the experiment. An alternative method based on dimensional analysis provided a straightforward way to relate C 1 to bed stress and sediment parameters. This new relationship based on dimensional analysis explained about 77% of the variability in C a1 derived from the ABS data.

California↗
USGS2011

North American Jurassic APW: The current dilemma

Geologists are in a quandary over the correct interpretation of paleomagnetic data for Jurassic rocks of the North American plate. Conflicting reference paleopoles and alternate methods of constructing apparent polar wander (APW) paths have led to a controversy regarding the configuration of Jurassic APW for cratonic North America. These differences have been recently disputed at meetings (see Eos , Spring Meeting Supplement, April 7, 1992, p. 94) and in an exchange of letters in the Journal of Geophysical Research-Solid Earth [ Butler et al ., 1992; Van Fossen and Kent , 1992a]. At stake is important information concerning the Jurassic paleogeography of North America and the whole of Pangea, as well as the nature and driving mechanism of plate movements. In addition, cratonic paleopoles provide a reference frame for measuring relative displacements of tectonostratigraphic terranes.

Eos Science News↗
USGS2001

Three-dimenstional crustal velocity structure beneath the Strait of Georgia, British Columbia

The Strait of Georgia is a topographic depression straddling the boundary between the Insular and Coast belts in southwestern British Columbia. Two shallow earthquakes located within the strait ( M = 4.6 in 1997 and M = 5.0 in 1975) and felt throughout the Vancouver area illustrate the seismic potential of this region. As part of the 1998 Seismic Hazards Investigation of Puget Sound (SHIPS) experiment, seismic instruments were placed in and around the Strait of Georgia to record shots from a marine source within the strait. We apply a tomographic inversion procedure to first-arrival traveltime data to derive a minimum-structure 3-D P -wave velocity model for the upper crust to about 13 km depth. We also present a 2-D velocity model for a profile orientated across the Strait of Georgia derived using a minimum-parameter traveltime inversion approach. This paper represents the first detailed look at crustal velocity variations within the major Cretaceous to Cenozoic Georgia Basin, which underlies the Strait of Georgia. The 3-D velocity model clearly delineates the structure of the Georgia Basin. Taking the 6 km s −1 isovelocity contour to represent the top of the underlying basement, the basin thickens from between 2 and 4 km in the northwestern half of the strait to between 8 and 9 km at the southeastern end of the study region. Basin velocities in the northeastern half are 4.5–6 km s −1 and primarily represent the Upper Cretaceous Nanaimo Group. Velocities to the south are lower (3–6 km s −1 ) because of the additional presence of the overlying Tertiary Huntingdon Formation and more recent sediments, including glacial and modern Fraser River deposits. In contrast to the relatively smoothly varying velocity structure of the basin, velocities of the basement rocks, which comprise primarily Palaeozoic to Jurassic rocks of the Wrangellia Terrane and possibly Jurassic to mid-Cretaceous granitic rocks of the Coast Belt, show significantly more structure, probably an indication of the varying basement rock lithologies. The 2-D velocity model more clearly reveals the velocity layering associated with the recent sediments, Huntingdon Formation and Nanaimo Group of the southern Georgia Basin, as well as the underlying basement. We interpret lateral variation in sub-basin velocities of the 2-D model as a transition from Wrangellian to Coast Belt basement rocks. The effect of the narrow, onshore–offshore recording geometry of the seismic experiment on model resolution was tested to allow a critical assessment of the validity of the 3-D velocity model. Lateral resolution throughout the model to a depth of 3–5 km below the top of the basement is generally 10–20 km.

British Columbia↗
USGS2001

Seismic tomography shows that upwelling beneath Iceland is confined to the upper mantle

We report the results of the highest-resolution teleseismic tomography study yet performed of the upper mantle beneath Iceland. The experiment used data gathered by the Iceland Hotspot Project, which operated a 35-station network of continuously recording, digital, broad-band seismometers over all of Iceland 1996-1998. The structure of the upper mantle was determined using the ACH damped least-squares method and involved 42 stations, 3159 P-wave, and 1338 S-wave arrival times, including the phases P, pP, sP, PP, SP, PcP, PKIKP, pPKIKP, S, sS, SS, SKS and Sdiff. Artefacts, both perceptual and parametric, were minimized by well-tested smoothing techniques involving layer thinning and offset-and-averaging. Resolution is good beneath most of Iceland from ??? 60 km depth to a maximum of ??? 450 km depth and beneath the Tjornes Fracture Zone and near-shore parts of the Reykjanes ridge. The results reveal a coherent, negative wave-speed anomaly with a diameter of 200-250 km and anomalies in P-wave speed, Vp, as strong as -2.7 per cent and in S-wave speed, Vs, as strong as -4.9 per cent. The anomaly extends from the surface to the limit of good resolution at ??? 450 km depth. In the upper ??? 250 km it is centred beneath the eastern part of the Middle Volcanic Zone, coincident with the centre of the ??? 100 mGal Bouguer gravity low over Iceland, and a lower crustal low-velocity zone identified by receiver functions. This is probably the true centre of the Iceland hotspot. In the upper ??? 200 km, the low-wave-speed body extends along the Reykjanes ridge but is sharply truncated beneath the Tjornes Fracture Zone. This suggests that material may flow unimpeded along the Reykjanes ridge from beneath Iceland but is blocked beneath the Tjornes Fracture Zone. The magnitudes of the Vp, Vs and Vp/Vs anomalies cannot be explained by elevated temperature alone, but favour a model of maximum temperature anomalies <200 K, along with up to ??? 2 per cent of partial melt in the depth range ??? 100-300 km beneath east-central Iceland. The anomalous body is approximately cylindrical in the top 250 km but tabular in shape at greater depth, elongated north-south and generally underlying the spreading plate boundary. Such a morphological change and its relationship to surface rift zones are predicted to occur in convective upwellings driven by basal heating, passive upwelling in response to plate separation and lateral temperature gradients. Although we cannot resolve structure deeper than ??? 450 km, and do not detect a bottom to the anomaly, these models suggest that it extends no deeper than the mantle transition zone. Such models thus suggest a shallow origin for the Iceland hotspot rather than a deep mantle plume, and imply that the hotspot has been located on the spreading ridge in the centre of the north Atlantic for its entire history, and is not fixed relative to other Atlantic hotspots. The results are consistent with recent, regional full-thickness mantle tomography and whole-mantle tomography images that show a strong, low-wave-speed anomaly beneath the Iceland region that is confined to the upper mantle and thus do not require a plume in the lower mantle. Seismic and geochemical observations that are interpreted as indicating a lower mantle, or core-mantle boundary origin for the North Atlantic Igneous Province and the Iceland hotspot should be re-examined to consider whether they are consistent with upper mantle processes.

Geophysical Journal International↗
USGS1998

Three-dimensional seismic structure and moment tensors of non-double-couple earthquakes at the Hengill-Grensdalur volcanic complex, Iceland

The volcanic and geothermal areas of Iceland are rich sources of non-double-couple (non-DC) earthquakes. A state-of-the-art digital seismometer network deployed at the Hengill–Grensdalur volcanic complex in 1991 recorded 4000 small earthquakes. We used the best recorded of these to determine 3-D V P and V P /V S structure tomographically and accurate earthquake moment tensors. The V P field is dominated by high seismic wave speed bodies interpreted as solidified intrusions. A widespread negative (−4 per cent) V P /V S anomaly in the upper 4km correlates with the geothermal field, but is too strong to be caused solely by the effect of temperature upon liquid water or the presence of vapour, and requires in addition mineralogical or lithological differences between the geothermal reservoir and its surroundings. These may be caused by geothermal alteration. Well-constrained moment tensors were obtained for 70 of the best-recorded events by applying linear programming methods to P - and S -wave polarities and amplitude ratios. About 25 per cent of the mechanisms are, within observational error, consistent with DC mechanisms consistent with shear faulting. The other 75 per cent have significantly non-DC mechanisms. Many have substantial explosive components, one has a substantial implosive component, and the deviatoric component of many is strongly non-DC. Many of the non-DC mechanisms are consistent, within observational error, with simultaneous tensile and shear faulting. However, the mechanisms occupy a continuum in source-type parameter space and probably at least one additional source process is occurring. This may be fluid flow into newly formed cracks, causing partial compensation of the volumetric component. Studying non-shear earthquakes such as these has great potential for improving our understanding of geothermal processes and earthquake source processes in general.

Geophysical Journal International↗
USGS1999

The thin hot plume beneath Iceland

We present the results of a seismological investigation of the frequency-dependent amplitude variations across Iceland using data from the HOTSPOT array currently deployed there. The array is composed of 30 broad-band PASSCAL instruments. We use the parameter t(*), defined in the usual manner from spectral ratios (Halderman and Davis 1991), to compare observed S-wave amplitude variations with those predicted due to both anelastic attenuation and diffraction effects. Four teleseismic events at a range of azimuths are used to measure t(*). A 2-D vertical cylindrical plume model with a Gaussian-shaped velocity anomaly is used to model the variations. That part of t(*) caused by attenuation was estimated by tracing a ray through IASP91, then superimposing our plume model velocity anomaly and calculating the path integral of 1/vQ. That part of t(*) caused by diffraction was estimated using a 2-D finite difference code to generate synthetic seismograms. The same spectral ratio technique used for the data was then used to extract a predicted t(*). The t(*) variations caused by anelastic attenuation are unable to account for the variations we observe, but those caused by diffraction do. We calculate the t(*) variations caused by diffraction for different plume models and obtain our best-fit plume, which exhibits good agreement between the observed and measured t(*). The best-fit plume model has a maximum S-velocity anomaly of - 12 per cent and falls to 1/e of its maximum at 100 km from the plume centre. This is narrower than previous estimates from seismic tomography, which are broadened and damped by the methods of tomography. This velocity model would suggest greater ray theoretical traveltime delays than observed. However, we find that for such a plume, wave-front healing effects at frequencies of 0.03-0.175 Hz (the frequency range used to pick S-wave arrivals) causes a 40 per cent reduction in traveltime delay, reducing the ray theoretical delay to that observed.

Geophysical Journal International↗
USGS1999

The complex frequencies of long-period seismic events as probes of fluid composition beneath volcanoes

Long-period (LP) events have been widely observed in relation to magmatic and hydrothermal activities in volcanic areas. LP waveforms characterized by their harmonic signature have been interpreted as oscillations of a fluid-filled resonator, and mixtures of liquid and gas in the form of bubbly liquids have been mainly assumed for the fluid. To investigate the characteristic properties of the resonator system, we analyse waveforms of LP events observed at four different volcanoes in Hawaii, Alaska, Colombia and Japan using a newly developed spectral method. This method allows an estimation of the complex frequencies of decaying sinusoids based on an autoregressive model. The results of our analysis show a wide variety of Q factors, ranging from tens to several hundred. We compare these complex frequencies with those predicted by the fluid-filled crack model for various mixtures of liquid, gas and ash. Although the oscillations of LP events with Q smaller than 50 can be explained by various combinations of liquids and gases, we find that ash-laden gases are required to explain long-lasting oscillations with Q larger than 100. The complex frequencies of LP events yield useful information on the types of fluids. Temporal and spatial variations of the complex frequencies can be used as probes of fluid composition beneath volcanoes.

Geophysical Journal International↗
USGS2000

Constraints on a plume in the mid-mantle beneath the Iceland region from seismic array data

Teleseismic P waves passing through low-wave-speed bodies in the mantle are refracted, causing anomalies in their propagation directions that can be measured by seismometer arrays. Waves from earthquakes in the eastern Pacific and western North America arriving at the NORSAR array in Norway and at seismic stations in Scotland pass beneath the Iceland region at depths of ~ 1000-2000 km. Waves arriving at NORSAR have anomalous arrival azimuths consistent with a low-wave-speed body at a depth of ~ 1500 km beneath the Iceland-Faeroe ridge with a maximum diameter of ~ 250 km and a maximum wave-speed contrast of ~ 1.5 per cent. This agrees well with whole-mantle tomography results, which image a low-wave-speed body at this location with a diameter of ~ 500 km and a wave-speed anomaly of ~ 0.5 per cent, bearing in mind that whole-mantle tomography, because of its limited resolution, broadens and weakens small anomalies. The observations cannot resolve the location of the body, and the anomaly could be caused in whole or in part by larger bodies farther away, for example by a body imaged beneath Greenland by whole-mantle tomography.

Geophysical Journal International↗
USGS2001

P-wave velocity structure of the uppermost mantle beneath Hawaii from traveltime tomography

We examine the P-wave velocity structure beneath the island of Hawaii using P-wave residuals from teleseismic earthquakes recorded by the Hawaiian Volcano Observatory seismic network. The station geometry and distribution of events makes it possible to image the velocity structure between ~ 40 and 100 km depth with a lateral resolution of ~ 15 km and a vertical resolution of ~ 30 km. For depths between 40 and 80 km, P-wave velocities are up to 5 per cent slower in a broad elongated region trending SE-NW that underlies the island between the two lines defined by the volcanic loci. No direct correlation between the magnitude of the lithospheric anomaly and the current level of volcanic activity is apparent, but the slow region is broadened at ~ 19.8??N and narrow beneath Kilauea. In the case of the occanic lithosphere beneath Hawaii, slow seismic velocities are likely to be related to magma transport from the top of the melting zone at the base of the lithosphere to the surface. Thermal modelling shows that the broad elongated low-velocity zone cannot be explained in terms of conductive heating by one primary conduit per volcano but that more complicated melt pathways must exist.

Geophysical Journal International↗
USGS2001

Solving the dynamic rupture problem with different numerical approaches and constitutive laws

We study the dynamic initiation, propagation and arrest of a 2-D in-plane shear rupture by solving the elastodynamic equation by using both a boundary integral equation method and a finite difference approach. For both methods we adopt different constitutive laws: a slip-weakening (SW) law, with constant weakening rate, and rate- and state-dependent friction laws (Dieterich-Ruina). Our numerical procedures allow the use of heterogeneous distributions of constitutive parameters along the fault for both formulations. We first compare the two solution methods with an SW law, emphasizing the required stability conditions to achieve a good resolution of the cohesive zone and to avoid artificial complexity in the solutions. Our modelling results show that the two methods provide very similar time histories of dynamic source parameters. We point out that, if a careful control of resolution and stability is performed, the two methods yield identical solutions. We have also compared the rupture evolution resulting from an SW and a rate- and state-dependent friction law. This comparison shows that despite the different constitutive formulations, a similar behaviour is simulated during the rupture propagation and arrest. We also observe a crack tip bifurcation and a jump in rupture velocity (approaching the P-wave speed) with the Dieterich-Ruina (DR) law. The rupture arrest at a barrier (high strength zone) and the barrier-healing mechanism are also reproduced by this law. However, this constitutive formulation allows the simulation of a more general and complex variety of rupture behaviours. By assuming different heterogeneous distributions of the initial constitutive parameters, we are able to model a barrier-healing as well as a self-healing process. This result suggests that if the heterogeneity of the constitutive parameters is taken into account, the different healing mechanisms can be simulated. We also study the nucleation phase duration Tn, defined as the time necessary for the crack to reach the half-length Ic. We compare the Tn values resulting from distinct simulations calculated using different constitutive laws and different sets of constitutive parameters. Our results confirm that the DR law provides a different description of the nucleation process than the SW law adopted in this study. We emphasize that the DR law yields a complete description of the rupture process, which includes the most prominent features of SW.

Geophysical Journal International↗