Relation between earthquake swarm activity and tides in the Noto region, Japan

The Noto region in Japan has been experiencing earthquake swarm activity since mid-2018, with repeated rises and falls in activity. Crustal deformation (expansion and uplift) observed there since the end of 2020 has been connected to crustal fluids. Studies in other regions have suggested that tides are related to earthquake swarm activities associated with crustal fluids. Therefore, we investigated whether tides are also involved in the Noto earthquake swarm activity. Our results suggest a tidal correlation only at greater depths in the southern part of the analyzed area (‘region Sd’). There, we inferred that an increase in pore fluid pressure caused by the inflow of deep fluids may have led to a decrease in fault fracture strength, making the local seismicity relatively susceptible to the effects of tidal forces. The significantly high value in region Sd of the scaling parameter b of the Gutenberg–Richter law (describing the earthquake magnitude–frequency distribution) and observed crustal deformation are consistent with this interpretation.


Introduction
In the Noto region of Japan, increased seismicity has persisted since around 2018, with the area of seismicity broadening in December 2020 and further still in July 2021 (Japan Meteorological Agency (JMA) 2023a).Additionally, crustal deformation (expansion and uplift) has been observed there since around the end of 2020 (Nishimura et al. 2023).The presence of crustal fluids has been suggested as a factor behind this swarm activity (Nakajima 2022;Amezawa et al. 2023).Elsewhere, the well-known Matsushiro earthquake swarm activity that started in 1965 was similarly attributed to crustal fluids (e.g., Yoshioka et al. 1970;Ohtake 1974;Yoshida et al. 2002), as well as to tides (e.g., Iwata and Nakanishi 1998).Miguelsanz et al. (2021Miguelsanz et al. ( , 2023) ) examined the behavior of magma in the Canary Islands based on the relationship between earthquake swarms in volcanic areas and tides.Han et al. (2018) examined the relationship between deep long-period earthquakes associated with magmatic and/or fluid behavior beneath Mt.St. Helens and tidal stresses, finding a correlation with compressional volumetric tidal stress.Further investigation of the relationship between seismic swarm activity and tides may therefore help elucidate the mechanism(s) behind seismic swarm activity.Therefore, in this study, we analyzed the relationship between earthquake swarm activity and tides in the Noto region.

