Search USGS⌕ Search

SEARCH · Search USGS

Results for “Journal of Seismology”

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 55 records · Page 3Linked to original sources

Hilbert-Huang transform analysis of dynamic and earthquake motion recordings

This study examines the rationale of Hilbert-Huang transform (HHT) for analyzing dynamic and earthquake motion recordings in studies of seismology and engineering. In particular, this paper first provides the fundamentals of the HHT method, which consist of the empirical mode decomposition (EMD) and the Hilbert spectral analysis. It then uses the HHT to analyze recordings of hypothetical and real wave motion, the results of which are compared with the results obtained by the Fourier data processing technique. The analysis of the two recordings indicates that the HHT method is able to extract some motion characteristics useful in studies of seismology and engineering, which might not be exposed effectively and efficiently by Fourier data processing technique. Specifically, the study indicates that the decomposed components in EMD of HHT, namely, the intrinsic mode function (IMF) components, contain observable, physical information inherent to the original data. It also shows that the grouped IMF components, namely, the EMD-based low- and high-frequency components, can faithfully capture low-frequency pulse-like as well as high-frequency wave signals. Finally, the study illustrates that the HHT-based Hilbert spectra are able to reveal the temporal-frequency energy distribution for motion recordings precisely and clearly.

Journal of Engineering Mechanics↗

An evaluation of the strong ground motion recorded during the May 1, 2003 Bingol Turkey, earthquake

An important record of ground motion from a M 6.4 earthquake occurring on May 1, 2003, at epicentral and fault distances of about 12 and 9 km, respectively, was obtained at a station near the city of Bingöl, Turkey. The maximum peak ground values of 0.55g and 36 cm/s are among the largest ground-motion amplitudes recorded in Turkey. From simulations and comparisons with ground motions from other earthquakes of comparable magnitude, we conclude that the ground motion over a range of frequencies is unusually high. Site response may be responsible for the elevated ground motion, as suggested from analysis of numerous aftershock recordings from the same station. The mainshock motions have some interesting seismological features, including ramps between the P- and S-wave that are probably due to near- and intermediate-field elastic motions and strong polarisation oriented at about 39 degrees to the fault (and therefore not in the fault-normal direction). Simulations of motions from an extended rupture explain these features. The N10E component shows a high-amplitude spectral acceleration at a period of 0.15 seconds resulting in a site specific design spectrum that significantly overestimates the actual strength and displacement demands of the record. The pulse signal in the N10E component affects the inelastic spectral displacement and increases the inelastic displacement demand with respect to elastic demand for very long periods.

Journal of Earthquake Engineering↗

Seismological analyses of the 2010 March 11, Pichilemu, Chile Mw 7.0 and Mw 6.9 coastal intraplate earthquakes

On 2010 March 11, a sequence of large, shallow continental crust earthquakes shook central Chile. Two normal faulting events with magnitudes around M w 7.0 and M w 6.9 occurred just 15 min apart, located near the town of Pichilemu. These kinds of large intraplate, inland crustal earthquakes are rare above the Chilean subduction zone, and it is important to better understand their relationship with the 2010 February 27, M w 8.8, Maule earthquake, which ruptured the adjacent megathrust plate boundary. We present a broad seismological analysis of these earthquakes by using both teleseismic and regional data. We compute seismic moment tensors for both events via a W-phase inversion, and test sensitivities to various inversion parameters in order to assess the stability of the solutions. The first event, at 14 hr 39 min GMT, is well constrained, displaying a fault plane with strike of N145°E, and a preferred dip angle of 55°SW, consistent with the trend of aftershock locations and other published results. Teleseismic finite-fault inversions for this event show a large slip zone along the southern part of the fault, correlating well with the reported spatial density of aftershocks. The second earthquake (14 hr 55 min GMT) appears to have ruptured a fault branching southward from the previous ruptured fault, within the hanging wall of the first event. Modelling seismograms at regional to teleseismic distances (Δ > 10°) is quite challenging because the observed seismic wave fields of both events overlap, increasing apparent complexity for the second earthquake. We perform both point- and extended-source inversions at regional and teleseismic distances, assessing model sensitivities resulting from variations in fault orientation, dimension, and hypocentre location. Results show that the focal mechanism for the second event features a steeper dip angle and a strike rotated slightly clockwise with respect to the previous event. This kind of geological fault configuration, with secondary rupture in the hanging wall of a large normal fault, is commonly observed in extensional geological regimes. We propose that both earthquakes form part of a typical normal fault diverging splay, where the secondary fault connects to the main fault at depth. To ascertain more information on the spatial and temporal details of slip for both events, we gathered near-fault seismological and geodetic data. Through forward modelling of near-fault synthetic seismograms we build a kinematic k −2 earthquake source model with spatially distributed slip on the fault that, to first-order, explains both coseismic static displacement GPS vectors and short-period seismometer observations at the closest sites. As expected, the results for the first event agree with the focal mechanism derived from teleseismic modelling, with a magnitude M w 6.97. Similarly, near-fault modelling for the second event suggests rupture along a normal fault, M w 6.90, characterized by a steeper dip angle (dip = 74°) and a strike clockwise rotated (strike = 155°) with respect to the previous event.

