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.

SourceResearch connectionFrom yearThrough year

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

At least 181 records · Page 10Linked to original sources

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 ﬁrst 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 deﬁning 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 ﬁts 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 conﬁrm 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 ﬂow, 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 conﬁrm 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 inﬂuence 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

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↗

USGS2012

Imaging with cross-hole seismoelectric tomography

We propose a cross-hole imaging approach based on seismoelectric conversions (SC) associated with the transmission of seismic waves from seismic sources located in a borehole to receivers (electrodes) located in a second borehole. The seismoelectric (seismic-to-electric) problem is solved using Biot theory coupled with a generalized Ohm's law with an electrokinetic streaming current contribution. The components of the displacement of the solid phase, the fluid pressure, and the electrical potential are solved using a finite element approach with Perfect Match Layer (PML) boundary conditions for the seismic waves and boundary conditions mimicking an infinite material for the electrostatic problem. We develop an inversion algorithm using the electrical disturbances recorded in the second borehole to localize the position of the heterogeneities responsible for the SC. Because of the ill-posed nature of the inverse problem (inherent to all potential-field problems), regularization is used to constrain the solution at each time in the SC-time window comprised between the time of the seismic shot and the time of the first arrival of the seismic waves in the second borehole. All the inverted volumetric current source densities are aggregated together to produce an image of the position of the heterogeneities between the two boreholes. Two simple synthetic case studies are presented to test this concept. The first case study corresponds to a vertical discontinuity between two homogeneous sub-domains. The second case study corresponds to a poroelastic inclusion (partially saturated by oil) embedded into an homogenous poroelastic formation. In both cases, the position of the heterogeneity is recovered using only the electrical disturbances associated with the SC. That said, a joint inversion of the seismic and seismoelectric data could improve these results.

Geophysical Journal International↗

USGS2012

Source parameters of microearthquakes on an interplate asperity off Kamaishi, NE Japan over two earthquake cycles

We have estimated the source parameters of interplate earthquakes in an earthquake cluster off Kamaishi, NE Japan over two cycles of M~ 4.9 repeating earthquakes. The M~ 4.9 earthquake sequence is composed of nine events that occurred since 1957 which have a strong periodicity (5.5 ± 0.7 yr) and constant size ( M 4.9 ± 0.2), probably due to stable sliding around the source area (asperity). Using P - and S -wave traveltime differentials estimated from waveform cross-spectra, three M~ 4.9 main shocks and 50 accompanying microearthquakes ( M 1.5–3.6) from 1995 to 2008 were precisely relocated. The source sizes, stress drops and slip amounts for earthquakes of M 2.4 or larger were also estimated from corner frequencies and seismic moments using simultaneous inversion of stacked spectral ratios. Relocation using the double-difference method shows that the slip area of the 2008 M~ 4.9 main shock is co-located with those of the 1995 and 2001 M~ 4.9 main shocks. Four groups of microearthquake clusters are located in and around the mainshock slip areas. Of these, two clusters are located at the deeper and shallower edge of the slip areas and most of these microearthquakes occurred repeatedly in the interseismic period. Two other clusters located near the centre of the mainshock source areas are not as active as the clusters near the edge. The occurrence of these earthquakes is limited to the latter half of the earthquake cycles of the M~ 4.9 main shock. Similar spatial and temporal features of microearthquake occurrence were seen for two other cycles before the 1995 M 5.0 and 1990 M 5.0 main shocks based on group identification by waveform similarities. Stress drops of microearthquakes are 3–11 MPa and are relatively constant within each group during the two earthquake cycles. The 2001 and 2008 M~ 4.9 earthquakes have larger stress drops of 41 and 27 MPa, respectively. These results show that the stress drop is probably determined by the fault properties and does not change much for earthquakes rupturing in the same area. The occurrence of microearthquakes in the interseismic period suggests the intrusion of aseismic slip, causing a loading of these patches. We also found that some earthquakes near the centre of the mainshock source area occurred just after the earthquakes at the deeper edge of the mainshock source area. These seismic activities probably indicate episodic aseismic slip migrating from the deeper regions in the mainshock asperity to its centre during interseismic periods. Comparison of the source parameters for the 2001 and 2008 main shocks shows that the seismic moments (1.04 x 10 16 Nm and 1.12 x 10 16 Nm for the 2008 and 2001 earthquakes, respectively) and source sizes (radius = 570 m and 540 m for the 2008 and 2001 earthquakes, respectively) are comparable. Based on careful phase identification and hypocentre relocation by constraining the hypocentres of other small earthquakes to their precisely located centroids, we found that the hypocentres of the 2001 and 2008 M~ 4.9 events are located in the southeastern part of the mainshock source area. This location does not correspond to either episodic slip area or hypocentres of small earthquakes that occurred during the earthquake cycle.

