The 2022 Goesan earthquake of the moment magnitude 3.8 along the buried fault in the central Korean Peninsula

On October 28, 2022, a moment magnitude (Mw) 3.8 earthquake occurred in Goesan, South Korea, typically characterized as a stable continental region. Herein, we analyze 42 earthquakes, including the Mw 3.8 earthquake, the largest foreshock (Mw 3.3), which preceded the mainshock by 17 s, and the largest aftershock (Mw 2.9). The primary aim of this study is to identify interactions among the seismic events. To this end, we utilized the permanent seismic networks with the closest station at 8.3 km from the epicenter, and the temporary network deployed eight hours after the mainshock’s occurrence. Relocation results delineate that the mainshock occurred at the southeastern tip of the hypocenter distribution of three foreshocks, trending west-northwest–east-southeast. The aftershocks form an overall spatially diffused seismic pattern that propagates toward both ends of the inferred lineament in the downdip direction. The rupture directivity of the mainshock, along with wave-form similarity across the mainshock and foreshocks, confirms the inferred geometry, corresponding well with the focal mechanisms of the mainshock and the largest foreshock. We demonstrate that the change in Coulomb failure stress (ΔCFS) by the largest fore-shock was positive where the mainshock occurred and that the mainshock generated ΔCFS capable of triggering the propagation of the aftershocks.


Introduction
On October 28, 2022, at 23:27:50 UTC, a moment magnitude (Mw) 3.8 earthquake occurred in the town of Goesan in the central region of South Korea (Fig. 1); there have been no reports of casualties or significant damage caused by the mainshock.The Korea Meteorological Administration (KMA) reported its local magnitude as 4.1.This region had Highlights • The 2022 Goesan earthquake (Mw 3.8) was preceded by the Mw 3.3 event in ~ 17s.• This study determined how one seismic event influenced the next event in stable continental regions.
• The Mw 3.3 foreshock generated enough stress to trigger the mainshock along the buried fault.
The Korean Peninsula is located in stable continental regions with a relatively uniform ambient stress field in South Korea (Fig. 1a).The direction of the maximum horizontal compressional stress axis in inland South Korea is east-northeast-west-southwest (Fig. 1a).The region of the epicenter was within the Okcheon Belt, a geological province in South Korea, approximately 1,000 km from the closest Nankai Trough in Japan (Fig. 1a and b).
Geological faults, including the Jogok Fault (JF), are located in the region (Fig. 1c).The JF is subparallel to the strike of the mainshock's nodal plane, with a maximum deviation of 10° (Fig. 1c).The main sense of slip across the JF is left-lateral (Park 2005), which is consistent with the focal mechanism of the mainshock (Fig. 1c).However, the lack of seismicity between the surface and subsurface at a focal depth of ~ 10 km limits further analysis.
Notably, the mainshock was preceded by a Mw 3.3 foreshock that occurred at 23:27:33 on October 28, 2022 (UTC), approximately 17 s before the mainshock (events 3 and 4 in Fig. 2).Permanent stations recorded the largest foreshock's S-wave coda interfering with the P wave of the mainshock from a distance greater than ~ 135 km (Fig. 3).The Korea Institute of Geoscience and Mineral Resources deployed four temporary broadband stations within 10 km from the epicenter of the Mw 3.8 mainshock (see Fig. 1c).These stations began operation at 8:00:00 on October 29, 2022 (UTC), approximately eight hours after the mainshock occurred (see the dashed line in Fig. 2a) and recorded the largest aftershock of Mw 2.9 (event 30 in Fig. 2a).
This study utilized data from both the temporary and permanent stations to enhance constraints for locating hypocenters, complementing previous studies on the Goesan earthquakes (Hong et al. 2023;Kim et al. 2023;Sheen et al. 2023).We delved into possible interactions among the foreshock, mainshock, and aftershocks, precisely relocating hypocenters and examining whether a coseismic static stress model can account for the associated successive seismicity.
Fig. 2 Chronology of the earthquake sequence.a Period from October 28, 2022 to November 6, 2022.Red and grey circles represent events included in the main fault and others, respectively.A dashed line indicates the time when four temporary stations were deployed.b A magnified plot for a two-hour window, including events 1-15, as marked with a green tick on the time axis of (a) This effort could contribute to a better understanding of earthquake triggering and interaction in stable continental regions.

