An asymptotic model for time-harmonic motion in fully-coupled linear thermoelastic orthorhombic materials is presented. The asymptotic approach takes advantage of the observation that the parameter expressing departure from the purely adiabatic regime is extremely small in practice. Consequently, the leading order bulk response turns out to be non-dissipative, and is governed by the usual equations of elastodynamics with adiabatic material constants. In the case of isothermal stress-free boundary conditions, it is shown that thermoelastic interaction is dominated by a thermoelastic boundary layer. Hence, effective boundary conditions may be constructed, which duly account for the influence of this boundary layer and successfully describe dispersion and dissipation of surface waves to leading order. As an illustration, in the special case of an isotropic half-space with free isothermal boundary conditions, we recover the asymptotic results by Chadwick and Windle (1964). Numerical comparison of the dispersion curves for surface waves in an orthorhombic half-space shows excellent agreement between the exact fully-coupled thermoelastic problem and the corresponding quasi-adiabatic approximation, even for relatively large wavenumbers.