Search USGS⌕ Search

SEARCH · Search USGS

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 217 records · Page 12Linked to original sources
USGS2007

A crustal seismic velocity model for the UK, Ireland and surrounding seas

A regional model of the 3-D variation in seismic P -wave velocity structure in the crust of NW Europe has been compiled from wide-angle reflection/refraction profiles. Along each 2-D profile a velocity–depth function has been digitised at 5 km intervals. These 1-D velocity functions were mapped into three dimensions using ordinary kriging with weights determined to minimise the difference between digitised and interpolated values. An analysis of variograms of the digitised data suggested a radial isotropic weighting scheme was most appropriate. Horizontal dimensions of the model cells are optimised at 40 × 40 km and the vertical dimension at 1 km. The resulting model provides a higher resolution image of the 3-D variation in seismic velocity structure of the UK, Ireland and surrounding areas than existing models. The construction of the model through kriging allows the uncertainty in the velocity structure to be assessed. This uncertainty indicates the high density of data required to confidently interpolate the crustal velocity structure, and shows that for this region the velocity is poorly constrained for large areas away from the input data.

Geophysical Journal International↗
USGS2008

Effects of topography and crustal heterogeneities on the source estimation of LP event at Kilauea volcano

The main goal of this study is to improve the modelling of the source mechanism associated with the generation of long period (LP) signals in volcanic areas. Our intent is to evaluate the effects that detailed structural features of the volcanic models play in the generation of LP signal and the consequent retrieval of LP source characteristics. In particular, effects associated with the presence of topography and crustal heterogeneities are here studied in detail. We focus our study on a LP event observed at Kilauea volcano, Hawaii, in 2001 May. A detailed analysis of this event and its source modelling is accompanied by a set of synthetic tests, which aim to evaluate the effects of topography and the presence of low velocity shallow layers in the source region. The forward problem of Green's function generation is solved numerically following a pseudo-spectral approach, assuming different 3-D models. The inversion is done in the frequency domain and the resulting source mechanism is represented by the sum of two time-dependent terms: a full moment tensor and a single force. Synthetic tests show how characteristic velocity structures, associated with shallow sources, may be partially responsible for the generation of the observed long-lasting ringing waveforms. When applying the inversion technique to Kilauea LP data set, inversions carried out for different crustal models led to very similar source geometries, indicating a subhorizontal cracks. On the other hand, the source time function and its duration are significantly different for different models. These results support the indication of a strong influence of crustal layering on the generation of the LP signal, while the assumption of homogeneous velocity model may bring to misleading results.

Geophysical Journal International↗
USGS2009

The postseismic response to the 2002 M 7.9 Denali Fault earthquake: Constraints from InSAR 2003-2005

InSAR is particularly sensitive to vertical displacements, which can be important in distinguishing between mechanisms responsible for the postseismic response to large earthquakes (afterslip, viscoelastic relaxation). We produce maps of the surface displacements resulting from the postseismic response to the 2002 Denali Fault earthquake, using data from the Canadian Radarsat-1 satellite from the periods summer 2003, summer 2004 and summer 2005. A peak-to-trough signal of amplitude 4 cm in the satellite line of sight was observed between summer 2003 and summer 2004. By the period between summer 2004 and summer 2005, the displacement rate had dropped below the threshold required for observation with InSAR over a single year. The InSAR observations show that the principal postseismic relaxation process acted at a depth of ∼50 km, equivalent to the top of the mantle. However, the observations are still incapable of distinguishing between distributed (viscoelastic relaxation) and localized (afterslip) deformation. The imposed coseismic stresses are highest in the lower crust and, assuming a Maxwell rheology, a viscosity ratio of at least 5 between lower crust and upper mantle is required to explain the contrast in behaviour. The lowest misfits are produced by mixed models of viscoelastic relaxation in the mantle and shallow afterslip in the upper crust. Profiles perpendicular to the fault show significant asymmetry, which is consistent with differences in rheological structure across the fault.

Geophysical Journal International↗
USGS2009

Developing framework to constrain the geometry of the seismic rupture plane on subduction interfaces a priori - A probabilistic approach

