We report the analysis of the statistics of the phase fluctuations in the coda of earthquakes recorded during a temporary experiment deployed at Pinyon Flats Observatory, California. The observed distributions of the first, second and third derivatives of the phase in the seismic coda exhibit universal power-law decays whose exponents agree accurately with circular Gaussian statistics. The correlation function of the spatial phase derivative is measured and used to estimate the mean free path of Rayleigh waves.PACS numbers: 46.65.+g, 91.30.Ab, 46.40.Cd In the short-period band (> 1 Hz) , ballistic arrivals of seismic waves are often masked by scattered waves due to small-scale heterogeneities in the lithosphere. The scattered elastic waves form the pronounced tail of seismograms known as the seismic coda [1,2]. Even when scattering is prominent, it is still possible to define the phase of the seismic record by introducing the complex analytic signal ψ(t, r) = A(t, r)e iφ(t,r) , with A the amplitude and φ the phase. In the past, many studies have focused on the modeling of the mean field intensity I(t) = A(t) 2 [see 3, for review]. The goal of the present paper is to study the statistics of the phase field in the coda. In the coda, the measured displacements result from the superposition of many partial waves which have propagated along different paths between the source and the receiver. Each path consists of a sequence of scattering events that affect the phase of the corresponding partial wave in a random way. For narrow-band signals, the phase field can therefore be written as φ(t, r) = ωt + δφ(t, r), where ω is the central frequency, and δφ denotes the random fluctuations. The trivial cyclic phase ωt cancels when a spatial phase difference is considered between two neighbouring points. Spatially resolved measurements are facilitated by dense arrays of seismometers that have been set up occasionally. We note that the phase of coda waves has not been given much attention so far. The advantage of phase is that it is not affected by the earthquake magnitude, and that it contains pure information on scattering, not blurred by absorption effects. For the statistical analysis of amplitude and phase fluctuations of direct arrivals, we refer the reader to e.g. Zheng and Wu [4].