It is a requirement of nuclear energy and research facilities to conduct continuous and comprehensive atmospheric monitoring in order to better forecast public or environmental exposure to routine or accidental releases of radioactive substances to the atmosphere. A key aspect of such monitoring programs is the assessment of the atmospheric mixing state (or “stability”). Whether these facilities are in dense urban areas, or surrounded by heavily vegetated exclusion zones, local roughness heterogeneity can hamper attempts to accurately categorise stability by conventional meteorological techniques. Based on an analysis of 8 months of hourly climatology and atmospheric radon observations from a 60 m tower at the IFIN-HH nuclear research facility (Bucharest, Romania), we develop and apply a continuous (i.e. not categorical) radon-based scheme for the classification of the nocturnal atmospheric stability state. We demonstrate the superior performance of the radon-based technique to Pasquill-Gifford or bulk Richardson number stability typing at this site where heterogeneous roughness elements reach to 15 m a.g.l. Under stable nocturnal conditions the Pasquill-Gifford scheme overestimates the atmosphere’s capacity to dilute pollutants with near-surface sources by 20% compared to the radon-based scheme. Under these conditions, near-surface wind speeds drop well below 1 m s-1 and nocturnal mixing depths vary from ~25 m to less than 10 m a.g.l. Climatological parameters are characterised by season and 4 arbitrarily-defined nocturnal stability categories. Benchmarks (based on 10/50/90th percentile distributions) of 30-60 m wind and temperature gradients are devised for each stability category for evaluation of model performance. Lastly, nocturnal radon-derived effective mixing depth estimates constrained by tower observations are used to better-constrain the seasonal variability in the Bucharest regional radon flux: 13 mBq m-2 s-1 (winter), 18 mBq m-2 s-1 (summer)
