The National Science Foundation has announced a plan to establish a Deep Underground Science and Engineering Laboratory (DUSEL) for interdisciplinary research in physics, geosciences, biosciences and engineering. The proposed laboratory will extend to a depth of about 2200 meters and will consist of research facilities for long term study. To date, eight sites in North America have been proposed to host DUSEL. One of these sites, known as Kimballton, is located near Butt Mountain in Giles County in southwestern Virginia. Two seismic lines were acquired along the top of Butt Mountain in June of 2004 to support the ongoing integrated site characterization effort by the Kimballton Science Team. Both lines, approximately 3 km in length, are standard multifold seismic reflection data aimed at imaging faults, thrust sheets, and repeated sections of Paleozoic rocks in the vicinity of the proposed Kimballton site. Crooked line geometry, irregular geophone spacing, ground roll, and poor impedance contrasts between juxtapositioned rock units were challenges in processing the data. Non-standard processing techniques included the use of travel time tomography to accurately constrain near surface velocities, the use of 2D median filters to remove ground roll, and stacking only offsets exceeding 500 m. Interpretation of seismic data supports a triplicated stratigraphic section caused by the stacking of the the St. Clair and Narrows thrust sheets. The St. Clair and Narrows faults are interpreted as shear zones within ductile units of the Martinsburg Formation. 3D travel time tomography was used to build a near surface velocity model of Lines 1 and 2 for the purposes of imaging near surface structure and constraining the extent of topographic lineaments, which are interpreted as bedrock joint systems. Interpretation of the velocity models suggests that the broadly folded strata of the Butt Mountain synclinorium dip gently to the east along the hinge surface. The surface extrapolation of the Lookout Rock fault and the intersection of topographic lineaments with the seismic lines are expressed as low velocity zones that extend to depths of 150 m. This may be related to accelerated weathering along jointed rock surfaces. Results of this study have already been incorporated into the NSF proposal submitted by the Kimballton Science Team (http://www.phys.vt.edu/~kimballton/s2p/b2.pdf). / Master of Science

Dense array seismic techniques can be applied to multiple types of seismic data to understand regional tectonic processes via analysis of crustal velocity structure, imaging reflection surfaces, and calculating high-resolution hypocenter locations. The two regions presented here include an intraplate seismogenic fault zone in Virginia and a steep cratonic margin in eastern Oregon and Idaho. The intraplate seismicity study in Virginia consisted of using 201 short-period vertical-component seismographs, which recorded events as low as magnitude -2 during a period of 12 days. Dense array analysis revealed almost no variation in the seismic velocity within the hypocentral zone, indicating that the aftershock zone is confined to a single crystalline-rock terrane. The 1-2 km wide cloud of hypocenters is characterized by a 29° strike and 53° dip consistent with the focal mechanism of the main shock. A 5° bend along strike and a shallower dip angle below 6 km points toward a more complex concave shaped fault zone. The seismic study in Idaho and Oregon was centered on the inversion of controlled-source wide-angle reflection and refraction seismic P- and S-wave traveltimes to determine a seismic velocity model of the crust beneath this part of the U.S. Cordillera. We imaged a narrow, steep velocity boundary within the crust that juxtaposes the Blue Mountains accreted terranes and the North American craton at the western Idaho shear zone. We found a 7 km offset in Moho depth, separating crust with different seismic velocities and Poisson's ratios. The crust beneath the Blue Mountains terranes is consistent with an intermediate lithology dominated by diorite. In the lower crust there is evidence of magmatic underplating which is consistent with the location of the feeder system of the Columbia River Basalts. The cratonic crust east of the WISZ is thicker and characterized by a felsic composition dominated by granite through most of the crust, with an intermediate composition layer in the lower crust. This sharp lithologic and rheologic boundary strongly influenced subsequent deformation and magmatic events in the region. / Ph. D.

Travel times from a refraction and wide-angle reflection seismic survey across the Coast Plutonic Complex and Stikine terrane of British Columbia were inverted to derive two dimensional P and S-wave seismic velocity models of the crust and uppermost mantle. A felsic upper crust and a felsic to intermediate middle crust are observed in both the batholith complex and the accreted Stikine island arc terrane. The P and S wave models demonstrate a high-velocity (P 7.0 km/s, S 3.8 km/s) layer in the lower crust beneath the youngest (late Cretaceous to Eocene) portion of the continental arc complex. In contrast, the lower crust under the Stikine terrane has lower velocities consistent with amphibolite or other hydrated mafic rocks. The Moho is at ~35 km depth under the Stikine terrane, deepens to ~38 km beneath the youngest portion of the arc, then shallows towards the coast. The high velocity zone under the younger portion of the Coast Plutonic Complex has a 1.81 Vp/Vs ratio and is interpreted to have a bulk composition of mafic garnet granulite. This garnet granulite and large volumes of granodiorite-dominated melt were created by arc dehydration melting of amphibolite (or hydrated gabbro) in the pre-existing lower crust Reverse time migration method was applied to image aftershocks recorded by a dense array deployed after the 2011 Virginia earthquake. Events as tiny as magnitude -2 were successfully imaged as point sources. The propagation of energy release as a function of time and space was observed for events larger than magnitude 2.5. Spatial resolution of the images was ~200 m, which synthetic data tests show was primarily limited by the temporal sampling rate. Improved temporal and spatial sampling could produce images with sharper resolution. / Ph. D.