Earthquake detection and relocation
We detected small-magnitude earthquakes using the template-matching method based on the list of 30 events provided by the KMA.From October 28 to November 10, 2022, vertical components of seismic waveforms were band-pass filtered at 2-14 Hz for time window durations of − 1 s to 3 s in terms of P wave at three stations, CHJ3, DKSA, and CGAA (Fig. 1c).We used waveforms with a mean correlation coefficient greater than 0.3.
Subsequently, P and S wave arrival times in the vertical and transverse components, respectively, were measured manually at 11 stations including the four temporary stations (Fig. 1c), to enhance relocation accuracy and decrease the minimum magnitude of locatable events.Manual measurement was Fig. 3 Interference between the Mw 3.3 foreshock (event 3) and Mw 3.8 mainshock (event 4). a Vertical record sections of event 3 and event 4. Their difference in origin time is 17 s (Table 1).The trace is filtered in 2-8 Hz by acausal Butterworth of four poles.P and S wave arrivals are guided by solid and dashed lines, respectively.Faint red and blue lines are events 3 and 4, respectively.The S waves of event 3 and the P waves of event 4 cross at ~ 135 km.The amplitude of the trace is normalized.Names of the stations are on the right, and the peak-to-peak amplitudes in μm/s are within parentheses.The azimuth of the stations ranges from 90° to 170°.b Map of the station used in (a).A box in (b) indicates the region of Fig. 1c Vol.: (0123456789) employed to prevent cyclic-skipping and undesirable effects resulting from varying durations (Bachura and Fischer 2019).The initial hypocenters were determined using bayesloc (Myers et al. 2007) with a 1-D velocity model for the Korean Peninsula based on the model described by Kim et al. (2022) originating from Kim and Kim (1983) (Table S1).
The double-difference method was applied to enhance the relative hypocenters based on the two sets of differential times (DT) for the P and S waves (Waldhauser and Ellsworth 2000).The first set comprised differential travel times of event pairs from P and S arrival picks (DT 1 ) while the second set was based on waveform cross-correlation (DT 2 ) after the waveform was interpolated to 1,000 samples per second.In the waveform cross-correlation, we used a short window from − 0.015 s to 0.015 s in terms of P or S arrival picked manually.We only used the DT 2 if 01s to prevent cyclic-skipping and unintended effects caused by various durations.The 1,921 DT 1 helped to determine the overall distribution of the relative hypocenters; the 1,591 DT 2 finetuned the hypocenters.
Finally, the hypocenters were shifted by the absolute hypocenter of the largest aftershock (event 30), presumed to have the highest accuracy as it was recorded by the temporary stations.We estimated the uncertainty of the relative hypocenters through a bootstrapping approach.A total of 370 pairs were considered (Fig. S1 in the supplementary materials), with 70% (i.e., 259 pairs) randomly selected without duplication.The perturbed relocated hypocenters were shifted by the unperturbed absolute hypocenter of event 30, and this process was repeated 1,000 times.

Source spectrum fitting
We estimated the moment magnitude (Mw) and corner frequency (f C ) of earthquakes by fitting the S wave source spectrum with the Brune model (Brune 1970(Brune , 1971)).The geometrical spreading and intrinsic attenuation were corrected to displacement source spectra, and the best-fitting model for the spectra was determined through a grid-searching scheme with an adopted range of domains and spacing (see details in Text S1 in the supplementary materials).The searching domain was set near the initial values of Mw and f C with the Snoke method (Snoke 1987; see details in Supplementary Text S2).
For the mainshock (event 4), we used data from 8 stations in distance of 116-325 km (see the list of stations in Table S2 in the supplementary materials).For the largest aftershock (event 30), we used more stations, up to 17, in distance of 39-397 km, because it was free from coda interferences caused by preceding large-sized events.However, for the largest foreshock (event 3), we had no other choice but to use the closest single station, CHJ3, at a distance of 8 km, due to the S coda interfering with the P wave of the mainshock (Fig. 3).This single measurement limited the reliability of the estimated f C for event 3.
We employed the single-station approach using CHJ3 to estimate moment magnitudes for 40 earthquakes, which were events 1 to 3, 5 to 29, and 31 to 42.However, for event 30, we utilized 17 stations owing to the large amplitudes of the event, making it the largest aftershock.

