We develop a theoretical model to study (dense) two-dimensional gravity current flow in a laterally extensive porous medium experiencing leakage through a discrete fissure situated along this boundary at some finite distance from the injection point. Our model, which derives from the depth-averaged mass and buoyancy equations in conjunction with Darcy's law, considers dispersive mixing between the gravity current and the surrounding ambient by allowing two different gravity current phases. Thus do we define a bulk phase consisting of fluid whose density is close to that of the source fluid and a dispersed phase consisting of fluid whose density is close to that of the ambient. We characterize the degree of dispersion by estimating, as a function of time, the buoyancy of the dispersed phase and the separation distance between the bulk nose and the dispersed nose. On this basis, it can be shown that the amount of dispersion depends on the flow conditions upstream of the fissure, the fissure permeability and the vertical and horizontal extents of the fissure. We also show that dispersion is larger when the gravity current propagates along an inclined barrier rather than along a horizontal barrier. Model predictions are fitted against numerical simulations. The simulations in question are performed using COMSOL and consider different inclination angles and fissure and upstream flow conditions. Our study is motivated by processes related to underground $\mathrm {H}_2$ storage e.g. an irrecoverable loss of $\mathrm {H}_2$ when it is injected into the cushion gas saturating an otherwise depleted natural gas reservoir.