In recent years, sum–product estimates in Euclidean space and finite fields have received great attention. They can often be interpreted in terms of Erdős type incidence problems involving the distribution of distances, dot products, areas, and so on, which have been studied quite extensively by way of combinatorial and Fourier analytic techniques. We use both kinds of techniques to obtain sharp or near-sharp results on the distribution of volumes (as examples of d -linear homogeneous forms) determined by sufficiently large subsets of vector spaces over finite fields and the associated arithmetic expressions. Arithmetic–combinatorial techniques turn out to be optimal for dimension d ≥ 4 to this end, while for d = 3 they have failed to provide us with a result that follows from the analysis of exponential sums. To obtain the latter result we prove a relatively straightforward function version of an incidence results for points and planes previously established in [D. Hart, A. Iosevich, Sums and products in finite fields: An integral geometric viewpoint, in: Radon Transforms, Geometry, and Wavelets, Contemp. Math. 464 (2008); D. Hart, A. Iosevich, D. Koh, M. Rudnev, Averages over hyperplanes, sum–product theory in vector spaces over finite fields and the Erdős–Falconer distance conjecture, arXiv:math/0711.4427, preprint 2007]. More specifically, we prove that if E = A × ⋯ × A is a product set in F q d , d ≥ 4 , the d -dimensional vector space over a finite field F q , such that the size | E | of E exceeds q d 2 (i.e. the size of the generating set A exceeds q ) then the set of volumes of d -dimensional parallelepipeds determined by E covers F q . This result is sharp as can be seen by taking A = F p , a prime sub-field of its quadratic extension F q , with q = p 2 . For in three dimensions, however, we are able to establish the same result only if | E | ≳ q 15 8 (i.e., | A | ≥ C q 5 8 , for some C ; in fact, the q 15 8 bound can be justified for a slightly wider class of “Cartesian product-like” sets), and this uses Fourier methods. Yet we do prove a weaker near-optimal result in three dimensions: that the set of volumes generated by a product set E = A × A × A covers a positive proportion of F q if | E | > q 3 2 (so | A | > q ). Besides, without any assumptions on the structure of E , we show that in three dimensions the set of volumes covers a positive proportion of F q if | E | ≥ C q 2 , which is again sharp up to the constant C , as taking E to be a 2-plane through the origin shows.