Pichilemu↗

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↗

The 1992 M=7 Cape Mendocino, California, earthquake: Coseismic deformation at the south end of the Cascadia megathrust

We invert geodetic measurements of coseismic surface displacements to determine a dislocation model for the April 25, 1992, M = 7 Cape Mendocino, California, earthquake. The orientation of the model slip vector, which nearly parallels North America-Juan de Fuca relative plate convergence, and the location and orientation of the model fault relative to the offshore Cascadia megathrust, suggest that the 1992 Cape Mendocino earthquake is the first well-recorded event to relieve strain associated with the Cascadia subduction zone. We use data from three geodetic techniques: (1) the horizontal and vertical displacements of 13 monuments surveyed with the Global Positioning System, corrected for observed horizontal interseismic strain accumulation, (2) 88 section-elevation differences between leveling monuments, and (3) the uplift of 12 coastal sites observed from the die-off of intertidal marine organisms. Maximum observed displacements are 0.4 m of horizontal movement and 1.5 m of uplift along the coast. We use Monte Carlo techniques to estimate an optimal uniform slip rectangular fault geometry and its uncertainties. The optimal model using all the data resolves 4.9 m of slip on a 14 by 15 km fault that dips 28° SE. The fault extends from 1.5 to 8.7 km in depth and the main-shock hypocenter is close to the downdip projection of the fault. The shallowly dipping fault plane is consistent with the observed aftershock locations, and the estimated geodetic moment is 3.1 × 10 19 N m, 70% of the seismic moment. Other models that exclude leveling data collected in 1935 and 1942 are more consistent with seismological estimates of the fault geometry. If the earthquake is characteristic for this segment, the estimated horizontal slip vector compared with plate convergence rates suggests a recurrence interval of 140 years, with a 95% confidence range of 100–670 years. The coseismic uplift occurred in a region that also has high Quaternary uplift rates determined from marine terrace studies. If repeated ruptures of this southernmost segment of the Cascadia megathrust are responsible for the Quaternary uplift, a comparison of the coseismic uplift with coastal uplift rates suggests a recurrence interval of 200–400 years. Thus comparing horizontal and vertical coseismic to long-term deformation suggests a recurrence interval of about 100–300 years for M = 7 events at the south end of the Cascadia megathrust.

Journal of Geophysical Research B: Solid Earth↗

Missing great earthquakes

The occurrence of three earthquakes with moment magnitude (Mw) greater than 8.8 and six earthquakes larger than Mw 8.5, since 2004, has raised interest in the long-term global rate of great earthquakes. Past studies have focused on the analysis of earthquakes since 1900, which roughly marks the start of the instrumental era in seismology. Before this time, the catalog is less complete and magnitude estimates are more uncertain. Yet substantial information is available for earthquakes before 1900, and the catalog of historical events is being used increasingly to improve hazard assessment. Here I consider the catalog of historical earthquakes and show that approximately half of all Mw ≥ 8.5 earthquakes are likely missing or underestimated in the 19th century. I further present a reconsideration of the felt effects of the 8 February 1843, Lesser Antilles earthquake, including a first thorough assessment of felt reports from the United States, and show it is an example of a known historical earthquake that was significantly larger than initially estimated. The results suggest that incorporation of best available catalogs of historical earthquakes will likely lead to a significant underestimation of seismic hazard and/or the maximum possible magnitude in many regions, including parts of the Caribbean.

Journal of Geophysical Research B: Solid Earth↗

The proliferation of induced seismicity in the Permian Basin, Texas

