Geology topics

A. Villasenor

Publications and source records attributed to A. Villasenor.

9 recordsLinked to original sources

USGS2011

Toward a comprehensive catalog of global historical seismicity

The United States Geological Survey (USGS) and the Colorado School of Mines (CSM) have initiated a project to locate more accurately all earthquakes recorded by instruments during the period 1900 to 1963. Seismicity for this period (hereafter referred to as historical seismicity, following Lee et al . [1988]) is still poorly understood, even for basic parameters such as earthquake locations (see Figure 1). In some cases this is the result of inherent limitations in the distribution, response characteristics, and timing of the instruments. However, locations for most of the pre-1964 earthquakes are poorly determined simply because modern data analysis techniques have yet to be applied to the available arrival-time observations, which are mainly preserved as printed bulletins and not in a computer-ready digital format. The arduous task of hand-entering these data has prevented the systematic analysis and relocation of historical seismicity.

Eos, Transactions, American Geophysical Union↗

USGS2009

A model for the termination of the Ryukyu subduction zone against Taiwan: A junction of collision, subduction/separation, and subduction boundaries

The NW moving Philippine Sea plate (PSP) collides with the Eurasian plate (EUP) in the vicinity of Taiwan, and at the same time, it subducts toward the north along SW Ryukyu. The Ryukyu subduction zone terminates against eastern Taiwan. While the Ryukyu Trench is a linear bathym??trie low about 100 km east of Taiwan, closer to Taiwan, it cannot be clearly identified bathymetrically owing to the deformation related to the collision, making the location of the intersection of the Ryukyu with Taiwan difficult to decipher. We propose a model for this complex of boundaries on the basis of seismicity and 3-D velocity structures. In this model the intersection is placed at the latitude of about 23.7??N, placing the northern part of the Coastal Range on EUP. As PSP gets deeper along the subduction zone it collides with EUP on the Taiwan side only where they are in direct contact. Thus, the Eurasian plate on the Taiwan side is being pushed and compressed by the NW moving Philippine Sea plate, at increasing depth toward the north. Offshore of northeastern Taiwan the wedge-shaped EUP on top of the Ryukyu subducting plate is connected to the EUP on the Ryukyu side and coupled to the NW moving PSP by friction at the plate interface. The two sides of the EUP above the western end of the subduction zone are not subjected to the same forces, and a difference in motions can be expected. The deformation of Taiwan as revealed by continuous GPS measurements, geodetic movement along the east coast of Taiwan, and the formation of the Hoping Basin can be understood in terms of the proposed model. Copyright 2009 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

USGS2003

Surface wave tomography of North America and the Caribbean using global and regional broad-band networks: Phase velocity maps and limitations of ray theory

We present phase velocity maps of fundamental mode Rayleigh waves across the North American and Caribbean plates. Our data set consists of 1846 waveforms from 172 events recorded at 91 broad-band stations operating in North America. We compute phase velocity maps in four narrow period bands between 50 and 150 s using a non-linear waveform inversion method that solves for phase velocity perturbations relative to a reference Earth model (PREM). Our results show a strong velocity contrast between high velocities beneath the stable North American craton, and lower velocities in the tectonically active western margin, in agreement with other regional and global surface wave tomography studies. We perform detailed comparisons with global model results, which display good agreement between phase velocity maps in the location and amplitude of the anomalies. However, forward modelling shows that regional maps are more accurate for predicting waveforms. In addition, at long periods, the amplitude of the velocity anomalies imaged in our regional phase velocity maps is three time larger than in global phase velocity models. This amplitude factor is necessary to explain the data accurately, showing that regional models provide a better image of velocity structures. Synthetic tests show that the raypath coverage used in this study enables one to resolve velocity features of the order of 800-1000 km. However, only larger length-scale features are observed in the phase velocity maps. The limitation in resolution of our maps can be attributed to the wave propagation theory used in the inversion. Ray theory does not account for off-great-circle ray propagation effects, such as ray bending or scattering. For wavelengths less than 1000 km, scattering effects are significant and may need to be considered.

Geophysical Journal International↗

USGS2000

Models of downdip frictional coupling for the Cascadia Megathrust

We have developed models of downdip frictional coupling along two transects across the Cascadia subduction zone in northern Washington and central Oregon. The models involve complicated downdip coupling profiles. Although not unique, our models closely predict available GPS shortening rates and vertical uplift data. We are able to explain relatively low horizontal shortening rates along the Washington coast and small vertical uplift rates in central Oregon. Our models depart from previous models by inclusion of a deeply coupled region assumed to be related to mafic rocks in both the upper and lower plates of the thrust.

Oregon, Washington↗

USGS1999

Three-dimensional velocity structure of the Kilauea caldera, Hawaii

High-resolution velocity models (0.5 km resolution) of the Kilauea caldera region are obtained by the tomographic inversion of both P - and S -wave arrival times. Data are from the permanent Hawaiian Volcano Observatory (HVO) seismic network, a broadband seismic network, and a temporary array of stations centered on the southern boundary of the caldera. A low-velocity P -wave anomaly is imaged centered on the southeastern edge of the caldera, with a velocity contrast of about 10% and a volume of 27 km³. The V P /V S model mimics the spatial extent of the P-wave anomaly, but is partitioned into two discrete anomalous volumes centered on the southern boundary of the caldera and on the upper east rift of the volcano. The corresponding Poisson's ratio in these zones is high (ν=0.25–0.32) which is consistent with a densely-cracked, hot volume which may contain partial melt. The large-scale features of the models are consistent with results obtained from an earlier, larger-scale (2 km resolution) tomographic image of Kilauea Volcano based on HVO network data.