Focal mechanism determination
We utilized the Time Domain Moment Tensor inversion tool (TDMT) (Dreger 2003;Minson and Dreger 2008) to derive moment tensor solutions for the Mw 3.8 mainshock (event 4) and Mw 3.3 foreshock (event 3), based on the velocity model of Kim et al. (2011) (Fig. S2 in the supplementary materials).However, for the foreshock (event 3), we only applied data from three stations (CHS, SND, and TJN) within limited time windows because the S wave of the foreshock and P wave of the mainshock interfered at distances > 135 km (see Fig. 3).For the mainshock (event 4), we applied data from 11 stations spanning a distance of 72-306 km.The applied passbands for the mainshock and foreshock were 0.04-0.10Hz and 0.09-0.20 Hz, respectively.The passband for the mainshock was tuned by trial and error to minimize the interference effect of coda from the largest foreshocks.The centroid depth was assessed every 1 km between 1 and 17 km.(Fig. S2 in the supplementary materials).The focal mechanism of events 3 and 30 were determined based on P wave polarity in the vertical component, using focmec (Snoke 2003) (Fig. S3).

Rupture directivity estimation
We used unilateral and bilateral rupture models to assess an azimuthal dependency of f C .The reciprocal of the f C is the apparent source duration Δt (Brune 1970(Brune , 1971)).In the unilateral model, the rupture propagated toward only one side laterally.In the bilateral model, the rupture propagated toward both sides with equal length and rupture velocity.The apparent source duration, Δt, can be modeled using Eqs.
The AZ and STK stand for the azimuth to a station and the strike of the fault, respectively.The AZ is determined from the locations of the event and the station.Equations ( 1) and ( 2) assume that the fault has a lateral rupture direction with no vertical component.The A describes the fault's lateral width, L, with the relation A = L∕ , where represents the shear wave velocity (Cesca et al. 2011).Similarly, the B is defined as B = t r + L∕V R , where t r is the rise time, and V R is the rupture velocity.
More general setups for rupture directivity, such as asymmetric bilateral rupture or involving the vertical component of the rupture direction, can be found in Cesca et al. (2011) and Park and Ishii (2015).However, we applied Eqs. ( 1) and (2) to our analysis, simplifying the fault rupture to consist solely of lateral directions.We aimed to determine the positive unknowns A and B through least-square linear inversion, using the four possible STK values.These values include opposing directions of a nodal plane.If either A or B becomes negative, we attempted another value for the STK.

Evaluation of change in the Coulomb failure stress
We also evaluated whether the changes in coseismic static stress of one earthquake induce subsequent earthquakes based on the concept of the change in the Coulomb Failure Stress (ΔCFS) (Lin and Stein (1) Ree et al. 2021;Kim et al. 2022).A positive ΔCFS indicates that the stress change causes the following event to have the assumed strike, dip, and rake.Conversely, a negative ΔCFS means that the stress change prevents the following event from occurring.The region with the negative ΔCFS is typically called the "stress shadow" (Freed 2005).We calculated the second-order tensor stress fields changed by the Mw 3.3 foreshock (event 3) and the Mw 3.8 mainshock (event 4) using the Coulomb 3 program (Lin and Stein 2004).
For the mainshock (event 4), we assigned the "source fault" mechanism with the strike, dip, and rake values of 105°, 87°, and 13°, respectively (listed in Table 2), obtained from moment tensor inversion.These parameters were also applied to the source fault for the largest foreshock (event 3), instead of employing its own moment tensor solution.The moment tensor inversion for event 3 employed only three stations (TJN, SND, and CHS) and consequently provided a solution with low variance reduction (54.8%) and a relatively high nondouble-couple component (34%) (see the solution shown in Fig. S2).
The fault size was simplified as a square with widths of 0.18 km for event 3 (Mw 3.3) and 0.35 km for event 4 (Mw 3.8).The latter was based on f C estimated in the source spectrum fitting.The f C can be converted to the radius r 0 of a circular rupture model using 0.21 ∕f C , where β is the shear wave velocity (Madariaga 1976).We set β as 3.5 km/s.The size of square was determined to match the area of the corresponding circle.Changing the shape of rupture model was intended to simplify the calculation of Coulomb 3, which does not carry a physical meaning to make a change in interpretation.For event 3, as we were unable to obtain a reasonable f C , the rupture area (r 0 of 0.1 km) was derived from the assumption that the earthquake has a stress drop of 30 MPa.This value is equal to the rounded-off stress drop of event 4 determined using its f C .
The final procedure involved conversion from a stress tensor to a scalar, ΔCFS, in terms of the targeted focal mechanism defined by strike, dip, and rake, designated the "receiver fault."We set Young's modulus, Poisson ratio, and effective frictional coefficient as 80 GPa, 0.25, and 0.4, respectively.

