Co-Seismic Ionospheric Disturbance with Alaska Strike-Slip Mw7.9 Earthquake on 23 January 2018 Monitored by GPS

: The Mw7.9 Alaska earthquake at 09:31:40 UTC on 23 January 2018 occurred as the result of strike slip faulting within the shallow lithosphere of the Paciﬁc plate. Global positioning system (GPS) data were used to calculate the slant total electron contents above the epicenter. The singular spectrum analysis (SSA) method was used to extract detailed ionospheric disturbance information, and to monitor the co-seismic ionospheric disturbances (CIDs) of the Alaska earthquake. The results show that the near-ﬁeld CIDs were detected 8–12 min after the main shock, and the typical compression-rarefaction wave (N-shaped wave) appeared. The ionospheric disturbances propagate to the southwest at a horizontal velocity of 2.61 km/s within 500 km from the epicenter. The maximum amplitude of CIDs appears about 0.16 TECU (1TECU = 10 16 el m − 2 ) near the epicenter, and gradually decreases with the location of sub-ionospheric points (SIPs) far away from the epicenter. The attenuation rate of amplitude slows down as the distance between the SIPs and the epicenter increases. The direction of the CIDs caused by strike-slip faults may be affected by the horizontal direction of fault slip. The propagation characteristics of the ionospheric disturbance in the Alaska earthquake may be related to the complex conditions of focal mechanisms and fault location.


Introduction
The deformation of the land and seafloor surface during earthquakes will cause oscillations in the total electron contents (TECs) of the ionosphere. Various detection techniques have been used to study the ionospheric anomalies since the ionospheric disturbance was first detected after the 1964 Alaska Mw9.2 earthquake in the United States [1]. In the early days, traditional monitoring techniques such as high-frequency Doppler sounders, ionosonde, and over-the-horizon radar were used to observe ionospheric disturbances [2][3][4][5].
As the global navigation satellite system (GNSS) was used for ionospheric anomaly detection, it breaks through the space-time limitations of traditional seismic monitoring techniques, and realizes direct observation of the total electron content of ionosphere in the seismic region for the analysis of ionospheric disturbance propagation features and generation mechanisms [6][7][8][9][10][11][12][13].
The TEC oscillations in the ionosphere caused by earthquakes were called co-seismic ionospheric disturbances (CIDs) [14]. A CID has a variety of propagation characteristics (such as propagation velocity, waveform and period) [15][16][17][18][19][20][21][22][23]. These propagation characteristics were different due to different focal mechanisms and fault conditions. Using global positioning system (GPS) data from the GPS Earth Observation Network, it was found that the CID propagation characteristics for the 2011 Mw9.1 Tohoku earthquake could be caused by acoustic waves generated from Rayleigh waves, acoustic waves from the epicenter, and atmospheric gravity waves from the epicenter [24,25]. However, the CID propagation velocity and range caused by the three types of wave were different. The acoustic waves generated by focal rupture usually propagate in the near field (within 1000 km) at the velocity of sound, while the acoustic waves excited by Rayleigh waves can propagate far away in a short time [26]. The CIDs caused by different disturbance sources will overlap near the epicenter and gradually separate in the far field [11]. The waveform of CIDs will be affected by co-seismic crustal movement. The CID caused by three major earthquakes in the Kuril Arc area combined with the focal mechanism found that the waveform polarity of CID is related to the polarity of vertical co-seismic crustal movement [11]. In addition, the direction was also an important feature of CID propagation [14]. The direction of the CID made it easier to follow the direction of the magnetic field in the seismic area [13] and was also affected by the rupture zone [20].
Strike-slip earthquakes mainly produce horizontal fault movement, and the vertical displacement is smaller than that of thrust faulting [6,19]. Vertical ground motion plays a major role in the formation of CID [20,27]. There were differences in CID propagation characteristics caused by strike-slip faulting and thrust faulting [28,29]. Among CIDs caused by three strike-slip earthquakes that occurred near Indonesia, the amplitude of strike-slip faults was smaller than that of earthquakes with large vertical displacements [30]. The CID waveform caused by the two strike-slip earthquakes discovered in Turkey were all N-shaped waves and propagate at a velocity greater than the velocity of sound [31]. Seismological parameters such as the horizontal displacement of the earthquake, the size of the fault, the sliding amplitude and the magnitude of the earthquake will also be crucial with the CID caused by the strike-slip faults, in the case of controlling the vertical displacement of the co-seismic crust [28]. Although the ionospheric disturbance caused by the strike-slip earthquake was found, the earthquake details and CID propagation characteristics were not combined.
This study monitors the ionospheric disturbance after the Mw7.9 Alaska earthquake on 23 January 2018 and calculated the slant total electron content (STEC) on the basis of the GPS data provided by continuously operating reference stations (CORS) in Alaska. Detailed ionospheric disturbance signals will be extracted through singular spectrum analysis (SSA), and the CID phenomenon caused by the strike-slip earthquake was analyzed. The CID propagation characteristics (including amplitude, waveform and velocity) caused by the Alaska earthquake will be discovered. The propagation direction and amplitude attenuation rate of CIDs will be further analyzed in detail, combined with the fault condition of the Alaska earthquake. The connection between the focal mechanism of the Alaska earthquake and the ionospheric response will be revealed, which will help to monitor the short-period variations in the ionosphere after the earthquake.

