  We study the small-mass asymptotic behavior of so-called angular integrals,
appearing in phase-space calculations in perturbative quantum field theory. For
this purpose we utilize the strategy of expansion by regions, which is a
universal method both for multiloop Feynman integrals and various parametric
integrals. To apply the technique to angular integrals, we convert them into
suitable parametric integral representations, which are accessible to existing
automation tools. We use the the code \texttt{asy.m} to reveal regions
contributing to the asymptotic expansion of angular integrals. To evaluate the
contributions of these regions in an epsilon expansion we apply the method of
Mellin-Barnes representation. Our approach is checked against existing results
on angular integrals revealing a connection between contributing regions and
angular integrals constructed from an algebraic decomposition. We explicitly
calculate the previously unknown asymptotics for angular integrals with three
and four denominators and formulate a conjecture for the leading asymptotics
and the pole part for a general number of denominators and masses.
