Deployment of temporary seismic stations after the 2011 Mineral, Virginia (USA), earthquake produced a well-recorded aftershock sequence. The majority of aftershocks are in a tabular cluster that delineates the previously unknown Quail fault zone. Quail fault zone aftershocks range from ~3 to 8 km in depth and are in a 1-km-thick zone striking ~036° and dipping ~50°SE, consistent with a 028°, 50°SE main-shock nodal plane having mostly reverse slip. This cluster extends ~10 km along strike. The Quail fault zone projects to the surface in gneiss of the Ordovician Chopawamsic Formation just southeast of the Ordovician–Silurian Ellisville Granodiorite pluton tail. The following three clusters of shallow (<3 km) aftershocks illuminate other faults. (1) An elongate cluster of early aftershocks, ~10 km east of the Quail fault zone, extends 8 km from Fredericks Hall, strikes ~035°–039°, and appears to be roughly vertical. The Fredericks Hall fault may be a strand or splay of the older Lakeside fault zone, which to the south spans a width of several kilometers. (2) A cluster of later aftershocks ~3 km northeast of Cuckoo delineates a fault near the eastern contact of the Ordovician Quantico Formation. (3) An elongate cluster of late aftershocks ~1 km northwest of the Quail fault zone aftershock cluster delineates the northwest fault (described herein), which is temporally distinct, dips more steeply, and has a more northeastward strike. Some aftershock-illuminated faults coincide with preexisting units or structures evident from radiometric anomalies, suggesting tectonic inheritance or reactivation.

Characterizing geologic features associated with major earthquakes provides insights into mechanisms contributing to fault slip and assists evaluation of seismic hazard. We use high-resolution airborne geophysical data combined with ground sample measurements to image subsurface geologic features associated with the 2011 moment magnitude (M w) 5.8 central Virginia (USA) intraplate earthquake and its aftershocks. Geologic mapping and magnetic data analyses suggest that the earthquake occurred near a complex juncture of geologic contacts. These contacts also intersect a >60-km-long linear gravity gradient. Distal aftershocks occurred in tight, ~1-km-wide clusters near other obliquely oriented contacts that intersect gravity gradients, in contrast to more linearly distributed seismicity observed at other seismic zones. These data and corresponding models suggest that local density contrasts (manifested as gravity gradients) modified the nearby stress regime in a manner favoring failure. However, along those gradients seismic activity is localized near structural complexities, suggesting a significant contribution from variations in associated rock characteristics such as rheological weakness and/or rock permeability, which may be enhanced in those areas. Regional magnetic data show a broader bend in geologic structures within the Central Virginia seismic zone, suggesting that seismic activity may also be enhanced in other nearby areas with locally increased rheological weaknesses and/or rock permeability. In contrast, away from the M w 5.8 epicenter, geophysical lineaments are nearly continuous for tens of kilometers, especially toward the northeast. Continuity of associated geologic structures probably contributed to efficient propagation of seismic energy in that direction, consistent with moderate to high levels of damage from Louisa County to Washington, D.C., and neighboring communities.

We use magnetic susceptibility and remanent magnetization measurements of the Eyreville and Cape Charles cores in combination with new and previously collected magnetic field data in order to constrain structural features within the inner basin of the Chesapeake Bay impact structure. The Eyreville core shows the first evidence of several-hundred-meter-thick basement-derived megablocks that have been transported possibly kilometers from their pre-impact location. The magnetic anomaly map of the structure exhibits numerous short-wavelength (<2 km) variations that indicate the presence of magnetic sources within the crater fill. With core magnetic properties and seismic reflection and refraction results as constraints, forward models of the magnetic field show that these sources may represent basement-derived megablocks that are a few hundred meters thick or melt bodies that are a few dozen meters thick. Larger-scale magnetic field properties suggest that these bodies overlie deeper, pre-impact basement contacts between materials with different magnetic properties such as gneiss and schist or gneiss and granite. The distribution of the short-wavelength magnetic anomalies in combination with observations of small-scale (1–2 mGal) gravity field variations suggest that basement-derived megablocks are preferentially distributed on the eastern side of the inner crater, not far from the Eyreville core, at depths of around 1–2 km. A scenario where additional basement-derived blocks between 2 and 3 km depth are distributed throughout the inner basin—and are composed of more magnetic materials, such as granite and schist, toward the east over a large-scale magnetic anomaly high and less magnetic materials, such as gneiss, toward the west where the magnetic anomaly is lower—provides a good model fit to the observed magnetic anomalies in a manner that is consistent with both gravity and seismic-refraction data.