Moment magnitude
We denote moment magnitudes that were estimated through the source spectrum fitting as Mw B ('B' stands for Brune) and those estimated through the TDMT inversion as M W TDMT .The 42 earthquakes have Mw B ranging from 0.6 to 3.8 (Table 1).The Mw B of the mainshock (event 4) was same as the Mw TDMT of it (Mw 3.8).The Mw B and Mw TDMT of the largest foreshock (event 3) had a slight difference by 0.1.(Mw B 3.2 in Table 1; Mw TDMT 3.3 in Table 2).We preferred the Mw TDMT for event 3 because it was constrained by a larger number of stations, that is 3, compared to the Mw B , which was based on the single station (CHJ3).The Mw B of the largest aftershock (event 30) was 2.9.See Fig. S2 in the supplementary materials for the detailed results of the TDMT inversion for events 3 and 4. Also, see Fig. S4-S6 in the supplementary materials for source spectrum fittings of events 3, 4, and 30, respectively.The asymptotic value of the source spectrum model to 0 Hz directly refers to the seismic moment.

Corner frequency (f C ) and directivity
We observed azimuthal dependency in the reciprocal of f C across 8 stations for event 4 and 17 stations for event 30 (Fig. 4).The f C is known to be affected by the directivity of rupture propagation.The f −1 C of event 4 were distributed having the maximum at 116° in azimuth (refer to station CHS for event 4 in Table S2; Fig. 4a).That is, the f −1 C exhibited a singlelobe pattern in azimuth, implying a unilateral rupture that propagated toward the west-northwest (Fig. 4a).On the other hand, the f −1 C of event 30 had two minima at 12° and 202° in azimuth (refer to stations KSA and BGD, respectively, in Table S2; Fig. 4b).The f −1 C had a two-lobe pattern in azimuth with the two estimated maxima (Fig. 4b), which indicates that the rupture had elongated apparent durations toward both of the directions, the west-northwest-east-southeast.This describes a bilateral directivity along the fault strike.It is noteworthy that the estimated directivity constrains the fault plane of the two nodal planes derived from the focal mechanism, not depending on the spatial distribution of relocated hypocenters.
For event 4 (mainshock), the representative f C was 3.61 Hz, and for event 30, it was 9.0 Hz.These values were obtained as the reciprocal of the constant term (B −1 ) from the best-fitting sinusoidal curves in Fig. 4a and b, respectively.The r 0 of the mainshock, obtained from the representative f C , was 0.20 km.The corresponding stress drop Δσ was calculated with Δ = 7M 0 ∕16r 0 3 , where M 0 is the seismic moment (Eshelby 1957).The Δσ of the mainshock was 29 MPa.Similarly, the r 0 and Δσ of event 30 were 0.08 km and 20 MPa, respectively.

Fault attributable to the mainshock
The focal mechanism solutions of the Mw 3.3 foreshock (event 3) and Mw 3.8 mainshock (event 4) had coherent nodal planes, which are west-northwest-east-southeast vertical strike-slip faults with a dominant double-couple composition of the mainshock by 97% and the foreshock by 66% (Figs.5a and S2; Table 2).Due to the small number of stations, we were unable to interpret the high non-double-couple composition of event 3 (Fig. S2 in the supplementary materials).Relocated hypocenters of events 1, 2, and 3 and the mainshock (event 4) constitute a plane with a horizontal direction that is consistent with the west-northwest-east-southeast nodal plane in the focal mechanisms of event 3 and the mainshock (Fig. 5a).Also, this nodal plane was consistent with the directivity observed in f C to west-northwest (285°; Fig. 4a).Lateral and vertical extension of this plane covered smaller events.All the events within this extended plane were labeled with 'm' standing for 'main fault' in Table 1.They were marked as red color in Fig. 5 to distinguish the other events which were not included in the main fault.The westnorthwest-east-southeast trend was also confirmed in previous studies (Hong et al. 2023;Kim et al. 2023;Sheen et al. 2023).Double-couple component of the focal mechanism for event 3 was similar to that constrained with P wave polarity (Table 2).
Of the two nodal planes, the largest aftershock (event 30) could be attributed to a west-northwest-east-southeast plane.This plane is subparallel to the main fault (Figs. 5 and S6), and the rupture directivity of event 30 (Fig. 4b) also aligns the fault plane trending west-northwest-east-southeast.The dip of the aftershock (68°) was slightly gentler than that of the mainshock (87°; Table 2).The P-axis of the faulting mechanisms observed in Fig. 5a was similar to the direction of the maximum horizontal principal stress of the crust in South Korea suggested by Soh et al. (2018).