The Permian Basin has a long history of induced earthquakes, but the seismicity rates have increased dramatically over the past two decades and included a M W 5.0 likely induced by wastewater disposal (WD) in March 2020. A detailed characterization of the proliferation of seismicity in the Permian Basin throughout this time period is needed for improving the scientific understanding of the mechanisms responsible and for mitigating future seismic hazard. Due to a sparse regional seismic network before the advent of Texas Seismological Network in 2017, we characterize seismicity using the 10-station TXAR array that is 100s of km away from most of the seismicity, with the objective of improving upon the substantial contributions from previous work. By exploiting the nature of waveform similarity, we detect events with template matching, performing a quantitative analysis of spatially varying detection capabilities throughout the study area. From an initial catalog of 10,753 events, we identify 45,009 earthquakes and 10,208 quarry blasts. Using our catalog of earthquakes, we improve epicentral locations, compare relative magnitude techniques, and associate earthquakes to WD or hydraulic stimulations. We further use our earthquake catalog to investigate the relationship between seismicity and human activities near the city of Pecos, Texas. Through a comparison of our earthquake catalog with industrial records, we determine that the vast majority seismicity near Pecos, Texas, since 2000 is likely induced by an increase of WD at wells injecting at depths greater than 1.5 km.

Texas↗

Ground-motion residuals, path effects, and crustal properties: A pilot study in southern California

To improve models of ground motion estimation and probabilistic seismic hazard analyses, the engineering seismology field is moving toward developing fully nonergodic ground motion models, models specific for individual source‐to‐site paths. Previous work on this topic has examined systematic variations in ground‐motion along particular paths (from either recorded or simulated earthquake data) and has not included physical properties of the path. We present here a framework to include physical path properties, by seeking correlations between ground motion amplitudes along specific paths and crustal properties, specifically seismic velocity and anelastic attenuation, along that path. Using a large data set of small‐magnitude earthquakes recorded in Southern California, we find a correlation between the gradient of seismic S wave velocity and the path term residual, after accounting for an average geometric spreading and anelastic attenuation, indicating that heterogeneity in crustal velocity primarily controls the path‐specific attenuation. Even in aseismic regions, details of path‐specific ground motion prediction equations can be developed from crustal structure and property data.

California↗

Estimating tectonic history through basin simulation-enhanced seismic inversion: Geoinformatics for sedimentary basins

A data assimilation approach is demonstrated whereby seismic inversion is both automated and enhanced using a comprehensive numerical sedimentary basin simulator to study the physics and chemistry of sedimentary basin processes in response to geothermal gradient in much greater detail than previously attempted. The approach not only reduces costs by integrating the basin analysis and seismic inversion activities to understand the sedimentary basin evolution with respect to geodynamic parameters-but the technique also has the potential for serving as a geoinfomatics platform for understanding various physical and chemical processes operating at different scales within a sedimentary basin. Tectonic history has a first-order effect on the physical and chemical processes that govern the evolution of sedimentary basins. We demonstrate how such tectonic parameters may be estimated by minimizing the difference between observed seismic reflection data and synthetic ones constructed from the output of a reaction, transport, mechanical (RTM) basin model. We demonstrate the method by reconstructing the geothermal gradient. As thermal history strongly affects the rate of RTM processes operating in a sedimentary basin, variations in geothermal gradient history alter the present-day fluid pressure, effective stress, porosity, fracture statistics and hydrocarbon distribution. All these properties, in turn, affect the mechanical wave velocity and sediment density profiles for a sedimentary basin. The present-day state of the sedimentary basin is imaged by reflection seismology data to a high degree of resolution, but it does not give any indication of the processes that contributed to the evolution of the basin or causes for heterogeneities within the basin that are being imaged. Using texture and fluid properties predicted by our Basin RTM simulator, we generate synthetic seismograms. Linear correlation using power spectra as an error measure and an efficient quadratic optimization technique are found to be most effective in determining the optimal value of the tectonic parameters. Preliminary 1-D studies indicate that one can determine the geothermal gradient even in the presence of observation and numerical uncertainties. The algorithm succeeds even when the synthetic data has detailed information only in a limited depth interval and has a different dominant frequency in the synthetic and observed seismograms. The methodology presented here even works when the basin input data contains only 75 per cent of the stratigraphic layering information compared with the actual basin in a limited depth interval.

Geophysical Journal International↗

Constraints on the viscosity of the continental crust and mantle from GPS measurements and postseismic deformation models in western Mongolia