Data
The earthquake swarm activity in the Noto region can be spatially divided into four clusters in the northern, southern, eastern, and western areas (regions N, S, E, and W, respectively) (Fig. 1a).Furthermore, in region S, activity can be divided two depth clusters: Ss (hypocenters shallower than 14 km) and Sd (deeper than 14 km) (Fig. 1b).Thus, we examined five clusters of earthquake swarm activity: regions N, Ss, Sd, E, and W. From the JMA unified hypocenter catalog, we extracted 6,202 high-precision hypocenters of events with magnitude M ≥ 1.3 (Fig. 1c) and flags K, k, or A (where K indicates a fully reviewed hypocenter, k a simply reviewed hypocenter, and A an automatically processed hypocenter; see JMA 2023b for details) occurring between 1 January 2018 and 31 December 2022 in the area within 37.41-37.) and at hypocenter depths of less than 20 km.
Several methods for estimating the threshold magnitude M th have been proposed (Woessner and Wiemer  (3528, 1521, 630, 203, and 320 in regions N, E, W, Ss, and Sd, respectively).Symbol colors in c-e are the same as in b 2005): For example, the MAXC method (Wiemer and Katsumata 1999) define magnitude with the maximum curvature of the frequency-magnitude distribution as M th , and this generally matches the magnitude bin with the highest frequency of events in the non-cumulative frequency-magnitude distribution (Woessner and Wiemer 2005).Another method, the EMR method (Woessner and Wiemer 2005), estimates M th based on the Gutenberg-Richter (G-R) law (Gutenberg and Richter 1944) using the entire data set, including the range of magnitudes reported incompletely.On average, M th determined by the EMR method corresponds to M th determined by the MAXC method plus 0.2.In each region analyzed in this study, the magnitude of bin with the highest frequency of events (that is, corresponding to M th of the MAXC method) ranges from 0.6 to 0.9, corresponding to 0.8 to 1.1 of M th of the EMR method.We also visually confirmed the shape of the frequency-magnitude distribution, and to be on the safe side, adopted 1.3 as M th common to each region (Fig. 1c).
During this period, the number of earthquakes first began to increase in region Ss in mid-2018, region Sd became seismically active in December 2020, and the activity then expanded to regions W, N, and E in 2021 (Fig. 1d, e).The largest earthquake ( M 5.4) occurred in region E on 19 June 2022 and was produced a maximum seismic intensity of 6 Lower on the JMA Seismic Intensity Scale (https:// www.jma.go.jp/ jma/ en/ Activ ities/ intta ble.html).We note that on 5 May 2023 (after our analysis period), a larger earthquake of M 6.5 (maximum seismic intensity 6 Upper) occurred in region E, and the subsequent expansion of the active region in the ocean area has blurred the boundary between regions N and E. Added on February 8, 2024: after the paper was submitted (December 11, 2023), on January 1, 2024, M 7.6 earthquake (maximum seismic intensity 7) occurred with the rupture starting point at 37.5°N near the boundary between regions N and E.
Calculation of theoretical tidal responses requires the hypocentral location, occurrence time, and fault parameters.For the hypocentral location and occurrence time, we used catalog information.For the fault parameters, we assumed that reverse faults with a northeasterly strike and southeasterly dip predominated (i.e., we assigned to all events a strike of 45°, a dip of 45°, and a rake of 90°).This is consistent with the regional tectonic stress field (Terakawa and Matsu'ura 2022), the focal mechanisms of major earthquakes, and the hypocentral distribution (Amezawa et al. 2023).

Theoretical tidal response
We calculated theoretical tidal responses on the faults using the "TidalStrain.2"Fortran software (Hirose et al. 2019).We investigated changes in the following 12 components as tidal indices: volume strain ΔV at the hypocenter, shear stress Δτ, normal stress Δσ, and the Coulomb failure function �CFF = �τ + µ ′ �σ for assumed apparent friction coefficients µ ′ of 0.1, 0.2,…, 0.9 (denoted ΔCFF (0.1) , ΔCFF (0.2) ,…, ΔCFF (0.9) , respectively) on the assumed fault plane.If a highly anisotropic fault hosted fluid-filled fractures aligned along the fault zone, µ′ = µ 0 (1 − B) (Houston 2015), where µ 0 is the static friction coefficient and B is Skempton's coefficient, which indicates the ratio of the change in pore fluid pressure to the change in confining pressure acting on the rock under undrained conditions.Accordingly, if the response of pore pressure to Δσ is strong (i.e., B ∼ 1 ) at the fault plane, µ ′ becomes small and the contribution of Δτ to ΔCFF becomes dominant.In the case of ΔV and Δσ, we defined expansion/dilatation as positive and contraction/compression as negative, and we defined Δτ and ΔCFF as positive when they promoted fault slip.We assigned phase angles of − 180° and + 180° to the local tidal minima before and after an event, respectively, and 0° to the tidal maximum between these two local minima.The phase angle ψ at the earthquake occurrence time was estimated by linear interpolation over the time interval between − 180º and 0° or between 0° and + 180° (Fig. 2).The tidal level L at the earthquake occurrence time was defined as the strain or stress value with a positive or negative sign relative to the zero line.

Sensitivity to tidal phase angle: Schuster test
The possibility of earthquakes occurring at a particular tidal phase angle is commonly evaluated by using the p-value (Schuster 1897), calculated as follows: where N is the total number of earthquakes, ψ i is the phase angle at the occurrence time of the i th earthquake.An approximation of Eq. ( 1) is sufficient when N is larger than 10 (Heaton 1975).The p-value (which ranges from 0 to 1) represents the significance level for rejecting the null hypothesis that earthquakes occur randomly with respect to the tidal phase angle, such that the confidence that the null hypothesis is rejected is greater when the p-value is smaller.In general, p-values of 0.05 (Tanaka et al. 2002) are used to judge tidal correlations.
If the occurrence of an earthquake is influenced by tides, there should be a peak in tidal phase angle histograms near a phase that promotes fault slip (for example, ΔCFF of 0°).Therefore, we fit a sinusoidal curve �(ψ) to the frequency distribution of tidal phase angles as follows: (1) (2) where A 0 is the mean relative frequency, fixed at 8.333% (= 1/12 of 100%) for the adopted phase angle bin width of 30°, and A 1 and ψ p are the amplitude and phase, respec- tively, of the least-squares-fitted curve.ψ p becomes sig- nificant only when p is small because larger p-values suggest lower correlations with tides.In addition, when a peak appears at times when fault slip is suppressed, it is physically reasonable to infer a spurious correlation.

Sensitivity to tidal level: Houston plot
We also investigated the dependence of an earthquake on the tidal level (amplitude).An example is shown in Fig. 3.As a first step, we investigated the case that the tidal level is divided into two bins, positive or negative (Fig. 3b), because if the number of earthquakes in each bin is small and the number of bins is large, the estimation error becomes large and meaningful results cannot be obtained.The tidal sensitivity in this case can be simply defined as the absolute value of the difference between the N obs /N exp ratios in the two bins (corre- sponding to the slope of the blue line in Fig. 3b), where N obs is the number of earthquakes observed during the period when the tidal force takes the value within the strain/stress bin (i.e., positive or negative), and N exp is the expected number of earthquakes in the same period if there were no correlation to tides (here  4)) to the N obs /N exp values in the six bins.In this example, tidal sensitivity α = 0.5 the tidal amplitude sampled at 15-min intervals during 4 days before and after each earthquake).To test whether a correlation was statistically significant, we used a Chisquare test at a significance level of 5% to test the null hypothesis that the event occurrence rate and the tidal occurrence rate (background distribution) belong to the same population.
If the first step suggested a correlation between earthquakes and tides, we proceeded to the second step.To evaluate whether the likelihood of an earthquake increases as the tidal level increases, we used the following formula (Houston 2015) for the case that the tidal level is subdivided into six bins (Fig. 3c): where S is the intermediate value in a strain or stress bin, and α represents tidal sensitivity.If the ratio N obs /N exp (blue diamonds in Fig. 3c) of the frequency distribution of tidal values at the time of an earthquake's occurrence is proportional to the tidal stress, then a relationship with the tide is indicated.A Chi-square test with a significance level of 5% was used to test the statistical significance of the correlation.In this study, we refer to plots such as those illustrated in Fig. 3b, c as 'Houston plots' .