A key step in many earthquake source inversions requires knowledge of the geometry of the fault surface on which the earthquake occurred. Our knowledge of this surface is often uncertain, however, and as a result fault geometry misinterpretation can map into significant error in the final temporal and spatial slip patterns of these inversions. Relying solely on an initial hypocentre and CMT mechanism can be problematic when establishing rupture characteristics needed for rapid tsunami and ground shaking estimates. Here, we attempt to improve the quality of fast finite-fault inversion results by combining several independent and complementary data sets to more accurately constrain the geometry of the seismic rupture plane of subducting slabs. Unlike previous analyses aimed at defining the general form of the plate interface, we require mechanisms and locations of the seismicity considered in our inversions to be consistent with their occurrence on the plate interface, by limiting events to those with well-constrained depths and with CMT solutions indicative of shallow-dip thrust faulting. We construct probability density functions about each location based on formal assumptions of their depth uncertainty and use these constraints to solve for the ‘most-likely’ fault plane. Examples are shown for the trench in the source region of the M w 8.6 Southern Sumatra earthquake of March 2005, and for the Northern Chile Trench in the source region of the November 2007 Antofagasta earthquake. We also show examples using only the historic catalogues in regions without recent great earthquakes, such as the Japan and Kamchatka Trenches. In most cases, this method produces a fault plane that is more consistent with all of the data available than is the plane implied by the initial hypocentre and CMT mechanism. Using the aggregated data sets, we have developed an algorithm to rapidly determine more accurate initial fault plane geometries for source inversions of future earthquakes.

Geophysical Journal International↗
USGS2009

Validation of the rupture properties of the 2001 Kunlun, China (Ms = 8.1), earthquake from seismological and geological observations

We determine the finite-fault slip distribution of the 2001 Kunlun earthquake ( M s = 8.1) by inverting teleseismic waveforms, as constrained by geological and remote sensing field observations. The spatial slip distribution along the 400-km-long fault was divided into five segments in accordance with geological observations. Forward modelling of regional surface waves was performed to estimate the variation of the speed of rupture propagation during faulting. For our modelling, the regional 1-D velocity structure was carefully constructed for each of six regional seismic stations using three events with magnitudes of 5.1–5.4 distributed along the ruptured portion of the Kunlun fault. Our result shows that the average rupture velocity is about 3.6 km s −1 , consistent with teleseismic long period wave modelling. The initial rupture was almost purely strike-slip with a rupture velocity of 1.9 km s −1 , increasing to 3.5 km s −1 in the second fault segment, and reaching a rupture velocity of about 6 km s −1 in the third segment and the fourth segment, where the maximum surface offset, with a broad fault zone, was observed. The rupture velocity decelerated to a value of 3.3 km s −1 in the fifth and final segment. Coseismic slip on the fault was concentrated between the surface and a depth of about 10 km. We infer that significant variations in rupture velocity and the observed fault segmentation are indicative of variations in strength along the interface of the Kunlun fault, as well as variations in fault geometry.

Qinghai Province↗
USGS2009

Shallow seismic structure of Kunlun fault zone in northern Tibetan Plateau, China: Implications for the 2001 M s8.1 Kunlun earthquake

The shallow seismic velocity structure of the Kunlun fault zone (KLFZ) was jointly deduced from seismic refraction profiling and the records of trapped waves that were excited by five explosions. The data were collected after the 2001 Kunlun M s8.1 earthquake in the northern Tibetan Plateau. Seismic phases for the in-line record sections (26 records up to a distance of 15 km) along the fault zone were analysed, and 1-D P - and S -wave velocity models of shallow crust within the fault zone were determined by using the seismic refraction method. Sixteen seismic stations were deployed along the off-line profile perpendicular to the fault zone. Fault-zone trapped waves appear clearly on the record sections, which were simulated with a 3-D finite difference algorithm. Quantitative analysis of the correlation coefficients of the synthetic and observed trapped waveforms indicates that the Kunlun fault-zone width is 300 m, and S -wave quality factor Q within the fault zone is 15. Significantly, S -wave velocities within the fault zone are reduced by 30–45 per cent from surrounding rocks to a depth of at least 1–2 km, while P -wave velocities are reduced by 7–20 per cent. A fault-zone with such P - and S -low velocities is an indication of high fluid pressure because Vs is affected more than Vp . The low-velocity and low- Q zone in the KLFZ model is the effect of multiple ruptures along the fault trace of the 2001 M s8.1 Kunlun earthquake.

Tibetan Plateau↗
USGS2009

A grid-doubling finite-element technique for calculating dynamic three-dimensional spontaneous rupture on an earthquake fault

