Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Geophysics Journal International”

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
USGS1999

Frequency-magnitude statistics and spatial correlation dimensions of earthquakes at Long Valley caldera, California

Intense earthquake swarms at Long Valley caldera in late 1997 and early 1998 occurred on two contrasting structures. The first is defined by the intersection of a north-northwesterly array of faults with the southern margin of the resurgent dome, and is a zone of hydrothermal upwelling. Seismic activity there was characterized by high b-values and relatively low values of D, the spatial fractal dimension of hypocentres. The second structure is the pre-existing South Moat fault, which has generated large-magnitude seismic activity in the past. Seismicity on this structure was characterized by low b-values and relatively high D. These observations are consistent with low-magnitude, clustered earthquakes on the first structure, and higher-magnitude, diffuse earthquakes on the second structure. The first structure is probably an immature fault zone, fractured on a small scale and lacking a well-developed fault plane. The second zone represents a mature fault with an extensive, coherent fault plane.

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

Remarks on the travelling wave decomposition

In elastic wave propagation on a spherically symmetric earth model, a normal mode sum is converted into a sum of equivalent travelling waves by means of a travelling wave decomposition (TWD). For two decades, seismologists have assumed that each travelling wave in the TWD is associated with only real phase velocities, that is, no evanescent waves travel on a spherically symmetric earth model. In this paper, this assumption is proven false. By including a countably infinite set of waves travelling as evanescent waves, several conceptual difficulties confronting the TWD are resolved.

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

Tomographic inversion of P-wave velocity and Q structures beneath the Kirishima volcanic complex, Southern Japan, based on finite difference calculations of complex traveltimes

We estimate the P-wave velocity and attenuation structures beneath the Kirishima volcanic complex, southern Japan, by inverting the complex traveltimes (arrival times and pulse widths) of waveform data obtained during an active seismic experiment conducted in 1994. In this experiment, six 200-250 kg shots were recorded at 163 temporary seismic stations deployed on the volcanic complex. We use first-arrival times for the shots, which were hand-measured interactively. The waveform data are Fourier transformed into the frequency domain and analysed using a new method based on autoregressive modelling of complex decaying oscillations in the frequency domain to determine pulse widths for the first-arrival phases. A non-linear inversion method is used to invert 893 first-arrival times and 325 pulse widths to estimate the velocity and attenuation structures of the volcanic complex. Wavefronts for the inversion are calculated with a finite difference method based on the Eikonal equation, which is well suited to estimating the complex traveltimes for the structures of the Kirishima volcano complex, where large structural heterogeneities are expected. The attenuation structure is derived using ray paths derived from the velocity structure. We obtain 3-D velocity and attenuation structures down to 1.5 and 0.5 km below sea level, respectively. High-velocity pipe-like structures with correspondingly low attenuation are found under the summit craters. These pipe-like structures are interpreted as remnant conduits of solidified magma. No evidence of a shallow magma chamber is visible in the tomographic images.

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↗
USGS2002

Geodetic imaging: Reservoir monitoring using satellite interferometry

Fluid fluxes within subsurface reservoirs give rise to surface displacements, particularly over periods of a year or more. Observations of such deformation provide a powerful tool for mapping fluid migration within the Earth, providing new insights into reservoir dynamics. In this paper we use Interferometric Synthetic Aperture Radar (InSAR) range changes to infer subsurface fluid volume strain at the Coso geothermal field. Furthermore, we conduct a complete model assessment, using an iterative approach to compute model parameter resolution and covariance matrices. The method is a generalization of a Lanczos-based technique which allows us to include fairly general regularization, such as roughness penalties. We find that we can resolve quite detailed lateral variations in volume strain both within the reservoir depth range (0.4-2.5 km) and below the geothermal production zone (2.5-5.0 km). The fractional volume change in all three layers of the model exceeds the estimated model parameter uncertainly by a factor of two or more. In the reservoir depth interval (0.4-2.5 km), the predominant volume change is associated with northerly and westerly oriented faults and their intersections. However, below the geothermal production zone proper [the depth range 2.5-5.0 km], there is the suggestion that both north- and northeast-trending faults may act as conduits for fluid flow.

