In a subsoil bioremediation intervention air or oxygen is injected in the polluted region and then a model for unsaturated porous media it is required, based on the theory of the dynamics of multiphase fluids in porous media. In order to optmize the costs of the intervention it is useful to consider the gas as compressible and this fact introduces nonlinearity in the mathematical model. The physical problem is described by a system of equations and the unknowns are: pollutant; bacteria concentration; oxygen saturation and oxygen pressure.

We investigate the long-time properties of the two-dimensional inviscid Boussinesq equations near a stably stratified Couette flow, for an initial Gevrey perturbation of size ?. Under the classical Miles-Howard stability condition on the Richardson number, we prove that the system experiences a shear-buoyancy instability: the density variation and velocity undergo an O(t-1/2) inviscid damping while the vorticity and density gradient grow as O(t1/2). The result holds at least until the natural, nonlinear timescale t??-2.