Previous studies indicate that liner defects can significantly accelerate contaminant migration in landfill, shortening the migration time by orders of magnitude. Although the effect of defects upon contaminant migration and barrier capacity of liner system has been investigated based on several theoretical models, they are always limited to one-dimensional scopes. Here a two-dimensional (2-D) model was established for contaminant transport through composite liner to the soil layer, with geomembrane (GM) defects arbitrarily distributed along the horizontal direction, and then Laplace transform, Fourier transform, and boundary transform method were utilized to derive its semi-analytical solution, which was verified by experimental data and numerical results. The obtained solution can take nonuniformly distributed defects into account and demonstrates the 2-D contaminant concentration and flux profiles in the spatiotemporal domain. Computational results of the solution are employed to illustrate the roles of material coefficients and modeling sizes upon the variations of contaminant concentration and flux. Results support that contaminant migration is expedited by the GM defect in landfill, and geosynthetics clay liner can reduce the pollution level caused by GM failure. Moreover, the contaminant flux will not be significantly influenced by the defect positions, which, however, leads to non-negligible impact on concentration profiles.