The temperature dependence of the symmetry energy for isotopic chains of even-even Ni, Sn, and Pb nuclei is investigated in the framework of the local density approximation (LDA). The Skyrme energy density functional with two Skyrme-class effective interactions, SkM* and SLy4, is used in the calculations. The temperaturedependent proton and neutron densities are calculated through the HFBTHO code that solves the nuclear SkyrmeHartree-Fock-Bogoliubov problem by using the cylindrical transformed deformed harmonic-oscillator basis. In addition, two other density distributions of 208 Pb, namely the Fermi-type density determined within the extended Thomas-Fermi (TF) method and symmetrized-Fermi local density obtained within the rigorous density functional approach, are used. The kinetic energy densities are calculated either by the HFBTHO code or, for a comparison, by the extended TF method up to second order in temperature (with T 2 term). Alternative ways to calculate the symmetry energy coefficient within the LDA are proposed. The results for the thermal evolution of the symmetry energy coefficient in the interval T = 0-4 MeV show that its values decrease with temperature. The temperature dependence of the neutron and proton root-mean-square radii and corresponding neutron skin thickness is also investigated, showing that the effect of temperature leads mainly to a substantial increase of the neutron radii and skins, especially in the more neutron-rich nuclei, a feature that may have consequences on astrophysical processes and neutron stars.