The quality and accuracy of remote sensing instruments continue to increase, allowing geoscientists to perform various quantitative retrieval applications to observe the geophysical variables of land, atmosphere, ocean etc. The explosive growth of time-series RS data over large-scales poses great challenges on managing, processing and interpreting RS ''Big Data''. To explore these time series RS data efficiently, in this paper, we design and implement a high performance framework to address the time-consuming time-series quantitative retrieval issue on a GPU cluster, taking the aerosol optical depth (AOD) retrieval from satellite images as a study case. The presented framework exploits the multi-level parallelism for time-series quantitative RS retrieval to promote efficiency. At the coarse-grained level of parallelism, the AOD time series retrieval is represented as multi-directed acyclic graph (DAG) workflows and scheduled based on a list-based heuristic algorithm, heterogeneous earliest finish time (HEFT), taking the idle slot and priorities of retrieval jobs into account. At the fine-grained level, the parallel strategies for the major remote sensing image processing algorithms divided into three categories, i.e. the point or pixel-based operations, the local operations and the global or irregular operations have been summarized. The parallel framework was implemented with message passing interface (MPI) and compute unified device architecture (CUDA), and experimental results with the AOD retrieval case verify the effectiveness of presented framework.