Earthquake Overview
The United States Geological Survey reported that a Mw7.9 earthquake occurred at approximately 280 km southeast of Kodiak Island in the Gulf of Alaska at 09:31:40 UTC on 23 January 2018. The epicenter was (56.004 • N, 149.116 • W), and the focal depth was 14.1 km. The earthquake occurred on a steeply dipping fault striking either west-southwest (left lateral) or north-northwest (right lateral) caused by strike-slip motion in the shallow lithosphere of the Pacific Plate.
The horizontal component of the seismic displacement of the Alaska seismic waveform was significantly larger than the vertical component [32]. The co-seismic displacement of the Alaska in the southwestern region of the epicenter has an upward trend from the co-seismic displacement fields released by the Navstar University Association (UNAVCO). The 2018 Mw 7.9 of Alaska Earthquake caused multiple fault ruptures, which ruptured the Pacific plate seaward of the Alaska subduction zone [33]. The largest slip of the Alaska earthquake occurred on the right-lateral south-southeast-striking (sliding amount 2.27 m) and the left-lateral northeast-striking (sliding amount 2.42 m) fault slipping segments [34].

Global Positioning System (GPS) Data
In this study, the STECs within 1 h around the seismic area of the Alaska earthquake, which are derived from the CORS GPS data at the sampling interval of 15 s, were used as the basis for detecting the CIDs after the Alaska earthquake. The CORS stations in Alaska were largely dense, providing sufficient data for subsequent research. The geographical distribution of approximately 81 CORS GPS tracking stations selected in the study is shown in Figure 1.
(UNAVCO). The 2018 Mw 7.9 of Alaska Earthquake caused multiple fault ruptures, which ruptured the Pacific plate seaward of the Alaska subduction zone [33]. The largest slip of the Alaska earthquake occurred on the right-lateral south-southeast-striking (sliding amount 2.27 m) and the left-lateral northeast-striking (sliding amount 2.42 m) fault slipping segments [34].

Global Positioning System (GPS) Data
In this study, the STECs within 1 h around the seismic area of the Alaska earthquake, which are derived from the CORS GPS data at the sampling interval of 15 s, were used as the basis for detecting the CIDs after the Alaska earthquake. The CORS stations in Alaska were largely dense, providing sufficient data for subsequent research. The geographical distribution of approximately 81 CORS GPS tracking stations selected in the study is shown in Figure 1.

Solar Geomagnetic Data
The ionosphere was affected by solar and geomagnetic activity [35][36][37][38]. Solar and geomagnetic activity should be detected in the period before and after the earthquake, before extracting the ionospheric anomaly information of the Alaska earthquake. This study uses the geomagnetic index (Dst) data provided by the Geomagnetic Data Center in Kyoto, Japan, and the solar radiation index F10.7 data provided by the Space Environment Forecast Center. The Dst index was derived from the H component of the geomagnetic field and reflects the geomagnetic activity, which was normally in the range of -30~30 nT during the period of low geomagnetic activity, and the temporal resolution of the data was 1 h [39]. The solar radiation index F10.7 reflects the solar activity, which was usually less than 100 SFU (1SFU = 10 −22 W m −2 Hz −1 ) during the period of low solar activity, and the temporal resolution for the data was 1 day [35].

Methodology
The relative change in STEC is sufficient for CID detection [39]. The equation for calculating the relative STEC using the GPS phase observation data is expressed as:

Solar Geomagnetic Data
The ionosphere was affected by solar and geomagnetic activity [35][36][37][38]. Solar and geomagnetic activity should be detected in the period before and after the earthquake, before extracting the ionospheric anomaly information of the Alaska earthquake. This study uses the geomagnetic index (Dst) data provided by the Geomagnetic Data Center in Kyoto, Japan, and the solar radiation index F10.7 data provided by the Space Environment Forecast Center. The Dst index was derived from the H component of the geomagnetic field and reflects the geomagnetic activity, which was normally in the range of −30~30 nT during the period of low geomagnetic activity, and the temporal resolution of the data was 1 h [39]. The solar radiation index F10.7 reflects the solar activity, which was usually less than 100 SFU (1SFU = 10 −22 W m −2 Hz −1 ) during the period of low solar activity, and the temporal resolution for the data was 1 day [35].

Methodology
The relative change in STEC is sufficient for CID detection [39]. The equation for calculating the relative STEC using the GPS phase observation data is expressed as: where f 1 and f 2 are the GPS carrier frequencies, L 1 and L 2 are dual-frequency carrier phase measurements, λ 1 and λ 2 are the carrier wavelengths, N 1 and N 2 are the integer ambiguities of phase observations, and ε is the observation noise obtained from the GPS dual-frequency phase with sampling rate of 15 s. We set the cut-off elevation angle of the satellite for GPS data to 15 • to calculate STEC [27]. SSA is one method for spectral decomposition [39], and can extract non-linear trends and identify periodic signals from time series, without being constrained by the assumption of sine waves [40]. Extracting seismic ionospheric disturbance signals from the relative STEC is affected by the trend and period of the TEC background caused by the orbital motion of GPS satellites and the temporal and spatial changes of the ionosphere [41,42]. SSA extracts more reliable information from data with unknown physical properties by extracting linear trend terms, identifying periods, reducing noise processing, and reconstructing time series [43]. The SSA extended empirical orthogonal function is used to extract the principal component features from the noise-containing STEC time series. The weighted correlation analysis (w-correlation) method is used to select the appropriate eigenvectors for reconstruction [41], and the remaining sequence after subtracting the reconstructed sequence is used as the STEC disturbance signals. The detailed ionospheric disturbance signals extraction steps were as follows.

Building the Hysteresis Matrix
The STECs at seismic period (UTC 9:00-10:00 on 23 January 2018) obtained the line of sight (LOS) between the station and GPS satellites were used as the original time series. The time series length is set to 240 under the sampling frequency of 15 s. Given a time series {x} of x 1 , x 2 , · · · , x N , a delay matrix X is constructed as: where N is the length of the time series, and M is the nested dimension. The window length is set to one-third of the STEC data length, which is 80. We can obtain:

Decomposing Singular Values
Given matrix S = XX T , where X T is the transpose matrix of X, let λ 1 , · · · , λ M be the eigenvalue of S and U 1 , · · · , U M is the feature vector corresponding to λ 1 Delay matrix X can be expressed as: Elementary matrix X i = √ λU i V T i , rank(X i ) = 1, which has the same matrix structure as X. U and V are X left and right singular vectors, λ i as the top r contribution rate of X i .

of 16
As shown in Equation (2), we can know that the contribution rate of composite matrix X I corresponding to I can be expressed as:

Establishing Diagonal Average
A w-correlation method is used to analyze the correlation between reconstruction components (RCs) and grouped RCs with the same signal characteristics [44]. Assuming the obtained RC is Y i , the w-correlation of any two RC can be expressed as: When the w-correlation of the two RCs is close to one, indicating that the two RCs are more correlated. The matrix obtained by grouping is transformed into a series of new time series with a length of 240. A time series of length 240 is defined as the RCs, and the original sequence is the sum of all RCs. Z = X I k (z 1 , z 2 , · · · z N ) is defined as the sequence obtained by Z diagonal averaging, and let M * = min(M, K), The equation for diagonal averaging can be expressed as:

Obtaining the Disturbance Signal
The l RCs and corresponding feature vector λ i (1 ≤ I ≤ l) obtained through SVD and eigenvalues λ k of the RCs are arranged in accordance with size λ 1 ≥ λ 2 ≥ · · · λ l ≥ 0. The feature vector corresponding to the magnitude of the eigenvalue represents the magnitude of the change trend of the TEC signal. The cumulative contribution rate of the eigenvalues of the first five RCs obtained using Equation (3) exceeds 95%. The first five large eigenvalues are intercepted and sorted, and the corresponding sum of χ k is found and reconstructed to fully reflect the overall characteristics of the original TEC sequence, which can be expressed as: The reconstructed STEC time series main component STEC main containing the ionospheric short-period signal is used as a high-precision background reference value to extract the disturbance signal.
STEC main is subtracted from the original STEC time series to obtain the preliminary ionospheric STEC disturbance signal: 2.4.6. Eliminating Noise The preliminary STEC disturbance signal obtained is then processed by the above steps 1-5, and the component with small contribution rate is eliminated as noise. Fifteen The CID amplitude caused by strike-slip earthquake is small [29][30][31]. The detailed information of CID propagation characteristics caused by strike-slip earthquakes is difficult to extract. The disturbance information of the time series extracted by filtering will leak the signal, causing signal loss when extracting the CID with small amplitude. We use SSA to process the time series to retain more details of the disturbance signal and to analyze the details of CID features caused by the strike-slip earthquake. Figure 2 shows the ionospheric disturbance information along LOS between the station av07 and GPS satellite 05 extracted by the SSA, band-pass filter (cutoffs of 2 mHz and 8 mHz) [28,45], and high-pass filter (cutoff of 2 mHz) [11,14]. The features of ionospheric disturbance extracted by SSA and band-pass filter were more detailed, and the signal processed by high-pass filter causes greater changes. The ionospheric disturbance peak amplitude extracted by SSA is larger than that extracted by band-pass filter and high-pass filter, and the signal amplitude extracted by the band-pass filter is very small. The detailed ionospheric disturbance information extracted by SSA helps to analyze the details of CID propagation characteristics. main nary ionospheric STEC disturbance signal: The preliminary STEC disturbance signal obtained is then processed by the above steps 1-5, and the component with small contribution rate is eliminated as noise. Fifteen eigenvectors with larger eigenvalues are extracted to construct a time series as the extracted STEC disturbance signal.
The CID amplitude caused by strike-slip earthquake is small [29][30][31]. The detailed information of CID propagation characteristics caused by strike-slip earthquakes is difficult to extract. The disturbance information of the time series extracted by filtering will leak the signal, causing signal loss when extracting the CID with small amplitude. We use SSA to process the time series to retain more details of the disturbance signal and to analyze the details of CID features caused by the strike-slip earthquake. Figure 2 shows the ionospheric disturbance information along LOS between the station av07 and GPS satellite 05 extracted by the SSA, band-pass filter (cutoffs of 2 mHz and 8 mHz) [28,45], and high-pass filter (cutoff of 2 mHz) [11,14]. The features of ionospheric disturbance extracted by SSA and band-pass filter were more detailed, and the signal processed by high-pass filter causes greater changes. The ionospheric disturbance peak amplitude extracted by SSA is larger than that extracted by band-pass filter and high-pass filter, and the signal amplitude extracted by the band-pass filter is very small. The detailed ionospheric disturbance information extracted by SSA helps to analyze the details of CID propagation characteristics.

Solar-Geomagnetic Activity Analysis
The changes in the geomagnetic D st index and solar radiation F10.7 for 10 days (18-27 January 2018), including the earthquake day. The solar radiation index F10.7 is less than 80 SFU, indicating that the solar activity is low during this period. The variation range of the geomagnetic index D st is |D st | < 25, indicating that the geomagnetic field is relatively quiet. The solar-geomagnetic activity was low on the earthquake day, and would not cause strong ionospheric disturbances.

