International audienceWe determine the statistical distribution of the co-polarized phase difference of fields scattered from a stack of two two-dimensional rough interfaces in the incidence plane. The electromagnetic fields are represented by Rayleigh expansions and a perturbation method is used to solve the boundary value problem and to determine the first-order scattering amplitudes. For slightly rough interfaces with infinite length and Gaussian height distributions, we show that the probability density function is only a function of two parameters. For a sand layer on a granite surface in backscattering configurations, we study the influence of the incidence angle, the layer thickness, the cross-spectral density and the wave frequency upon both parameters of the probability law