<p>Onset of thermal and haline convection was studied separately by Lapwood (1948) and Wooding (1956) in theoretical models using analytical methods. They established that the buoyancy force caused by difference in temperature (&#916;<em>T</em>) or concentration (&#916;<em>c</em>) can induce natural convection depending on the model properties (e.g. geometry, permeability, etc.). In the course of further numerical simulations, the thermal (<em>Ra<sub>T</sub></em>) and the haline Rayleigh number (<em>Ra<sub>H</sub></em>) proved itself useful to characterise the type, the intensity and the form of the natural convection (e.g. Diersch and Kolditz, 2002). The main purpose of our study was to examine numerically the combined effect of temperature- and salinity-driven natural convection in a two-dimensional homogeneous porous medium.</p><p>Two-dimensional finite element base model was set up in agreement with the Elder problem (Wooding, 1956) in order to verify the numerical calculation. First, it was established that (1) the critical Rayleigh numbers are mathematically equivalent in the two separated cases (<em>Ra</em><sub><em>Tcr</em></sub>=<em>Ra<sub>Hcr</sub></em>=4&#960;<sup>2</sup>), and (2) time-dependent thermal or haline natural convections evolve, when the Rayleigh number lies within the range of 300&#8211;600. Numerical simulations were accomplished to investigate the interaction of the temperature- and salinity-driven natural convection. Non-dimensional thermal expansion and haline concentration were increased from &#945;&#916;<em>T</em>=0.01 to 1 and from &#946;&#916;<em>c</em>=10<sup>-5</sup> to 10<sup>-3</sup>, respectively, while the variation of the Darcy flux, the temperature, the concentration, the Nusselt and the Sherwood numbers was computed. The main points of this study were that (1) how the onset of the thermohaline convection is facilitated by the positive interaction of the thermal and haline effects (<em>Ra<sub>THcr</sub></em>); (2) under what conditions time-dependent flow evolves in the theoretical models; (3) whether a new non-dimensional number can be defined instead of the two separated Rayleigh numbers in order to characterise the behaviour of the thermonaline convection. These simulations draw attention to the importance of understanding the physical background of thermohaline convection, for instance, at the margin of confined and unconfined carbonate systems (e.g. Buda Thermal Karst), or in the case of groundwater flow induced by water pumping/injection of deep geothermal power plants.</p><p>The project was supported by the &#218;NKP-19-3 and &#218;NKP-19-4 New National Excellence Program of the Ministry for Innovation and Technology, the Hungarian Research Fund (K 129279) and the J&#225;nos Bolyai Scholarship of the Hungarian Academy of Science. This research is a part of a project that has received funding from the European Union&#8217;s Horizon 2020 research and innovation program under grant agreement No 810980.</p><p>References:</p><p>Diersch, H.-J.G., Kolditz, O. (2002), Variable-density flow and transport in porous media: approaches and challenges. Advances in Water Resources, 25, 899-944.</p><p>Lapwood, E.R. (1948), Convection of a fluid in a porous medium. Mathematical Proceedings of the Cambridge Philosophical Society, 44, 508-521.</p><p>Wooding, R.A. (1956), Steady state free thermal convection of liquid in a saturated permeable medium. Journal of Fluid Mechanics, 2, 273-285.</p>