Accurate Land surface temperature (LST) retrievals from sensors aboard orbiting satellites are dependent on the corresponding atmospheric correction, especially in the Thermal InfraRed (TIR) spectral domain (8-14 µm). To remove the atmospheric effects from at-sensor measured radiance in the TIR range it is needed to characterize the atmosphere by means of three specific variables; the upwelling path and the hemispherical downwelling radiances plus the atmospheric transmissivity. Those variables can be derived from the previous knowledge of vertical atmospheric profiles of air temperature and relative humidity at different geo-potential heights and pressures. In this work, the above mentioned atmospheric variables were analyzed for three specific weather station site located in Spain, at three different altitudes. These variables were calculated with atmospheric profiles retrieved from three different sources; The National Centers for Environmental Prediction (NCEP) web-tool atmospheric profiles calculator, the MODIS (MOD07) product and the radiosoundings available on the web of the University of Wyoming (WYO), which are launched by the Agencia Estatal de Meteorologia (AEMET), in the particular case of Spain. Atmospheric profiles from 2010 to 2013 were obtained to carry out the present study. Results from comparison of these three different atmospheric profiles show that the NCEP profiles characterize the atmosphere in a better manner than MOD07. Average results values of the three MODIS spectral bands 29, 31 and 32 show a BIAS of 0.06 Wm-2µm-1sr-1 and RMSE of ±0.2 Wm-2µm-1sr-1 for upwelling radiance, a BIAS of 0.13 Wm-2µm-1sr-1 and RMSE of ±0.3 Wm-2µm-1sr-1 for the donwelling radiance and a BIAS of -0.008 and RMSE of ±0.03 for the atmospheric transmissivity. In terms of simulated LST, these errors yield a deviation of ±0.9 K when applying a single-channel method.