Study on Low-Frequency Abnormal Signal and Structural Characteristics of 2015 Azuoqi Ms5.8 Earthquake

: In this article, the phenomenon of low-frequency abnormal signals before earthquakes, which reﬂects the three elements of earthquakes and the beneath structure change information, is discussed. Based on the data recorded at the Shizuishan (SZS), Wuhai (WUH) and Dongshenmiao seismic stations around the epicenter of the Ms5.8 earthquake in Azuoqi, Inner Mongolia, in 2015, the low-frequency abnormal signal from the seismic waves before this earthquake is extracted. At the same time, the autocorrelation method is used to extract the reﬂected waves of the main interface from teleseismic events recorded by the seismic array in the epicenter area, and then the change information from the beneath structure is obtained. It is explained in time and space that the low-frequency abnormal signal before the main earthquake, extracted from the continuous waveform, is directly related to the change in the underground structure near the epicenter, and it can be determined that the wave propagation direction f the crustal stress before the earthquake is from south to north, and it continues to accumulate near the epicenter until the main earthquake occurs.


Introduction
A large number of experimental studies show that rocks under high stress will experience the development process of "micro fracture initiation (micro scale and micro energy)-fracture propagation, penetration (small-scale and small energy fracture)-local rock instability (mesoscale and medium energy impact)-rock mass structure failure and large-scale instability (large-scale and large-energy release) [1]".
One significance of extracting the low-frequency seismic signals is to explore the abnormal information in the above process. After the huge energy generated by the source is released, most of the high-frequency energy is absorbed by the medium or continuously weakened in the form of aftershocks, and a small part of the low-frequency energy has the characteristics of long wave length and slow attenuation, which is not easily absorbed and can be transmitted to thousands of kilometers away. Materials with strong high-frequency absorption will have positive prospects in seismic fortification [2,3].
The low-frequency signal from seismic waves also belongs to a kind of underground pulsation information, which may be the external manifestation of slow earthquakes (also known as intermittent deep tremors, slow slip events, low-frequency earthquakes, etc.) [4,5]. Compared with the conventional seismic radiation energy, the energy radiated from slow earthquakes is less than that from conventional earthquakes [6]. This micro fracture until rock mass instability and macro failure play an important role in the process of seismic nucleation [7,8]. The micro fracture behavior is different from that of rigid structural media [9]. The mechanism explanation of slow earthquakes may be a slow slip. This mechanism explanation has been proved by the rock fracturing in the laboratory [10,11], that is, the stable slip and stick slip phenomena before a rock fracture in the laboratory may correspond to the stable slip and stick slip phenomena near the subduction zone of a seismic fault. The occurrence of these phenomena mainly depend on the critical stiffness coefficient of the fault [6]. When the fault is unstable or conditionally stable, it leads to the phenomenon of a slow slip. However, the mechanical model of moderately strong earthquakes in nature is very complex and may be triggered by the joint action of many factors [12]. The occurrence of a slow slip event may not necessarily produce a destructive earthquake, but if the slow slip event occurs in the locking area deep within a plate or fault, resulting in a continuous increase in stress accumulation in the locking area, a strong earthquake may be triggered. Extracting the low-frequency abnormal signals from the continuous seismic waveforms is an attempt to find the phenomenon of a slow slip caused by the continuous stress accumulation in this locked area. It is certain that the strongest earthquakes in history showed a locked state underground before the main earthquake: for example, the Wenchuan earthquake of magnitude 8 that occurred in the Chinese mainland on 12 May 2008, and the epicentral area of the Inner Mongolia Left Banner 5.8 earthquake on 15 April 2015. By studying the low-frequency abnormal signals from seismic waves before and after a main earthquake, we can identify the characteristics of a slow slip and of the energy release from the stress accumulation in the locking area [7,13]. Therefore, it is very necessary to make a retrospective study of historical moderate and strong earthquakes. This study can first determine the location of the locking area. If a slow earthquake occurs in the locked section of the fault, when the stress accumulates to the critical value, the locked section of the fault will be unlocked. At the moment of unlocking, the fault friction will suddenly decrease and produce rapid sliding or fracturing, which explains the hypothesis that a slow earthquake in the locked area will trigger the occurrence of a normal earthquake. If slow earthquakes occur in areas with nonlocking or developed fractures, the occurrence of slow earthquakes is only a process of normal stress release equilibrium, which can explain why there are no strong earthquakes in some areas even if there are slow earthquakes. In this way, in order to judge the impact of slow earthquakes or the low-frequency anomaly information about seismic waves on a region, it is necessary to first determine whether the region is a stress locking area. On this premise, the monitoring of slow earthquakes and the low-frequency anomaly information will have more practical significance.
In order to prove that the occurrence of the low-frequency abnormal signal in field seismic waves originates from the change in underground structure, while extracting the low-frequency signal from continuous seismic waveforms, this paper extracts the underground structure information from distant seismic waveforms according to the same frequency period and reflects the change characteristics of an underground structure in different time periods.