Total Electron Content (TEC) Anomalies Following the Alaska Earthquake
Based on the GPS data from CORS stations in Alaska, the disturbances of the ionosphere over the epicenter and surrounding areas during 9:00-10:00 UTC on 23 January 2018 was extracted. The height of the ionospheric shell model was given the height of the hmF2 which was calculated as 300 km by IRI2016 model [8]. Then, the location information of CIDs calculated by the ionospheric shell model was used to analyze the CID propagation characteristics. The coordinates of ionospheric pierce point, which were the intersection of line from the satellite receiver to the thin layer, were calculated, and the epicenter distances of its corresponding sub-ionospheric point (SIP) were obtained. Figure 3 shows the information along LOSs between the stations near the epicenter and all the GPS satellites (02, 05, 25, 29, and 31) used for the analysis covered the ionospheric area around the epicenter during 9:00-10:00 UTC.
sphere over the epicenter and surrounding areas during 9:00-10:00 UTC on 23 January 2018 was extracted. The height of the ionospheric shell model was given the height of the hmF2 which was calculated as 300km by IRI2016 model [8]. Then, the location information of CIDs calculated by the ionospheric shell model was used to analyze the CID propagation characteristics. The coordinates of ionospheric pierce point, which were the intersection of line from the satellite receiver to the thin layer, were calculated, and the epicenter distances of its corresponding sub-ionospheric point (SIP) were obtained. Figure 3 shows the information along LOSs between the stations near the epicenter and all the GPS satellites (02, 05, 25, 29, and 31) used for the analysis covered the ionospheric area around the epicenter during 9:00-10:00 UTC. We combine the location and time information of ionospheric anomalies recorded to determine the CID associated with the earthquake. In the information along LOS between the stations and GPS satellites (05 and 25) during 9:00-10:00 UTC on the earthquake day, CID was found above the epicenter and southwest. Figure 4 shows the time series related to the earthquake show typical "N"-shaped disturbances, appearing 8-12 We combine the location and time information of ionospheric anomalies recorded to determine the CID associated with the earthquake. In the information along LOS between the stations and GPS satellites (05 and 25) during 9:00-10:00 UTC on the earthquake day, CID was found above the epicenter and southwest. Figure 4 shows the time series related to the earthquake show typical "N"-shaped disturbances, appearing 8-12 min after the shock with amplitudes of 0.16 TECU (1TECU = 10 16 el m −2 ). The time series along LOS between the station av07 and the GPS satellite 05 is taken as an example for detailed analysis of TEC response. The disturbances began at 9:38 UTC, and the first peak occurred at 9:44 UTC with amplitude of 0.14 TECU, and the distance between SIP and the epicenter was 168 km. After a few minutes, a second peak occurred at 9:48 UTC with amplitude of 0.05 TECU, and the distance between SIP and the epicenter was 162 km. The disturbance lasted for approximately 15 min, and the ionosphere gradually recovered after 9:53 UTC. The CID propagation characteristics were similar to the 2006 Mw8.2 Kuril earthquake in northeast Japan, showing typical compression-rarefaction wave (N-shaped waves) caused by co-seismic vertical crustal movements [11]. The typical waveform of CID may be related to the co-seismic vertical ground motion in the southwest of the epicenter of the Alaska earthquake [34]. The co-seismic vertical crustal movement caused by a strike-slip fault was smaller than the thrust earthquake of the same magnitude [23,29]. Although the co-seismic vertical displacement amplitude of strike-slip faults is small, it can still affect the CID waveform. earthquake in northeast Japan, showing typical compression-rarefaction wave (N-shaped waves) caused by co-seismic vertical crustal movements [11]. The typical waveform of CID may be related to the co-seismic vertical ground motion in the southwest of the epicenter of the Alaska earthquake [34]. The co-seismic vertical crustal movement caused by a strike-slip fault was smaller than the thrust earthquake of the same magnitude [23,29]. Although the co-seismic vertical displacement amplitude of strike-slip faults is small, it can still affect the CID waveform.

Velocity of Co-Seismic Ionospheric Disturbance (CID)
When the earthquake's fault was located on the seafloor, the near-field CID was induced by acoustic gravity waves (AGWs) and internal gravity waves (IGWs) [11,14,31,46]. The observed frequency of ionospheric disturbances caused by gravity waves is 1-2 mHz, and the frequency of ionospheric disturbances caused by acoustic waves is 3-5 mHz [24,26]. The Pacific Tsunami Warning Center issued a tsunami warning after the Alaska earthquake, using frequency to identify ionospheric disturbances caused by IGWs and AGWs [10,21,47].

