As a part of the series of Source Physics Experiments (SPE) conducted on the Nevada National Security Site in southern Nevada, we have developed a local‐to‐regional scale seismic velocity model of the site and surrounding area. Accurate earth models are critical for modeling sources like the SPE to investigate the role of earth structure on the propagation and scattering of seismic waves. We combine seismic body waves, surface waves, and gravity data in a joint inversion procedure to solve for the optimal 3D seismic compressional and shear‐wave velocity structures and earthquake locations subject to model smoothness constraints. Earthquakes, which are relocated as part of the inversion, provide P‐ and S‐body‐wave absolute and differential travel times. Active source experiments in the region augment this dataset with P‐body‐wave absolute times and surface‐wave dispersion data. Dense ground‐based gravity observations and surface‐wave dispersion derived from ambient noise in the region fill in many areas where body‐wave data are sparse. In general, the top 1–2 km of the surface is relatively poorly sampled by the body waves alone. However, the addition of gravity and surface waves to the body‐wave dataset greatly enhances structural resolvability in the near surface. We discuss the methodology we developed for simultaneous inversion of these disparate data types and briefly describe results of the inversion in the context of previous work in the region.