Gravity and magnetic anomalies used to delineate geologic features associated with earthquakes and aftershocks in the central Virginia seismic zone

Estimating seismic hazard in intraplate environments can be challenging partly because events are relatively rare and associated data thus limited. Additionally, in areas such as the central Virginia seismic zone, numerous pre-existing faults may or may not be candidates for modern tectonic activity, and other faults may not have been mapped. It is thus important to determine whether or not specific geologic features are associated with seismic events. Geophysical and geologic data collected in response to the Mw5.8 August 23, 2011 central Virginia earthquake provide excellent tools for this purpose. Portable seismographs deployed within days of the main shock showed a series of aftershocks mostly occurring at depths of 3-8 km along a southeast-dipping tabular zone ~10 km long, interpreted as the causative fault or fault zone. These instruments also recorded shallow (< 4 km) aftershocks clustered in several areas at distances of ~2-15 km from the main fault zone. We use new airborne geophysical surveys (gravity, magnetics, radiometrics, and LiDAR) to delineate the distribution of various surface and subsurface geologic features of interest in areas where the earthquake and aftershocks took place. The main (causative fault) aftershock cluster coincides with a linear, NE-trending gravity gradient (~ 2 mgal/km) that extends over 20 km in either direction from the Mw5.8 epicenter. Gravity modeling incorporating seismic estimates of Moho variations suggests the presence of a shallow low-density body overlying the main aftershock cluster, placing it within the upper 2-4 km of the main-fault hanging wall. The gravity, magnetic, and radiometric data also show a bend in generally NE-SW orientation of anomalies close to the Mw5.8 epicenter. Most shallow aftershock clusters occur near weaker short-wavelength gravity gradients of one to several km length. In several cases these gradients correspond to geologic contacts mapped at the surface. Along the gravity gradients, the aftershocks appear to cluster near areas with cross-cutting geologic features such as Jurassic diabase dikes. These associations suggest that local variations in rock density and/or rheology may have contributed to modifications of local stress regimes in a manner encouraging localized seismicity associated with the Mw5.8 event and its aftershocks. Such associations are comparable to results of previous studies recognizing correspondences between seismicity and features such as intrusive bodies and failed rifts in the New Madrid seismic zone and elsewhere. To explore whether similar correspondences may have occurred in the past, we use regional gravity and magnetic data to consider possible relations between historical earthquakes and comparable geologic features elsewhere in the central Virginia seismic zone.