Changes in the Coulomb failure stress
We examined the spatiotemporal development of the Goesan earthquakes in terms of the coseismic static stress transfer using ΔCFSs.Figure 6 presents the ΔCFSs for the Mw 3.3 foreshock (event 3) and Mw 3.8 mainshock (event 4).The ΔCFS of the mainshock accumulated from the foreshock.The receiver fault was set to the strike, dip, and rake of the mainshock to assess the tendency for planar extension along the "main fault." The general patterns in ΔCFS include the presence of four positive lobes in the horizontal slice at a depth of 12.8 km (i.e., the focal depth of the mainshock; Fig. 6a; Table 1).Of these, two lobes were in the direction of both strikes, while the others were perpendicular.In detail, the mainshock (event 4) is within the positive zone induced by the largest foreshock (event 3) (Fig. 6a), indicating that the coseismic stress of the foreshock induced the mainshock.This remained valid when we considered the uncertainty of the location and the 3-D geometry in the vertical cross-section (Fig. 6b Table 2 Source properties of events 3, 4, and 30 The preferred nodal plane is provided first for the strike, dip, and rake.Focal mechanisms for events 3 and 4 were determined using Time Domain Moment Tensor inversion (TDMT, Dreger 2003;Minson and Dreger 2008; the first and third rows).The focal mechanisms of events 3 and 30 were determined based on P wave polarity (Snoke 2003; the second and fourth rows).DC and CLVD indicate components of double couple component and compensated linear vector dipole, respectively   and c).Most of the perturbed epicenters by bootstrapping were included in the positive zone (Fig. 6b); the perturbed depths did not enter the negative zone in the cross-section (Fig. 6c).Similarly, Fig. 6d-f shows the horizontal and vertical slices of the ΔCFS with the perturbed locations of the aftershocks of the main fault (Table 1).The aftershocks were located in the positive zone in the vertical cross-section (Fig. 6f).