We present a new finite-element technique for calculating dynamic 3-D spontaneous rupture on an earthquake fault, which can reduce the required computational resources by a factor of six or more, without loss of accuracy. The grid-doubling technique employs small cells in a thin layer surrounding the fault. The remainder of the modelling volume is filled with larger cells, typically two or four times as large as the small cells. In the resulting non-conforming mesh, an interpolation method is used to join the thin layer of smaller cells to the volume of larger cells. Grid-doubling is effective because spontaneous rupture calculations typically require higher spatial resolution on and near the fault than elsewhere in the model volume. The technique can be applied to non-planar faults by morphing, or smoothly distorting, the entire mesh to produce the desired 3-D fault geometry. Using our FaultMod finite-element software, we have tested grid-doubling with both slip-weakening and rate-and-state friction laws, by running the SCEC/USGS 3-D dynamic rupture benchmark problems. We have also applied it to a model of the Hayward fault, Northern California, which uses realistic fault geometry and rock properties. FaultMod implements fault slip using common nodes, which represent motion common to both sides of the fault, and differential nodes, which represent motion of one side of the fault relative to the other side. We describe how to modify the traction-at-split-nodes method to work with common and differential nodes, using an implicit time stepping algorithm.

Geophysical Journal International↗
USGS2010

Rayleigh-wave phase-velocity maps and three-dimensional shear velocity structure of the western US from local non-plane surface wave tomography

We utilize two-and-three-quarter years of vertical-component recordings made by the Transportable Array (TA) component of Earthscope to constrain three-dimensional (3-D) seismic shear wave velocity structure in the upper 200 km of the western United States. Single-taper spectral estimation is used to compile measurements of complex spectral amplitudes from 44 317 seismograms generated by 123 teleseismic events. In the first step employed to determine the Rayleigh-wave phase-velocity structure, we implement a new tomographic method, which is simpler and more robust than scattering-based methods (e.g. multi-plane surface wave tomography). The TA is effectively implemented as a large number of local arrays by defining a horizontal Gaussian smoothing distance that weights observations near a given target point. The complex spectral-amplitude measurements are interpreted with the spherical Helmholtz equation using local observations about a succession of target points, resulting in Rayleigh-wave phase-velocity maps at periods over the range of 18–125 s. The derived maps depend on the form of local fits to the Helmholtz equation, which generally involve the nonplane-wave solutions of Friederich et al. In a second step, the phase-velocity maps are used to derive 3-D shear velocity structure. The 3-D velocity images confirm details witnessed in prior body-wave and surface-wave studies and reveal new structures, including a deep (>100 km deep) high-velocity lineament, of width ∼200 km, stretching from the southern Great Valley to northern Utah that may be a relic of plate subduction or, alternatively, either a remnant of the Mojave Precambrian Province or a mantle downwelling. Mantle seismic velocity is highly correlated with heat flow, Holocene volcanism, elastic plate thickness and seismicity. This suggests that shallow mantle structure provides the heat source for associated magmatism, as well as thinning of the thermal lithosphere, leading to relatively high stress concentration. Our images also confirm the presence of high-velocity mantle at 100 km depth beneath areas of suspected mantle delamination (southern Sierra Nevada; Grande Ronde uplift), low velocity mantle underlying active rift zones, and high velocity mantle associated with the subducting Juan de Fuca plate. Structure established during the Proterozoic appears to exert a lasting influence on subsequent volcanism and tectonism up to the Present.

Geophysical Journal International↗
USGS2010

Viscoelastic-cycle model of interseismic deformation in the northwestern United States

We apply a viscoelastic cycle model to a compilation of GPS velocity fields in order to address the kinematics of deformation in the northwestern United States. A viscoelastic cycle model accounts for time-dependent deformation following large crustal earthquakes and is an alternative to block models for explaining the interseismic crustal velocity field. Building on the approach taken in Pollitz et al., we construct a deformation model for the entire western United States-based on combined fault slip and distributed deformation-and focus on the implications for the Mendocino triple junction (MTJ), Cascadia megathrust, and western Washington. We find significant partitioning between strike-slip and dip-slip motion near the MTJ as the tectonic environment shifts from northwest-directed shear along the San Andreas fault system to east-west convergence along the Juan de Fuca Plate. By better accounting for the budget of aseismic and seismic slip along the Cascadia subduction interface in conjunction with an assumed rheology, we revise a previous model of slip for the M~ 9 1700 Cascadia earthquake. In western Washington, we infer slip rates on a number of strike-slip and dip-slip faults that accommodate northward convergence of the Oregon Coast block and northwestward convergence of the Juan de Fuca Plate. Lateral variations in first order mechanical properties (e.g. mantle viscosity, vertically averaged rigidity) explain, to a large extent, crustal strain that cannot be rationalized with cyclic deformation on a laterally homogeneous viscoelastic structure. Our analysis also shows that present crustal deformation measurements, particularly with the addition of the Plate Boundary Observatory, can constrain such lateral variations.

