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 199 records · Page 11Linked to original sources
USGS2010

Time-dependent seismic tomography

Of methods for measuring temporal changes in seismic-wave speeds in the Earth, seismic tomography is among those that offer the highest spatial resolution. 3-D tomographic methods are commonly applied in this context by inverting seismic wave arrival time data sets from different epochs independently and assuming that differences in the derived structures represent real temporal variations. This assumption is dangerous because the results of independent inversions would differ even if the structure in the Earth did not change, due to observational errors and differences in the seismic ray distributions. The latter effect may be especially severe when data sets include earthquake swarms or aftershock sequences, and may produce the appearance of correlation between structural changes and seismicity when the wave speeds are actually temporally invariant. A better approach, which makes it possible to assess what changes are truly required by the data, is to invert multiple data sets simultaneously, minimizing the difference between models for different epochs as well as the rms arrival-time residuals. This problem leads, in the case of two epochs, to a system of normal equations whose order is twice as great as for a single epoch. The direct solution of this system would require twice as much memory and four times as much computational effort as would independent inversions. We present an algorithm, tomo4d, that takes advantage of the structure and sparseness of the system to obtain the solution with essentially no more effort than independent inversions require. No claim to original US government works Journal compilation ?? 2010 RAS.

Geophysical Journal International↗
USGS2010

The relationship between noise correlation and the Green's function in the presence of degeneracy and the absence of equipartition

Recent derivations have shown that when noise in a physical system has its energy equipartitioned into the modes of the system, there is a convenient relationship between the cross correlation of time-series recorded at two points and the Green's function of the system. Here, we show that even when energy is not fully equipartitioned and modes are allowed to be degenerate, a similar (though less general) property holds for equations with wave equation structure. This property can be used to understand why certain seismic noise correlation measurements are successful despite known degeneracy and lack of equipartition on the Earth. No claim to original US government works Journal compilation ?? 2010 RAS.

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

Lower crustal relaxation beneath the Tibetan Plateau and Qaidam Basin following the 2001 Kokoxili earthquake

In 2001 November a magnitude 7.8 earthquake ruptured a 400 km long portion of the Kunlun fault, northeastern Tibet. In this study, we analyse over five years of post-seismic geodetic data and interpret the observed surface deformation in terms of stress relaxation in the thick Tibetan lower crust. We model GPS time-series (first year) and InSAR line of sight measurements (years two to five) and infer that the most likely mechanism of post-seismic stress relaxation is time-dependent distributed creep of viscoelastic material in the lower crust. Since a single relaxation time is not sufficient to model the observed deformation, viscous flow is modelled by a lower crustal Burgers rheology, which has two material relaxation times. The optimum model has a transient viscosity 9 × 10 17 Pa s, steady-state viscosity 1 × 10 19 Pa s and a ratio of long term to Maxwell shear modulus of 2:3. This model gives a good fit to GPS stations south of the Kunlun Fault, while displacements at stations north of the fault are over-predicted. We attribute this asymmetry in the GPS residual to lateral heterogeneity in rheological structure across the southern margin of the Qaidam Basin, with thinner crust/higher viscosities beneath the basin than beneath the Tibetan Plateau. Deep afterslip localized in a shear zone beneath the fault rupture gives a reasonable match to the observed InSAR data, but the slip model does not fit the earlier GPS data well. We conclude that while some localized afterslip likely occurred during the early post-seismic phase, the bulk of the observed deformation signal is due to viscous flow in the lower crust. To investigate regional variability in rheological structure, we also analyse post-seismic displacements following the 1997 Manyi earthquake that occurred 250 km west of the Kokoxili rupture. We find that viscoelastic properties are the same as for the Kokoxili area except for the transient viscosity, which is 5 × 10 17 Pa s. The viscosities estimated for the Manyi and Kokoxili areas are consistent with constraints obtained from other earthquakes in the northwest and south central parts of the Tibetan Plateau.

Tibetan Plateau and Qaidam Basin↗
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 &plusmn; 0.7 yr) and constant size ( M 4.9 &plusmn; 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&ndash;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&ndash;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&ndash;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↗
USGS2012

Extension of the spatial autocorrelation (SPAC) method to mixed-component correlations of surface waves

Using ambient seismic noise for imaging subsurface structure dates back to the development of the spatial autocorrelation (SPAC) method in the 1950s. We present a theoretical analysis of the SPAC method for multicomponent recordings of surface waves to determine the complete 3 × 3 matrix of correlations between all pairs of three-component motions, called the correlation matrix. In the case of isotropic incidence, when either Rayleigh or Love waves arrive from all directions with equal power, the only non-zero off-diagonal terms in the matrix are the vertical–radial (ZR) and radial–vertical (RZ) correlations in the presence of Rayleigh waves. Such combinations were not considered in the development of the SPAC method. The method originally addressed the vertical–vertical (ZZ), RR and TT correlations, hence the name spatial autocorrelation. The theoretical expressions we derive for the ZR and RZ correlations offer additional ways to measure Rayleigh wave dispersion within the SPAC framework. Expanding on the results for isotropic incidence, we derive the complete correlation matrix in the case of generally anisotropic incidence. We show that the ZR and RZ correlations have advantageous properties in the presence of an out-of-plane directional wavefield compared to ZZ and RR correlations. We apply the results for mixed-component correlations to a data set from Akutan Volcano, Alaska and find consistent estimates of Rayleigh wave phase velocity from ZR compared to ZZ correlations. This work together with the recently discovered connections between the SPAC method and time-domain correlations of ambient noise provide further insights into the retrieval of surface wave Green’s functions from seismic noise.

Geophysical Journal International↗