The mass transport velocity in a thin layer of muddy fluid is studied theoretically. The mud motion is driven by a periodic pressure load on the free surface, and the mud is described by a power-law model. Based on the key assumptions of the shallowness and the small deformation, a perturbation analysis is conducted up to the second order to find the mean Eulerian velocity in an Eulerian coordinate system. The numerical iteration method is adopted to solve these non-linear equations of the leading order. From the numerical results, both the first-order flow fields and the second-order mass transport velocities are examined. The verifications are made by comparing the numerical results with experimental results in the literature, and a good agreement is confirmed. The mass transport velocity in a thin layer of muddy fluid is studied theoretically. The mud motion is driven by a periodic pressure load on the free surface, and the mud is described by a power-law model. Based on the key assumptions of the shallowness and the small deformation, a perturbation analysis is conducted up to the second order to find the mean Eulerian velocity in an Eulerian coordinate system. The numerical iteration method is adopted to solve these non-linear equations of the leading order. both the first-order flow fields and the second-order mass transport velocities are examined the numerical results with experimental results in the literature, and a good agreement is confirmed.