Geophysical Journal International↗
USGS2010

Locations and magnitudes of historical earthquakes in the Sierra of Ecuador (1587–1996)

The whole territory of Ecuador is exposed to seismic hazard. Great earthquakes can occur in the subduction zone (e.g. Esmeraldas, 1906, M w 8.8), whereas lower magnitude but shallower and potentially more destructive earthquakes can occur in the highlands. This study focuses on the historical crustal earthquakes of the Andean Cordillera. Several large cities are located in the Interandean Valley, among them Quito, the capital (∼2.5 millions inhabitants). A total population of ∼6 millions inhabitants currently live in the highlands, raising the seismic risk. At present, precise instrumental data for the Ecuadorian territory is not available for periods earlier than 1990 (beginning date of the revised instrumental Ecuadorian seismic catalogue); therefore historical data are of utmost importance for assessing seismic hazard. In this study, the Bakun & Wentworth method is applied in order to determine magnitudes, locations, and associated uncertainties for historical earthquakes of the Sierra over the period 1587–1976. An intensity-magnitude equation is derived from the four most reliable instrumental earthquakes ( M w between 5.3 and 7.1). Intensity data available per historical earthquake vary between 10 (Quito, 1587, Intensity ≥VI) and 117 (Riobamba, 1797, Intensity ≥III). The bootstrap resampling technique is coupled to the B&W method for deriving geographical confidence contours for the intensity centre depending on the data set of each earthquake, as well as confidence intervals for the magnitude. The extension of the area delineating the intensity centre location at the 67 per cent confidence level (±1σ) depends on the amount of intensity data, on their internal coherence, on the number of intensity degrees available, and on their spatial distribution. Special attention is dedicated to the few earthquakes described by intensities reaching IX, X and XI degrees. Twenty-five events are studied, and nineteen new epicentral locations are obtained, yielding equivalent moment magnitudes between 5.0 and 7.6. Large earthquakes seem to be related to strike slip faults between the North Andean Block and stable South America to the east, while moderate earthquakes ( M w ≤ 6) seem to be associated with to thrust faults located on the western internal slopes of the Interandean Valley.

Geophysical Journal International↗
USGS2010

Historical seismograms for unravelling a mysterious earthquake: The 1907 Sumatra Earthquake

History of instrumental seismology is short. Seismograms are available only for a little more than 100 years; high-quality seismograms are available only for the last 50 years and the seismological database is very limited in time. To extend the database, seismograms of old events are of vital importance. Many unusual earthquakes are known to have occurred, but their seismological characteristics are poorly known. The 1907 Sumatra earthquake is one of them (1907 January 4, M = 7.6). Gutenberg and Richter located this event in the outer-rise area of the Sunda arc. This earthquake is known to be anomalous because of its extensive tsunami, which is disproportionate of its magnitude. The tsunami affected the coastal areas over 950 km along the Sumatran coast. We investigated this earthquake using the historical seismograms we could collect from several seismological observatories. We examined the P -wave arrival times listed in the Strassburg Bulletin (1912) and other station bulletins. The scatter of the Observed−Computed traveltime residuals ranges from –30 to 30 s, too large to locate the event accurately. The uncertainty of the epicentre estimated from an S-P grid-search relocation study is at least 1° (∼110 km). We interpreted the Omori seismograms from Osaka, Mizusawa and Tokyo, and the Wiechert seismograms from Göttingen and Uppsala by comparing them with the seismograms simulated from modern broad-band seismograms of the 2002, 2008 and two 2010 Sumatra earthquakes which occurred near the 1907 earthquake. From the amplitude of Rayleigh waves recorded on the Omori seismograms we conclude that the magnitude of the 1907 earthquake at about 30 to 40 s is about 7.8 (i.e. 7.5 to 8.0). The SH waveforms recorded on the Göttingen and Uppsala seismograms suggest that the 1907 earthquake is a thrust earthquake at a shallow depth around 30 km. The most likely scenario is that the 1907 earthquake initiated on the subduction interface, and slowly ruptured up-dip into the shallow sediments and caused the extensive tsunami. Although their quantity and quality are limited, historical seismograms provide key quantitative information about old events that cannot be obtained otherwise. This underscores the importance of preserving historical seismograms.