Extraction of Underground Structure
Formation wave impedance is restricted by two factors: lithology and physical properties. Lithology refers to the characteristics of a rock structure and its components, while physical properties refer to the characteristics of rock porosity, elasticity and fluid in the aquifer. Since the reflection coefficient can reflect the geological elements in the seismic records, we should extract the reflection coefficient series from the seismic records to obtain the geological information, and then use the obtained reflection coefficient series to convert the seismic waveform section into the wave impedance section, so as to comprehensively and accurately display most of the geological information contained in the seismic records [14,15]. The following figures (Figures 1 and 2) show the relationship between the seismic wavelet, reflection coefficient and wave impedance. Appl. Sci. 2022, 12, x FOR PEER REVIEW 3 of 16  For non-interference seismic traces, the results of a convolution in the seismic wavelet and reflection coefficient series are shown in Formula (1): where ) (t ω is the seismic wavelet, representing the reflection coefficient series and seismic trace. According to Formula (1), the geological information contained in the seismic record will only exist in the only two structural element carriers in the seismic wavelet or reflection coefficient [15,16]. It is impossible to separate every pure instantaneous reflection wavelet from seismic records to study the changes between the top and bottom reflection wavelets in a single target layer to obtain some geological information about the layer. In other words, it is not easy to obtain underground formation information from the seismic wavelet. Therefore, we can only start from another angle: the reflection coefficient series ( Figure 2). The reflection coefficient series depends on the wave impedance difference between the upper and lower media at each formation interface, that is, it is determined by the wave impedance in the formation (the product of the wave velocity and density) [17].
The reflection coefficient of a single interface is expressed by Equation (2):  For non-interference seismic traces, the results of a convolution in the seismic wavelet and reflection coefficient series are shown in Formula (1): where ) (t ω is the seismic wavelet, representing the reflection coefficient series and seismic trace. According to Formula (1), the geological information contained in the seismic record will only exist in the only two structural element carriers in the seismic wavelet or reflection coefficient [15,16]. It is impossible to separate every pure instantaneous reflection wavelet from seismic records to study the changes between the top and bottom reflection wavelets in a single target layer to obtain some geological information about the layer. In other words, it is not easy to obtain underground formation information from the seismic wavelet. Therefore, we can only start from another angle: the reflection coefficient series ( Figure 2). The reflection coefficient series depends on the wave impedance difference between the upper and lower media at each formation interface, that is, it is determined by the wave impedance in the formation (the product of the wave velocity and density) [17].
The reflection coefficient of a single interface is expressed by Equation (2): For non-interference seismic traces, the results of a convolution in the seismic wavelet and reflection coefficient series are shown in Formula (1): where ω(t) is the seismic wavelet, representing the reflection coefficient series and seismic trace. According to Formula (1), the geological information contained in the seismic record will only exist in the only two structural element carriers in the seismic wavelet or reflection coefficient [15,16]. It is impossible to separate every pure instantaneous reflection wavelet from seismic records to study the changes between the top and bottom reflection wavelets in a single target layer to obtain some geological information about the layer. In other words, it is not easy to obtain underground formation information from the seismic wavelet. Therefore, we can only start from another angle: the reflection coefficient series ( Figure 2). The reflection coefficient series depends on the wave impedance difference between the upper and lower media at each formation interface, that is, it is determined by the wave impedance in the formation (the product of the wave velocity and density) [17]. The reflection coefficient of a single interface is expressed by Equation (2):

Data Processing
Two types of data were used for the study data. One is the continuous seismic waveform data with a one-hour window length, where the sampling rate is 0.01 Hz. This type of data is used to extract the maximum value of a spectral envelope. The other is the data between 10 s before and 100 s after the arrival of teleseismic P wave, that is, the direct wave and its coda data, which is used to extract the information from the underground structure. In fact, the two different types of data are best recorded from the same dense array. For this article, the continuous seismic wave data are recorded by three fixed stations (DSM, WUH, SZS) around the epicenter, while the teleseismic data are recorded by the China Science Array [18].

Spectrum Envelope
Intuitively, the continuous waveform recorded by a seismometer with a window length of 1 h is in a basically stable form when no obvious events such as earthquakes or blasting are recorded, and the maximum value of the corresponding spectrum envelope amplitude fluctuates in a basically horizontal state [19]. Therefore, two key features are extracted when processing the continuous waveform. One is to obtain the spectrum from the 1 h continuous seismic waveform through a fast Fourier transform, and the other is to find the maximum amplitude of the spectrum envelope.
{x i } is the continuous waveform data sequence recorded in 1 h, which is the amplitude spectrum. The curve formed by connecting the highest lines of the different frequency amplitudes from the obtained spectrum is called a spectrum envelope, as shown in Figures 3 and 4 [20].
array. For this article, the continuous seismic wave data are record stations (DSM, WUH, SZS) around the epicenter, while the teleseismi by the China Science Array [18].

Spectrum Envelope
Intuitively, the continuous waveform recorded by a seismomet length of 1 h is in a basically stable form when no obvious events such blasting are recorded, and the maximum value of the corresponding amplitude fluctuates in a basically horizontal state [19]. Therefore, tw extracted when processing the continuous waveform. One is to obtain the 1 h continuous seismic waveform through a fast Fourier transform find the maximum amplitude of the spectrum envelope.

Autocorrelation Extraction of Reflected Wave
The reflected P wave can be obtained by autocorrelation accor component of the continuous waveform recorded by the fixed station lation amount from this method is too large when extracting the ref resulting in low efficiency. In this article, the P wave, S wave and coda w mic events are used to extract the reflected wave, and the calculation improved. The core idea is that the reflection seismogram with a surf face receiver can be obtained from the autocorrelation of the transmit deep source and the receiver at the same position [22], that is, the au vertically incident transmitted wave is equivalent to the reflected wav sence of this idea is to extract reflected waves from the multiples bet and a discontinuity using autocorrelation [23]. In this way, the trans the teleseismic P wave and its coda is transformed into a reflection pro reflection coefficient of the layered media can be studied [24]. The pro a station superposition autocorrelation map to extract the crustal ref can improve the consistency of the crustal reflection coefficient and certain extent. Compared with an active source exploration, autoco technology, as passive source detection technology, makes up for the active source exploration in cost and depth [25]. Compared with the function, the autocorrelation method can obtain deeper underground The teleseismic data recorded by a seismometer includes the s transmission path information from a source to discontinuity, reflectio information of multiples from a discontinuity to a free surface, in whi mation is negligible, so a non-interference seismic channel teleseismi pressed in the frequency domain by Equation (6):

Autocorrelation Extraction of Reflected Wave
The reflected P wave can be obtained by autocorrelation according to the vertical component of the continuous waveform recorded by the fixed station [21], but the calculation amount from this method is too large when extracting the reflected wave signal, resulting in low efficiency. In this article, the P wave, S wave and coda wave of the teleseismic events are used to extract the reflected wave, and the calculation efficiency is greatly improved. The core idea is that the reflection seismogram with a surface source and surface receiver can be obtained from the autocorrelation of the transmitted waves from the deep source and the receiver at the same position [22], that is, the autocorrelation of the vertically incident transmitted wave is equivalent to the reflected wave. The physical essence of this idea is to extract reflected waves from the multiples between a free surface and a discontinuity using autocorrelation [23]. In this way, the transmission problem of the teleseismic P wave and its coda is transformed into a reflection problem, and then the reflection coefficient of the layered media can be studied [24]. The process of constructing a station superposition autocorrelation map to extract the crustal reflection information can improve the consistency of the crustal reflection coefficient and suppress noise to a certain extent. Compared with an active source exploration, autocorrelation extraction technology, as passive source detection technology, makes up for the shortcomings of an active source exploration in cost and depth [25]. Compared with the traditional receiver function, the autocorrelation method can obtain deeper underground information [26,27].
The teleseismic data recorded by a seismometer includes the source information, transmission path information from a source to discontinuity, reflection and transmission information of multiples from a discontinuity to a free surface, in which the source information is negligible, so a non-interference seismic channel teleseismic record can be expressed in the frequency domain by Equation (6): From the relationship between autocorrelation and convolution, it can be seen that after autocorrelation: When the teleseismic P wave is vertically incident to the horizontal layer as the transmission wave, the source information as the point source can be ignored, and the Appl. Sci. 2022, 12, 259 6 of 15 relationship between the transmission response and reflection response can be expressed by Equation (8): The autocorrelation method is used to obtain the seismic reflection wave with the station as the virtual source and receiving point. The processing flow mainly includes four steps: pretrending, autocorrelation, moveout correction and stacked, as shown in Figure 5.
The autocorrelation method is used to obtain the seismic reflection wave with the station as the virtual source and receiving point. The processing flow mainly includes four steps: pretrending, autocorrelation, moveout correction and stacked, as shown in Figure  5. Because this method uses the P wave of the teleseismic seismic events recorded by the vertical component and the seismic phases, such as global seismic phases PKP and PKIKP and their coda, to construct the reflected wave, the most critical step is NMO correction, which is used to eliminate the travel time difference of the seismic reflection phase, caused by the incident slowness, and correct the time difference of the rays to make them incident vertically [28].
According to the Taylor formula, Equation (9) can be obtained: The specific determination of the moveout correction is as follows: it is assumed that the interface depth is 0 h , the velocity of the overlying medium is 0 v and the two-way time of the vertical reflection is  (9), so we can find the moveout correlation with the incident slowness p is

Low-Frequency Abnormal Signal of Seismic Wave
The maximum value of the amplitude frequency envelope is tracked for the continuous seismic waveform data recorded at the Wuhai (WUH), Shizuishan (SZS) and Dongshengmiao (DSM) seismic stations from 0:00 on April 1 to 0:00 on 16 April 2015. The results show that the maximum values of the amplitude frequency envelopes of the three stations have obvious morphological differences (Figures 6-8). On the basis of the relatively gentle changes since 1 April 2015, the curve shape of WUH, closest to the epicenter of the Azuoqi Ms5.8 earthquake, began to fluctuate from 6:00 on 10 April 2015. It is considered that this morphology from WUH corresponds to the theory of fracture nucleation, which may be the manifestation of sub-instability. The result from DSM is relatively stable and the same sharp pulse signal as that from WUH before 6 April appears only at some individual times. Moreover, the intensity of this sharp pulse signal weakened after the continuous curve fluctuation at WUH on 6 April. It is speculated that in the process of the rock fracture and nucleation under WUH, the rock fissures are squeezed due to the increasing stress, so the low-frequency signal is blocked and cannot be better transmitted to Because this method uses the P wave of the teleseismic seismic events recorded by the vertical component and the seismic phases, such as global seismic phases PKP and PKIKP and their coda, to construct the reflected wave, the most critical step is NMO correction, which is used to eliminate the travel time difference of the seismic reflection phase, caused by the incident slowness, and correct the time difference of the rays to make them incident vertically [28].
According to the Taylor formula, Equation (9) can be obtained: The specific determination of the moveout correction is as follows: it is assumed that the interface depth is h 0 , the velocity of the overlying medium is v 0 and the two-way time of the vertical reflection is τ 0 = 2h 0 /v 0 when the source is located on the surface. The reflection time in τ − P domain can be expressed as Equation (10). Let x = −v 2 0 p 2 and m = 1 2 in Equation (9), so we can find the moveout correlation with the incident slowness p is − 1 2 v 2 0 p 2 .

Low-Frequency Abnormal Signal of Seismic Wave
The maximum value of the amplitude frequency envelope is tracked for the continuous seismic waveform data recorded at the Wuhai (WUH), Shizuishan (SZS) and Dongshengmiao (DSM) seismic stations from 0:00 on 1 April to 0:00 on 16 April 2015. The results show that the maximum values of the amplitude frequency envelopes of the three stations have obvious morphological differences (Figures 6-8). On the basis of the relatively gentle changes since 1 April 2015, the curve shape of WUH, closest to the epicenter of the Azuoqi Ms5.8 earthquake, began to fluctuate from 6:00 on 10 April 2015. It is considered that this morphology from WUH corresponds to the theory of fracture nucleation, which may be the manifestation of sub-instability. The result from DSM is relatively stable and the same sharp pulse signal as that from WUH before 6 April appears only at some individual times. Moreover, the intensity of this sharp pulse signal weakened after the continuous curve fluctuation at WUH on 6 April. It is speculated that in the process of the rock fracture and nucleation under WUH, the rock fissures are squeezed due to the increasing stress, so the low-frequency signal is blocked and cannot be better transmitted to the area below DSM. The difference in the low-frequency signals tracked by the three stations may be related to the epicentral distance. Assuming that a region is equipped with enough dense seismometers, the epicenter location of moderate and strong earthquakes can be predicted through this difference feature.
the area below DSM. The difference in the low-frequency signals tracked by the three stations may be related to the epicentral distance. Assuming that a region is equipped with enough dense seismometers, the epicenter location of moderate and strong earthquakes can be predicted through this difference feature.    the area below DSM. The difference in the low-frequency signals tracked by the three stations may be related to the epicentral distance. Assuming that a region is equipped with enough dense seismometers, the epicenter location of moderate and strong earthquakes can be predicted through this difference feature.    the area below DSM. The difference in the low-frequency signals tracked by the three stations may be related to the epicentral distance. Assuming that a region is equipped with enough dense seismometers, the epicenter location of moderate and strong earthquakes can be predicted through this difference feature.

Feature Extraction of Stratigraphic Information
According to the existing seismogenic mechanisms such as slow earthquakes and micro fracture, we try further work to explore whether the medium under the monitoring station has also changed accordingly. When detecting the change in the underground media under the three stations, we use seismic waveform data less than 5 Hz. The data come from the teleseismic data recorded by the Chinese seismological array in western Inner Mongolia (Figure 9). The analysis data are selected as teleseismic events within the range of 30-90 • from the array in three different time periods: October 2013-October 2014, June 2014-April 2015 and April 2015-October 2016 ( Figure 10). Select the teleseismic event waveform P wave and its coda ([−10,100] s waveform data), and then perform autocorrelation calculation after band-pass filtering on the selected data at [0.8,2] Hz, [0. 5,4] Hz and [1.5] Hz. Finally, stack the autocorrelation results of the teleseismic events with different slowness. The results show that the structural characteristics recorded by the same station during the three different time periods have obvious changes, and the depth information recorded by the stations close to the source on the same survey line has obvious changes (Figures 11-14). In addition, the change degree of the underground structure at the different seismic stations within the same survey is also different. For example, 15673 has the most significant change in Line 4 (L4), which has the most prominent change from June 2014 to April 2015 ( Figure 11). The most significant change in survey Line 6 (L6) is 64052, which has shown obvious underground structure changes from October 2013 to July 2014 ( Figure 12). The most significant change in survey Line 7 (L7) is 15670, which has shown obvious underground structure changes from April 2015 to October 2016 after the Ms5.8 earthquake in Azuoqi, Inner Mongolia, in 2015 ( Figure 13). The most significant change in survey Line 9 (L9) is 15597, which has shown obvious underground structure changes from June 2014 to April 2015 ( Figure 14).

Feature Extraction of Stratigraphic Information
According to the existing seismogenic mechanisms such as slow earthquakes and micro fracture, we try further work to explore whether the medium under the monitoring station has also changed accordingly. When detecting the change in the underground media under the three stations, we use seismic waveform data less than 5 Hz. The data come from the teleseismic data recorded by the Chinese seismological array in western Inner Mongolia (Figure 9). The analysis data are selected as teleseismic events within the range of 30-90° from the array in three different time periods: October 2013-October 2014, June 2014-April 2015 and April 2015-October 2016 ( Figure 10). Select the teleseismic event waveform P wave and its coda ([−10,100] s waveform data), and then perform autocorrelation calculation after band-pass filtering on the selected data at [0. 8,2] Hz, [0.5,4] Hz and [1.5] Hz. Finally, stack the autocorrelation results of the teleseismic events with different slowness. The results show that the structural characteristics recorded by the same station during the three different time periods have obvious changes, and the depth information recorded by the stations close to the source on the same survey line has obvious changes (Figures 11-14). In addition, the change degree of the underground structure at the different seismic stations within the same survey is also different. For example, 15673 has the most significant change in Line 4 (L4), which has the most prominent change from June 2014 to April 2015 ( Figure 11). The most significant change in survey Line 6 (L6) is 64052, which has shown obvious underground structure changes from October 2013 to July 2014 ( Figure 12). The most significant change in survey Line 7 (L7) is 15670, which has shown obvious underground structure changes from April 2015 to October 2016 after the Ms5.8 earthquake in Azuoqi, Inner Mongolia, in 2015 ( Figure 13). The most significant change in survey Line 9 (L9) is 15597, which has shown obvious underground structure changes from June 2014 to April 2015 ( Figure 14).

Relationship between Low-Frequency Signal and Beneath Structure
As described earlier, the morphological characteristics of the low-frequency signals recorded at DSM, WUH and SZS before the Ms5.8 earthquake in Azuoqi, Inner Mongolia, in 2015, are completely different, and the low-frequency abnormal signals lasting about

Relationship between Low-Frequency Signal and Beneath Structure
As described earlier, the morphological characteristics of the low-frequency signals recorded at DSM, WUH and SZS before the Ms5.8 earthquake in Azuoqi, Inner Mongolia, in 2015, are completely different, and the low-frequency abnormal signals lasting about 100 h before the main earthquake are extracted from the data recorded at WUH (Figures 6-8). Does the low-frequency abnormal signal form from the three seismic stations reflect the state and process of rock micro fracture below the epicenter? The different low-frequency anomaly patterns from the three seismic stations also reflect that the unstable change direction of the stress field is from south to north, that is, the unstable micro fracture from the stress is transmitted from below SZS to below WUH and then blocked under WUH. Until the main earthquake of the Ms5.8 earthquake occurs, the energy continues to be transmitted north to DSM. Therefore, we compared the positions of the temporary stations around the epicenter with DSM, WUH and SZS and found that the erection positions of the stations have the following characteristics: 64052 in L6 is equivalent to SZS, 15597 in L9 is equivalent to WUH, and 15670 in L7 is equivalent to DSM. The underground structure change characteristics of 64052, 15597 and represent the underground structure change characteristics of SZS, WUH and DSM, respectively. The results show that the underground structure of SZS (64052), located about 60 km from the epicenter of the Ms5.8 earthquake in Azuoqi, Inner Mongolia, in 2015, has changed from October 2013 to July 2014. Until the Ms5.8 earthquake in Azuoqi, Inner Mongolia, its underground morphological characteristics were constantly changing, which was consistent with the characteristics of the low-frequency signal and extreme instability of Shizuishan station (Figures 15 and 16).   The underground structure morphology of WUH (15597), which is located 42 km from the epicenter of the 2015 Azuoqi Ms5.8 earthquake in Inner Mongolia, did not change as significantly as that of DSM (15670) from October 2013 to July 2014. It only showed significant changes from June 2014 to April 2015. After the earthquake, the medium morphology recovered to a state similar to that of DSM (15670). The SZS appeared earlier than the WUH, in terms of the time when the instability of the underground structure began (Figures 17 and 18).   This feature shows that the stress transmission direction is from south to north, which is consistent with the northeast compression of the Qinghai Tibet block on the Ordos block in the southwest, and that the stress transmission direction reflected by the maximum tracking results from the seismic wave spectrum envelope (Figures 19 and 20).

Conclusions and Discussion
In the research on the underground structure of the epicenter of the 2015 Azuoqi earthquake with Ms5.8 in Inner Mongolia, on the one hand, it is studied according to the survey line profile, while on the other hand, it is studied one by one according to the stations. Using the data from different time periods before and after the main earthquake, it can be clearly seen that the underground structure has changed significantly and that the change characteristics of the underground structure under the stations at different locations from the epicenter before the main earthquake are different. This difference corresponds well to the low-frequency abnormal signal before the earthquake. The details are summarized as follows: In the three seismic stations, WUH, SZS and DSM, near the Ms5.8 earthquake in 2015, the low-frequency abnormal signals extracted from the recorded seismic waveforms showed completely different forms before the earthquake, and the station closest to the epicenter, WUH, showed unstable fluctuation forms about 100 h before the earthquake. This may be the manifestation of the nucleation process in the metastable state of the source.
In the process of rock fracture and nucleation beneath WUH, due to the increase in stress concentration, the rock fissures are squeezed so that the low-frequency signal is blocked and cannot be better transmitted to the position below DSM. The low-frequency signal recorded by SZS has always been in a sharp pulse state, indicating the development of underground rock fractures, which is not conducive to stress accumulation and may be the transition zone of the energy transmission process.

Conclusions and Discussion
In the research on the underground structure of the epicenter of the 2015 Azuoqi earthquake with Ms5.8 in Inner Mongolia, on the one hand, it is studied according to the survey line profile, while on the other hand, it is studied one by one according to the stations. Using the data from different time periods before and after the main earthquake, it can be clearly seen that the underground structure has changed significantly and that the change characteristics of the underground structure under the stations at different locations from the epicenter before the main earthquake are different. This difference corresponds well to the low-frequency abnormal signal before the earthquake. The details are summarized as follows: In the three seismic stations, WUH, SZS and DSM, near the Ms5.8 earthquake in 2015, the low-frequency abnormal signals extracted from the recorded seismic waveforms showed completely different forms before the earthquake, and the station closest to the epicenter, WUH, showed unstable fluctuation forms about 100 h before the earthquake. This may be the manifestation of the nucleation process in the metastable state of the source.
In the process of rock fracture and nucleation beneath WUH, due to the increase in stress concentration, the rock fissures are squeezed so that the low-frequency signal is blocked and cannot be better transmitted to the position below DSM. The low-frequency signal recorded by SZS has always been in a sharp pulse state, indicating the development of underground rock fractures, which is not conducive to stress accumulation and may be the transition zone of the energy transmission process.

Conclusions and Discussion
In the research on the underground structure of the epicenter of the 2015 Azuoqi earthquake with Ms5.8 in Inner Mongolia, on the one hand, it is studied according to the survey line profile, while on the other hand, it is studied one by one according to the stations. Using the data from different time periods before and after the main earthquake, it can be clearly seen that the underground structure has changed significantly and that the change characteristics of the underground structure under the stations at different locations from the epicenter before the main earthquake are different. This difference corresponds well to the low-frequency abnormal signal before the earthquake. The details are summarized as follows: In the three seismic stations, WUH, SZS and DSM, near the Ms5.8 earthquake in 2015, the low-frequency abnormal signals extracted from the recorded seismic waveforms showed completely different forms before the earthquake, and the station closest to the epicenter, WUH, showed unstable fluctuation forms about 100 h before the earthquake. This may be the manifestation of the nucleation process in the metastable state of the source.
In the process of rock fracture and nucleation beneath WUH, due to the increase in stress concentration, the rock fissures are squeezed so that the low-frequency signal is blocked and cannot be better transmitted to the position below DSM. The low-frequency signal recorded by SZS has always been in a sharp pulse state, indicating the development of underground rock fractures, which is not conducive to stress accumulation and may be the transition zone of the energy transmission process.
According to the low-frequency abnormal signal forms recorded by SZS, WUH and DSM, it is judged that the underground stress transmission direction is from south to north and accumulates below WUH until the main earthquake occurs. From the change time of the underground structure, the results show that 64052 equivalent to SZS is the station with the earliest underground structure change, followed by 15597 equivalent to WUH, and, finally, 15670 corresponding to DSM, which is consistent with judging the stress change direction of an underground structure according to the form of low-frequency abnormal signal.
Through the retrospective study of the 2015 Azuoqi Ms5.8 earthquake in Inner Mongolia, it is found that the underground structures near the epicenter changed significantly before the main earthquake, and this change even appeared 1-2 years before the earthquake. The low-frequency signal of a seismic wave has not been studied with a completely matched time period on the time scale, so it is impossible to state whether it changed during a longer time before the main earthquake. In particular, the fluctuation of the low-frequency abnormal signal before the earthquake at WUH may be the sub-instability state before instability and may be the extraterrestrial manifestation of a stick slip process in the sub-instability state.