The lens of freshwater in a tropical island––2d withdrawal

Steady solutions are found for the shape of the interface between salt and freshwater beneath a tropical island. A Green’s function approach is used and proves to be much more robust than spectral methods. Computations of the surface shape during withdrawal through a line sink are described and the maximum pumping rate before critical breakthrough of the salt water is found for a range of different sink locations. The effects of density ratio, withdrawal depth and horizontal location are also calculated and found to be significant. Decreasing the density ratio between the layers, or increasing withdrawal height is shown to make a significant difference to the critical flow rate. Further, it is shown that for realistic densities, the island width can be a factor for quite substantial widths, whereas for higher, unrealistic density differences an increase in width quickly becomes irrelevant to the interface shape.