Geophysical Journal International↗
USGS2011

Slip rate and slip magnitudes of past earthquakes along the Bogd left-lateral strike-slip fault (Mongolia)

We carried out morphotectonic studies along the left-lateral strike-slip Bogd Fault, the principal structure involved in the Gobi-Altay earthquake of 1957 December 4 (published magnitudes range from 7.8 to 8.3). The Bogd Fault is 260 km long and can be subdivided into five main geometric segments, based on variation in strike direction. West to East these segments are, respectively: the West Ih Bogd (WIB), The North Ih Bogd (NIB), the West Ih Bogd (WIB), the West Baga Bogd (WBB) and the East Baga Bogd (EBB) segments. Morphological analysis of offset streams, ridges and alluvial fans—particularly well preserved in the arid environment of the Gobi region—allows evaluation of late Quaternary slip rates along the different faults segments. In this paper, we measure slip rates over the past 200 ka at four sites distributed across the three western segments of the Bogd Fault. Our results show that the left-lateral slip rate is ∼1 mm yr –1 along the WIB and EIB segments and ∼0.5 mm yr –1 along the NIB segment. These variations are consistent with the restraining bend geometry of the Bogd Fault. Our study also provides additional estimates of the horizontal offset associated with the 1957 earthquake along the western part of the Bogd rupture, complementing previously published studies. We show that the mean horizontal offset associated with the 1957 earthquake decreases progressively from 5.2 m in the west to 2.0 m in the east, reflecting the progressive change of kinematic style from pure left-lateral strike-slip faulting to left-lateral-reverse faulting. Along the three western segments, we measure cumulative displacements that are multiples of the 1957 coseismic offset, which may be consistent with a characteristic slip. Moreover, using these data, we re-estimate the moment magnitude of the Gobi-Altay earthquake at M w 7.78–7.95. Combining our slip rate estimates and the slip distribution per event we also determined a mean recurrence interval of ∼2500–5200 yr for past earthquakes along the different segments of the western Bogd Fault. This suggests that the three western segments of the Bogd Fault and the Gurvan Bulag thrust fault (a reverse fault bounding the southern side of the Ih Bogd range that ruptured during the 1957 earthquake) have similar average recurrence times, and therefore may have ruptured together in previous earthquakes as they did in 1957. These results suggest that the western part of the Bogd Fault system, including the Gurvan Bulag thrust fault, usually behaves in a ‘characteristic earthquake’ mode.

Bogd Fault↗
USGS2011

Evidence from lava flows for complex polarity transitions: The new composite Steens Mountain reversal record

Geomagnetic polarity transitions may be significantly more complex than are currently depicted in many sedimentary and lava-flow records. By splicing together paleomagnetic results from earlier studies at Steens Mountain with those from three newly studied sections of Oregon Plateau flood basalts at Catlow Peak and Poker Jim Ridge 70–90 km to the southeast and west, respectively, we provide support for this interpretation with the most detailed account of a magnetic field reversal yet observed in volcanic rocks. Forty-five new distinguishable transitional (T) directions together with 30 earlier ones reveal a much more complex and detailed record of the 16.7 Ma reversed (R)-to-normal (N) polarity transition that marks the end of Chron C5Cr. Compared to the earlier R-T-N-T-N reversal record, the new record can be described as R-T-N-T-N-T-R-T-N. The composite record confirms earlier features, adds new west and up directions and an entire large N-T-R-T segment to the path, and fills in directions on the path between earlier directional jumps. Persistent virtual geomagnetic pole (VGP) clusters and separate VGPs have a preference for previously described longitudinal bands from transition study compilations, which suggests the presence of features at the core–mantle boundary that influence the flow of core fluid and distribution of magnetic flux. Overall the record is consistent with the generalization that VGP paths vary greatly from reversal to reversal and depend on the location of the observer. Rates of secular variation confirm that the flows comprising these sections were erupted rapidly, with maximum rates estimated to be 85–120 m ka −1 at Catlow and 130–195 m ka −1 at Poker Jim South. Paleomagnetic poles from other studies are combined with 32 non-transitional poles found here to give a clockwise rotation of the Oregon Plateau of 11.4°± 5.6° with respect to the younger Columbia River Basalt Group flows to the north and 14.5°± 4.6° with respect to cratonic North America (95 per cent confidence interval).

