We present 126 eclipsing binary candidates among 4683 Catalina Sky Surveys (CSS) detached and semi-detached eclipsing binary systems (EBs) showing cyclic or quadratic period variations over a 12 yr time span. By using inverse Gaussian profiles of the eclipses coupled with a Markov chain Monte Carlo procedure, times of minima (ToM) were calculated and diagrams with eclipse timing variations (ETVs) were constructed. Numerical tests were performed, involving synthetic EBs with period variations generated by the PHOEBE 2.0 engine and actual data for EBs with well-known period variations from the literature, to verify that the calculation of ToM variations for our CSS systems is reliable. A total of 63 out of the 126 EBs show likely cyclic ETVs, while the remainder present quadratic behaviour instead. Periods, amplitudes, period change rates, and associated errors were determined by using sinusoidal and parabolic models. 12 out of the 63 EBs (19 per cent) that appear to exhibit periodic ETVs are low-mass candidates. Additionally, four out of 126 also have maximum quadrature light variations. The possibility that the cyclic variations are caused by the light traveltime effect due to the presence of a tertiary companion is investigated. The possible nature of the quadratic ETVs is also discussed.

Cone search capability for table J/MNRAS/503/2979/table1 (Linear ephemeris (epoch HJD and period Pbin) for 126 CSS EBs and mean error of the calculated ToM)