From 1980 to 2021, earthquakes have caused more than 846,000 casualties and about US$ 661 billion in economic losses. At present, there are no reliable earthquake precursors to generate alerts. Currently, the link between earthquakes and Total Electron Content (TEC) variations measured by Global Navigation Satellite Systems (GNSS) monitoring ground stations has been studied. However, GNSS ground-based monitoring stations are irregularly disseminated around the globe with significant gaps, particularly in the ocean's regions. In this study, we analyze ionospheric intensity scintillation anomalies as potential proxies of earthquakes. NASA CYGNSS GNSS-R (Reflectometry) acquired by CYGNSS/DDMI (Delay Doppler Mapping Instrument) from 2017 to 2021 have been used to calculate and analyze the anomalies in the S4 parameter over ocean areas affected by earthquakes. More than 30,000 ocean earthquakes within ±40° in latitude and with a magnitude larger than M4 have been examined.The daily planetary geomagnetic index Kp was used to discard data that may be disturbed due to space weather conditions. Additionally, a daily sea wind speed mask was used to eliminate sea states that impact the signal reflectivity. The Standard Deviation (STD) and the Interquartile (IQT) time series methods have been used to detect these S4 anomalies. The confusion matrix, the Receiver Operating Curve (ROC), and some other Figures of Merit (FoM) have been used for the first time to evaluate and improve the performance of the prediction parameters and identify the optimum configuration to be used as a potential proxy of earthquake occurrence. As a result, a small, but detectable positive S4 anomaly was detected between 1 and 7 days before the earthquakes under study.