The long-term impact of aerosol on ocean in the subarctic Pacific was investigated based on satellite remote sensing data from MODIS (2003-2018). The spatial, temporal distributions of aerosol optical depth (AOD) and marine net primary productivity (NPP) and their correlations in the area of interest divided into four areas every 10° in longitude in the subarctic Pacific (40° N～50° N, 160° E～160° W) were analyzed and compared. The results show that AOD and NPP in each ares vary seasonally, with AOD peaking in July and NPP peaking in August or September, and both of them have an obvious decreasing trend from west to east. According to the delay correlation analysis of the monthly-averaged long time series data, it is found that AOD in each area has a strong correlation with NPP with a lag of 1-2 months, and the highest Pearson correlation coefficients of the four aeras from west to east are 0.75, 0.84, 0.79 and 0.74, respectively. Through the classification of aerosol according to the AOD peak in July, it is found that dust is an important aerosol type affecting NPP in this area. And it is also found that there is no obvious change trend from west to east for sea surface temperature (SST) and photosynthetically available radiation (PAR) in each area, which effectively excluded their influence on the spatial variation trend of NPP from west to east.