We use GPS measurements and models of postseismic deformation caused by seven M6.8 to 8.4 earthquakes that occurred in the past 100 years in Mongolia to assess the viscosity of the lower crust and upper mantle. We find an upper mantle viscosity between 1 ?? 1018 and 4 ?? 1018 Pa s. The presence of such a weak mantle is consistent with results from independent seismological and petrological studies that show an abnormally hot upper mantle beneath Mongolia. The viscosity of the lower crust is less well constrained, but a weak lower crust (3 ?? 1016 to 2 ?? 1017 Pa s) is preferred by the data. Using our best fit upper mantle and lower crust viscosities, we find that the postseismic effects of viscoelastic relaxation on present-day horizontal GPS velocities are small (<2 mm yr-1) but still persist 100 years after the 1905, M8.4, Bolnay earthquake. This study shows that the GPS velocity field in the Baikal-Mongolia area can be modeled as the sum of (1) a rigid translation and rotation of the whole network, (2) a 3-5 mm yr-1 simple shear velocity gradient between the Siberian platform to the north and northern China to the south, and (3) the contribution of postseismic deformation, mostly caused by the 1905 Bolnay-Tsetserleg sequence and by the smaller, but more recent, 1957 Bogd earthquake. Copyright 2003 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

Afterslip and creep in the rate-dependent framework: Joint inversion of borehole strain and GNSS displacements for the Mw 7.1 Ridgecrest earthquake

The elusive transition toward afterslip following an earthquake is challenging to capture with typical data resolution limits. A dense geodetic network recorded the Mw 7.1 Ridgecrest earthquake, including 16 Global Navigation Satellite System (GNSS) stations and 3 borehole strainmeters (BSM). The sub-nanostrain precision and sub-second sampling rate of BSMs bridges a gap between conventional seismologic and geodetic methods, exemplified by atypical postseismic shear strain reversals observed at nearfield (<2 km) station B921 that remain unexplained. We jointly invert GNSS displacements and BSM strains for coseismic and postseismic slip spanning hours to months over 7 independent periods. Cosiesmically, our model resolves the largest slip magnitudes of up to 6.6 m on the mainshock rupture plane, with similar patterns to other inferred slip distributions. The foreshock fault appears to slip coincidently with mainshock, revealing potential asperities activated during the preceding Mw 6.4 event. Postseismically, the best-fitting models adhere to mechanical rate-and-state expectations of logarithmically decaying slip adjacent to the coseismic rupture terminus, and where deep rheologic conditions favor creep. Most spatial variation occurs in the early postseismic timeframe (<1–2 weeks), with evidence for regional rheologic control and static stress dependence. Triggered creep on the neighboring Garlock Fault unexpectedly persists for >178 days—further highlighting the importance of fault networks in postseismic stress redistribution, critical to assessing future hazard.

Journal of Geophysics Research↗

Use of a hybrid computer in engineering-seismology research

A hybrid computer is an important tool in the seismological research conducted by the U.S. Geological Survey in support of the Energy Research and Development Administration nuclear explosion testing program at the Nevada Test Site and the U.S. Geological Survey Earthquake Hazard Reduction Program. The hybrid computer system, which employs both digital and analog computational techniques, facilitates efficient seismic data processing. Standard data processing operations include: (1) preview of dubbed magnetic tapes of data; (2) correction of data for instrument response; (3) derivation of displacement and acceleration time histories from velocity recordings; (4) extraction of peak-amplitude data; (5) digitization of time histories; (6) rotation of instrumental axes; (7) derivation of response spectra; and (8) derivation of relative transfer functions between recording sites. Catalog of time histories and response spectra of ground motion from nuclear explosions and earthquakes that have been processed by the hybrid computer are used in the Earthquake Hazard Research Program to evaluate the effects of source, propagation path, and site effects on recorded ground motion; to assess seismic risk; to predict system response; and to solve system design problems.

Journal of Research of the U.S. Geological Survey↗

Seismicity within a propagating ice shelf rift: the relationship between icequake locations and ice shelf structure

Iceberg calving is a dominant mass loss mechanism for Antarctic ice shelves, second only to basal melting. An important known process involved in calving is the initiation and propagation of through-penetrating fractures called rifts; however, the mechanisms controlling rift propagation remain poorly understood. To investigate the mechanics of ice-shelf rifting, we analyzed seismicity associated with a propagating rift tip on the Amery Ice Shelf, using data collected during the Austral summers of 2004-2007. We investigated seismicity associated with fracture propagation using a suite of passive seismological techniques including icequake locations, back projection, and moment tensor inversion. We confirm previous results that show that seismicity is characterized by periods of relative quiescence punctuated by swarms of intense seismicity of one to three hours. However, even during periods of quiescence, we find significant seismic deformation around the rift tip. Moment tensors, calculated for a subset of the largest icequakes (M W > -2.0) located near the rift tip, show steeply dipping fault planes, horizontal or shallowly plunging stress orientations, and often have a significant volumetric component. They also reveal that much of the observed seismicity is limited to the upper 50 m of the ice shelf. This suggests a complex system of deformation that involves the propagating rift, the region behind the rift tip, and a system of rift-transverse crevasses. Small-scale variations in the mechanical structure of the ice shelf, especially rift-transverse crevasses and accreted marine ice, play an important role in modulating the rate and location of seismicity associated with propagating ice shelf rifts.

