temps = WeatherData["KHOU", "MeanTemperature", {{1975}, {2007}, "Year"}];
fit = LinearModelFit[{#1[[1]], #2} & @@@ temps, {1, year}, year];
DateListPlot[{temps, {{#}, fit[#]} & /@ temps[[All, 1, 1]]}, Joined -> {False, True}, PlotStyle -> Thick, PlotRange -> All]