High-resolution modelling of atmospheric dispersion of dense gas using TWODEE-2.1: application to the 1986 Lake Nyos limnic eruption

Abstract. Atmospheric dispersal of a gas denser than air can threat the environment and surrounding communities if the terrain and meteorological conditions favour its accumulation in topographic depressions, thereby reaching toxic concentration levels. Numerical modelling of atmospheric gas dispersion constitutes a useful tool for gas hazard assessment studies, essential for planning risk mitigation actions. In complex terrains, microscale winds and local orographic features can have a strong influence on the gas cloud behaviour, potentially leading to inaccurate results if not captured by coarser-scale modelling. We introduce a methodology for microscale wind field characterisation based on transfer functions that couple a mesoscale numerical weather prediction model with a microscale computational fluid dynamics (CFD) model for the atmospheric boundary layer. The resulting time-dependent high-resolution microscale wind field is used as input for a shallow-layer gas dispersal model (TWODEE-2.1) to simulate the time evolution of CO₂ gas concentration at different heights above the terrain. The strategy is applied to review simulations of the 1986 Lake Nyos event in Cameroon, where a huge CO₂ cloud released by a limnic eruption spread downslopes from the lake, suffocating thousands of people and animals across the Nyos and adjacent secondary valleys. Besides several new features introduced in the new version of the gas dispersal code (TWODEE-2.1), we have also implemented a novel impact criterion based on the percentage of human fatalities depending on CO₂ concentration and exposure time. New model results are quantitatively validated using the reported percentage of fatalities at several locations. The comparison with previous simulations that assumed coarser-scale steady winds and topography illustrates the importance of high-resolution modelling in complex terrains.