- Full paper
- Open access
- Published:
Characteristics of secondary-ruptured faults in the Aso Caldera triggered by the 2016 Mw 7.0 Kumamoto earthquake
Earth, Planets and Space volume 72, Article number: 175 (2020)
Abstract
The 16 April 2016 Mw 7.0 Kumamoto earthquake caused prominent fault displacements and crustal deformation, not only around the main rupture faults but also around numerous secondary-ruptured faults. The physics and characteristics of such secondary faulting have not yet been studied in detail. We investigated a set of two secondary faults that appeared at the timing of the Mw 7.0 quake in the Aso Caldera by mainly using synthetic aperture radar interferometry and fault slip modeling. The two faults were found to be associated with surface displacement offsets of several centimeters or more, in the oblique sense of right-lateral and vertical motion. Fault slip inversions found that the slip was dominantly in normal sense with smaller contribution from the right-lateral component. The deeper limit of the slips was estimated to be around 1.3 km, which may coincide with the boundary between the superficial sediment layer and the basement rock. The shallowness of the slip and the difference in the dip angles of the main secondary fault and the Mw 7.0 seismogenic fault suggest separation of the two fault systems, although the fault strike and sense of motions were similar. The amount of slip on the two secondary faults was larger than that expected from the scaling law derived from seismogenic faults, which may indicate the difference in the physics of seismogenic and secondary faultings.
Introduction
The 16 April 2016 Mw 7.0 Kumamoto earthquake, the largest event of the 2016 Kumamoto earthquake sequence, caused prominent fault displacements and crustal deformation. Complex rupture patterns associated with the earthquake were revealed by InSAR (e.g., Fujiwara et al. 2016; Ozawa et al. 2016), SAR pixel offset analysis (e.g., Himematsu and Furuya 2016; Yue et al. 2017), differential Lidar (Scott et al. 2018), and field surveys (Goto et al. 2017; Kumahara and Research Group of Inter-University 2016; Shirahama et al. 2016). Most studies on fault slip estimation using geodetic and seismic data indicated that the maximum fault slip of the earthquake was around 5 m or larger (e.g., Himematsu and Furuya 2016; Kobayashi et al. 2017; Kubo et al. 2016; Tanaka et al. 2019; Yue et al. 2017). Multiple faults having different dip and rake angles were required to explain the geodetic and seismic data. Such complex faulting has been revealed by modern instruments and analysis techniques also for other crustal earthquakes of similar size (e.g., Elliott et al. 2012; Hamling et al. 2017; Oskin et al. 2012; Wei et al. 2011).
What makes the Kumamoto earthquake unique was numerous triggered secondary faultings captured by InSAR (Fujiwara et al. 2016; Ozawa et al. 2016). We define secondary-ruptured faults, or secondary faults, as those (1) that were triggered to slip in response to the main seismogenic faulting; (2) that are significantly smaller in their slip dimensions than those of the main faulting, and (3) whose continuity from the main seismogenic fault is not discernible. Fujiwara et al. (2016) found approximately 230 secondary surface faultings, some along previously mapped fault traces, and many others along unmapped locations. Since the deformation around each of the secondary faults did not extend over a broad area, Ozawa et al. (2016) speculated that the fault slips must have been limited to shallow depths.
In this study, we focus on a set of secondary-ruptured faults in the Miyaji district located in the northern floor of Aso caldera. Ishimura et al. (2017) reported that vertical and dextral displacements of less than 10 cm were found in the field along the faults in Miyaji. They also excavated a pit across a surface break and found a nearly vertical fault structure. As stated earlier, the secondary faults do not seem to be connected to the main source fault of the Kumamoto earthquake. If so, how do these secondary faults extend to deeper depths? What is the relation between the rupture depth and the structural properties of the rocks forming the shallow part of the crust? The purpose of this paper is to estimate the fault geometry and slip distribution of the secondary faults in the Aso caldera by inverting InSAR data observed from three different directions, in order to obtain insights on the nature of secondary faults.
This paper pairs with Ishimura et al. (submitted to Earth, Planets and Space), who present and discuss the results of a trench survey conducted on the same secondary faults, along with some geomorphological discussion taking into account the detailed deformation field captured by InSAR. The set of two studies enables us to deduce insights not only on the ruptures that occurred in 2016, but also on past events including the long-term slip rates.
InSAR analysis and features of the displacement field
We investigate two faults located in the Miyaji district, Aso City, Kumamoto Prefecture, within the Aso caldera (Fig. 1). The targeted faults run in the direction of N50\(^\circ\) E, roughly parallel to the main rupture fault of the Mw 7.0 Kumamoto earthquake. Hereafter, we call these faults as Miyaji faults. No aftershocks were located around the Miyaji faults (Fig. 1), which is a piece of supporting evidence that the faults are not connected to the main rupture fault of the Kumamoto earthquake.
In order to capture displacements around the faults from different directions, three sets of ALOS-2 SAR images obtained along three different line-of-sight (LOS) directions were used to form interferograms (Table 1). According to the local residents, the faulting in Miyaji occurred simultaneously with the Mw 7.0 rupture. In principle, measurements from three independent directions enable us to obtain the three-dimensional displacement field, but the north-south component has relatively large uncertainty (Wright et al. 2004; Morishita et al. 2016). This limitation is inherent to InSAR, originating from the fact that the LOS horizontal directions are close to east or west (Table 1). Standard deviation of 1 cm in each interferograms leads to standard deviation of 1.3 cm, 6.5 cm, and 0.8 cm for east–west, north–south, and vertical components, respectively, under the configuration of the data used in this study.
We conducted a 2-pass interferometry with RINC software (Ozawa et al. 2016). We multi-looked the images by 4 and 7 in the range and azimuth directions, respectively, which resulted in the spatial resolution of approximately 15–20 m in the two directions. Digital ellipsoidal height model developed by M. Tobita (Geospatial Information Authority of Japan, GSI) and T. Ozawa (National Research Institute for Earth Science and Disaster Resilience of Japan) based on the digital elevation model of GSI was used for geocoding and computation of topographic phase. We used an adaptive spectral filter (Baran et al. 2003), followed by median filter (\(5 \times 5\)) and mean filter (\(3 \times 3\)) to reduce noise. The resulting interferograms show high coherence in the target area (Fig. 2). While similar fringe patterns were obtained for Paths 126 and 23 (obliquely from ENE and ESE, respectively), the interferogram of Path 130 (obliquely from WSW) exhibits considerably different fringe pattern.
For phase unwrapping of the interferograms (conversion from phase to displacements), SNAPHU (Chen and Zebker 2000) was used. We subsequently converted the relative displacements mapped by the interferograms to absolute (referenced) displacements by using the GNSS displacements as a reference. For that purpose, we used the coseismic GNSS displacements in the SAR line-of-sight (LOS) direction, observed at two stations Aso (960703) and Takamori (960704) of the GNSS Earth Observation Network System (GEONET) operated by GSI, which were the nearest from our target area (Fig. 2a–c). The offset for each interferogram was estimated in such a way that the misfit between the corrected InSAR displacements and GNSS displacements was minimized at the GNSS sites. Namely, we estimated the offset l for each interferogram by
where \(u_i^{\text{InSAR}}\) and \(u_i^{\text{GNSS}}\) denote the InSAR and GNSS displacements in the LOS direction at the ith GNSS station location, respectively. For each unwrapped interferograms, the estimated offset was subtracted to produce referenced LOS displacements (Fig. 3).
We decomposed the three referenced LOS displacements along different directions into eastward, northward and upward components by (Wright et al. 2004)
where \(\mathbf {u}\) is the column vector containing the vector components of displacement \((u_{\text{east}},u_{\text{north}},u_{\text{up}})^\text{T}\), \(\mathbf {P}\) is formed by a stack LOS vectors (unit row vectors pointing from the ground to the satellite in a local east, north, and up directions) for the Path 126, Path 130 and Path 23 acquisitions,
and the vector \(\mathbf {r} = (r_{126},r_{130},r_{23})^\text{T}\) is composed of the referenced InSAR displacements (Fig. 3). Although we used different LOS vectors for each pixel taking into account the small differences in the pointing directions to the satellite, the LOS vectors at the center of the Miyaji fault location are shown in Table 2. To better extract the displacement features related to the faultings, the obtained horizontal displacements were further converted to a horizontal component along N50\(^\circ\) E, approximately parallel to the strike of the Miyaji faults, by
The displacement field along the N50\(^\circ\) E horizontal direction (Fig. 4a) shows localized displacements around the Miyaji faults, superimposed by broad regional right-lateral movement that mainly resulted from the rupture of the main fault of the Mw 7.0 Kumamoto earthquake. The vertical displacements (Fig. 4b), on the other hand, exhibits north-side-down vertical displacements across two strands of faults. More closely looking, each of two strands of faults is composed of multiple smaller faults located en echelon (see Ishimura et al. (submitted to Earth, Planets and Space) for detailed interpretation). The N50\(^\circ\) E horizontal component is noisier than the vertical one because of the relatively large uncertainty in the NS displacements.
The profiles of displacements across the faults show clear displacement offsets at the location of the faults (Fig. 5). The offsets approximately amount to 5 cm and 8 cm, for the northern and southern faults, respectively, in both vertical and right-lateral profiles. On the northern fault, the maximum vertical and dextral surface slip measured from the InSAR displacements were 8 cm and 19 cm, respectively, whereas those for the southern fault were 12 cm and 19 cm, respectively (Ishimura et al., submitted to Earth, Planets and Space).
Fault slip inversion
Method
Our main interests on the secondary faulting can be summarized by two questions: (1) what are the depths of the slipped portion of the faults? (2) How are the two secondary faults related to each other? To answer these questions, we estimated the fault slip models of the two secondary faults using the inversion procedure developed and conducted in previous studies (Fukushima et al. 2013, 2018; Ghayournajarkar and Fukushima 2020). For more detailed explanation on the inversion procedure, see Fukushima et al. (2013).
For preparation of the inversion dataset, we first removed the regional long-wavelength displacements caused mainly by the fault rupture of the Kumamoto earthquake from each of the unwrapped interferograms (Fig. 6a). The regional displacements were modeled by a two-dimensional cubic polynomial function expressed by:
where \(a_1, \ldots a_{10}\) are constants estimated by the least-squares method, and x and y are longitude and latitude. This estimation was performed in the areas shown in Fig. 6. The resulting two-dimensional polynomial function (Fig. 6b) was subtracted from the original (Fig. 6a) to extract the displacements due to the rupture of Miyaji faults (Fig. 6c). An alternative approach for removing the long-wavelength displacements would be to model the coseismic displacements of the main rupture of the Kumamoto earthquake and subtract them, but this approach requires caution because complete removal of the long-wavelength displacements is not guaranteed. Our approach has an advantage that all the long-wavelength variations can be eliminated including long-wavelength noise due to orbital inaccuracies and atmospheric delay.
The detrended interferograms show localized displacements around the Miyaji faults and little displacements far from them except for localized depressions seemingly unrelated to the faultings, indicating the adequacy of the detrending processing. After removal of the long-wavelength displacements, the data were subsampled with concentric grids (Fig. 6d).
We used triangular dislocation elements to express curved fault surfaces, and faults were discretized with an average interval of 200 m. We solved for the geometry and slip distribution on the northern and southern secondary faults simultaneously. We assumed that, for each of the northern or southern fault, the small en echelon structures mentioned above were connected to a single fault at depth, and approximated the top of the fault to be a smooth curve close to a straight line. The bottom side of the fault was assumed to be a straight line parallel to the top side. The lengths of the faults were taken large enough to cover the extent where deformation was observed in the interferograms (black curves in Fig. 6). The displacement fields (e.g., Fig. 4) seem to indicate smaller amount of faulting extending beyond the assumed fault location, to the northeast from the southern Miyaji fault en echelon, and perhaps to the southwest of it as well, but we did not attempt to model these displacements because they are small and the displacement gradients across the fault that are important for constraining the depth of the faulting were hardly discernible. The widths of the faults were taken large enough so that the slipped area does not reach the bottom of the assumed faults.
We solved for the amount of slip on the triangular elements, dip and rake angles of the two faults, and the weight of the smoothing of the slip distribution. Slips on the fault elements have linear relationship to the observed data, whereas the other parameters have nonlinear relationship to the data. This is hence a mixed linear–nonlinear inverse problem.
The likelihood function p was defined as follows (Fukuda and Johnson 2008):
where C is an arbitrary constant, \(\mathbf {d}\) is the data vector (length N), \(\mathbf {s}\) is a model vector composed of linear parameters (slip on the triangular elements, length M), \(\sigma ^2\) is the variance of the data which was set to \(1.0 \times 10^{-4}\) \(\text {m}^2\) based on the calculation in the far-field, \(\alpha ^2\) is the weight of the smoothing matrix \(\mathbf {L(m)}\), \(\mathbf {m}\) is a model vector composed of the dip and rake angles of the two faults, \(\mathbf {G(m)}\) is the design matrix relating the slip on the fault elements to displacements on the ground surface, and \(\Sigma _{d}^{-1}\) is the inverse of the data covariance matrix.
Since \(\mathbf {L(m)}\) and \(\mathbf {G(m)}\) are functions of the geometrical model parameters, these matrices are fixed once the values of \(\mathbf {m}\) are fixed. For the design matrix \(\mathbf {G(m)}\), we used solutions of angular dislocations in a homogeneous elastic half-space (Thomas 1993) with an assumption of Poisson’s ratio of 0.25. We adopted the smoothing matrix defined by Maerten et al. (2005), where the operator concerning the ith element and the three elements in contact of it is expressed by:
Here, \(h_{ij}\) is the distance between the center of the ith and jth elements, \(s_i\) and \(s_j\) are the ith and jth elements in \(\mathbf {s}\), and \(l_i=\sum _{j=1}^{M}h_{ij}\), where M is the number of elements, is the sum of the distance between element centers.
We assumed that the inverse of the data covariance matrix \(\Sigma _{d}^{-1}\) was an identity matrix, meaning that we assumed no spatial correlation of noise. Trials with data correlation taken into account led to models that had slip deeper than a few kilometers and large residuals in the vicinity of the faults. This problem occurred probably because taking the data correlation into account is equivalent to placing less weight into short-scale variations in data; that is, the inversion algorithm attempted to search for models that fit with relatively small amplitude of large wavelength noise variations rather than the short-wavelength signals around the faults.
The nonlinear parameters were solved using particle swarm optimisation (PSO) (Eberhart and Kennedy 1995). PSO is a Monte Carlo search algorithm and it evaluates the likelihood function for many sets of nonlinear parameter values. For each set, the slip distribution (values of slip on the fault elements) on each fault was solved with a damped least-squares method with non-negativity constraints. The misfit of the least-squares inversion was used to evaluate the likelihood function.
As for the model uncertainties, neighborhood approximation and Monte Carlo integration of the likelihood function (Eq. 6) was performed to form the probability density function (Sambridge 1999), and the standard deviation of the slip distribution was calculated from the diagonal elements of the model covariance matrix. The 95% confidence intervals of the nonlinear parameters were calculated by discarding the top and bottom 2.5% of the one-dimensional marginal probability density functions (Fukushima et al. 2013).
Results
The convergence of the nonlinear parameters (Fig. 7 top) and the one-dimensional marginal probability distribution of the parameters (Fig. 7 bottom) show that wide ranges of model parameters were searched and the search converged toward one single set of model parameters, leaving a clear peak in each marginal probability distribution. The dip angle of the northern fault was less constrained compared with the southern one (Table 3), but the difference in the confidence intervals indicate that the northern fault is dipping steeper than the southern one. The dip angles of the northern and southern faults were 63\(^\circ\) and 42\(^\circ\), respectively, for the maximum likelihood model. As for the rake angles, similar values of around 240\(^\circ\) were estimated for the two faults, indicating dominantly normal faulting with smaller amount of right-lateral slip.
The geometry and slip distributions of the faults of the maximum likelihood model and the uncertainty of slip are shown in Figs. 8 and 9. Larger dip angle for the northern fault means that the two faults are converged at depth, forming splay fault branching. It is worth pointing out that the convergence depth coincides with the deeper limit of the slipped area of the southern fault. One should note, however, that our results do not warrant the connection of the two faults at depth because the shallow slip on the northern fault does not extend down to the depth where two fault surfaces are converged. This set of fault slip models well explains the data (Fig. 10).
The estimated slip extends from the surface to \(\sim\)0.5 km with a peak slip of 11 cm, and to \(\sim\)1.3 km with a peak slip of 19 cm, for the northern and southern faults, respectively. For the northern fault, the maximum standard deviation was 3.7 cm, whereas the standard deviation for most of the southern fault elements shallower than 1.5 km was less than 5 cm (Fig. 9). The estimated slip at most area along the faults can be considered real, given the small amount of standard deviation of slip.
The moment magnitudes for the northern and southern faults amount to 3.52 and 4.36, respectively. Here, we assumed the rigidity of 6.2 GPa calculated from the density of 1900 kg/\(\text {m}^3\) (Komazawa 1995) and the shear velocity of 1.8 km/s (Abe et al. 2010). As far as the authors are aware, the smallest moment magnitude of surface-rupturing earthquakes found so far in the world is Mw 4.7 of the 2007 Katanning earthquake in Australia (Dawson et al. 2008; King et al. 2019), highlighting the peculiarity of the secondary ruptures.
We calculated the Coulomb stress changes caused by the Kumamoto mainshock on the Miyaji faults using the Coulomb software (version 3.3.01) (Lin and Stein 2004; Toda et al. 2005) and the mainshock slip model of Ozawa et al. (2016). In this calculation, we assumed a homogeneous half-space with Poisson’s ratio 0.25 and rigidity of 6.2 GPa, effective fault friction coefficient of 0.4, and the rake directions same as the estimated rake angles (227\(^\circ\) and 233\(^\circ\) for the northern and southern Miyaji faults, respectively). The Coulomb stress changes on the two faults were positive and amounted up to 32 kPa and 85 kPa for the northern and southern faults, respectively. This result indicates that the secondary faultings could have been statically triggered by the mainshock.
Discussion
Relation to the subsurface structure
The superficial structure beneath the northern part of the Aso caldera floor, also called the Aso Valley, consists of sediments deposited in the caldera lake underlain by the pumice flow from the Kuju volcano (Matsumoto and Fujimoto 1969). Further down, basement rock of granite has been identified deeper than \(\sim\)500 m near the northern rim of the caldera (Matsumoto and Fujimoto 1969).
Beneath the location of the Miyaji faults, the data of gravity (Komazawa 1995; Miyakawa et al. 2016), seismic velocity structure (Sudo and Kong 2001; Abe et al. 2010), and electrical resistivity (Handa et al. 1998; Hase et al. 2005; Hata et al. 2016, 2018; Matsushima et al. 2020) altogether indicate a low-density (porous having high water content) superficial layer of 1 to 2 km, which can be attributed to the lake deposits and pumice flow.
The consistency between the thickness of the superficial sediment layer and the depth extent of the estimated slip (\(\sim\)1.3 km) indicates that the fault ruptured the sediment layer. Although our results cannot exclude the possibility that the fault extends further down, it is more likely that Miyaji faults are “root-less” and cause minor ruptures within the superficial layer, considering the fact that the slip of the penultimate event was almost the same as the 2016 one (Ishimura et al., submitted to Earth, Planets and Space).
The optimal dip angle of the northern and southern faults were 63\(^{\circ }\) and 42\(^{\circ }\), respectively. Since the dip angles found in the trench were nearly vertical (Ishimura et al., submitted to Earth, Planets and Space), it is implied that the Miyaji faults are listric. The dip angle of the main southern fault is smaller than the Futagawa fault responsible for the Mw 7.0 Kumamoto main rupture, estimated in a range of 65\(^{\circ }\)–79\(^{\circ }\) (Kubo et al. 2016; Himematsu and Furuya 2016; Ozawa et al. 2016; Yue et al. 2017). This supports the lack of continuation of Miyaji faults to the main Futagawa fault structure.
Considering the fact that numerous secondary faultings of the Kumamoto earthquake appeared not only in the Aso caldera, but also outside of it, we speculate that the volcanic products (lava and pyroclastic material), commonly constituting the superficial layer in the region, facilitate secondary faultings. Around the Unzen volcano and the Beppu–Haneyama fault zone (10–20 km from the Kuju volcano), both located in the Beppu–Shimabara graben in Kyushu, there are areas with densely clustered active faults of a few to several kilometers (Fig. 11a). Similarly, the Tanna fault connected to the Hakone volcano in central Japan is surrounded by a cluster of minor faults of similar lengths (Fig. 11b). The latter resembles the tectonic setting of the Kumamoto area, having a major fault (Futagawa fault) surrounded by minor faults, and hence it is possible that the minor faults in the Tanna fault area are activated by secondary faulting.
Effect of Earth’s layered structure
The deformation field in the crust due to faulting is dependent on the properties of the surrounding medium. It has been known that ignoring the effect of layering biases the estimated fault depth to be shallower than the actual depth, when the fault is located beneath a compliant superficial layer (Cattin et al. 1999; Savage 1998; Segall 2010). In our case, the faults are inferred to be located in a compliant superficial layer overlain by stiffer bedrock. We performed a synthetic test to evaluate the degree of bias in the depth and dip angle estimated under the simplified assumption of homogeneous half-space.
As the actual fault slip, we assumed a uniform slip of 1 m on a rectangular area reaching the ground surface with length of 2 km, dip angle of 45\(^{\circ }\), and rake angle of 240\(^{\circ }\). We tested with different widths, hence with different bottom depths of the faulting.
We assumed that the Earth is formed by a compliant layer of 2 km thickness having rigidity of 6.2 GPa and a half-space beneath it having rigidity of 33 GPa; the latter was given based on the shear velocity of 3.5 km/s (Abe et al. 2010) and the density of 2700 kg/\(\text {m}^3\) (Komazawa 1995) (Fig. 12a).
The observation points were assumed to be distributed in a concentric grid, similar to the data we used for the inversion. The eastern, northern, and vertical displacements at the observation points were calculated by the EDGRN/EDCMP package (Wang et al. 2003) (Fig. 12b–d).
The three-dimensional displacements were inverted to estimate the fault width and dip angle using PSO for each case of actual fault width. The results revealed that the effect of layering was minor for our studied case (Fig. 13). When the actual width was smaller than or equal to 1.0 km, the estimated width was nearly equal to or slightly smaller than the actual width. When the actual width was larger than 1.0 km, the estimated width was slightly larger than the actual width. The maximum degree of overestimation was 7.5% (for the case where the actual width was 1.5 km). As for the dip angle, non-systematic deviation of less than 2.5° from the actual dip angle was observed when the actual width was smaller than 0.5 km, which probably resulted from a relatively small number of points on which the deformation was observed. Systematic underestimation of approximately 1° or less was observed when the actual width was larger than or equal to 0.5 km.
From the synthetic test, we conclude, for the purpose of this research, that the bias was limited and that the assumption of homogeneous half-space and hence the discussion on the rupture depth made in the previous section were valid.
Fault scaling
Wells and Coppersmith (1994) derived empirical relationships among magnitude, rupture length and width, rupture area, and average displacement, from 244 historical events in the world. Although the relationships are basically for seismogenic faults, here we compared them with our findings on the Miyaji faults to shed light on the differences in the seismogenic and secondary faultings. Wells and Coppersmith (1994) gave, for strike–slip and normal faults,
and
where \(\text {MD}_x\) and \(\text {SRL}_x\) (x = s for strike–slip and x = n for normal slip) correspond to the maximum surface slip (m) and surface rupture length (km), respectively. We assumed here that the relationships were applicable to the Miyaji faults even though their parameter values lie around the smaller limit of the data range used to derive the empirical relationships of Eqs. (8) and (9) (minimum of \(\text {MD}_x\) and \(\text {SRL}_x\) for both strike–slip and normal faults used by Wells and Coppersmith (1994) were 2 cm and 3.8 km, respectively).
As the surface fault lengths, we used 1.5 km and 3.2 km for the northern and southern faults, each of which correspond to the extent where slip was estimated on the surface elements by the inversion (Fig. 8). Substituting these values to the equations, the expected amount of strike–slip would be 3.2 cm and 7.9 cm, respectively, and that of normal slip would be 1.9 cm and 6.1 cm, respectively. The amount of maximum dextral surface slip measured from the InSAR displacements on the northern and southern faults was, on the other hand, 8 cm and 19 cm, respectively, and that of normal slip was 12 cm and 19 cm, respectively (Ishimura et al., submitted to Earth, Planets and Space). The large deviations of the observed values from the relationship of Eqs. (8) and (9) may indicate that secondary faults tend to slip more than the average of seismogenic faulting. Whether this tendency is a ubiquitous characteristic is an intriguing question to be investigated in future studies.
Conclusions
The Mw 7.0 Kumamoto earthquake triggered numerous secondary faultings. We obtained detailed fault-parallel and vertical displacement fields around a set of two secondary-ruptured faults in Miyaji located in the Aso caldera by conducting InSAR analysis using the ALOS-2 SAR data. The results revealed clear dextral and vertical surface displacement offsets across the two faults.
By conducting fault slip inversions, we further found that the northern and southern faults ruptured down to \(\sim\)0.5 km and \(\sim\)1.3 km, with the optimal dip angles of 65\(^{\circ }\) and 42\(^{\circ }\), respectively. The two faults seemed to form a splay structure whose connection depth lay at \(\sim\)1.3 km. The deeper limit of the rupture may correspond to the base of the superficial sediment layer. The shallowness of the ruptures as well as the difference in the dip angles of the main southern fault and the steeper dip angle of the Futagawa fault, the main seismogenic fault responsible for the Mw 7.0 quake, suggest that the secondary faults are separated from the seismogenic fault system.
Our analysis suggested that the slip on the Miyaji secondary faults was larger than expected from the scaling law proposed for seismogenic faults (Wells and Coppersmith 1994). Systematic studies of secondary faulting, being detected for the Kumamoto earthquake case (Fujiwara et al. 2016; Ozawa et al. 2016) and other earthquakes by InSAR and differential Lidar techniques, will be able to elucidate the characteristics of secondary faults and provide insights into the physics of faulting.
Data availability statement
The datasets used and/or analyzed during the current study are available from the authors upon request.
Abbreviations
- ALOS-2:
-
Advanced Land Observing Satellite 2
- InSAR:
-
Interferometric Synthetic Aperture Radar
- LOS:
-
Line of sight
- GEONET:
-
GNSS Earth Observation Network System
- GSI:
-
Geospatial Information Authority of Japan
- PSO:
-
Particle Swarm Optimisation
- SAR:
-
Synthetic Aperture Radar
References
Abe Y, Ohkura T, Shibutani T, Hirahara K, Kato M (2010) Crustal structure beneath Aso Caldera, Southwest Japan, as derived from receiver function analysis. J Volcanol Geotherm Res 195(1):1–12. https://doi.org/10.1016/j.jvolgeores.2010.05.011
Baran I, Stewart MP, Kampes BM, Perski Z, Lilly P (2003) A modification to the Goldstein radar interferogram filter. IEEE Trans Geosci Remote Sens 41(9):2114–2118. https://doi.org/10.1109/Tgrs.2003.817212
Cattin R, Briole P, Lyon-Caen H, Bernard P, Pinettes P (1999) Effects of superficial layers on coseismic displacements for a dip-slip fault and geophysical implications. Geophys J Int 137(1):149–158. https://doi.org/10.1046/j.1365-246x.1999.00779.x
Chen CW, Zebker HA (2000) Network approaches to two-dimensional phase unwrapping: intractability and two new algorithms. J Opt Soc Am A 17(3):401–414. https://doi.org/10.1364/Josaa.17.000401
Dawson J, Cummins P, Tregoning P, Leonard M (2008) Shallow intraplate earthquakes in Western Australia observed by interferometric synthetic aperture radar. J Geophys Res Solid Earth. https://doi.org/10.1029/2008jb005807
Eberhart RC, Kennedy J (1995) A new optimizer using particle swarm theory. In: Proceedings of the sixth international symposium on micromachine and human science, Nagoya, Japan, pp 39–43
Elliott JR, Nissen EK, England PC, Jackson JA, Lamb S, Li Z, Oehlers M, Parsons B (2012) Slip in the 2010–2011 Canterbury earthquakes, New Zealand. J Geophys Res Solid Earth. https://doi.org/10.1029/2011jb008868
Fujiwara S, Yarai H, Kobayashi T, Morishita Y, Nakano T, Miyahara B, Nakai H, Miura Y, Ueshiba H, Kakiage Y, Une H (2016) Small-displacement linear surface ruptures of the 2016 Kumamoto earthquake sequence detected by ALOS-2 SAR interferometry. Earth Planets Space. https://doi.org/10.1186/s40623-016-0534-x
Fukuda J, Johnson KM (2008) A fully Bayesian inversion for spatial distribution of fault slip with objective smoothing. Bull Seismol Soc Am 98(3):1128–1146. https://doi.org/10.1785/0120070194
Fukushima Y, Takada Y, Hashimoto M (2013) Complex ruptures of the 11 April 2011 Mw 6.6 Iwaki earthquake triggered by the 11 March 2011 Mw 9.0 Tohoku earthquake, Japan. Bull Seismol Soc Am 103(2B):1572–1583. https://doi.org/10.1785/0120120140
Fukushima Y, Toda S, Miura S, Ishimura D, Fukuda J, Demachi T, Tachibana K (2018) Extremely early recurrence of intraplate fault rupture following the Tohoku-Oki earthquake. Nat Geosci 11(10):777–781. https://doi.org/10.1038/s41561-018-0201-x
Ghayournajarkar N, Fukushima Y (2020) Determination of the dipping direction of a blind reverse fault from InSAR: case study on the 2017 Sefid Sang earthquake, northeastern Iran. Earth Planets Space 72(1):64. https://doi.org/10.1186/s40623-020-01190-6
Goto H, Tsutsumi H, Toda S, Kumahara Y (2017) Geomorphic features of surface ruptures associated with the 2016 Kumamoto earthquake in and around the downtown of Kumamoto City, and implications on triggered slip along active faults. Earth Planets Space. https://doi.org/10.1186/s40623-017-0603-9
Hamling IJ, Hreinsdottir S, Clark K, Elliott J, Liang CR, Fielding E, Litchfield N, Villamor P, Wallace L, Wright TJ, D’Anastasio E, Bannister S, Burbidge D, Denys P, Gentle P, Howarth J, Mueller C, Palmer N, Pearson C, Power W, Barnes P, Barrell DJA, Van Dissen R, Langridge R, Little T, Nicol A, Pettinga J, Rowland J, Stirling M (2017) Complex multifault rupture during the 2016 M-w 7.8 Kaikoura earthquake, New Zealand. Science. https://doi.org/10.1126/science.aam7194
Handa S, Suzuki A, Tanaka Y (1998) The electrical resistivity structure of Aso caldera, Japan. Bull Volcanol Soc Jpn 43:15–23
Hase H, Hashimoto T, Sakanaka S, Kanda W, Tanaka Y (2005) Hydrothermal system beneath Aso volcano as inferred from self-potential mapping and resistivity structure. J Volcanol Geotherm Res 143(4):259–277. https://doi.org/10.1016/j.jvolgeores.2004.12.005
Hata M, Takakura S, Matsushima N, Hashimoto T, Utsugi M (2016) Crustal magma pathway beneath Aso caldera inferred from three-dimensional electrical resistivity structure. Geophys Res Lett 43(20):10,720–10,727. https://doi.org/10.1002/2016gl070315
Hata M, Matsushima N, Takakura S, Utsugi M, Hashimoto T, Uyeshima M (2018) Three-dimensional electrical resistivity modeling to elucidate the crustal magma supply system beneath Aso caldera, Japan. J Geophys Res Solid Earth 123(8):6334–6346. https://doi.org/10.1029/2018jb015951
Himematsu Y, Furuya M (2016) Fault source model for the 2016 Kumamoto earthquake sequence based on ALOS-2/PALSAR-2 pixel-offset data: evidence for dynamic slip partitioning. Earth Planets Space. https://doi.org/10.1186/s40623-016-0545-7
Ishimura D, Toda S, Ichihara T, Takahashi N, Konno A, Sato H (2017) A study on surface ruptures around Miyaji, Aso City, Kumamoto Prefecture, associated with the 2016 Kumamoto earthquake sequence and upward slip tapering on pit excavation walls. Act Fault Res 47:9–16 (in Japanese with English abstract)
King TR, Quigley M, Clark D (2019) Surface-rupturing historical earthquakes in Australia and their environmental effects: new insights from re-analyses of observational data. Geosciences. https://doi.org/10.3390/geosciences9100408
Kobayashi H, Koketsu K, Miyake H (2017) Rupture processes of the 2016 Kumamoto earthquake sequence: causes for extreme ground motions. Geophys Res Lett 44(12):6002–6010. https://doi.org/10.1002/2017gl073857
Komazawa M (1995) Gravimetric analysis of Aso Volcano and its interpretation. J Geod Soc Jpn 41(1):17–45. https://doi.org/10.11366/sokuchi1954.41.17
Kubo H, Suzuki W, Aoi S, Sekiguchi H (2016) Source rupture processes of the 2016 Kumamoto, Japan, earthquakes estimated from strong-motion waveforms. Earth Planets Space. https://doi.org/10.1186/s40623-016-0536-8
Kumahara Y, Research Group of Inter-University (2016) The characteristics of surface rupture associated by the 2016 Kumamoto earthquake, programme and abstracts JSAF 2016 fall meeting, Hosei University, Tokyo, 29–30 October 2016 (in Japanese)
Lin J, Stein RS (2004) Stress triggering in thrust and subduction earthquakes and stress interaction between the southern San Andreas and nearby thrust and strike-slip faults. J Geophys Res Solid Earth. https://doi.org/10.1029/2003jb002607
Maerten F, Resor P, Pollard D, Maerten L (2005) Inverting for slip on three-dimensional fault surfaces using angular dislocations. Bull Seismol Soc Am 95(5):1654–1665. https://doi.org/10.1785/0120030181
Matsumoto H, Fujimoto Y (1969) Noticeable result of a drilling on the northeastern floor in Aso caldera. Bull Volcanol Soc Jpn 14:1–7 (in Japanese with English abstract)
Matsushima N, Utsugi M, Takakura S, Yamasaki T, Hata M, Hashimoto T, Uyeshima M (2020) Magmatic-hydrothermal system of Aso volcano, Japan, inferred from electrical resistivity structures. Earth Planets Space 72(1):57. https://doi.org/10.1186/s40623-020-01180-8
Miyakawa A, Sumita T, Okubo Y, Okuwaki R, Otsubo M, Uesawa S, Yagi Y (2016) Volcanic magma reservoir imaged as a low-density body beneath Aso volcano that terminated the 2016 Kumamoto earthquake rupture. Earth Planets Space. https://doi.org/10.1186/s40623-016-0582-2
Morishita Y, Kobayashi T, Yarai H (2016) Three-dimensional deformation mapping of a dike intrusion event in Sakurajima in 2015 by exploiting the right- and left-looking ALOS-2 InSAR. Geophys Res Lett 43(9):4197–4204. https://doi.org/10.1002/2016gl068293
Nakata T, Imaizumi T (2002) Digital active fault map of Japan. University of Tokyo Press, Tokyo (in Japanese with abstract in English)
Oskin ME, Arrowsmith JR, Corona AH, Elliott AJ, Fletcher JM, Fielding EJ, Gold PO, Garcia JJG, Hudnut KW, Liu-Zeng J, Teran OJ (2012) Near-field deformation from the El Mayor-Cucapah Earthquake revealed by differential LIDAR. Science 335(6069):702–705. https://doi.org/10.1126/science.1213778
Ozawa T, Fujita E, Ueda H (2016) Crustal deformation associated with the 2016 Kumamoto earthquake and its effect on the magma system of Aso volcano. Earth Planets Space. https://doi.org/10.1186/s40623-016-0563-5
Sambridge M (1999) Geophysical inversion with a neighbourhood algorithm-II. Appraising the ensemble. Geophys J Int 138(3):727–746. https://doi.org/10.1046/j.1365-246x.1999.00900.x
Savage JC (1998) Displacement field for an edge dislocation in a layered half-space. J Geophys Res Solid Earth 103(B2):2439–2446. https://doi.org/10.1029/97jb02562
Scott CP, Arrowsmith JR, Nissen E, Lajoie L, Maruyama T, Chiba T (2018) The M7 2016 Kumamoto, Japan, earthquake: 3-D deformation along the fault and within the damage zone constrained from differential lidar topography. J Geophys Res Solid Earth 123(7):6138–6155. https://doi.org/10.1029/2018jb015581
Segall P (2010) Earthquake and volcano deformation. Princeton University Press, New Jersey
Shirahama Y, Yoshimi M, Awata Y, Maruyama T, Azuma T, Miyashita Y, Mori H, Imanishi K, Takeda N, Ochi T, Otsubo M, Asahina D, Miyakawa A (2016) Characteristics of the surface ruptures associated with the 2016 Kumamoto earthquake sequence, Central Kyushu. Earth Planets Space, Japan. https://doi.org/10.1186/s40623-016-0559-1
Sudo Y, Kong LSL (2001) Three-dimensional seismic velocity structure beneath Aso Volcano, Kyushu, Japan. Bull Volcanol 63(5):326–344. https://doi.org/10.1007/s004450100145
Tanaka Y, Ohta Y, Miyazaki S (2019) Real-time coseismic slip estimation via the GNSS carrier phase to fault slip approach: a case study of the 2016 Kumamoto earthquake. Geophys Res Lett 46(3):1367–1374. https://doi.org/10.1029/2018gl080741
Thomas AL (1993) Poly3D: a three-dimensional, polygonal element, displacement discontinuity boundary element computer program with applications to fractures, faults, and cavities in the Earth’s crust. MS thesis, Stanford University
Toda S, Stein RS, Richards-Dinger K, Bozkurt SB (2005) Forecasting the evolution of seismicity in southern California: animations built on earthquake stress transfer. J Geophys Res Solid Earth. https://doi.org/10.1029/2004jb003415
Wang RJ, Martin FL, Roth F (2003) Computation of deformation induced by earthquakes in a multi-layered elastic crust– FORTRAN programs EDGRN/EDCMP. Comput Geosci 29(2):195–207. https://doi.org/10.1016/S0098-3004(02)00111-5
Wei SJ, Fielding E, Leprince S, Sladen A, Avouac JP, Helmberger D, Hauksson E, Chu RS, Simons M, Hudnut K, Herring T, Briggs R (2011) Superficial simplicity of the 2010 El Mayor-Cucapah earthquake of Baja California in Mexico. Nat Geosci 4(9):615–618. https://doi.org/10.1038/ngeo1213
Wells DL, Coppersmith KJ (1994) New empirical relationships among magnitude, rupture length, rupture width, rupture area, and surface displacement. Bull Seismol Soc Am 84(4):974–1002
Wright TJ, Parsons BE, Lu Z (2004) Toward mapping surface deformation in three dimensions using InSAR. Geophys Res Lett. https://doi.org/10.1029/2003gl018827
Yue H, Ross ZE, Liang CR, Michel S, Fattahi H, Fielding E, Moore A, Liu Z, Jia B (2017) The 2016 Kumamoto Mw = 7.0 earthquake: a significant event in a fault-volcano system. J Geophys Res Solid Earth 122(11):9166–9183. https://doi.org/10.1002/2017jb014525
Acknowledgements
We used the PALSAR and PALSAR-2 data shared among the PALSAR Interferometry Consortium to Study our Evolving Land Surface (PIXEL). The data were provided by the Japan Aerospace Exploration Agency (JAXA) under a cooperative research contract with the Earthquake Research Institute of the University of Tokyo. We are grateful to Dr. Taku Ozawa who provided the RINC interferometry processing software, to Dr. Mehdi Nikkhoo for helps with stress calculation, and to the Geospatial Authority of Japan for the GNSS data. Discussion with Dr. Shinji Toda helped to improve the quality of the paper. This study was partially supported by the Japan Society for the Promotion of Sciences KAKENHI Grant Number JP18K03795 and JP17H04730.
Funding
This study was partially supported by the Japan Society for the Promotion of Sciences KAKENHI Grant Numbers JP18K03795 and JP17H04730.
Author information
Authors and Affiliations
Contributions
YF led the preparation of the paper and conducted InSAR analysis and modeling. DI conducted field survey and discussed the geological interpretation of the modeled results. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Fukushima, Y., Ishimura, D. Characteristics of secondary-ruptured faults in the Aso Caldera triggered by the 2016 Mw 7.0 Kumamoto earthquake. Earth Planets Space 72, 175 (2020). https://doi.org/10.1186/s40623-020-01306-y
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s40623-020-01306-y