The statistical analysis of impact of the top 50 solar flares of X-class (1997-2023) on the global seismic activity, as well as on the earthquake preparation zones located in illuminated part of the globe and in the area of 5000 km radius around the subsolar point was carried out. It is shown by a method of epoch superposition that for all cases the increase of seismicity is observed, especially in the region around the subsolar point (up to 38%) during 10 days after the solar flare in comparison with preceding 10 days. The case study of aftershock sequence of strong M=9.1 earthquake (Su-matra-Andaman Islands, 26.12.2004) after the solar flare of X7.2 class (20.01.2005) demonstrated that the number of aftershocks with magnitude M≥2.5 increases more than 20 times after the solar flare with a delay of 7 days. For the case of the Darfield earthquake (M=7.1, 04.09.2010, New Zealand) it was shown that strong solar flares of class X and M probably triggered two strong aftershocks (M>6) with the same delay of 6 days on the Port Hills fault, which is the most sensitive to external electromagnetic impact from point of view of the fault electrical conductivity and orientation. Based on the obtained results the possible application of natural electromagnetic triggering of earthquakes is discussed for a short-term earthquake prediction using confidently recorded strong external electromagnetic triggering impacts on the specific earthquake preparation zones, as well as ionospheric perturbations due to aerosol emission from the earthquake sources recorded by satellites.