  We prove a theorem describing the limiting fine-scale statistics of orbits of
a point in hyperbolic space under the action of a discrete subgroup. Similar
results have been proved only in the lattice case, with two recent
infinite-volume exceptions by Zhang for Apollonian circle packings and certain
Schottky groups. Our results hold for general Zariski dense, non-elementary,
geometrically finite subgroups in any dimension. Unlike in the lattice case,
orbits of geometrically finite subgroups do not necessarily equidistribute on
the whole boundary of hyperbolic space. But rather, they may equidistribute on
a fractal subset. Understanding the behaviour of these orbits near the boundary
is central to Patterson-Sullivan theory and much further work. Our theorem
characterizes the higher order spatial statistics and thus addresses a very
natural question. As a motivating example our work applies to sphere packings
(in any dimension) which are invariant under the action of such discrete
subgroups. At the end of the paper we show how this statistical
characterization can be used to prove convergence of moments and to write down
the limiting formula for the two-point correlation function and nearest
neighbor distribution. Moreover we establish an formula for the 2 dimensional
limiting gap distribution (and cumulative gap distribution) which was not known
previously even in the lattice case.