Earthquakes following the static stress increase
We confirmed that the change in the coseismic static stress of the Mw 3.3 foreshock explained the occurrence of the mainshock (Fig. 6a).The occurrence of aftershocks was well understood based on the   1).e Modified version from (d), like (b), with the perturbed locations and the iso-contour of 0 MPa.f Vertical crosssection of (d) Vol:. ( 1234567890) cross-sectional view of the ΔCFS in the framework of the coseismic static stress (Fig. 6f).Events 7, 9, 17, 18, 26, and 31 were located near the edge of the rupture zone where the ΔCFS was highest (up to 0.15 MPa; Fig. 6f); although the upper limit near the rupture zone is sensitive to heterogeneous slip distribution, this was not determined in this study.Events 13, 25, and 42 were located relatively far from the rupture zone of the mainshock, extending the seismic zone to ~ 1.7 km horizontally and ~ 0.8 km vertically.
Notably, this discussion presumes that the aftershocks had a similar focal mechanism as the mainshock.
4.2 The Mw 2.9 aftershock in the stress-shadow zone The occurrence of the largest aftershock (event 30) in the negative ΔCFS zone is not simply explained in the ΔCFS (Fig. S7 in the supplementary materials).The applied stress tensor was the same as that for Fig. 6, however, the receiver fault was changed to the strike, dip, and rake of event 30.The epicenter of event 30 was at positive ΔCFS at the time of the largest foreshock (Fig. S7a in the supplementary materials); however, became negative at the time of the mainshock (Fig. S7b in the supplementary materials).Unlike other aftershocks in Fig. 6, event 30 does not conform to the 3-D geometry, mainly due to its depth of 12.89 km, a value close to the mainshock depth of 12.83 km (Table 1).However, we noticed that another study suggested that event 30 (event '25' in their article) was located in the positive ΔCFS zone (Kim et al. 2023).The crucial distinction is in the depth of event 30, which differs from that of event 4 in their study by 0.6 km (hypocentral) or − 1.5 km (centroid), while our results do not show this discrepancy.We also suggested alternative contributions to the occurrence of event 30: smaller events near event 30, such as event 16, could indicate a possible interaction between events 4 and 30 (Fig. 5).However, it is important to note that our understanding is limited, and further studies are needed to explore and validate these relationships.An oversimplified uniform rupture could be involved.Also, secondary mechanisms resulting from time-dependent changes in the frictional state can also affect the stress field (Dieterich 1994;Kilb et al. 2000;Freed 2005  be applied to cross-check the relocation outcomes (Schaff and Richards 2004;Kim et al. 2010Kim et al. , 2022;;Son et al. 2015Son et al. , 2020;;Sheen 2021).We assessed waveform similarity by computing cross-correlation coefficients (CCs) for the high-pass-filtered at 1 Hz in a 0-6 s window from the time of origin to include P wave arrival time and S coda.We used waveforms recorded by the north-south component of station CHJ3, positioned with an azimuth of 93°, nearly parallel to the transverse direction of ten events (1-6, 14, 16, 17, and 30).The selected events encompassed the three of foreshocks (events 1-3), relatively largemagnitude events (3-6, 14, 16, 17, and 30).
The CCs were grouped into two clusters: C1 comprised events 1-4, and 30, while C2 comprised events 5, 6, 14, 16, and 17 (Fig. 8a).Although the CC of event 6 was low and events 5 and 6 had negative polarity, we included them in C2 due to shared monochromatic features observed in each spectrum, similar to those present in the spectra of other events (14, 16, and 17) in C2 (Fig. 8b and c).The low CC of event 6 was due to a slightly low-frequency peak at 27 Hz (Fig. 8c).The monochromatic feature with an average peak at ~ 29 Hz, after the S wave arrival (Fig. 8b and  c) could result from a trapped mode in a fault zone (Li et al. 1990;Ben-Zion 1998;Li and Vernon 2001;Ross and Ben-Zion 2015).However, the exceptions observed in events 17 and 30, along with their opposite positions toward each group, are contradictory to this possibility.
Instead, we observed the spatial coherence among hypocenters in C1 and C2, except for events 17 and 30 in C2 and C1, respectively (see also Fig. 5).Events 1-4, constituting the main fault, were included in C1, whereas events 5, 6, 14, and 16 from C2 were not associated with the main fault.It is also noteworthy that C1 includes the foreshock-to-mainshock sequence (4 of 5), while C2 predominantly comprises aftershocks.This observation may suggest the presence of structural or frictional heterogeneity at a local scale, indicating different faulting properties for each cluster within a relatively confined region around the main fault.The temporal distinction between foreshocks and aftershocks could further emphasize a possible mainshock-induced coseismic changes in the stress field of the confined region.This hypothesis needs additional waveform similarity analyses with more stations; however, we were not able to observe high-frequency waveform similarity at stations DKSA and CGAA, located relatively far from the epicenter to CHJ3, due to attenuated waveforms.

Coefficient A and stress drop (Δσ)
The A in Eq. ( 1) is described as A = L∕ where L is the fault's lateral width (Cesca et al. 2011).When we substituted A and β with 0.078 s and 3.5 km/s, respectively, L became 0.27 km, which was similar to the assumed square-shaped fault's width of 0.35 km in the calculation of the ΔCFS.
The Δσ of 29 MPa for the mainshock (event 4) was larger than the average in the Korean Peninsula 1-10 MPa for Mw larger than 3.5 of Rhee and Sheen (2016).Also, the Δσ was slightly larger than those of the 2016 Mw 5.5 Gyeongju earthquake, 11 MPa of Chai et al. (2020) based on the code wave, 13 MPa of Son et al. (2018) based on the Snoke method of S wave (Snoke 1987), and 23 MPa of Uchide and Song (2018) based on the finite fault inversion with the empirical Green's function method.Uchide and Song (2018) reported that the local Δσ on the fault was up to 62 MPa.Caution must be applied when comparing the absolute values of stress drops that were estimated using different methods (Shearer et al. 2019;Abercrombie 2021;Calderoni and Abercrombie 2023).

Conclusions
The Mw 3.8 Goesan earthquake and its fore-and aftershocks represent a case with well-determined seismic properties owing to the dense national seismic networks and their geographical location in the central region of South Korea.Moreover, the rapid deployment of the temporary stations to provide azimuthal coverage within 10 km and one atop the mainshock epicenter improved the accuracy of the foreshock, mainshock, and aftershock locations.Indeed, we shifted the hypocenters in terms of the coordinate of the largest aftershock that was determined through full utilization of the aftershock array.
The focal mechanisms of the mainshock and its largest foreshock (Mw 3.3; event 3) and aftershock (Mw 2.9; event 30) indicate a left-lateral strike-slip, of which the P-axis corresponds with the maximum horizontal axis of the current tectonic stress.The west-northwest-east-southeast trending strikes of the faulting mechanisms agree with the spatial distribution of the relocated hypocenters and the rupture directivity inferred from the azimuthal variation in corner frequency estimates obtained through source spectrum fitting.Based on the source spectrum fitting, we estimated the fault attributable to the mainshock.Its area of approximately 0.13 km 2 and the stress drop of 29 MPa.
The seismic properties identified the interaction between the mainshock and foreshock that preceded the mainshock by 17 s, as well as the development of most of the aftershock sequence.The mainshock occurred at the edge of the rupture zone, where the ΔCFS was positive.This indicates that the foreshock triggered the mainshock.Following the mainshock, the seismic zone propagated outside of the rupture zone.The cross-sectional view of the ΔCFS indicated that the aftershock distribution scattered horizontally and propagated downward from the rupture zone of the mainshock.

Fig. 1
Fig. 1 Distribution of the Goesan earthquakes and tectonic setting.a Maximum horizontal crustal stress with the tectonic setting around the Korean Peninsula.Stress is obtained from the World Stress Map (Heidbach et al. 2018).Plate boundary is from Bird (2003).Larger and smaller dashed boxes indicate regions of (b) and (c), respectively.b Stations used in this study.Inset shows the geological province, Okcheon Belt (OB), including the epicenters.The box indicates the region of (c).c Magnified map including the stations used to determine

Fig. 4
Fig. 4 Directivity of the mainshock and the largest aftershock.A reciprocal of f C is plotted as a function of azimuth (AZ).a Mainshock (event 4) with unilateral rupture direction at 285° (west-northwest).b Largest aftershock (event 30) with bilateral rupture directions at 111° and 291° (westnorthwest-east-southeast).The equations of the fitting

Fig. 6
Fig. 6 Change in the Coulomb Failure Stress (ΔCFS) of the Mw 3.3 foreshock (event 3) and Mw 3.8 mainshock (event 4).The strike, dip, and rake of the receiver fault are 105°, 87°, and 13°, respectively.a ΔCFS of the Mw 3.3 foreshock.b Modified version from (a).Perturbed locations of the Mw 3.3 foreshock on the ΔCFS simplified by an iso-contour of 0 MPa.c Vertical cross-section of (a).d ΔCFS of the Mw 3.8 main- implications from fore-and aftershocks Foreshocks hold significant interest due to their potential role in triggering a mainshock; thus, it is crucial to comprehend the spatiotemporal features associated with their occurrence.Tape et al. (2018) found a notable observation that the high-frequency precursor in 2-8 Hz and the signal in the very low frequency band (0.02-0.05 Hz) are coincident within tens of seconds before the P arrival of the earthquake in central Alaska.However, we are unable to observe the high frequency precursor, following the same process asTape et al. (2018).Traces filtered in the 2-8 Hz range and their log-scaled envelope in the top and bottom of Fig.7did not display any anomalous signal.The inspection for the wave in the very low frequency band is limited because of the short length of the window before the P arrival of the mainshock.Waveform similarity can provide insights into the spatial feature of hypocenter distribution and

Fig. 7
Fig. 7Traces filtered in a high-frequency band.The station CHJ3's vertical trace is filtered in 2-8 Hz at the top.P3 and S3 indicate the P and S wave arrival times of the Mw 3.3 foreshock (event 3); P4 and S4 are for the Mw 3.8 mainshock (event 4).The trace at the bottom is a log-scaled envelope from the top image to show the possible existence of a small amplitude precursor, which is absent in this figure

Fig. 8
Fig. 8 Waveform similarity.a Cross-correlation (CC) matrix between the waveforms of station CHJ3 for the events labeled in Fig. 1.The events belonging to the main fault are marked by an asterisk.Cluster C1 comprises events 1-4, and 30, and C2 events 5, 6, 14, 16, and 17. b Transverse component aligned by the S arrival in a window from − 0.5 to 0.5 s.The waveforms