Given two CM elliptic curves over a number field and a natural number m, we establish a polynomial lower bound (in terms of m) for the number of rational primes p such that the reductions of these elliptic curves modulo a prime above p are m-isogenous. The proof relies on higher Green functions and theorems of Gross-Zagier and Gross-Kohnen-Zagier. A crucial observation is that the Fourier coefficients of incoherent Eisenstein series can be approximated by those of coherent Eisenstein series of increasing level. Another key ingredient is an explicit upper bound for the Petersson norm of an arbitrary elliptic modular form in terms of finitely many of its Fourier coefficients at the cusp infinity, which is a result of independent interest.
Assing et al. (Fri,) studied this question.