Declustering method
When an earthquake occurs, earthquakes such as aftershocks are likely to be induced in the surrounding area as a result of stress redistribution.When the tidal correlation of seismic activity is investigated, if such aftershocks are not excluded (i.e., declustered) in advance, apparently high tidal correlations tend to appear (Nagata et al. 2022).In general, declustering is a process that extracts independent events by excluding earthquakes that occurred closely in space and time.In this study, we analyzed data obtained by using a general declustering method (e.g., De Natale and Zollo 1986).
We applied this declustering method as follows: in each region, a series of earthquakes in which the time interval between successive earthquakes was within the threshold time t th was considered to be one group, and all but the first earthquake in each group were excluded.To obtain the declustered catalog, we performed a Kolmogorov-Smirnov test with the significance level set to 5% using the frequency distribution of the time intervals between consecutive earthquakes and the theoretical time interval for a stationary Poisson process.We finally adopted a minimum value of t th of 1, 2, 3,4,5,6,12,24,48,72,96, or 120 h if the null hypothesis "the time series after declustering is a stationary Poisson process" could not be rejected.
(4) N obs (�S)/N exp (�S) = e α�S , When the test period was set to the entire study period (1 January 2018 to 31 December 2022), there was no value of t th in any region for which the null hypothesis could not be rejected (note that the null hypothesis was rejected when the test period was the entire study period because the predominant occurrence interval was significantly different between the period of normal activity and the swarm period).Therefore, we limited the test period to the swarm period in each region (rightmost column of Table 1), defined as the period after a distinct change in slope was observed in the cumulative seismicity trend (Fig. 1e).To obtain the declustered catalog, we then determined the t th value for which the null hypothesis could not be rejected in each region (Table 1).Note that when we applied declustering using the determined t th values (Table 1) to data from other periods in each region, we were able to confirm by a Kolmogorov-Smirnov test that the null hypothesis could not be rejected except in region N.Because in region N, however, t th was 48 h (that is, the interval between earthquakes in the declustered catalog was longer than the durations of the principal tidal cycles of ~ 12 and 24 h), we inferred that the declustered catalog obtained in region N would not affect the evaluation of tidal correlation.Therefore, in subsequent analyses, we used declustered data for the entire period for each region (Table 1).

