A 3D Tomographic Model of Asia Based on Pn and P Travel Times from GT Events

Increasingly, nuclear explosion monitoring is focusing on detection, location, and identification of small events recorded at regional distances. Because Earth structure is highly variable on regional scales, locating events accurately at these distances requires the use of region-specific models to provide accurate travel times. Improved results have been achieved with composites of 1D models and with approximate 3D models with simplified upper mantle structures, but both approaches introduce non-physical boundaries that are problematic for operational monitoring use. Ultimately, what is needed is a true, seamless 3D model of the Earth. Towards that goal, we have developed a 3D tomographic model of the P velocity of the crust and mantle for the Asian continent. Our model is derived by an iterative least squares travel time inversion of more than one million Pn and teleseismic P picks from some 35,000 events recorded at 4,000+ stations. We invert for P velocities from the top of the crust to the core mantle boundary, along with source and receiver static time terms to account for the effects of event mislocation and unaccounted for fine-scale structure near the receiver. Because large portions of the model are under-constrained, we apply spatially varying damping, which constrains the inversion to update the starting model only where good data coverage is available. Our starting crustal model is taken from the a priori crust and upper mantle model of Asia developed through National Nuclear Security Administration laboratory collaboration, which is based on various global and regional studies, and we substantially increase the damping in the crust to discourage changes from this model. Our starting mantle model is AK135. To simplify the inversion, we fix the depths of the major mantle discontinuities (Moho, 410 km, 660 km). 3D rays are calculated using an implementation of the Um and Thurber ray pseudo-bending approach, with full enforcement of Snell's Law in 3D at the major discontinuities. Due to the highly non-linear nature of our ray tracer, we are forced to substantially damp the inversion in order to converge on a reasonable model. We apply both horizontal and vertical regularization to produce smooth models with velocity feature scale lengths that are consistent with established conventions for mantle velocity structure. To investigate the importance of using true 3D rays for the inversion, as opposed to proxy rays through a reference model, we compare our model and ray paths with the model and ray paths resulting from inverting the same data set using rays traced through a 1D reference model. Finally, we validate the model by performing several inversions with random portions of the data set omitted and then testing the predictive capability of the model against those portions compared with AK135. We test the location performance of the model by relocating the GT events using our model and using AK135. Sandia is a multiprogram laboratory operated by Sandia Corporation, a Lockheed Martin Company, for the United States Department of Energy's National Nuclear Security Administration under Contract DE-AC04- 94AL85000.