Journal of Geophysical Research F: Earth Surface↗

Hypothesis: Many earthquakes in the central and southeastern United States are causally related to mafic intrusive bodies

Assessment of earthquake hazards in the United States is based largely upon knowledge acquired in the seismically active parts of the western United States. Earthquakes in the central and southeastern United States are seismologically and geologically very anomalous, however, compared with those in the western United States. For example, shallow earthquakes of magnitude 7-f in the western United States in general have tectonic surface faulting associated with them. No tectonic surface faulting, however, has been found to be associated with the New Madrid, Mo., earthquakes of 1811 and 1812, several of which had estimated magnitudes of 7-f. The lack of surface faulting is, in part, responsible for the fact that the New Madrid earthquakes have not been related with certainty to any particular geologic structural element, which makes assessment of the seismic hazards associated with possible similar future earthquakes very difficult. Geologic and seismologic data from the central and southeastern United States, however, suggest a spatial correlation between mafic intrusive rocks and seismic activity. An additional supportive correlation is that the trend of many dikes coincides with the trend of nodal planes in many fault-plane solutions. These correlations differ greatly from the direct correlation of earthquakes with faults in the western United States. Consideration of different concepts of earthquake sources in the central and southeastern. United States seems necessary- The preferred concept is that local high-stress concentrations, which may result in earthquakes, occur at or close to the contact of mafic intrusives and their felsic host rocks. Ancient rift zones may have been the primary control of the location of the mafic intrusives. Implications of a causal relation between mafic intrusives and seismicity are that more useful seismo-tectonic maps may be possible and that the source dimensions of the earthquakes would be small. Because of low attenuation, however, effects of earthquakes in the central and southeastern United States, even if relatively small because of smaller source dimensions, must be considered hazardous.

Journal of Research of the U.S. Geological Survey↗

Icelandic-type crust

Numerous seismic studies, in particular using receiver functions and explosion seismology, have provided a detailed picture of the structure and thickness of the crust beneath the Iceland transverse ridge. We review the results and propose a structural model that is consistent with all the observations. The upper crust is typically 7 ?? 1 km thick, heterogeneous and has high velocity gradients. The lower crust is typically 15-30 ?? 5 km thick and begins where the velocity gradient decreases radically. This generally occurs at the V p ??? 6.5 km s-1 level. A low-velocity zone ??? 10 000 km2 in area and up to ??? 15 km thick occupies the lower crust beneath central Iceland, and may represent a submerged, trapped oceanic microplate. The crust-mantle boundary is a transition zone ???5 ?? 3 km thick throughout which V p increases progressively from ???7.2 to ???8.0 km s-1. It may be gradational or a zone of alternating high- and low-velocity layers. There is no seismic evidence for melt or exceptionally high temperatures in or near this zone. Isostasy indicates that the density contrast between the lower crust and the mantle is only ???90 kg m-3 compared with ???300 kg m-3 for normal oceanic crust, indicating compositional anomalies that are as yet not understood. The seismological crust is ???30 km thick beneath the Greenland-Iceland and Iceland-Faeroe ridges, and eastern Iceland, ???20 km beneath western Iceland, and ???40 km thick beneath central Iceland. This pattern is not what is predicted for an eastward-migrating plume. Low attenuation and normal V p/V s ratios in the lower crust beneath central and southwestern Iceland, and normal uppermost mantle velocities in general, suggest that the crust and uppermost mantle are subsolidus and cooler than at equivalent depths beneath the East Pacific Rise. Seismic data from Iceland have historically been interpreted both in terms of thin-hot and thick-cold crust models, both of which have been cited as supporting the plume hypothesis. This suggests that the plume model for Iceland is an a priori assumption rather than a hypothesis subject to testing. The long-extinct Ontong-Java Plateau, northwest India and Parana??, Brazil large igneous provinces, beneath which mantle plumes are not expected are all underlain by mantle low-velocity bodies similar to that beneath Iceland. A plume interpretation for the mantle anomaly beneath Iceland is thus not required.

