An unsaturated flow and non-passive transport model for water-soluble organic compounds has been implemented in cylindrical coordinates, with a top boundary condition that accounts for different zones of the surface that can be under infiltration or volatilization independently. We simulated two-dimensional infiltration of aqueous mixtures of methanol from a disk source, its redistribution and volatilization in both homogeneous and heterogeneous soils. Simulations showed that the incoming composition significantly affects volumetric liquid content and concentration profiles, as well as the fraction of infiltrated mass of methanol that is released to the atmosphere. Concentration-dependent viscosity had the major impact on the liquid flow. The differences in volumetric liquid content and normalized concentration of methanol became more pronounced during transport through a soil composed of a clay lens embedded within a main matrix of sandy clay loam texture. Dispersion in the liquid-phase was the predominant transport mechanism when dispersivity at saturation was set to 7.8 cm. However, for dispersivity of 1.0 cm, changes in composition led to changes in surface tension inducing significantly higher liquid flow. In this case, liquid-phase advection was the most active transport mechanism for homogeneous soils and highly concentrated infiltrating mixtures.