Geophysical Journal International↗
USGS2002

Seismic anisotropy and mantle creep in young orogens

Seismic anisotropy provides evidence for the physical state and tectonic evolution of the lithosphere. We discuss the origin of anisotropy at various depths, and relate it to tectonic stress, geotherms and rheology. The anisotropy of the uppermost mantle is controlled by the orthorhombic mineral olivine, and may result from ductile deformation, dynamic recrystallization or annealing. Anisotropy beneath young orogens has been measured for the seismic phase Pn that propagates in the uppermost mantle. This anisotropy is interpreted as being caused by deformation during the most recent thermotectonic event, and thus provides information on the process of mountain building. Whereas tectonic stress and many structural features in the upper crust are usually orientated perpendicular to the structural axis of mountain belts, Pn anisotropy is aligned parallel to the structural axis. We interpret this to indicate mountain-parallel ductile (i.e. creeping) deformation in the uppermost mantle that is a consequence of mountain-perpendicular compressive stresses. The preferred orientation of the fast axes of some anisotropic minerals, such as olivine, is known to be in the creep direction, a consequence of the anisotropy of strength and viscosity of orientated minerals. In order to explain the anisotropy of the mantle beneath young orogens we extend the concept of crustal 'escape' (or 'extrusion') tectonics to the uppermost mantle. We present rheological model calculations to support this hypothesis. Mountain-perpendicular horizontal stress (determined in the upper crust) and mountain-parallel seismic anisotropy (in the uppermost mantle) require a zone of ductile decoupling in the middle or lower crust of young mountain belts. Examples for stress and mountain-parallel Pn anisotropy are given for Tibet, the Alpine chains, and young mountain ranges in the Americas. Finally, we suggest a simple model for initiating mountain parallel creep.

Geophysical Journal International↗
USGS2002

Identifying elements of the plumbing system beneath Kilauea Volcano, Hawaii, from the source locations of very-long-period signals

We analyzed 16 seismic events recorded by the Hawaiian broad-band seismic network at Kilauca Volcano during the period September 9-26, 1999. Two distinct types of event are identified based on their spectral content, very-long-period (VLP) waveform, amplitude decay pattern and particle motion. We locate the VLP signals with a method based on analyses of semblance and particle motion. Different source regions are identified for the two event types. One source region is located at depths of ~1 km beneath the northeast edge of the Halemaumau pit crater. A second region is located at depths of ~8 km below the northwest quadrant of Kilauea caldera. Our study represents the first time that such deep sources have been identified in VLP data at Kilauea. This discovery opens the possibility of obtaining a detailed image of the location and geometry of the magma plumbing system beneath this volcano based on source locations and moment tensor inversions of VLP signals recorded by a permanent, large-aperture broad-band network.

Geophysical Journal International↗
USGS2002

Crustal structure beneath western and eastern Iceland from surface waves and receiver functions

We determine the crustal structures beneath 14 broad-band seismic stations, deployed in western, eastern, central and southern Iceland, using surface wave dispersion curves and receiver functions. We implement a method to invert receiver functions using constraints obtained from genetic algorithm inversion of surface waves. Our final models satisfy both data sets. The thickness of the upper crust, as defined by the velocity horizon Vs = 3.7 km s-1, is fairly uniform at ???6.5-9 km beneath the Tertiary intraplate areas of western and eastern Iceland, and unusually thick at 11 km beneath station HOT22 in the far south of Iceland. The depth to the base of the lower crust, as defined by the velocity horizon Vs = 4.1 km s-1 is ???20-26 km in western Iceland and ???27-33 km in eastern Iceland. These results agree with those of explosion profiles that detect a thinner crust beneath western Iceland than beneath eastern Iceland. An earlier report of a substantial low-velocity zone beneath the Middle Volcanic Zone in the lower crust is confirmed by a similar observation beneath an additional station there. As was found in previous receiver function studies, the most reliable feature of the results is the clear division into an upper sequence that is a few kilometres thick where velocity gradients are high, and a lower, thicker sequence where velocity gradients are low. The transition to typical mantle velocities is variable, and may range from being very gradational to being relatively sharp and clear. A clear Moho, by any definition, is rarely seen, and there is thus uncertainty in estimates of the thickness of the crust in many areas. Although a great deal of seismic data are now available constraining the structures of the crust and upper mantle beneath Iceland, their geological nature is not well understood.