Geophysical Journal International↗
USGS2011

A distribution-based parameterization for improved tomographic imaging of solute plumes

Difference geophysical tomography (e.g. radar, resistivity and seismic) is used increasingly for imaging fluid flow and mass transport associated with natural and engineered hydrologic phenomena, including tracer experiments, in situ remediation and aquifer storage and recovery. Tomographic data are collected over time, inverted and differenced against a background image to produce ‘snapshots’ revealing changes to the system; these snapshots readily provide qualitative information on the location and morphology of plumes of injected tracer, remedial amendment or stored water. In principle, geometric moments (i.e. total mass, centres of mass, spread, etc.) calculated from difference tomograms can provide further quantitative insight into the rates of advection, dispersion and mass transfer; however, recent work has shown that moments calculated from tomograms are commonly biased, as they are strongly affected by the subjective choice of regularization criteria. Conventional approaches to regularization (Tikhonov) and parametrization (image pixels) result in tomograms which are subject to artefacts such as smearing or pixel estimates taking on the sign opposite to that expected for the plume under study. Here, we demonstrate a novel parametrization for imaging plumes associated with hydrologic phenomena. Capitalizing on the mathematical analogy between moment-based descriptors of plumes and the moment-based parameters of probability distributions, we design an inverse problem that (1) is overdetermined and computationally efficient because the image is described by only a few parameters, (2) produces tomograms consistent with expected plume behaviour (e.g. changes of one sign relative to the background image), (3) yields parameter estimates that are readily interpreted for plume morphology and offer direct insight into hydrologic processes and (4) requires comparatively few data to achieve reasonable model estimates. We demonstrate the approach in a series of numerical examples based on straight-ray difference-attenuation radar monitoring of the transport of an ionic tracer, and show that the methodology outlined here is particularly effective when limited data are available.

Geophysical Journal International↗
USGS2011

Pressure waves in a supersaturated bubbly magma

We study the interaction of acoustic pressure waves with an expanding bubbly magma. The expansion of magma is the result of bubble growth during or following magma decompression and leads to two competing processes that affect pressure waves. On the one hand, growth in vesicularity leads to increased damping and decreased wave amplitudes, and on the other hand, a decrease in the effective bulk modulus of the bubbly mixture reduces wave velocity, which in turn, reduces damping and may lead to wave amplification. The additional acoustic energy originates from the chemical energy released during bubble growth. We examine this phenomenon analytically to identify conditions under which amplification of pressure waves is possible. These conditions are further examined numerically to shed light on the frequency and phase dependencies in relation to the interaction of waves and growing bubbles. Amplification is possible at low frequencies and when the growth rate of bubbles reaches an optimum value for which the wave velocity decreases sufficiently to overcome the increased damping of the vesicular material. We examine two amplification phase-dependent effects: (1) a tensile-phase effect in which the inserted wave adds to the process of bubble growth, utilizing the energy associated with the gas overpressure in the bubble and therefore converting a large proportion of this energy into additional acoustic energy, and (2) a compressive-phase effect in which the pressure wave works against the growing bubbles and a large amount of its acoustic energy is dissipated during the first cycle, but later enough energy is gained to amplify the second cycle. These two effects provide additional new possible mechanisms for the amplification phase seen in Long-Period (LP) and Very-Long-Period (VLP) seismic signals originating in magma-filled cracks.

Geophysical Journal International↗
USGS2011

Time-lapse three-dimensional inversion of complex conductivity data using an active time constrained (ATC) approach

