Global Positioning System (GPS) measurements to study regional deformation were initiated in northern Cascadia in the late 1980s and early 1990s. On the basis of a decade of GPS data, we derive a crustal velocity field for NW Washington-SW British Columbia. The permanent and campaign GPS velocities are defined with respect to North America in the ITRF2000 reference frame. Velocity uncertainties are estimated using a model of time series noise spectra. This new velocity field is the basis for interpretation of the tectonics of the northern Cascadia subduction system. GPS velocities are interpreted in terms of interseismic loading of the megathrust using different coupling models. Our data confirm that the upper part of the megathrust is nearly fully locked. An exponential model for the downdip transition zone gives slightly better agreement with the data compared to the common linear transition. The landward decrease of forearc strain loading is smaller than predicted by any of the current subduction interseismic models. This could be a consequence of a small (0--3 mm/yr) long-term motion of the southern Vancouver Island forearc, with respect to North America, or of a concentration of interseismic strain across the elastically weaker Cascadia volcanic arc. In northern Vancouver Island, our velocity field supports the existence of an independent Explorer microplate currently underthrusting underneath North America, at least up to Brooks Peninsula. Further north, GPS velocities indicate transient and/or permanent deformation of northernmost Vancouver Island related to the interaction with the Explorer microplate and possibly with the Queen Charlotte transform margin.