Vertical mixing of the nocturnal stable boundary layer (SBL) over a complex land surface is investigated for a range of stabilities, using a decoupling index (0 < D rb < 1) based on the 2-50 m bulk gradient of the ubiquitous natural trace gas radon-222. The relationship between D rb and the bulk Richardson number (R ib ) exhibits three broad regions: (1) a well-mixed region (D rb ≈ 0.05) in weakly stable conditions (R ib < 0.03); (2) a steeply increasing region (0.05 < D rb < 0.9) for "transitional" stabilities (0.03 < R ib < 1); and (3) a decoupled region (D rb ≈ 0.9-1.0) in very stable conditions (R ib > 1). D rb exhibits a large variability within individual R ib bins, however, due to a range of competing processes influencing bulk mixing under different conditions. To explore these processes in R ib -D rb space, we perform a bivariate analysis of the bulk thermodynamic gradients, various indicators of external influences, and key turbulence quantities at 10 and 50 m. Strong and consistent patterns are found, and five distinct regions in R ib -D rb space are identified and associated with archetypal stable boundary-layer regimes. Results demonstrate that the introduction of a scalar decoupling index yields valuable information about turbulent mixing in the SBL that cannot be gained directly from a single bulk thermodynamic stability parameter. A significant part of the high variability observed in turbulence statistics during very stable conditions is attributable to changes in the degree of decoupling of the SBL from the residual layer above. When examined in R ib -D rb space, it is seen that very different turbulence regimes can occur for the same value of R ib , depending on the particular combination of values for the bulk temperature gradient and wind shear, together with external factors. Extremely low turbulent variances and fluxes are found at 50 m height when R ib > 1 and D rb ≈ 1 (fully decoupled). These "quiescent" cases tend to occur when geostrophic forcing is very weak and subsidence is present, but are not associated with the largest bulk temperature gradients. Humidity and net radiation data indicate the presence of low cloud, patchy fog or dew, any of which may aid decoupling in these cases by preventing temperature gradients from increasing sufficiently to favour gravity wave activity. The largest temperature gradients in our dataset are actually associated with smaller values of the decoupling index (D rb < 0.7), indicating the presence of mixing. Strong evidence is seen from enhanced turbulence levels, fluxes and submeso activity at 50 m, as well as high temperature variances and heat flux intermittencies at 10 m, suggesting this region of the R ib -D rb distribution can be identified as a top-down mixing regime. This may indicate an important role for gravity waves and other wave-like phenomena in providing the energy required for sporadic mixing at this complex terrain site.