In seismic processing, velocity fields are commonly represented on finely sampled Cartesian grids. Attractive alternatives are unstructured grids such as meshes composed of triangles or tetrahedra. Meshes provide a space-filling framework that enables editing of velocity models while facilitating numerical tasks such as seismic modeling and inversion. In this paper, we introduce an automated process to generate meshes of subsurface velocity structures for highly resolved velocity fields without providing additional external constraints such as horizons and faults. Our analysis shows that these new meshes can represent both smooth and discontinuous velocity profiles accurately and with less computer memory than grids.