Sensitivity to tidal phase angle
Tidal phase angle-frequency histograms for each region and each tidal index are shown in Fig. 4. We obtained statistically significant tidal correlations in regions Ss and Sd (i.e., p < 0.05 in the Schuster test).However, in region Ss, the peak phase angle ψ p is approximately + 90°; by definition, slip begins to be

Table 1 Declustering results
Number of events for which the null hypothesis that 'declustered seismic activity is a stationary Poisson process' could not be rejected at the significance level of 5% after the declustering method was applied to events during the swarm period (i.e., declustering was successful).The numbers in parentheses, which show the number of earthquakes remaining after the declustering method was applied to the entire period (1 January 2018 to 31 December 2022), were used in this analysis suppressed at this phase angle.According to the Dieterich model (Dieterich 1994), which is based on the rateand state-dependent friction constitutive law, the peak of seismicity rate may appear delayed from the peak of the stressing rate.For example, in Nakata et al. (2008), which investigated the relationship between deep tectonic tremors and tides in eastern Shikoku, a peak difference of about 1-2 h was observed.Regarding region Ss, there was a delay of about 180°, or about 6 h on average, from the peak of stressing rate (tidal phase angle of − 90°) to the peak of the number of earthquakes (90°).
We preliminarily input a 12-h periodic sinusoidal wave that simulates the tidal stressing rate into the Dieterich model and performed calculations by changing several parameter values, but the maximum delay was only 3 h.Therefore, tidal correlations for region Ss are questionable.In region Sd, we obtained significant tidal correlations for ΔV, Δτ, ΔCFF (0.1) , and ΔCFF (0.2) .The peak phase angle ψ p of ΔV of ~ + 140° indicates that earth- quakes tended to occur during volumetric contraction.The peak phase angles ψ p of Δτ, ΔCFF (0.1) , and ΔCFF (0.2) , all ~ − 40°, indicate that earthquakes tended to occur at times when tidal stress accelerated fault slip.
Figure 5 shows the relationship of the tidal phase angles between pairs of the four tidal indices in region Sd for which a tidal correlation was indicated.Almost no phase difference exists between Δτ and ΔCFF (0.1) or ΔCFF (0.2) (Fig. 5e, f ) or between ΔCFF (0.1) and ΔCFF (0.2) (Fig. 5c), whereas the phase difference between ΔV and the other indices was ~ 180° (Fig. 5a, b, d).The Schuster test is sensitive to bias in the phase angle distribution but returns similar p-values for samples with a definite phase differ- ence.Therefore, these four tidal indices are inferred to show the same degree of tidal correlation.
Because the p-value of ΔCFF (0.1) was the smallest among these four indices, ΔCFF (0.1) may be the tidal index that controls earthquake swarm activity.However, because the actual fault shape may differ from the assumed fault parameters, the dominance of this index cannot be determined here.Previous studies (Fig. 7c of Japan Meteorological Agency 2023a and Fig. 1 of Amezawa et al. 2023) show focal mechanism solutions (P-wave first-motion polarity solution, CMT solution, and NIED MT solution) for the Noto earthquake swarm activity, and show that the reverse-fault event with fault parameters assumed in this study are dominant.However, the events for which focal mechanisms have been estimated are mainly in regions N and E, and for events with M < 3.2 , even the P-wave first-motion polarity solution has not been determined (note that focal mechanisms for all events with M ≥ 3.2 are not determined).Therefore, it is unclear whether focal mechanisms of small earthquakes are consistent with the assumed fault parameters.Nonetheless, because a tidal correlation is indicated even for ΔV, which does not depend on fault parameters, at minimum it can be inferred that the earthquake swarm activity in region Sd is influenced by tides.