Geophysical Journal International↗
USGS2002

Strike-slip earthquakes in the oceanic lithosphere: Observations of exceptionally high apparent stress

The radiated energies, E S , and seismic moments, M 0 , for 942 globally distributed earthquakes that occurred between 1987 to 1998 are examined to find the earthquakes with the highest apparent stresses (τ a =μ E S / M 0 , where μ is the modulus of rigidity). The globally averaged τ a for shallow earthquakes in all tectonic environments and seismic regions is 0.3 MPa. However, the subset of 49 earthquakes with the highest apparent stresses (τ a greater than about 5.0 MPa) is dominated almost exclusively by strike-slip earthquakes that occur in oceanic environments. These earthquakes are all located in the depth range 7–29 km in the upper mantle of the young oceanic lithosphere. Many of these events occur near plate-boundary triple junctions where there appear to be high rates of intraplate deformation. Indeed, the small rapidly deforming Gorda Plate accounts for 10 of the 49 high-τ a events. The depth distribution of τ a , which shows peak values somewhat greater than 25 MPa in the depth range 20–25 km, suggests that upper bounds on this parameter are a result of the strength of the oceanic lithosphere. A recently proposed envelope for apparent stress, derived by taking 6 per cent of the strength inferred from laboratory experiments for young (less than 30 Ma) deforming oceanic lithosphere, agrees well with the upper-bound envelope of apparent stresses over the depth range 5–30 km. The corresponding depth-dependent shear strength for young oceanic lithosphere attains a peak value of about 575 MPa at a depth of 21 km and then diminishes rapidly as the depth increases. In addition to their high apparent stresses, which suggest that the strength of the young oceanic lithosphere is highest in the depth range 10–30 km, our set of high-τ a earthquakes show other features that constrain the nature of the forces that cause interplate motion. First, our set of events is divided roughly equally between intraplate and transform faulting with similar depth distributions of τ a for the two types. Secondly, many of the intraplate events have focal mechanisms with the T -axes that are normal to the nearest ridge crest or subduction zone and P -axes that are normal to the proximate transform fault. These observations suggest that forces associated with the reorganization of plate boundaries play an important role in causing high-τ a earthquakes inside oceanic plates. Extant transform boundaries may be misaligned with current plate motion. To accommodate current plate motion, the pre-existing plate boundaries would have to be subjected to large horizontal transform push forces. A notable example of this is the triple junction near which the second large aftershock of the 1992 April Cape Mendocino, California, sequence occurred. Alternatively, subduction zone resistance may be enhanced by the collision of a buoyant lithosphere, a process that also markedly increases the horizontal stress. A notable example of this is the Aleutian Trench near which large events occurred in the Gulf of Alaska in late 1987 and the 1998 March Balleny Sea M = 8.2 earthquake within the Antarctic Plate.

Geophysical Journal International↗
USGS2003

Gaussian statistics for palaeomagnetic vectors