Kamaishi↗

USGS2012

W phase source inversion for moderate to large earthquakes (1990-2010)

Rapid characterization of the earthquake source and of its effects is a growing field of interest. Until recently, it still took several hours to determine the first-order attributes of a great earthquake (e.g. M w ≥ 7.5), even in a well-instrumented region. The main limiting factors were data saturation, the interference of different phases and the time duration and spatial extent of the source rupture. To accelerate centroid moment tensor (CMT) determinations, we have developed a source inversion algorithm based on modelling of the W phase, a very long period phase (100–1000 s) arriving at the same time as the P wave. The purpose of this work is to finely tune and validate the algorithm for large-to-moderate-sized earthquakes using three components of W phase ground motion at teleseismic distances. To that end, the point source parameters of all M w ≥ 6.5 earthquakes that occurred between 1990 and 2010 (815 events) are determined using Federation of Digital Seismograph Networks, Global Seismographic Network broad-band stations and STS1 global virtual networks of the Incorporated Research Institutions for Seismology Data Management Center. For each event, a preliminary magnitude obtained from W phase amplitudes is used to estimate the initial moment rate function half duration and to define the corner frequencies of the passband filter that will be applied to the waveforms. Starting from these initial parameters, the seismic moment tensor is calculated using a preliminary location as a first approximation of the centroid. A full CMT inversion is then conducted for centroid timing and location determination. Comparisons with Harvard and Global CMT solutions highlight the robustness of W phase CMT solutions at teleseismic distances. The differences in M w rarely exceed 0.2 and the source mechanisms are very similar to one another. Difficulties arise when a target earthquake is shortly (e.g. within 10 hr) preceded by another large earthquake, which disturbs the waveforms of the target event. To deal with such difficult situations, we remove the perturbation caused by earlier disturbing events by subtracting the corresponding synthetics from the data. The CMT parameters for the disturbed event can then be retrieved using the residual seismograms. We also explore the feasibility of obtaining source parameters of smaller earthquakes in the range 6.0 ≤M w < 6.5. Results suggest that the W phase inversion can be implemented reliably for the majority of earthquakes of M w = 6 or larger.

Geophysical Journal International↗

USGS2012

Hybrid analysis of multiaxis electromagnetic data for discrimination of munitions and explosives of concern

