The border between Georgia and South Carolina has a moderate amount of seismicity typical of the Piedmont Province of the eastern United States and greater than most other intraplate regions. Historical records suggest on average a M-w 4.5 earthquake every 50 yr in the region of the J. Strom Thurmond Reservoir, which is located on the border between Georgia and South Carolina. The M-w 4.1 earthquake on 15 February 2014 near Edgefield, South Carolina, was one of the largest events in this region recorded by nearby modern seismometers, providing an opportunity to study its source properties and aftershock productivity. Using the waveforms of the M-w 4.1 mainshock and the only cataloged M w 3.0 aftershock as templates, we apply a matched-filter technique to search for additional events between 8 and 22 February 2014. The resulting six new detections are further employed as new templates to scan for more events. Repeating the waveform-matching method with new templates yields 13 additional events, for a total of 19 previously unidentified events with magnitude 0.06 and larger. The low number of events suggests that this sequence is deficient in aftershock production, as compared with expected aftershock productivities for other mainshocks of similar magnitudes. Hypocentral depths of the M-w 4.1 mainshock and M-w 3.0 aftershock are estimated by examining the differential time between a depth phase called sPL and P-wave arrivals, as well as by modeling the depth phase of body waves at shorter periods. The best-fitting depths for both events are around 3-4 km. The obtained stress drops for the M-w 4.1 mainshock and M-w 3.0 aftershock are 3.75 and 4.44 MPa, respectively. The corresponding updated moment magnitude for the aftershock is 2.91.

Between 1965 and 2003, the Carnegie Institution of Washington's Department of Terrestrial Magnetism operated a continuous network of nine broadband seismographs with locations in South America, Japan, Iceland, Papua New Guinea, and Washington, D.C. The Carnegie seismographs designed in the 1960s by Selwyn Sacks were among the earliest broadband instruments, sensing between at least 30 s and similar to 30 Hz. Given the scarcity of historic seismic data of comparable bandwidth and dynamic range prior to the widespread shift to force-feedback instruments and digital recording around the mid-1980s, this dataset is still of high scientific value today.

Processes related to eruptions at arc volcanoes are linked by structures that transect the entire crust. Imaging the mid- to lower-crustal portions (here, similar to 5-15km and>15km respectively) of these magmatic systems where intermediate storage may occur has been a longstanding challenge. Tomography, local seismic source studies, geodetic, and geochemical constraints, are typically most sensitive to shallow (<5km) storage and/or have insufficient resolution at these depths. Geophysical methods are even further limited at frequently-erupting volcanoes where well-developed trans-crustal magmatic systems are likely to exist, due to a lack of deep seismicity. Here we show direct evidence for mid-crustal magma storage beneath the frequently erupting Cleveland volcano, Alaska, using a novel application of seismic receiver functions. We use P-s scattered waves from the Moho as virtual sources to investigate S-wave velocities between the Moho and the surface. Our forward modeling approach allows us to provide direct constraints on the geometry of low velocity regions beneath volcanoes despite having a comparatively sparse seismic network. Our results show clear evidence of mid-crustal magma storage beneath the depths of located volcanic seismicity. Future work using similar approaches will enable an unprecedented comparative examination of magmatic systems beneath sparsely instrumented volcanoes globally.

The nature and cause of deep earthquakes remain enduring unknowns in the field of seismology. We present new models of thermal structures of subducted slabs traced to mantle transition zone depths that permit a detailed comparison between slab pressure/temperature (P/T) paths and hydrated/carbonated mineral phase relations. We find a remarkable correlation between slabs capable of transporting water to transition zone depths in dense hydrous magnesium silicates with slabs that produce seismicity below similar to 300-km depth, primarily between 500 and 700 km. This depth range also coincides with the P/T conditions at which oceanic crustal lithologies in cold slabs are predicted to intersect the carbonate-bearing basalt solidus to produce carbonatitic melts. Both forms of fluid evolution are well represented by sublithospheric diamonds whose inclusions record the existence of melts, fluids, or supercritical liquids derived from hydrated or carbonate-bearing slabs at depths (similar to 300-700 km) generally coincident with deep-focus earthquakes. We propose that the hydrous and carbonated fluids released from subducted slabs at these depths lead to fluid-triggered seismicity, fluid migration, diamond precipitation, and inclusion crystallization. Deep focus earthquake hypocenters could track the general region of deep fluid release, migration, and diamond formation in the mantle. The thermal modeling of slabs in the mantle and the correlation between sublithospheric diamonds, deep focus earthquakes, and slabs at depth demonstrate a deep subduction pathway to the mantle transition zone for carbon and volatiles that bypasses shallower decarbonation and dehydration processes.

Seismic tomography of the crust is an essential tool for studying the three-dimensional structure of magmatic plumbing systems feeding active volcanoes, but it is often limited in resolution by the absence of deep local seismicity. Teleseismic receiver functions can be used to illuminate local structural variations, but typically do not account for the effects of three-dimensional velocity heterogeneities. Here we harness the complementary strengths of both techniques by processing Ps-P delay times derived from teleseismic receiver functions in a tomographic S wave inversion. Using our inversion technique, we produce the first tomographic crustal velocity model beneath Cleveland Volcano, identifying a vertically extensive high V-P/V-S anomaly beneath the volcano that likely signifies a middle-to-lower crustal magma reservoir. The observation is the first of its kind in the central Aleutians, illustrating the potential of our technique to advance our understanding of crustal magmatic systems without broad seismic networks or distributed local seismicity.

