A linearized baroclinic, spectral-in-time tidal inverse model has been developed for assimilation of surface currents from coast-based high-frequency (HF) radars. Representer functions obtained as a part of the generalized inverse solution show that for superinertial flows information from the surface velocity measurements propagates to depth along wave characteristics, allowing internal tidal flows to be mapped throughout the water column. Application of the inverse model to a 38 km ϫ 57 km domain off the mid-Oregon coast, where data from two HF radar systems are available, provides a uniquely detailed picture of spatial and temporal variability of the M 2 internal tide in a coastal environment. Most baroclinic signal contained in the data comes from outside the computational domain, and so data assimilation (DA) is used to restore baroclinic currents at the open boundary (OB). Experiments with synthetic data demonstrate that the choice of the error covariance for the OB condition affects model performance. A covariance consistent with assumed dynamics is obtained by nesting, using representers computed in a larger domain. Harmonic analysis of currents from HF radars and an acoustic Doppler profiler (ADP) mooring off Oregon for May-July 1998 reveals substantial intermittence of the internal tide, both in amplitude and phase. Assimilation of the surface current measurements captures the temporal variability and improves the ADP/solution rms difference. Despite significant temporal variability, persistent features are found for the studied period; for instance, the dominant direction of baroclinic wave phase and energy propagation is always from the northwest. At the surface, baroclinic surface tidal currents (deviations from the depth-averaged current) can be 10 cm s Ϫ1 , 2 times as large as the depth-averaged current. Barotropic-to-baroclinic energy conversion is generally weak within the model domain over the shelf but reaches 5 mW m Ϫ2 at times over the slopes of Stonewall Bank.