We present analytical expressions for the Fourier analogue of the cosmic microwave background three-point and four-point correlation functions, the spatial bispectrum and trispectrum, of the Ostriker-Vishniac effect in the linear and mildly nonlinear regime. Through this systematic study, we illustrate a technique to tackle the calculation of such statistics making use of the effects of its small-angle and vectorlike properties through the Limber approximation. Finally we discuss its configuration dependence and detectability in the context of Gaussian theories for the currently favored flat LambdaCDM cosmology.