Two-dimensional (2D) ZrS2 monolayer (ML) has emerged as a promising candidate for thermoelectric (TE) device applications due to its high TE figure of merit, which is mainly contributed by its inherently low lattice thermal conductivity. However, the inadequate understanding of the lattice anharmonicity and its effect on the phonon dispersion and lattice thermal conductivity can be misleading in revealing the actual significance of ZrS2 ML in thermoelectrics. Here we investigated the temperature-dependent phonon dispersions and corresponding thermodynamic parameters along with the lattice thermal conductivity of ZrS2 ML by including the lattice anharmonicity based on the nonperturbative self-consistent phonon (SCP) theory, which was recently implemented and well tested with the experimental results. The higherorder (quartic) force constants were extracted by using an efficient compressive sensing lattice dynamics technique, which estimates the necessary data based on the emerging machine learning program. The phonon frequency of low-energy optical modes is enhanced upon including the quartic anharmonicity. The lattice thermal conductivity is found to be increased (by 21% at 300 K) within the SCP approach, which is due to the relatively low value of phonon linewidth contributed by the anharmonic frequency renormalization. The low energy phonons (below 160 cm -1 ) of ZrS2 ML are found to have major contributions to the lattice thermal conductivity.