[1] A method for 2.5D viscoacoustic waveform tomography that can be applied to generate 2D models of velocity and attenuation from inversion of refraction waveforms on land seismic reflection data acquired along crooked roads is developed. It is particularly useful for typical crustal reflection surveys. First-arrival travel time tomography is applied using a 3D method, but with constraints on the intermediate 3D velocity model; the result is the starting model for the next step. A 2.5D frequency-domain full-waveform inversion stage parameterizes 3D geometry in the seismic source and receiver arrays, with the assumption that the velocity and attenuation models are homogeneous in the out-of-plane direction. This approach results in superior results compared to a strictly 2D approach when the acquisition line is crooked, with a moderate increase in computational cost. A case study using data acquired in the Nechako Basin in south-central British Columbia, Canada, exemplifies and validates the procedure. The velocity model derived from 2.5D waveform tomography is compared with that from a previous study in which 2D waveform tomography was applied to the same data set and with results from 3D travel time tomography. The resolution and accuracy of the velocity model from 2.5D waveform tomography are demonstrated to be greater than those from travel time or 2D waveform tomography. A model of viscoacoustic attenuation, which was not possible in the 2D case, is also generated. These models are interpreted jointly to highlight features of geological interest, such as a sedimentary basin, basement rocks, and faults, from surface to about 3 km depth.