Sensitivity to tidal level
Figure 6 shows Houston plots for the case of two bins, positive and negative.The tidal sensitivity, that is, the degree of deviation from the background distribution (indicated by the slopes of the blue lines in Fig. 6), is greater in region Sd than in the other regions.In region Sd, significant tidal correlations are indicated for ΔV, Δτ, and ΔCFF (0.1-0.5) (red frames in Fig. 6).Although the Chi-square value for ΔV was largest among the seven indices, it is not possible to determine which tidal index controls earthquake swarm activity (see "Sensitivity to tidal phase angle").No relationship with tides is suggested for earthquakes in the other regions.
Figure 7 shows the relationships between tidal phase angle and tidal level, as indicated by tidal indices ΔV, Δτ, and ΔCFF (0.1, 0.2) in region Sd.Tidal levels are generally positive for phase angles from − 90° to + 90° and generally negative for other phase angles.Therefore, the tidal components suggesting a tidal correlation in the Houston plots with two bins and those suggested by the Schuster test results ("Sensitivity to tidal phase angle") are mostly the same.Because tidal levels were not necessarily negative when the tidal phase angle was outside of the − 90° to + 90° range, the Schuster test results did not suggest any correlation with ΔCFF (0.3-0.5) .Next, to gain a more detailed understanding of the trends in region Sd, we examined Houston plots for the case of 6 bins (Fig. 8).These frequency distributions, which show the same tendencies as those with two bins (Fig. 6), also suggest a correlation with tides.However, in no case could the null hypothesis be rejected at the 5% significance level based on the Chi-square test result.
To summarize, although it cannot be said that an earthquake is more likely to occur when the tidal level is larger, the Houston plots in region Sd indicate a tidal correlation.It is not possible to determine, however, which tidal Fig. 6 Houston plots with two bins, positive and negative, for each tidal index in each region.See Fig. 3b for how to read the figure.For each region, the number of earthquakes N EQ meeting the threshold magnitude M th of 1.3 is shown along with N BG which is the number sampled at 15-min intervals during four days before and four days after each earthquake, N BG = 769 × N EQ .The first value in each panel is the absolute value of the difference between the N obs /N exp ratios in the two bins (corresponding to the slope of the blue line); the slope becomes greater as the event occurrence distribution deviates from the background distribution.The second value in each panel is the Chi-square value.In panels showing a Chi-square value of 3.8415 or higher (framed in red), the null hypothesis is rejected at a significance level of 5% index is dominantly associated with the triggering of earthquakes.