Geophysical Research Letters↗

USGS1998

Seismic image of the Mount Spurr magmatic system

The three-dimenisonal P-wave velocity structure of Moutn Spurr is determined to depths of 10km by tomographic inversion of 3754 first-arriving P-wave times from local earthquakes recorded by a permanent network of 11 seismographs. Results show a prominent low-velocity zone extending from the surface to 3-4km below sea level beneath the southeastern flank of Crater Peak, spatially coincident with a geothermal system. P-wave velocities in this low-velocity zone are approximately 20% slower than those in the shallow crystalline basement rocks. Beneath Crater Peak an approximately 3km-wide zone of relative low velocities correlates with a near-vertical band of seismicity, suggestive of a magmatic conduit. No large low-velocity zone indicative of a magma chamber occurs within the upper 10km of the crust. These observations are consistent with petrologic and geochemical studies suggesting that Crater Peak magmas originate in the lower crust or upper mantle and have s short residence time in the shallow crust. Earthquakes relocaetd using the three-dimensional velocity structure correlate well with surface geology and other geophysical observations; thus, they provide additional constraints on the kinematics of the Mount Spurr magmatic system.

Bulletin of Volcanology↗

USGS1998

Three-dimensional P-wave velocity structure of Mt. Etna, Italy

The three-dimensional P -wave velocity structure of Mt. Etna is determined to depths of 15 km by tomographic inversion of first arrival times from local earthquakes recorded by a network of 29 permanent and temporary seismographs. Results show a near-vertical low-velocity zone that extends from beneath the central craters to a depth of 10 km. This low-velocity region is coincident with a band of steeply-dipping seismicity, suggesting a magmatic conduit that feeds the summit eruptions. The most prominent structure is an approximately 8-km-diameter high-velocity body located between 2 and 12 km depth below the southeast flank of the volcano. This high-velocity body is interpreted as a remnant mafic intrusion that is an important structural feature influencing both volcanism and east flank slope stability and faulting.

Geophysical Research Letters↗

USGS1998

Tectonic controls on magmatism in the Geysers-Clear Lake region: Evidence from new geophysical models

In order to study magmatism and geothermal systems in The Geysers–Clear Lake region, we developed a detailed three-dimensional tomographic velocity model based on local earthquakes. This high-resolution model resolves the velocity structure of the crust in the region to depths of approximately 12 km. The most significant velocity contrasts in The Geysers–Clear Lake region occur in the steam production area, where high velocities are associated with a Quaternary granitic pluton, and in the Mount Hannah region, where low velocities occur in a 5-km-thick section of Mesozoic argillites. In addition, a more regional tomographic model was developed using traveltimes from earthquakes covering most of northern California. This regional model sampled the whole crust, but at a lower resolution than the local model. The regional model outlines low velocities at depths of 8–12 km in The Geysers–Clear Lake area, which extend eastward to the Coast Range thrust. These low velocities are inferred to be related to unmetamorphosed Mesozoic sedimentary rocks. In addition, the regional velocity model indicates high velocities in the lower crust beneath the Clear Lake volcanic field, which we interpret to be associated with mafic underplating. No large silicic magma chamber is noted in either the local or regional tomographic models. A three-dimensional gravity model also has been developed in the area of the tomographic imaging. Our gravity model demonstrates that all density contrasts can be accounted for in the upper 5–7 km of the crust. Two-dimensional magnetotelluric models of data from a regional, east-west profile indicate high resistivities associated with the granitic pluton in The Geysers production area and low resistivities in the low-velocity section of Mesozoic argillites near Mount Hannah. No indication of midcrustal magma bodies is present in the magnetotelluric data. On the basis of heat flow and geologic evidence, Holocene intrusive activity is thought to have occurred near the Northwest Geysers, Mount Hannah, Sulphur Bank Mine, and perhaps other areas. The geophysical data provide no conclusive evidence for such activity, but the detailed velocity model is suggestive of intrusive activity near Mount Hannah similar to that in the “felsite” of The Geysers production area. The geophysical models, seismicity patterns, distribution of volcanic vents, heat flow, and other data indicate that small, young intrusive bodies that were injected along a northeast trend from The Geysers to Clear Lake probably control the thermal regime.

Geological Society of America Bulletin↗

USGS1998

Estimation of the intrinsic absorption and scattering attenuation in Northeastern Venezuela (Southeastern Caribbean) using coda waves

Northeastern Venezuela has been studied in terms of coda wave attenuation using seismograms from local earthquakes recorded by a temporary short-period seismic network. The studied area has been separated into two subregions in order to investigate lateral variations in the attenuation parameters. Coda-Q-1 (Q(c)-1) has been obtained using the single-scattering theory. The contribution of the intrinsic absorption (Q(i)-1) and scattering (Q(s)-1) to total attenuation (Q(t)-1) has been estimated by means of a multiple lapse time window method, based on the hypothesis of multiple isotropic scattering with uniform distribution of scatterers. Results show significant spatial variations of attenuation: the estimates for intermediate depth events and for shallow events present major differences. This fact may be related to different tectonic characteristics that may be due to the presence of the Lesser Antilles subduction zone, because the intermediate depth seismic zone may be coincident with the southern continuation of the subducting slab under the arc.

Pure and Applied Geophysics↗