Estimation of daily photosynthetically active radiation (PAR) is of primary importance for monitoring the ocean primary production and the subsequent production of carbon by phytoplankton at global scale from remote sensing ocean color sensors. On the other hand, aerosol abundance and composition play a critical role in the modulation of PAR. In this study, an original algorithm, so-called OLCIPAR, is proposed for routinely determining the daily PAR from optical satellite sensors such as the OLCI sensor aboard Sentinel-3 (ESA). The OLCIPAR algorithm has been developed to overcome some of the limitations of the current existing methods. In particular, multiple scattering effects induced by the atmospheric layer are taken into account based on exact radiative transfer calculations. Another advantage of OLCIPAR method is to consider a great variety of aerosol models to better account for their optical variability as observed in real world conditions. The OLCIPAR algorithm was applied to the archive of MERIS data, whose sensor is similar to OLCI. The validation of the retrieved daily PAR was carried out based on comparison with the time series acquired by the BOUSSOLE oceanographic buoy moored in the Mediterranean Sea. Results show a regression slope of 1% and an accuracy within 10% which confirms the robustness of the algorithm. The comparison of OLCIPAR retrievals with the products routinely distributed by NASA shows that estimates of PAR differ by up to 20% in the subtropical Atlantic Ocean where important amounts of dust aerosols are present. The improvements brought by OLCIPAR method for deriving the daily PAR could thus permit to better assess the impact of aerosols on reduction of PAR with implications on the estimation of oceanic primary production.