Spatial distribution of the b-value of the Gutenberg-Richter law
The G-R law, logN = a − bM , describes the earthquake magnitude-frequency distribution, and its b-value is a scaling parameter.The b-value was estimated using the maximum likelihood method (Utsu 1965;Aki 1965), b = log 10 e/ M − M th , where M is the average mag- nitude.The uncertainty of b-value was calculated from σ = b/ √ N derived by Aki (1965).Table 2 shows the G-R b-values for the entire study area and within each region.For earthquakes in region Sd with hypocentral depths below 14 km, the b-value was 1.85, considerably larger than those in the other four regions (~ 1.1) as well as the median value for all of Japan (~ 0.9; Nagata et al. 2022).Over the entire analyzed region, the b-value tends to gradually decrease from a depth of 10 km to 14 km, and below that depth it increases (Fig. 9a).Spada et al. (2013) investigated depth changes in the G-R b-value for earthquakes occurring in continen- tal crust around the world, and found that the b-value decreases with depth to about 15 km (about 25 km in Switzerland), and then increases below that depth.They indicated that the depth dependence of the b-value is due to the strength profile (differential stress dependence) of the continental crust, and that the inflection point, i.e., the depth at which the b-value begins to increase, cor- responds to the brittle-ductile transition in the crust.Because the brittle-ductile transition depth beneath the Noto region of Japan is ~ 13 km (the D90 depth; Omuralieva et al. 2012), the depth dependence of the b-value obtained here is generally consistent with the differential stress dependence reported by Spada et al. (2013).However, the average b-value in Japan at 15-20 km depth, i.e., after the b-value begins to increase below the brittle- ductile transition, is only about 0.85 ± 0.05 (Spada et al. 2013); thus, the b-value of 2.0 ± 0.1 (Fig. 9a) in this depth range in the Noto region is uniquely high.
High b-values are often observed around magma cham- bers, for example, in the active seismic swarm area off Ito on the Izu Peninsula (Wyss et al. 1997) and directly beneath volcanoes in the Tohoku region (Wyss et al. 2001).Because areas around magma chambers are highly fractured and hot (Wyss et al. 1997) and/or under a low stress field due to high pore fluid pressure (Wyss et al. 2001), it has been interpreted that large fractures cannot occur around magma chambers; as a result, the slope of the earthquake magnitude frequency distribution is high (i.e., the G-R b-value is high).
In region Sd, observations of hypocenter migration (Amezawa et al. 2023) and crustal deformation (Nishimura et al. 2023) suggest that high-pressure fluids are supplied from depth, and a tomographic analysis has suggested that beneath the Noto Peninsula, 2) , e ΔCFF (0.3) , f ΔCFF (0.4) , and g ΔCFF (0.5) .When the phase angle was in the range of − 90° to + 90°, the tidal level was generally positive, and generally negative otherwise (white areas in each panel).However, there were some exceptional events that fall into the gray area fluids originating from the mantle are being supplied to the lower crust (Nakajima 2022).In addition to the tendency toward increased b-values at greater depths, high-pressure fluids may be responsible for the uniquely high b-values in region Sd.

Possible explanation for the high tidal correlation in region Sd
Taken together, our results suggest that the high tidal correlation in region Sd results from the following: 1.In region Sd, the deeper hypocenters and higher b -values than in the other regions suggest that the differential stress in region Sd is small, based on the previously reported depth distribution of b-values (Spada et al. 2013).

Fig. 1
Fig. 1 Study area and hypocenter distribution.a The four regions of earthquake swarm activity (regions N, S, E, and W) based on the epicentral distribution of earthquakes in the study area with flags of K, k, or A, hypocentral depths shallower than 20 km, and magnitudes of M 1.3 or higher occurring between 1 January 2018 and 31 December 2022 extracted from the JMA unified hypocenter catalog.The color scale shows hypocentral depth, and symbol size is proportional to M .The JMA CMT solution for the largest earthquake ( M 5.4, M w 5.1), which occurred in region E on 19 June 2022, is also shown.The outsets show the location of the study area (red) on wider area maps.b Hypocenters projected onto a north-south cross section.Region S is divided into shallow (Ss) and deep regions (Sd).Symbol colors correspond to regions N (blue), E (orange), W (green), Ss (black), and Sd (red).c Cumulative frequency-magnitude distribution.The vertical dashed line indicates the threshold magnitude M th = 1.3 .d Time series of hypocentral depths.e Time series of the cumulative number of earthquakes normalized by the total number of earthquakes in each region (3528, 1521, 630, 203, and 320 in regions N, E, W, Ss, and Sd, respectively).Symbol colors in c-e are the same as in b

