Development and sensitivity analysis of a model for assessing stratification and safety of Lake Nyos during artificial degassing

Schmid, Martin ; Lorke, Andreas ; Wüest, Alfred ; Halbwachs, Michel ; Tanyileke, Gregory

In: Ocean Dynamics, 2003, vol. 53, no. 3, p. 288-301

Zum persönliche Liste hinzufügen
    Summary
    To prevent the recurrence of a disastrous eruption of carbon dioxide (CO2) from Lake Nyos, a degassing plan has been set up for the lake. Since there are concerns that the degassing of the lake may reduce the stability of the density stratification, there is an urgent need for a simulation tool to predict the evolution of the lake stratification in different scenarios. This paper describes the development of a numerical model to predict the CO2 and dissolved solids concentrations, and the temperature structure as well as the stability of the water column of Lake Nyos. The model is tested with profiles of CO2 concentrations and temperature taken in the years 1986 to 1996. It reproduces well the general mixing patterns observed in the lake. However, the intensity of the mixing tends to be overestimated in the epilimnion and underestimated in the monimolimnion. The overestimation of the mixing depth in the epilimnion is caused either by the parameterization of the k-epsilon model, or by the uncertainty in the calculation of the surface heat fluxes. The simulated mixing depth is highly sensitive to the surface heat fluxes, and errors in the mixing depth propagate from one year to the following. A precise simulation of the mixolimnion deepening therefore requires high accuracy in the meteorological forcing and the parameterization of the heat fluxes. Neither the meteorological data nor the formulae for the calculation of the heat fluxes are available with the necessary precision. Consequently, it will be indispensable to consider different forcing scenarios in the safety analysis in order to obtain robust boundary conditions for safe degassing. The input of temperature and CO2 to the lake bottom can be adequately simulated for the years 1986 to 1996 with a constant sublacustrine source of 18 l s−1 with a CO2 concentration of 0.395 mol l−1 and a temperature of 26 °C. The results of this study indicate that the model needs to be calibrated with more detailed field data before using it for its final purpose: the prediction of the stability and the safety of Lake Nyos during the degassing process