Velocity of Co-Seismic Ionospheric Disturbance (CID)
When the earthquake's fault was located on the seafloor, the near-field CID was induced by acoustic gravity waves (AGWs) and internal gravity waves (IGWs) [11,14,31,46]. The observed frequency of ionospheric disturbances caused by gravity waves is 1-2 mHz, and the frequency of ionospheric disturbances caused by acoustic waves is 3-5 mHz [24,26]. The Pacific Tsunami Warning Center issued a tsunami warning after the Alaska earthquake, using frequency to identify ionospheric disturbances caused by IGWs and AGWs [10,21,47]. Figure 5 shows the sequence and the corresponding spectrum diagram obtained by the av07 station to the GPS satellite 05. The perturbation signal was transformed from the time domain to the frequency domain for spectral analysis using wavelet transform to determine the frequency of the time series. At the time when the abnormal value appears in the time series, the peak frequency distribution at the corresponding time in the spectrum diagram is evident, and the peak frequency of 4.8 mHz is consistent with the atmospheric resonance frequencies caused by the upward propagation of acoustic waves [24,26,48].
The horizontal propagation velocity of the CID can be calculated on the basis of the occurrence time of the CID and the epicenter distance of the SIP. Figure 6 shows the distance-time diagram of post-seismic TEC disturbances following the Alaska earthquake. A linear fit was performed using least squares in accordance with the occurrence time and the corresponding position of the CID peaks, and the slope was approximately equal to the CID propagation velocity after the earthquake. the av07 station to the GPS satellite 05. The perturbation signal was transformed from the time domain to the frequency domain for spectral analysis using wavelet transform to determine the frequency of the time series. At the time when the abnormal value appears in the time series, the peak frequency distribution at the corresponding time in the spectrum diagram is evident, and the peak frequency of 4.8 mHz is consistent with the atmospheric resonance frequencies caused by the upward propagation of acoustic waves [24,26,48]. The horizontal propagation velocity of the CID can be calculated on the basis of the occurrence time of the CID and the epicenter distance of the SIP. Figure 6 shows the distance-time diagram of post-seismic TEC disturbances following the Alaska earthquake. A linear fit was performed using least squares in accordance with the occurrence time and the corresponding position of the CID peaks, and the slope was approximately equal to the CID propagation velocity after the earthquake.  Figure 6 shows that CIDs propagate at a certain velocity, and the amplitude decreases with time and distance between peaks and the epicenter. The propagation velocity of CID is 2.61 ± 0.06 km/s, which is much higher than the propagation velocity of acoustic waves from the epicenter and slightly lower than the horizontal propagation  Figure 6 shows that CIDs propagate at a certain velocity, and the amplitude decreases with time and distance between peaks and the epicenter. The propagation velocity of CID is 2.61 ± 0.06 km/s, which is much higher than the propagation velocity of acoustic waves from the epicenter and slightly lower than the horizontal propagation velocity of Rayleigh waves after the earthquake [47]. The velocity of CID was similar to the velocity of near-field CID caused by Rayleigh waves, such as the 2.2 km/s in the 2004 Mw8.1 Macquarie Island earthquake [28] and the 2.22 km/s in the 2012 Mw7.8 Haida Gwaii earthquake [49,50]. The difference in velocities can be attributed to the elevation angle and direction of the GPS satellite corresponding to the LOS [19], which is more obvious near the epicenter.
The solar-terrestrial environment is low in the period. The frequency of CID conforms to the frequency of atmospheric resonance caused by the upward propagation of acoustic waves [24,26], so we can infer that the ionospheric disturbances during this period were caused by the acoustic waves. There were two types of CID disturbance caused by acoustic waves; that is, the acoustic waves generated from Rayleigh wave, and acoustic wave effect of the focal rupture from the ionospheric epicenter. The horizontal velocity of CID caused by Rayleigh waves is close to the propagation velocity (2-4 km/s) of Rayleigh waves that usually exist in the near field [10,11,28,51]. The surface Rayleigh wave group velocity from the bottom pressure in Alaska was estimated to be approximately 2.7 km/s in the southwest directions [45], which was close to the CID velocity. These results indicate that the CID of the Alaska earthquake may be caused by the acoustic waves generated by Rayleigh waves. Figure 7 shows the detected ionospheric anomalies and SIP trajectories. The CIDs are distributed in the southwest region within 500 km of the epicenter. The CID propagation of the Alaska earthquake showed evident directional differences, propagating from the epicenter to the southwest. The difference in CID propagation direction was mainly affected by two aspects. On the one hand, the CID propagation was affected by the earthquake rupture zone [20]. On the other hand, as the geomagnetic field makes the plasma tend to move freely along the field lines, the CID in the northern hemisphere has a trend spreading south [13]. was mainly affected by two aspects. On the one hand, the CID propagation was affected by the earthquake rupture zone [20]. On the other hand, as the geomagnetic field makes the plasma tend to move freely along the field lines, the CID in the northern hemisphere has a trend spreading south [13]. Strike-slip fault-dominated earthquakes have a smaller rupture magnitude than thrust earthquakes, and horizontal slip is the main source of displacement [28]. The de- Strike-slip fault-dominated earthquakes have a smaller rupture magnitude than thrust earthquakes, and horizontal slip is the main source of displacement [28]. The detailed horizontal slip segment of the Alaska earthquake was used to analyze the propagation characteristics of the CID. Figure 7 shows that the horizontal propagation direction CID of the Alaska earthquake was consistent with the west-southwest direction of maximum slip segment F1. Co-seismic displacement is particularly important to the generation of CID, and the maximum co-seismic displacement of the Alaska earthquake is caused by the slip segment F1 [34]. The CID spread to the south of the epicenter due to the effect of the geomagnetic field on the plasma [13]. These results indicate that the propagation of the CID in the southwest direction was not only affected by the geomagnetic field, but may also be related to the horizontal sliding direction of the fault. The CID caused by the other main sliding segment F3 was not detected during 9:30-10:00 UTC. This may be due to the fact that the CID has a tendency to spread southward due to the influence of the geomagnetic field and the direction of the sliding segment F3, while the distribution of SIPs in the south of the epicenter was lacking within half an hour after the Alaska earthquake. The CID direction of earthquakes dominated by strike-slip faults may be affected by the horizontal sliding direction of the fault. Figure 8a shows the location of the CID peaks. As the CID propagates from the epicenter to the southwest, the amplitude changes with the time and the distance between SIPs and the epicenter. The largest ionosphere disturbance was located near the epicenter as 0.16 TECU, the amplitude decreases during the propagation process, and the amplitude decreases to 0.05 TECU in the region far from the epicenter.

