<p>This study introduces a novel nonlinear mathematical model by considering three dynamic variables: densities of human population and industries, and concentration of atmospheric pollutants. The formulated model assumes that the industries are established proportional to the density of human population for their livelihood, posing serious health risks due to the emission of pollutants. Also, due to atmospheric pollutants, the government mandates the relocation of industries. The formulated model exhibits a unique, two or three-interior equilibria depending on parameter values, which further elicits the robust behavior and rich dynamics. Bifurcation analysis shows the emergence of transcritical, saddle-node, and supercritical Hopf bifurcations when the establishment rate of industries crosses a critical threshold. The model is further extended using a sustainable pollution mitigation strategy ‘water spraying’, revealing the results that when water is sprayed proportional to the atmospheric pollutants, the equilibrium value of concentration of pollutants decreases under a specific condition; this minimizes the negative effect of pollution on human population and terminates the limit cycle oscillations, which leads to gain stability along with the industrialization. Also, at the fixed establishment rate of industries, the stability region enhances on increasing the spraying rate of water. It is also obtained that to reduce the atmospheric pollutants, the natural depletion rate of pollutants as well as scavenging rate of pollutants due to sprayed water should be large; and if the relocation of industries is unattainable then the minimization of pollutant’s emission due to anthropogenic or industrial activities or both together also work effectively.</p>