ABSTRACT Reconstructing high-quality Normalized Difference Vegetation Index time series data is essential for ecological and agricultural applications but remains challenging in heavily cloudy areas. Fusing Sentinel SAR and optical data with deep learning could be helpful but is also challenging for stable models due to unstable SAR-NDVI relationships caused by imaging mechanism differences and environmental complexities. In this study, we developed a new Bidirectional Recurrent Imputation for Optical-SAR fusion (BRIOS) model to reconstruct high-quality Sentinel-2 NDVI time series data. BRIOS designs a two-layer recurrent architecture that integrates the autocorrelation of discrete, cloud-free NDVI observations into the model for a more stable SAR-NDVI relationship. Evaluating BRIOS against three baseline methods (GF-SG spatiotemporal fusion, Harmonic regression interpolation, and MCNN-Seq deep learning) across three full Sentinel-2 tiles in reconstructing 8-day NDVI time series, BRIOS consistently outperformed in scenarios of either random or continuously missing data, as evidenced by lower RMSE values (e.g. 0.075 for BRIOS vs. 0.108 for GF-SG vs. 0.143 for Harmonic regression vs. 0.303 for MCNN-Seq), better Edge index, and high linear correlation coefficients (R values up to 0.97). Further ablation experiments revealed that deep integration of NDVI autocorrelation features and SAR temporal change patterns has improved the stability and generalization of BRIOS. Discussions on the model's scalability across various cloud sizes and training dataset sizes affirm its practicality for broad-scale application in vegetation monitoring under challenging cloudy conditions.