With the aim of treating the statistics of palaeomagnetic directions and intensities jointly and consistently, we represent the mean and the variance of palaeomagnetic vectors, at a particular site and of a particular polarity, by a probability density function in a Cartesian three-space of orthogonal magnetic-field components consisting of a single (unimoda) non-zero mean, spherically-symmetrical (isotropic) Gaussian function. For palaeomagnetic data of mixed polarities, we consider a bimodal distribution consisting of a pair of such symmetrical Gaussian functions, with equal, but opposite, means and equal variances. For both the Gaussian and bi-Gaussian distributions, and in the spherical three-space of intensity, inclination, and declination, we obtain analytical expressions for the marginal density functions, the cumulative distributions, and the expected values and variances for each spherical coordinate (including the angle with respect to the axis of symmetry of the distributions). The mathematical expressions for the intensity and off-axis angle are closed-form and especially manageable, with the intensity distribution being Rayleigh-Rician. In the limit of small relative vectorial dispersion, the Gaussian (bi-Gaussian) directional distribution approaches a Fisher (Bingham) distribution and the intensity distribution approaches a normal distribution. In the opposite limit of large relative vectorial dispersion, the directional distributions approach a spherically-uniform distribution and the intensity distribution approaches a Maxwell distribution. We quantify biases in estimating the properties of the vector field resulting from the use of simple arithmetic averages, such as estimates of the intensity or the inclination of the mean vector, or the variances of these quantities. With the statistical framework developed here and using the maximum-likelihood method, which gives unbiased estimates in the limit of large data numbers, we demonstrate how to formulate the inverse problem, and how to estimate the mean and variance of the magnetic vector field, even when the data consist of mixed combinations of directions and intensities. We examine palaeomagnetic secular-variation data from Hawaii and Re??union, and although these two sites are on almost opposite latitudes, we find significant differences in the mean vector and differences in the local vectorial variances, with the Hawaiian data being particularly anisotropic. These observations are inconsistent with a description of the mean field as being a simple geocentric axial dipole and with secular variation being statistically symmetrical with respect to reflection through the equatorial plane. Finally, our analysis of palaeomagnetic acquisition data from the 1960 Kilauea flow in Hawaii and the Holocene Xitle flow in Mexico, is consistent with the widely held suspicion that directional data are more accurate than intensity data.

Geophysical Journal International↗
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↗
USGS2003

Shallow-velocity models at the Kilauea Volcano, Hawaii, determined from array analyses of tremor wavefields

The properties of the surface wavefield at Kilauea Volcano are analysed using data from small-aperture arrays of short-period seismometers deployed in and around the Kilauea caldera. Tremor recordings were obtained during two Japan-US cooperative experiments conducted in 1996 and 1997. The seismometers were deployed in three semi-circular arrays with apertures of 300, 300 and 400 m, and a linear array with length of 1680 m. Data are analysed using a spatio-temporal correlation technique well suited for the study of the stationary stochastic wavefields of Rayleigh and Love waves associated with volcanic activity and scattering sources distributed in and around the summit caldera. Spatial autocorrelation coefficients are obtained as a function of frequency and are inverted for the dispersion characteristics of Rayleigh and Love waves using a grid search that seeks phase velocities for which the L-2 norm between data and forward modelling operators is minimized. Within the caldera, the phase velocities of Rayleigh waves range from 1400 to 1800 m s -1 at 1 Hz down to 300-400 m s -1 at 10 Hz, and the phase velocities of Love waves range from 2600 to 400 m s -1 within the same frequency band. Outside the caldera, Rayleigh wave velocities range from 1800 to 1600 m s -1 at 1 Hz down to 260-360 m s -1 at 10 Hz, and Love wave velocities range from 600 to 150 m s -1 within the same frequency band. The dispersion curves are inverted for velocity structure beneath each array, assuming these dispersions represent the fundamental modes of Rayleigh and Love waves. The velocity structures observed at different array sites are consistent with results from a recent 3-D traveltime tomography of the caldera region, and point to a marked velocity discontinuity associated with the southern caldera boundary.

Geophysical Journal International↗
USGS2003

A shallow fault-zone structure illuminated by trapped waves in the Karadere-Duzce branch of the North Anatolian Fault, western Turkey

We discuss the subsurface structure of the Karadere-Duzce branch of the North Anatolian Fault based on analysis of a large seismic data set recorded by a local PASSCAL network in the 6 months following the Mw = 7.4 1999 Izmit earthquake. Seismograms observed at stations located in the immediate vicinity of the rupture zone show motion amplification and long-period oscillations in both P- and S-wave trains that do not exist in nearby off-fault stations. Examination of thousands of waveforms reveals that these characteristics are commonly generated by events that are well outside the fault zone. The anomalous features in fault-zone seismograms produced by events not necessarily in the fault may be referred to generally as fault-zone-related site effects. The oscillatory shear wave trains after the direct S arrival in these seismograms are analysed as trapped waves propagating in a low-velocity fault-zone layer. The time difference between the S arrival and trapped waves group does not grow systematically with increasing source-receiver separation along the fault. These observations imply that the trapping of seismic energy in the Karadere-Duzce rupture zone is generated by a shallow fault-zone layer. Traveltime analysis and synthetic waveform modelling indicate that the depth of the trapping structure is approximately 3-4 km. The synthetic waveform modelling indicates further that the shallow trapping structure has effective waveguide properties consisting of thickness of the order of 100 m, a velocity decrease relative to the surrounding rock of approximately 50 per cent and an S-wave quality factor of 10-15. The results are supported by large 2-D and 3-D parameter space studies and are compatible with recent analyses of trapped waves in a number of other faults and rupture zones. The inferred shallow trapping structure is likely to be a common structural element of fault zones and may correspond to the top part of a flower-type structure. The motion amplification associated with fault-zone-related site effects increases the seismic shaking hazard near fault-zone structures. The effect may be significant since the volume of sources capable of generating motion amplification in shallow trapping structures is large.

