QuikSCAT backscatter, AMSR-E and DMSP SSM/I radiance data have been used to derive sea ice motion for both the Arctic and Antarctic region using the wavelet analysis tracking method. All results from QuikSCAT, AMSR-E and SSM/I in the Arctic for fall/winter period are compatible with buoys and can then be merged by some data fusion methods to generate composite sea ice motion maps for more complete coverage. Furthermore, based on this merged data set, daily sea ice deformation (divergence and maximum shear) maps have been produced and show consistent spatial and temporal patterns. For summer ice, in order to focus the tracking templates in time, descending data are separated from ascending data during pre-processing. Due to the high resolution of AMSR-E 89 GHz data, sea ice motion maps from AMSR-E 89 GHz descending and ascending data are complementary each other, and sea ice drift in summer has been derived by merging results from descending with ascending data for composite daily sea ice motion maps. The general circulation pattern of the derived summer sea ice drift agrees with buoy data. In this study, principal component analysis of both the merged ice tracking result from satellite data and pressure field from buoy have also been examined for the relationship between the principal components and eigenvectors from these two data sets. While the result shows that principal components of modes 1 and 2 from two data sets are highly correlated, which confirms that wind forcing is a major factor driving the ice drift, it also reveals that other high energy modes are not highly correlated, which may be caused by coastal effects. Principal component analysis of Arctic sea ice motion during fall/winter period in different years shows that the reverse of dominant modes or patterns is related to the Arctic Oscillation.