Thermospheric mass density from the GOCE satellite for Sun-synchronous orbits between 83.5° S and 83.5° N normalized to 270 km during 2009–2013 has been used to develop our GOCE model at dawn/dusk local solar time sectors based on the empirical orthogonal function (EOF) method. We find that both amplitude and phase of the seasonal variations have strong latitudinal and solar activity dependences, and the annual asymmetry and effect of the Sun–Earth distance vary with latitude and solar activity.