The occurrence of intermediate depth seismicity (70-300 km) is commonly attributed to the dehydration of hydrous phases within the downgoing oceanic plate. While some water is incorporated into the oceanic crust at formation, a significant amount of water is introduced into the plate immediately before subduction along outer-rise faults. These faults have been shown to extend to depths of over 30 km and can channel water to depths of 20 km or more beneath the seafloor. However, the amount of water introduced into the oceanic mantle lithosphere, and the role of that water in the formation of intermediate depth seismicity, has been the topic of ongoing research. Here we compile evidence from areas where the subducted oceanic crust is likely thicker than the penetration depth of water into the downgoing plate. These regions comprise aseismic plateaus and ridges (hot spot tracks) that can be compared directly to adjacent segments of the oceanic plate where oceanic crust of normal thickness is subducted. Regions with thick oceanic crust show little to no seismicity at intermediate depths, whereas adjacent regions with normal oceanic crust (similar to 6-8 km thick) have significant seismicity at similar depths and distances from the trench. We hypothesize that intermediate depth earthquakes observed in regions with thinner oceanic crust are caused by mantle dehydration reactions that are not possible in regions where the oceanic mantle was never hydrated because the thickness of the oceanic crust exceeded the penetration depth of water into the plate. We compare our observations to phase diagrams of hydrous basalt and hydrated depleted peridotite to determine pressures and temperatures that would be consistent with our observations. These can provide valuable constraints, not only on the degree of hydration and dehydration in the downgoing plate, but also as ground-truth for thermal models of these regions, all of which have complex, three-dimensional, time-variant subduction geometries and thermal histories.

Shear-wave splitting observations of SKS and SKKS phases have been used widely to map azimuthal anisotropy, as caused by the occurrence of olivine, to constrain the dominant directions of upper mantle deformation. While SK(K)S splitting measurements at individual seismic stations are often averaged before interpretation, it is useful to consider additional information, for example, based on the variation of splitting parameters with azimuth due to the non-vertical incidence of core-phases. These constraints in theory enable a differentiation between various types of olivine and may allow us to infer otherwise poorly known upper mantle parameters such as stress, temperature, and water content. In this study, we predict the azimuthal variation of splitting parameters for A-, C-, and E-type olivine fabrics and match them with observations from the High Lava Plains, Northwestern Basin and Range, and Western Yellowstone Snake River Plain in the Pacific Northwest US. This helps to constrain the amount of water in the upper mantle in the back-arc of the Cascadia subduction zone, known for its consistent E-W oriented seismic anisotropy, and particularly large splitting delay times. Our investigation renders the C-type olivine mechanism improbable for this location; A- and E-type fabrics match the observations, although differentiating between them is difficult. However, the agreement of the amplitude of backazimuthal variation of the fast orientation, plus the potential to explain large splitting delay times, suggest the occurrence of E-type olivine, and thus the likely presence of a moderately hydrated upper mantle beneath Cascadia's back-arc.

Most flat-slab subduction regions are marked by an absence of arc volcanism, which is consistent with closure of the hot mantle wedge as the subducting plate flattens below the continent. Farther inland, low surface heat flow is observed, which is generally attributed to cooling of the continent by the underlying flat slab. However, modern flat slabs have only been in place for <20 Ma, and it is unclear whether there has been sufficient time for cooling to occur. We use numerical models to assess temporal variations in continental thermal structure during flat-slab subduction. Our models show that the flat slab leads to continental cooling on timescales of tens of millions of years. Cool slab temperatures must diffuse through the continental lithosphere, resulting in a delay between slab emplacement and surface cooling. Therefore, the timescales primarily depend on the flat-slab depth with shallower slabs resulting in shorter timescales. The magnitude of cooling increases for a shallow or long-lived flat slab, old subducting plate, and fast convergence rates. For regions with flat slabs at 45-70 km depth (e.g., Mexico and Peru), shallow continental cooling initiates 5-10 Ma after slab emplacement, and low surface heat flow in these regions is largely explained by the presence of the flat slab. However, for the Pampean region in Chile, with an similar to 100-km-deep slab, our models predict that conductive cooling has not yet affected the surface heat flow. The low heat flow observed requires additional processes such as advective cooling from the infiltration of fluids released through dehydration of the flat slab.

The eastern margin of North America has been shaped by a series of tectonic events including the Paleozoic Appalachian Orogeny and the breakup of Pangea during the Mesozoic. For the past similar to 200 Ma, eastern North America has been a passive continental margin; however, there is evidence in the Central Appalachian Mountains for post-rifting modification of lithospheric structure. This evidence includes two co-located pulses of magmatism that post-date the rifting event (at 152 and 47 Ma) along with low seismic velocities, high seismic attenuation, and high electrical conductivity in the upper mantle. Here, we synthesize and evaluate constraints on the lithospheric evolution of the Central Appalachian Mountains. These include tomographic imaging of seismic velocities, seismic and electrical conductivity imaging along the Mid-Atlantic Geophysical Integrative Collaboration array, gravity and heat flow measurements, geochemical and petrological examination of Jurassic and Eocene magmatic rocks, and estimates of erosion rates from geomorphological data. We discuss and evaluate a set of possible mechanisms for lithospheric loss and intraplate volcanism beneath the region. Taken together, recent observations provide compelling evidence for lithospheric loss beneath the Central Appalachians; while they cannot uniquely identify the processes associated with this loss, they narrow the range of plausible models, with important implications for our understanding of intraplate volcanism and the evolution of continental lithosphere. Our preferred models invoke a combination of (perhaps episodic) lithospheric loss via Rayleigh-Taylor instabilities and subsequent small-scale mantle flow in combination with shear-driven upwelling that maintains the region of thin lithosphere and causes partial melting in the asthenosphere.