Geophysical Journal International↗
USGS2003

The relationship between the instantaneous velocity field and the rate of moment release in the lithosphere

Instantaneous velocity gradients within the continental lithosphere are often related to the tectonic driving forces. This relationship is direct if the forces are secular, as for the case of loading of a locked section of a subduction interface by the downgoing plate. If the forces are static, as for the case of lateral variations in gravitational potential energy, then velocity gradients can be produced only if the lithosphere has, on average, zero strength. The static force model may be related to the long-term velocity field but not the instantaneous velocity field (typically measured geodetically over a period of several years) because over short time intervals the upper lithosphere behaves elastically. In order to describe both the short- and long-term behaviour of an (elastic) lithosphere-(viscoelastic) asthenosphere system in a self-consistent manner, I construct a deformation model termed the expected interseismic velocity (EIV) model. Assuming that the lithosphere is populated with faults that rupture continually, each with a definite mean recurrence time, and that the Earth is well approximated as a linear elastic-viscoelastic coupled system, I derive a simple relationship between the instantaneous velocity field and the average rate of moment release in the lithosphere. Examples with synthetic fault networks demonstrate that velocity gradients in actively deforming regions may to a large extent be the product of compounded viscoelastic relaxation from past earthquakes on hundreds of faults distributed over large ( ≥106 km2) areas.

Geophysical Journal International↗
USGS2003

A numerical simulation of magma motion, crustal deformation, and seismic radiation associated with volcanic eruptions

The finite difference method is used to calculate the magma dynamics, seismic radiation, and crustal deformation associated with a volcanic eruption. The model geometry consists of a cylindrical reservoir and narrow cylindrical conduit embedded in a homogeneous crust. We consider two models of eruption. In the first model, a lid caps the vent and the magma is overpressurized prior to the eruption. The eruption is triggered by the instantaneous removal of the lid, at which point the exit pressure becomes equal to the atmospheric pressure. In the second model, a plug at the reservoir outlet allows pressurization of only the magmatic fluid in the reservoir before the eruption. Magma transfer between the reservoir and conduit is triggered by the instantaneous removal of the plug, and the eruption occurs when the pressure at the conduit orifice exceeds the material strength of the lid capping the vent. In both models, magma dynamics are expressed by the equations of mass and momentum conservation in a compressible fluid, in which fluid expansion associated with depressurization is accounted for by a constitutive law relating pressure and density. Crustal motions are calculated from the equations of elastodynamics. The fluid and solid are dynamically coupled by applying the continuity of wall velocities and normal stresses across the conduit and reservoir boundaries. Free slip is allowed at the fluid-solid boundary. Both models predict the gradual depletion of the magma reservoir, which causes crustal deformation observed as a long-duration dilatational signal. Superimposed on this very-long-period (VLP) signal generated by mass transport are long-period (LP) oscillations of the magma reservoir and conduit excited by the acoustic resonance of the reservoir-conduit system during the eruption. The volume of the reservoir, vent size, and magma properties control the duration of VLP waves and dominant periods of LP oscillations. The second model predicts that when the magmatic fluid reaches the vent, a high-pressure pulse occurs at this location in accordance with the basic theory of compressible fluid dynamics. This abrupt pressure increase just beneath the vent is consistent with observed seismograms in which pulse-like Rayleigh waves excited by a shallow source are dominant. The strength of the lid plays an important role in the character of the seismograms and in defining the type of eruption observed.

Geophysical Journal International↗