The remediation of land containing munitions and explosives of concern, otherwise known as unexploded ordnance, is an ongoing problem facing the U.S. Department of Defense and similar agencies worldwide that have used or are transferring training ranges or munitions disposal areas to civilian control. The expense associated with cleanup of land previously used for military training and war provides impetus for research towards enhanced discrimination of buried unexploded ordnance. Towards reducing that expense, a multiaxis electromagnetic induction data collection and software system, called ALLTEM, was designed and tested with support from the U.S. Department of Defense Environmental Security Technology Certification Program. ALLTEM is an on-time time-domain system that uses a continuous triangle-wave excitation to measure the target-step response rather than traditional impulse response. The system cycles through three orthogonal transmitting loops and records a total of 19 different transmitting and receiving loop combinations with a nominal spatial data sampling interval of 20 cm. Recorded data are pre-processed and then used in a hybrid discrimination scheme involving both data-driven and numerical classification techniques. The data-driven classification scheme is accomplished in three steps. First, field observations are used to train a type of unsupervised artificial neural network, a self-organizing map (SOM). Second, the SOM is used to simultaneously estimate target parameters (depth, azimuth, inclination, item type and weight) by iterative minimization of the topographic error vectors. Third, the target classification is accomplished by evaluating histograms of the estimated parameters. The numerical classification scheme is also accomplished in three steps. First, the Biot–Savart law is used to model the primary magnetic fields from the transmitter coils and the secondary magnetic fields generated by currents induced in the target materials in the ground. Second, the target response is modelled by three orthogonal dipoles from prolate, oblate and triaxial ellipsoids with one long axis and two shorter axes. Each target consists of all three dipoles. Third, unknown target parameters are determined by comparing modelled to measured target responses. By comparing the rms error among the self-organizing map and numerical classification results, we achieved greater than 95 per cent detection and correct classification of the munitions and explosives of concern at the direct fire and indirect fire test areas at the UXO Standardized Test Site at the Aberdeen Proving Ground, Maryland in 2010.

Geophysical Journal International↗

USGS2012

Determination and uncertainty of moment tensors for microearthquakes at Okmok Volcano, Alaska

Efforts to determine general moment tensors (MTs) for microearthquakes in volcanic areas are often hampered by small seismic networks, which can lead to poorly constrained hypocentres and inadequate modelling of seismic velocity heterogeneity. In addition, noisy seismic signals can make it difficult to identify phase arrivals correctly for small magnitude events. However, small volcanic earthquakes can have source mechanisms that deviate from brittle double-couple shear failure due to magmatic and/or hydrothermal processes. Thus, determining reliable MTs in such conditions is a challenging but potentially rewarding pursuit. We pursued such a goal at Okmok Volcano, Alaska, which erupted recently in 1997 and in 2008. The Alaska Volcano Observatory operates a seismic network of 12 stations at Okmok and routinely catalogues recorded seismicity. Using these data, we have determined general MTs for seven microearthquakes recorded between 2004 and 2007 by inverting peak amplitude measurements of P and S phases. We computed Green's functions using precisely relocated hypocentres and a 3-D velocity model. We thoroughly assessed the quality of the solutions by computing formal uncertainty estimates, conducting a variety of synthetic and sensitivity tests, and by comparing the MTs to solutions obtained using alternative methods. The results show that MTs are sensitive to station distribution and errors in the data, velocity model and hypocentral parameters. Although each of the seven MTs contains a significant non-shear component, we judge several of the solutions to be unreliable. However, several reliable MTs are obtained for a group of previously identified repeating events, and are interpreted as compensated linear-vector dipole events.

Alaska↗

USGS2012

The 1909 Taipei earthquake: implication for seismic hazard in Taipei

The 1909 April 14 Taiwan earthquake caused significant damage in Taipei. Most of the information on this earthquake available until now is from the written reports on its macro-seismic effects and from seismic station bulletins. In view of the importance of this event for assessing the shaking hazard in the present-day Taipei, we collected historical seismograms and station bulletins of this event and investigated them in conjunction with other seismological data. We compared the observed seismograms with those from recent earthquakes in similar tectonic environments to characterize the 1909 earthquake. Despite the inevitably large uncertainties associated with old data, we conclude that the 1909 Taipei earthquake is a relatively deep (50–100 km) intraplate earthquake that occurred within the subducting Philippine Sea Plate beneath Taipei with an estimated M_W of 7 ± 0.3. Some intraplate events elsewhere in the world are enriched in high-frequency energy and the resulting ground motions can be very strong. Thus, despite its relatively large depth and a moderately large magnitude, it would be prudent to review the safety of the existing structures in Taipei against large intraplate earthquakes like the 1909 Taipei earthquake.

Taipei↗