The determination of aerosol's effect contains much uncertainty. During quantification of aerosols via remote sensing, surface reflectance error of 0.01 could bring Aerosol Optical Depth (AOD) error of 0.1. In order to avoid this, simultaneous retrieval of AOD and surface properties would be a promising method. In this paper, we present a novel analytic solution to atmospheric radiative transfer equation and utilize this solution to retrieve AOD and surface albedo simultaneously from bi-temporal geostationary FY-2 remote sensing data.