We consider the OðN þ 1Þ=OðNÞ nonlinear sigma model for large N as an effective theory for lowenergy QCD at finite temperature T, in the chiral limit. At T ¼ 0 this formulation provides a good description of scattering data in the scalar channel and dynamically generates the f 0 ð500Þ pole, the pole position lying within experimental determinations. Previous T ¼ 0 results with this model are updated using newer analysis of pion scattering data. We calculate the pion scattering amplitude at finite T and show that it exactly satisfies thermal unitarity, which had been assumed but not formally proven in previous works. We discuss the main differences with the T ¼ 0 result, and we show that one can define a proper renormalization scheme with T ¼ 0 counterterms such that the renormalized amplitude can be chosen to depend only on a few parameters. Next, we analyze the behavior of the f 0 ð500Þ pole at finite T, which is consistent with chiral symmetry restoration when the scalar susceptibility is saturated by the f 0 ð500Þ state, in a second-order transition scenario and in accordance with lattice and theoretical analysis.