Indiana University-Purdue University Indianapolis (IUPUI) / Phytoplankton size classes (pico-plankton, nano-plankton, and micro-plankton) provide information about pelagic ocean ecosystem structure, and their spatiotemporal variation is crucial in understanding ocean ecosystem structure and global carbon cycling. Remote sensing provides an efficient approach to estimate phytoplankton size compositions on global scale. In the first part of this thesis, a global sensitivity analysis method was used to determine factors mainly controlling the variations of remote sensing reflectance and inherent optical properties inverse algorithms. To achieve these purposes, average mass-specific coefficients of particles were first calculated through Mie theory, using particle size distributions and refractive indices as input; and then a synthesis remote sensing reflectance dataset was created using Hydrolight. Based on sensitivity analysis results, an algorithm for estimating phytoplankton size composition was proposed in the second part. This algorithm uses five types of spectral features: original and normalized remote sensing reflectance, two-band ratios, continuum removed spectra, and spectral curvatures. With the spectral features, phytoplankton size compositions were regressed using support vector machine. According to validation results, this algorithm performs well with simulated and satellite Sea-viewing Wide Field-of-view Sensor (SeaWiFS) and Moderate Resolution Imaging Spectroradiometer (MODIS), indicating that it is possible to estimate phytoplankton size compositions through satellite data on global scale.