Fig. 2
Fig. 2 Example of temporal change in tidal stress.Example of changes in tidal stress in relation to time (black curve) and the estimated phase angle ψ and tidal level L at the earthquake occurrence time (cross)

Fig. 3
Fig. 3 Example of Houston plot.a Theoretical tidal response (gray circles, sampling interval 15 min).Black circles indicate earthquake occurrences.b Example of a Houston plot in which the tidal stress is divided into two bins: positive or negative.The gray bars show the distribution of the expected background relative frequency N exp of tidal stress (relative frequency distribution of tidal stress during 1 day before and 1 day after the earthquake, sampling interval 15 min), and the bars outlined by thick black lines show the distribution of the relative frequency N obs of tidal stress at the time of the earthquake (bottom axis).The error bars show the 1σ error (assuming a total of 100 events), where σ = √ np i (1 − p i ) for the binomial distribution B(n, p i ) , n is the total number of earthquakes, and p i is the probability of random events in the i th bin obtained from the frequency distribution of theoretical tidal values.The blue diamonds, which show the ratio of N obs to N exp (top axis), are connected by a straight blue line.c Example of a Houston plot in which the tidal stress is divided into six bins.The orange line indicates the value obtained by maximum likelihood method (Yabe et al. 2015) of N obs (�S)/N exp (�S) = e α�S (Eq.(4)) to the N obs /N exp values in the six bins.In this example, tidal sensitivity α = 0.5

Fig. 4
Fig. 4 Tidal phase angle histograms for each tidal index in each of the five regions.For each region, the number of earthquakes N is shown at the top right (above panel c).The dashed line in each panel shows the average relative frequency, and the thick curve is the sine wave approximation with amplitude A 1 and peak phase angle ϕ p .Cases with significant correlations ( p < 0.05 ) are framed in red

Fig. 5
Fig. 5 Relationships of tidal phase angles between pairs of the four significant tidal indices in region Sd.Δτ and ΔCFF (0.1, 0.2) each show a phase difference of ~ 180° with respect to ΔV

Fig. 7
Fig. 7 Relationship between tidal phase angle and the seven dominant tidal level indices in region Sd. a ΔV, b Δτ, c ΔCFF (0.1) , d ΔCFF (0.2) , e ΔCFF (0.3) , f ΔCFF (0.4) , and g ΔCFF (0.5) .When the phase angle was in the range of − 90° to + 90°, the tidal level was generally positive, and generally negative otherwise (white areas in each panel).However, there were some exceptional events that fall into the gray area

Fig. 8
Fig.8Houston plots with six bins for each tidal index in region Sd.See Fig.3cfor how to read the figure.The horizontal axis shows the minimum-to-maximum range of each index.In no cases could the null hypothesis be rejected at the 5% significance level based on the Chi-square test.Bins with N obs = 0 have error bars of zero length.�α represents the uncertainty in the estimated parameter α and corresponds to the 95% confidence interval of a Gaussian distribution(Yabe et al. 2015)

Fig. 9
Fig. 9 Depth change of b-value.a Change of the b-value with depth in the Noto region.The vertical bars show the depth range used to estimate each b-value, and the horizontal bars show the 1σ error.The symbol colors correspond to b cumulative frequency-magnitude distributions at selected depths

Table 2 b
-value of the G-R law N is the total number of events