Induced polarization (more precisely the magnitude and phase of impedance of the subsurface) is measured using a network of electrodes located at the ground surface or in boreholes. This method yields important information related to the distribution of permeability and contaminants in the shallow subsurface. We propose a new time-lapse 3-D modelling and inversion algorithm to image the evolution of complex conductivity over time. We discretize the subsurface using hexahedron cells. Each cell is assigned a complex resistivity or conductivity value. Using the finite-element approach, we model the in-phase and out-of-phase (quadrature) electrical potentials on the 3-D grid, which are then transformed into apparent complex resistivity. Inhomogeneous Dirichlet boundary conditions are used at the boundary of the domain. The calculation of the Jacobian matrix is based on the principles of reciprocity. The goal of time-lapse inversion is to determine the change in the complex resistivity of each cell of the spatial grid as a function of time. Each model along the time axis is called a ‘reference space model’. This approach can be simplified into an inverse problem looking for the optimum of several reference space models using the approximation that the material properties vary linearly in time between two subsequent reference models. Regularizations in both space domain and time domain reduce inversion artefacts and improve the stability of the inversion problem. In addition, the use of the time-lapse equations allows the simultaneous inversion of data obtained at different times in just one inversion step (4-D inversion). The advantages of this new inversion algorithm are demonstrated on synthetic time-lapse data resulting from the simulation of a salt tracer test in a heterogeneous random material described by an anisotropic semi-variogram.

Geophysical Journal International↗
USGS2011

Assessing historical rate changes in global tsunami occurrence

The global catalogue of tsunami events is examined to determine if transient variations in tsunami rates are consistent with a Poisson process commonly assumed for tsunami hazard assessments. The primary data analyzed are tsunamis with maximum sizes >1m. The record of these tsunamis appears to be complete since approximately 1890. A secondary data set of tsunamis >0.1m is also analyzed that appears to be complete since approximately 1960. Various kernel density estimates used to determine the rate distribution with time indicate a prominent rate change in global tsunamis during the mid-1990s. Less prominent rate changes occur in the early- and mid-20th century. To determine whether these rate fluctuations are anomalous, the distribution of annual event numbers for the tsunami catalogue is compared to Poisson and negative binomial distributions, the latter of which includes the effects of temporal clustering. Compared to a Poisson distribution, the negative binomial distribution model provides a consistent fit to tsunami event numbers for the >1m data set, but the Poisson null hypothesis cannot be falsified for the shorter duration >0.1m data set. Temporal clustering of tsunami sources is also indicated by the distribution of interevent times for both data sets. Tsunami event clusters consist only of two to four events, in contrast to protracted sequences of earthquakes that make up foreshock-main shock-aftershock sequences. From past studies of seismicity, it is likely that there is a physical triggering mechanism responsible for events within the tsunami source 'mini-clusters'. In conclusion, prominent transient rate increases in the occurrence of global tsunamis appear to be caused by temporal grouping of geographically distinct mini-clusters, in addition to the random preferential location of global M >7 earthquakes along offshore fault zones.

Geophysical Journal International↗
USGS2011

Analysis of dispersion and attenuation of surface waves in poroelastic media in the exploration-seismic frequency band

We analyse dispersion and attenuation of surface waves at free surfaces of possible vacuum/poroelastic media: permeable-‘open pore’, impermeable-‘closed pore’ and partially permeable boundaries, which have not been previously reported in detail by researchers, under different surface-permeable, viscous-damping, elastic and fluid-flowing conditions. Our discussion is focused on their characteristics in the exploration-seismic frequency band (a few through 200 Hz) for near-surface applications. We find two surface-wave modes exist, R 1 waves for all conditions, and R 2 waves for closed-pore and partially permeable conditions. For R 1 waves, velocities disperse most under partially permeable conditions and least under the open-pore condition. High-coupling damping coefficients move the main dispersion frequency range to high frequencies. There is an f 1 frequency dependence as a constant- Q model for attenuation at high frequencies. R 1 waves for the open pore are most sensitive to elastic modulus variation, but least sensitive to tortuosities variation. R 1 waves for partially permeable surface radiate as non-physical waves (Im( k ) < 0) at low frequencies. For R 2 waves, velocities are slightly lower than the bulk slow P 2 waves. At low frequencies, both velocity and attenuation are diffusive of f 1/2 frequency dependence, as P 2 waves. It is found that for partially permeable surfaces, the attenuation displays - f 1 frequency dependence as frequency increasing. High surface permeability, low-coupling damping coefficients, low Poisson′s ratios, and low tortuosities increase the slope of the - f 1 dependence. When the attenuation coefficients reach 0, R 2 waves for partially permeable surface begin to radiate as non-physical waves.

Geophysical Journal International↗