The coda waves of 21 Parkfield earthquakes were recorded simultaneously at the surface and at a depth of 198 m. We characterize these codas by computing (1) the time dependence of the integral of squared particle velocity, (2) their frequency content, (3) time-frequency coda decay planes, (4) the frequency dependence of root mean square ellipticity, and (5) the difference in their normalized decay rates. These five measures reveal significant differences in the early portion of the uphole versus downhole codas. Our analysis solves for coda Q as an explicit function of both time and frequency. Average apparent coda Q for the 21 events observed both at the surface and downhole is Qa(f, &tgr;) =39(f/f0)0.43(&tgr;/&tgr;0)0.37, with the reference frequency and time defined at 1 Hz and 1 s, respectively. Coda Q exhibits a trend in its value over the 1-month period of the study. The trend correlates with a spatial variation in source locations and is not produced by temporal variations of the Q of the medium. The spectral ratio between the uphole signal and the downhole signal varies with time through the coda.

The early portion of the uphole coda experiences a pulse of low-frequency energy not seen in the downhole coda. To explain this observation, we propose a model of the coda that includes body wave to surface wave scattering near the site. In this model, the near-surface acts as both a filter and a scattering waveguide. The coda at the surface is the sum of a coda produced by scattering deep in the lithosphere and a coda produced by scattering in the near-surface. The scattering in the near-surface is made up of near-vertically traveling multiples and horizontally travelling surface waves. The guided waves are recorded in the surface coda but are absent from the coda recorded at depth. The time-frequency structure of the surface waves in the uphole coda is determined by the velocity, scattering, and attenuation characteristics of the site. By fitting the model to our data, we have found that the fractional conversion factor of body waves to surface waves travelling toward the receiver is 0.03 at the Vineyard Canyon site. The relatively rapid damping of these trapped waves can be modelled by assuming a quality factor of QT=50. The principles of this model are general enough to be of use at other locations. ¿ American Geophysical Union 1990