The dynamics in three-dimensional (3D) billiards leads, using a Poincaré section, to a four-dimensional map, which is challenging to visualize. By means of the recently introduced 3D phase-space slices, an intuitive representation of the organization of the mixed phase space with regular and chaotic dynamics is obtained. Of particular interest for applications are constraints to classical transport between different regions of phase space which manifest in the statistics of Poincaré recurrence times. For a 3D paraboloid billiard we observe a slow power-law decay caused by long-trapped trajectories, which we analyze in phase space and in frequency space. Consistent with previous results for 4D maps, we find that (i) trapping takes place close to regular structures outside the Arnold web, (ii) trapping is not due to a generalized island-around-island hierarchy, and (iii) the dynamics of sticky orbits is governed by resonance channels which extend far into the chaotic sea. We find clear signatures of partial transport barriers. Moreover, we visualize the geometry of stochastic layers in resonance channels explored by sticky orbits.