Change of CID Amplitudes
In the lithosphere beneath the oceans, the infrasonic acoustic waves which are generated by the vertical ground motion of seismic waves are emitted into the atmosphere and excite the acoustic waves to propagate upward to the ionosphere [9]. Due to the attenuation of atmospheric density, the ionospheric disturbance caused by acoustic waves will be amplified when acoustic waves propagate upward to the ionosphere [12,26,45]. The near-field CID amplitudes caused by the Alaska earthquake were extremely small, which was related to the small co-seismic vertical displacement caused by strike-slip faults [28,29]. The near-field CID amplitude caused by the Alaska earthquake decreases rapidly. The characteristics of rapid amplitude attenuation were similar to the near-field CID caused by Rayleigh waves in the Mw8.1 2004 Macquarie Island earthquake [28] and the 2012 Mw7.8 Haida Gwaii earthquake [50]. It may be due to the attenuation of seismic wave amplitude caused by fluid sliding in the cracks as seismic waves propagate on the surface of the medium on the seabed [52][53][54].
We use the relative value of the processed TEC to calculate the amplitude attenuation rate of the ionospheric disturbance. Figure 8b shows the relationship between the average amplitude of the CIDs and the time and distance between the SIP's location and the epicenter. The amplitude of the CID within 200 km from the epicenter decreases at a faster rate. From 200 to 500 km, the amplitudes attenuation rate becomes slower. These characteristics were similar to the long-distance propagation features of CID found in the 2011 Mw9.1 Tohoku earthquake [55][56][57]. We found that the amplitude attenuation rate of CID decreases with increasing distance [58], not only in the CID with a longer propagation distance, but also in the near-field CID. Figure 8a shows the location of the CID peaks. As the CID propagates from the epicenter to the southwest, the amplitude changes with the time and the distance between SIPs and the epicenter. The largest ionosphere disturbance was located near the epicenter as 0.16 TECU, the amplitude decreases during the propagation process, and the amplitude decreases to 0.05 TECU in the region far from the epicenter. In the lithosphere beneath the oceans, the infrasonic acoustic waves which are generated by the vertical ground motion of seismic waves are emitted into the atmosphere and excite the acoustic waves to propagate upward to the ionosphere [9]. Due to the attenuation of atmospheric density, the ionospheric disturbance caused by acoustic waves will be amplified when acoustic waves propagate upward to the ionosphere [12,26,45]. The near-field CID amplitudes caused by the Alaska earthquake were extremely small, which was related to the small co-seismic vertical displacement caused by strike-slip faults [28,29]. The near-field CID amplitude caused by the Alaska earthquake decreases

Determining the CID Source Location
The source location was determined by selecting a known CID and using ray tracing technique [19]. The ray tracing technique assumes that the seismic wave propagates to the height of the ionosphere in the vertical direction and then propagates around the horizontal direction at the velocity of the CID. From the given velocity model, the CID origin position was calculated from the CID position and velocity information. The search range of the CID source was given between 50 • -60 • N, 140 • -155 • W. At a given horizontal propagation velocity V H of the CID, an offset of 0.1 • × 0.1 • was used as the assumed CID source location, and the occurrence time of the CID source was expressed as [19]: where T AK is the observation time when the CID peak points (the first peak of each CID) appears, D VK is 300 km of the ionosphere height, V Z is the vertical velocity, and D HK is the horizontal distance from the peaks to the CID source. The CID points detected by the Alaska earthquake were used to calculate the standard deviation (STD) of the earthquake occurrence time of the grid points. The grid point with the smallest STD can be regarded as the CID source location. Figure 9 shows that the CID source position was 56.4 • N, 148.9 • W, which was approximately located in the northeast of the epicenter at 47 km. The calculation of the CID source shows that the CID source was close to the location of the left-lateral northeast-striking slip segment, and both were located in the northeast direction of the epicenter. The upwardly propagating sound waves are affected by the horizontal neutral wind in the ionosphere, causing the estimated position of the CID source to be different from the actual position of the fault-sliding segment [19]. . CID source location of the Alaska earthquake, the black triangle represents the CID source location, the red star represents the actual epicenter location, and color depth represents standard deviation (STD).