Geophysical Journal International↗

Slab1.0: A three-dimensional model of global subduction zone geometries

We describe and present a new model of global subduction zone geometries, called Slab1.0. An extension of previous efforts to constrain the two-dimensional non-planar geometry of subduction zones around the focus of large earthquakes, Slab1.0 describes the detailed, non-planar, three-dimensional geometry of approximately 85% of subduction zones worldwide. While the model focuses on the detailed form of each slab from their trenches through the seismogenic zone, where it combines data sets from active source and passive seismology, it also continues to the limits of their seismic extent in the upper-mid mantle, providing a uniform approach to the definition of the entire seismically active slab geometry. Examples are shown for two well-constrained global locations; models for many other regions are available and can be freely downloaded in several formats from our new Slab1.0 website, http://on.doi.gov/d9ARbS . We describe improvements in our two-dimensional geometry constraint inversion, including the use of ‘average’ active source seismic data profiles in the shallow trench regions where data are otherwise lacking, derived from the interpolation between other active source seismic data along-strike in the same subduction zone. We include several analyses of the uncertainty and robustness of our three-dimensional interpolation methods. In addition, we use the filtered, subduction-related earthquake data sets compiled to build Slab1.0 in a reassessment of previous analyses of the deep limit of the thrust interface seismogenic zone for all subduction zones included in our global model thus far, concluding that the width of these seismogenic zones is on average 30% larger than previous studies have suggested.

Journal of Geophysical Research B: Solid Earth↗

Getting the story straight: Common misconceptions around volcanoes and eruptions with case studies of recent events

The ability to rapidly disseminate scientific knowledge across the world has never been greater. This is paralleled by the ability of rumors to spread equally far and fast. Popular misconceptions can take root and, without correction, take on the appearance of facts. Commonly misunderstood topics include those with attention-grabbing phrasing (e.g. “Pacific Ring of Fire”, “supervolcanoes”, or “mega-tsunamis”), the intersection of volcanic research with high-profile science topics (climate or human health), aspirational science goals (eruption forecasts), and basic terminology (volcanic ash versus smoke, or magma “chambers”). In this paper, we describe eighteen misconceptions grouped in six themes (seismology, volcanoes and climate, systems and structures, volcanic hazards, timescales, photographs and footage) commonly encountered by volcanologists who are actively engaging with the public and media both during and between eruptions. Each misconception description is accompanied by takeaway points and summaries of recent research to assist the reader with the information necessary to describe the issues. Finally, we present seven case studies where incorrect information was widely shared and counteracted at well-known volcanoes and recent eruptions at Mount Agung (Bali, Indonesia), Mount Etna (Sicily, Italy), Kīlauea (Hawaii, U.S.A.), Tajogaite (La Palma, Spain), La Soufrière (St. Vincent), Taal Volcano (Philippines), and Yellowstone caldera (U.S.A.). This work is intended to help both volcanologists and non-specialists get the story straight, and to effectively communicate processes and scales of volcanism and volcanic hazards, impacts, and risks.

Journal of Applied Volcanology↗

Earthquake fracture energy inferred from kinematic rupture models on extended faults

We estimate fracture energy on extended faults for several recent earthquakes by retrieving dynamic traction evolution at each point on the fault plane from slip history imaged by inverting ground motion waveforms. We define the breakdown work (Wb) as the excess of work over some minimum traction level achieved during slip. Wb is equivalent to "seismological" fracture energy (G) in previous investigations. Our numerical approach uses slip velocity as a boundary condition on the fault. We employ a three-dimensional finite difference algorithm to compute the dynamic traction evolution in the time domain during the earthquake rupture. We estimate Wb by calculating the scalar product between dynamic traction and slip velocity vectors. This approach does not require specifying a constitutive law and assuming dynamic traction to be collinear with slip velocity. If these vectors are not collinear, the inferred breakdown work depends on the initial traction level. We show that breakdown work depends on the square of slip. The spatial distribution of breakdown work in a single earthquake is strongly correlated with the slip distribution. Breakdown work density and its integral over the fault, breakdown energy, scale with seismic moment according to a power law (with exponent 0.59 and 1.18, respectively). Our estimates of breakdown work range between 4 ?? 105 and 2 ?? 107 J/m2 for earthquakes having moment magnitudes between 5.6 and 7.2. We also compare our inferred values with geologic surface energies. This comparison might suggest that breakdown work for large earthquakes goes primarily into heat production. Copyright 2005 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