Principal modes of thermospheric density variability: Empirical orthogonal function analysis of CHAMP 2001–2008 data