Conclusions
Based on the TEC calculated from GPS data, the ionospheric disturbances after the Alaska earthquake were monitored. The detailed ionospheric disturbances were extracted with SSA from the STEC data. We analyzed the propagation characteristics of the Alaska earthquake dominated by strike-slip faults, including waveform, velocity and amplitude attenuation. By analyzing the propagation characteristics of CID, it is found that the ionospheric disturbance after the earthquake may be related to the complex conditions of focal mechanisms and fault location.
The CID of the Alaska earthquake show typical "N"-shaped disturbances, appearing 8-12 min after the shock with amplitudes of 0.16 TECU. The waveform of CID may be related to the co-seismic vertical displacement in the southwest of the epicenter of the Alaska earthquake. The spectral analysis showed that the center frequency at the time of anomaly was 4.8 mHz, which was consistent with the atmospheric resonance frequencies. The propagation velocity of CID was southwestward to the distance of approximately 500 km with the horizontal phase velocity of 2.61 ± 0.06 km/s. The CID of the Alaska earthquake may be caused by the acoustic waves generated from the Rayleigh wave. The propagation direction of the CID of the Alaska earthquake was consistent with the direction of the largest slip segment. The CID direction caused by strike-slip faults was not only affected by the rupture zone and the geomagnetic field, but also by the horizontal sliding direction of the fault. As the distance between SIPs and the epicenter increases, the amplitude attenuation and the attenuation rate becomes slower. This may Figure 9. CID source location of the Alaska earthquake, the black triangle represents the CID source location, the red star represents the actual epicenter location, and color depth represents standard deviation (STD).

Conclusions
Based on the TEC calculated from GPS data, the ionospheric disturbances after the Alaska earthquake were monitored. The detailed ionospheric disturbances were extracted with SSA from the STEC data. We analyzed the propagation characteristics of the Alaska earthquake dominated by strike-slip faults, including waveform, velocity and amplitude attenuation. By analyzing the propagation characteristics of CID, it is found that the ionospheric disturbance after the earthquake may be related to the complex conditions of focal mechanisms and fault location.
The CID of the Alaska earthquake show typical "N"-shaped disturbances, appearing 8-12 min after the shock with amplitudes of 0.16 TECU. The waveform of CID may be related to the co-seismic vertical displacement in the southwest of the epicenter of the Alaska earthquake. The spectral analysis showed that the center frequency at the time of anomaly was 4.8 mHz, which was consistent with the atmospheric resonance frequencies.
The propagation velocity of CID was southwestward to the distance of approximately 500 km with the horizontal phase velocity of 2.61 ± 0.06 km/s. The CID of the Alaska earthquake may be caused by the acoustic waves generated from the Rayleigh wave. The propagation direction of the CID of the Alaska earthquake was consistent with the direction of the largest slip segment. The CID direction caused by strike-slip faults was not only affected by the rupture zone and the geomagnetic field, but also by the horizontal sliding direction of the fault. As the distance between SIPs and the epicenter increases, the amplitude attenuation and the attenuation rate becomes slower. This may be related to the amplitude attenuation of seismic waves propagating on the seafloor. The location of the CID source was close to the left-lateral northeast-striking slip segment, and the propagation direction of CID may be related to the direction of sliding fault segment.
We analyzed the CID propagation characteristics of the Alaska earthquake and discovered the effect of fault location on ionospheric disturbance propagation. The rapid attenuation characteristics of the amplitude of the Alaska earthquake CID were similar to the 2011 Mw9.1 Tohoku earthquake. The CID caused by the strike-slip fault in the 2002 Mw7.9 Denali earthquake in Alaska, propagated to 3000 km with the small amplitude of 0.05 TECU. The difference in amplitude attenuation rate and range of CIDs caused by the 2018 Mw7.9 Alaska earthquake may be related to the coupling differences of the lithosphere with atmosphere in the different fault slip (inland or seabed).
This study found the connection between the detailed focus mechanism and the characteristics of CID propagation, and showed the potential for inferring the focus mechanism based on GPS data monitoring ionospheric disturbances over the epicenter. In the future, higher precision TEC data will help improve the accuracy of ionospheric disturbance monitoring. Moreover, monitoring the attenuation of the CID amplitude caused by the Alaska earthquake helps to understand the coupling mechanism between the lithosphere and the atmosphere.