En este estudio, se seleccionó un valor de ventana de variación de TEC basado en el método del espectro de potencia de wavelet para analizar la relación de acoplamiento sísmico-ionosférico. En el dominio de tiempo completo, se obtuvo una periodicidad de 27 días del espectro de potencia de wavelet que pasó la prueba de significancia del 95%. Se utilizó el método del rango intercuartílico deslizante para analizar terremotos superiores a Mw6.5 en Japón desde 2011 hasta 2022, excluyendo los efectos híbridos entre terremotos cercanos. Se utilizaron el número de manchas solares (SSN), el flujo de radio de 10.7 cm (F10.7), la irradiancia solar total (TSI), la velocidad del viento solar (Vsw), el índice de actividad geomagnética en la región ecuatorial (DST) y el índice de actividad geomagnética global (KP) como índices que representan la actividad solar y geomagnética. Después de eliminar la interferencia solar y geomagnética de los cambios en las anomalías ionosféricas utilizando el método del rango intercuartílico deslizante, se verificó que los cambios en las anomalías de TEC antes del terremoto fueron causados por el terremoto y se analizaron. El análisis estadístico de las anomalías del contenido total de electrones ionosféricos (TEC) mostró que la magnitud del terremoto estaba correlacionada positivamente con la amplitud de las anomalías de TEC, pero no de manera lineal. El tiempo de ocurrencia de las anomalías ionosféricas se retrasó hasta cierto punto con el aumento de la magnitud del terremoto. Además, los cambios anormales en el día 29 (15 de febrero de 2022) antes del 20º terremoto no se ajustaron a las reglas de investigación anteriores. Según el mecanismo de acoplamiento litosfera-atmósfera-ionosfera (LAIC) y los estudios del mapa ionosférico global (GIM), la anomalía de TEC fue consistente con la proyección vertical del epicentro con una regularidad obvia. Los resultados muestran que estas anomalías de TEC pueden estar relacionadas con terremotos.