Long-period ground motion simulation using centroid moment tensor inversion solutions based on the regional three-dimensional model in the Kanto region, Japan

We conducted centroid moment tensor (CMT) inversions of moderate (Mw 4.5–6.5) earthquakes in the Kanto region, Japan, using a local three-dimensional (3D) model. We then investigated the effects of our 3D CMT solutions on long-period ground motion simulations. Grid search CMT inversions were conducted using displacement seismograms for periods of 25–100 s. By comparing our 3D CMT solutions with those from the local one-dimensional (1D) catalog, we found that our 3D CMT inversion systematically provides magnitudes smaller than those in the 1D catalog. The Mw differences between 3D and 1D catalogs tend to be significant for earthquakes within the oceanic slab. By comparing ground motion simulations between 1D and 3D velocity models, we confirmed that observed Mw differences could be explained by differences in the rigidity structures around the source regions between 3D and 1D velocity models. The 3D velocity structures (especially oceanic crust and mantle) are important for estimating seismic moments in intraslab earthquakes, which are related to fault size estimation. A detailed discussion for intraslab seismicity can be conducted using the 3D CMT catalog. The seismic moments also directly affect the amplitudes of ground motions. The 3D CMT catalog allows us to directly conduct the precise forward and inverse modeling of long-period ground motion without adjusting source models, which have been typically applied in the cases using the 1D CMT catalog. We also conducted long-period ground motion simulations using our 3D CMT solutions to evaluate the reproducibility of long-period ground motions at stations within the Kanto Basin. The simulations of our 3D CMT solutions well-reproduced observed ground motions for periods longer than 10 s, even at stations within the Kanto Basin. The reproducibility of simulations was improved from those using solutions in the 1D catalog.


Introduction
Long-period ground motions with predominant periods of 3-20 s have often been observed in large sedimentary basins and offshore regions. The duration and amplitude of long-period ground motion is enhanced due to thick low-velocity sediments along the propagation path (e.g., Boore 1999; Furumura et al. 2001Furumura et al. , 2008Olsen et al. 2006;Day et al. 2008;Gomberg 2018;Kaneko et al. 2019;Takemura et al. 2019a). Thus, long-period ground motions can cause severe resonance and damage of large-scale man-made structures in the densely populated areas. The characteristics of long-period ground motions have been summarized in Koketsu and Miyake (2008). In the Takemura et al. Earth, Planets and Space (2021) 73:15 Kanto region, Japan, long-period ground motions with predominant periods of 5-10 s have frequently been observed during shallow moderate-to-large earthquakes (e.g., Kinoshita et al. 1992;Miyake and Koketsu 2005;Yoshimoto and Takemura 2014). The propagation feature of long-period ground motion is complicated due to the lateral variations of sedimentary structures within the Kanto Basin (e.g., Koketsu and Kikuchi 2000;Furumura and Hayakawa 2007;Mukai et al. 2018). Long-period ground motion and the structural properties within the Kanto Basin have been studied extensively using observed seismograms and geological surveys to contribute to disaster mitigation in the Tokyo metropolitan area (e.g., Koketsu et al. 2009;Takemura et al. 2015).
Recent advances in numerical simulation codes (e.g., Gokhberg and Fichtner 2016;Maeda et al. 2017) and local/regional three-dimensional (3D) velocity structure models (e.g., Koketsu et al. 2012;Kennett et al. 2013;Stephenson et al. 2017) have enabled the implementation of realistic 3D simulations of long-period ground motion (e.g., Komatitsch 2004;Iwaki et al. 2018;Wirth et al. 2019;Takemura et al. 2019a), and the estimation of structural properties (e.g., Tape et al. 2009;Gao and Shen 2014;Miyoshi et al. 2017). In forward and inverse modeling of long-period ground motion and structural properties along propagation paths, an assumption of a double-couple point source is usually assumed. The centroid moment tensor (CMT) solutions based on displacement for periods longer than 20 s are generally considered robust against structural heterogeneities, compared to first-motion solutions (e.g., Takemura et al. 2016). As such, one-dimensional (1D) velocity models are adopted in local/global CMT inversion systems (e.g., Kubo et al. 2002;Bernardi et al. 2004;Vallée et al. 2011;Ekström et al. 2012); these solutions are typically used in ground motion simulations. However, in regions with strong heterogeneities, such as thick sediments and subducting oceanic plates, focal mechanisms could be incorrectly estimated using conventional CMT methods with a 1D velocity model. To address this issue, the CMT inversion based on Green's functions using the local/regional 3D model has been developed in such regions (e.g., Lee et al. 2013;Hejrani et al. 2017;Okamoto et al. 2018;Takemura et al. 2018aTakemura et al. , 2018bTakemura et al. , b, 2020Wang and Zhan 2020;Hejrani and Tkalčić 2020). By using the 3D CMT results of moderate earthquakes along the Nankai Trough, Takemura et al. (2020) demonstrated that the differences in centroid depths and focal mechanisms between 1D and their 3D CMT solutions were significant for offshore earthquakes due to offshore heterogeneities. These differences could affect ground motion simulations (e.g., Takemura et al. 2019c). To achieve precise forward and inverse modeling of long-period ground motions in the Kanto region, where large sedimentary basin and two subducting plates exist, accurate CMT solutions is required.
In this study, we conduct CMT inversions of moderate earthquakes in the Kanto region based on the Green's function calculated using the local 3D model. We evaluate differences in source parameters between CMT solutions obtained from 1D and 3D earth model. We conduct ground motion simulations using 3D CMT solutions to discuss the effects of CMT solutions on longperiod ground motion modeling in the Kanto Basin. To accurately model phases and amplitudes of long-period ground motion, we demonstrate that the adjusted source model should be incorporated in the used 3D model.

Method
In this study, we used the F-net broadband seismograms of the target earthquakes. At each F-net station (filled triangles in Fig. 1), a broadband velocity seismometer (STS-1, STS-2, or STS-2.5) was installed. The health of sensors has been systematically monitored by the National Research Institute for Earth Science and Disaster Resilience (NIED; Okada et al. 2004;Kimura et al. 2015; National Research Institute for Earth Science and Disaster Resilience 2019). Our target earthquakes were shallow (≤ 50 km) earthquakes with moment magnitudes (Mw) between 4.5 and 6.5, listed in the F-net moment tensor (MT) catalog. The F-net 1 D velocity structure model has been used in the F-net MT catalog (Fukuyama et al. 1998;Kubo et al. 2002). Target earthquakes (focal mechanisms in Fig. 1) occurred within the area of assumed source grids (crosses in Fig. 2a) between April 2017 and March 2020. The data of the Metropolitan Seismic Observation Network (MeSO-net) were also available from the NIED website in the analyzed period (e.g., Kasahara et al. 2009;Sakai and Hirata 2009). As MeSOnet stations were densely deployed around the Tokyo metropolitan area in Japan (inverse triangles in Fig. 1), we also evaluated long-period ground motion in the Kanto Basin using earthquakes that occurred after April 2017.
CMT inversions were conducted using displacement seismograms for periods of 25-100 s. By using the Opensource Seismic Wave Propagation Code software "Open-SWPC" (Maeda et al. 2017), we numerically evaluated Green's functions in the 3D model. We used the Japan Integrated Velocity Structure Model version 1 (JIVSM; Koketsu et al. 2012) as the 3D velocity structure model for this study. The JIVSM was constructed by combine use of geological and geophysical surveys around Japan. The structures beneath the bedrock were modeled using regional-scale tomography and gravity data (e.g., Ryoki 1999;Matsubara et al. 2008;Matsubara and Obara 2011). Detailed procedure for model construction was described in Koketsu et al. (2009Koketsu et al. ( , 2012. The JIVSM has been widely used in the evaluation of ground motion, crustal deformation, and seismic monitoring across Japan (e.g., Guo et al. 2016;Miyazawa 2016;Agata 2020;Oba et al. 2020;Baba et al. 2020). The simulation model covered an area of 600 × 600 × 160 km 3 (blue dashed rectangle in Fig. 2a), which was discretized by grid intervals of 0.5 km in the horizontal direction and 0.2 km in the vertical direction. The physical parameters of each layer in the JIVSM are listed in Table 1. The minimum S-wave velocity in the solid column was assumed to be 1.5 km/s for the calculation of 3D Green's functions. Source grids were uniformly distributed at a horizontal interval of 0.1° and vertical interval of 2 km from depth 6 to 60 km. In calculation of Green's functions, the Küpper wavelet with a duration of 1 s was employed. The "Seismic Analysis Code" (SAC; Goldstein and Snoke 2005;Helffrich et al. 2013) was used to store simulated seismograms and conduct a part of signal processing. We discarded the source grids within seawater column in Green's function evaluation and stations with epicentral distances equal to or larger than 550 km. Approximately 13 000 000 Green's function SAC files from 49 279 source grids to 15 F-net stations were obtained by conducting 45 reciprocal calculations. Each reciprocal calculation required 247 GB of computer memory and a wall-clock time of 1 h by parallel computing with 432 cores of the computer system of Earthquake and Volcano Information Center at the Earthquake Research Institute, the University of Tokyo.
We choose F-net stations within epicentral distances between 100-400 km from the initial epicenter, obtained from the F-net MT catalog. A set of Green's functions at the source grids, which were located in ± 0.4° grids from the initial epicenter and distributed at depths between 6 and 60 km, were selected for grid search inversion. To apply band-pass filter and integration stably, we used 10-min F-net velocity seismograms from three minutes before the initial origin minutes. A 200 s time window for each CMT inversion was adopted. CMT inversions were conducted at every 1 s within ± 1 min from the origin minute to determine the centroid time. Time shifts, which adjust travel time differences between observed and synthetic seismograms at each station, were not permitted during grid search inversion, despite typical 1D MT routines such as the F-net MT system, enabling time shifts at each station between synthetic and observed seismograms. After CMT inversions, we obtained seismic moments and focal mechanisms at all possible spatial and temporal source grids. Then, to identify the optimal solution from all CMT solutions, we evaluated variance reductions (VRs) between the observed and synthetic displacement seismograms for periods of 25-100 s. If there was a perfect match between the observations and synthetics, the VR is 100%. The maximum VR solution was  (Koketsu et al. 2012) considered the optimal solution, providing the optimal centroid location, depth, time, focal mechanism, and seismic moment. Other technical details of CMT inversions using the 3D model and the evaluation of Green's functions are described in Takemura et al. (2020).
We conducted forward simulations of ground motion in the Kanto region to investigate the effects of CMT solutions on long-period (> 5 s) ground motion using CMT solutions based on the 3D model. A more realistic velocity model, including velocity layers with Vs slower than 1.5 km/s, should be used in simulations of longperiod ground motion. As such, the model in long-period ground motion simulations included all sedimentary layers of the JIVSM listed in Table 1. The model covered an area of 480 × 480 × 100 km 3 (red dotted rectangle in Fig. 2a), which was discretized by grid intervals of 0.2 km in the horizontal direction and 0.1 km in the vertical direction. Simulated seismograms were evaluated at the F-net and MeSO-net stations within the calculation region (red dotted rectangle in Fig. 2a). The source parameters for events a, b, and A-C are listed in Table 3. The centroid locations, times, and moment tensors from CMT solutions based on the 3D model were used for the listed earthquakes. The source time functions were characterized by the Küpper wavelet, and durations were expected by the empirical relationship between seismic moments and rupture durations (Kanamori and Brodsky 2004). Each simulation of long-period ground motion required 1190 GB computer memory and 4 h computation time; simulations were performed by parallel computing with 64 nodes of the Fujitsu PRIMERGY CX600M1/CX1640M1 (Oakforest-PAC) in the Information Technology Center, University of Tokyo. F-net MT solution of this earthquake was also a similar focal mechanism; however, its centroid depth was slightly (2 km) deeper than that of our solution (see numbers above focal spheres in Fig. 2b). The spatial variation of the VRs at each source grid is illustrated in Fig. 2c. The resolution against centroid time is shown in Additional file 1: Figure S1. Although the optimal depth was very close to the upper surface of the Pacific Plate, high (> 80%) VR solutions appeared at a wider depth range (36-56 km). These features are similar to CMT solutions in the Hyuga-nada region, southwest Japan ( Fig. 4 of Takemura et al. 2020). The centroid depth constraint is also not high for shallower crustal earthquake (Additional file 1: Figure S2). The synthetic seismograms of the optimal solution accurately reproduce the observations. We also showed CMT solution of this earthquake using periods of 10-100 s (Additional file 1: Figure S3). As similar to Hejrani and Tkalčić (2020), we confirmed that broader period analysis is important for constraining centroid depths. However, we employed the period band of 25-100 s for CMT inversions because the optimal solutions were not changed and the maximum VR was smaller than that for periods of 25-100 s.

Results: 3D CMT solutions of moderate earthquakes beneath Kanto region
We obtained 74 CMT solutions for shallow earthquakes (< 50 km) with an Mw of 4.2-6.3; Fig. 3a illustrates the spatial distribution of these CMT solutions. All parameters for these CMT solutions are available from the data repository site (https ://doi.org/10.5281/ zenod o.39268 84). Hereafter, we simply call CMT solutions based on the local 3D model as "3D CMT solutions. " Cross-sections of profiles A and B are also plotted at the bottom of Fig. 3a. Along profile A (around 37.0ºN), distributions of CMT solutions between 3D and 1D catalogs show different futures. However, unfortunately, detailed discussions are difficult because of the lack of seismicity studies by ocean bottom seismometers along profile A. At the north region of profile A (around 37.5ºN), Shinohara et al. (2005) detected earthquakes of magnitudes ranging from 0.8 to 4.2 from 8  Takemura et al. Earth, Planets and Space (2021) 73:15 August to 24 August 1997 using temporal ocean bottom seismometers. According to their results, the reverse faulting solutions off the Fukushima in our CMT catalog (blue focal spheres) could be considered as the lower seismicity in the double seismic zone within the subducting Pacific slab. In our analysis period, many shallow (~ 20-30 km) small earthquakes were listed in the JMA catalog, but no CMT solutions for moderate earthquakes were located in the upper seismicity in the double seismic zone. Along with profile B, many earthquakes occurred just below the upper surface of the Pacific Plate. To validate our CMT solutions, especially in epicenter locations and depths, we compared the hypocenter distributions based on temporal ocean bottom seismometers (Ito et al. 2017a, b), which covered the area around the profile B. This seismicity was also confirmed in their hypocenter distribution. Although we could not compare detailed hypocenter distribution studies in other offshore regions (e.g., along profile A, around 36.0ºN), we consider that the 3D CMT inversion works well in offshore region of the Kanto region. The percentages of isotropic, compensated linear vector dipole (CLVD), and double-couple (DC) components in both catalogs are illustrated in Additional file 1: Figure S4. Differences between our and F-net catalogs could be caused by differences in velocity model and moment tensor inversion settings. In the F-net routine system, the isotropic component is not considered.

Differences between 3D CMT and F-net MT catalogs
The F-net solutions of corresponding earthquakes are also plotted in Fig. 3b. Spatial distributions of both catalogs seem to be similar. To quantitatively evaluate differences between the 3D CMT in this study and the F-net MT catalogs, we calculated cross-correlation coefficients of P-wave radiation patterns (e.g., Kuge and Kawakatsu 1993;Helffrich 1997), depth, and Mw differences between the 3D CMT and F-net MT catalogs (Fig. 4). A large negative value (-0.6) of cross-correlation coefficients only appeared in a solution for an earthquake on February 23, 2019. Only 4 F-net stations (N.JIZF,N.KZKF,N.YMZF,and N.KSKF;see Fig. 2b) were applicable for the CMT inversion of this earthquake because of the low signal-to-noise ratio (SNR) for 25-100 s periods. Additionally, the VR of the 3D CMT solution was not high (~ 57%). With the exception of this event, differences in focal mechanisms and centroid depths were not significant compared to offshore earthquakes along the Nankai Trough (Fig. 8 of Takemura et al. 2020).
On the other hand, we found that the Mw values based on the 3D CMT were systematically smaller than those of the F-net MT catalog (Fig. 4c). The Mw values are very important for ground motion simulations because values of seismic moments are directly related to the amplitude of the simulated ground motion. Using the 3D CMT catalog along the Nankai Trough ; https ://doi.org/10.5281/zenod o.36741 61), we also evaluated

Table 1 Physical parameters of each layer in JIVSM
The air and seawater layers were treated as being the same, following Maeda et al. (2017). The P-wave velocity (V P ), S-wave velocity (V S ), density (ρ), rigidity (μ) and inelastic attenuation (Q P and Q S ) are listed the differences in the Mw between 3D CMT and F-net MT solutions. We found both larger and smaller Mw values compared to the F-net catalog in the Nankai region (Fig. 5). In the Kanto and Nankai regions, the differences in Mw for offshore earthquakes were larger than those of onshore earthquakes; these differences may be caused by 3D heterogeneities.
To investigate the cause of these Mw differences, we conducted ground motion simulations for Mw 4.36 earthquake on 10:02:45 (JST) November 17, 2017 (Event a) and Mw 4.35 earthquake on 19:56:05 (JST) on August 6, 2018 (Event b). Using the 3D CMT method, Events a and b were located just below the upper surface of the oceanic crust layer 2 and the boundary between oceanic crust layers 2 and 3 of the Pacific Plate, respectively. The Mw differences for events a and b were -0.31 and -0.25, respectively, and the estimated seismic moments of the 3D CMT solutions were approximately 35% and 42% of the F-net 1D solutions, respectively. We conducted simulations using the same source models and three different heterogeneous models; the JIVSM (Koketsu et al. 2012), the JIVSM without sediments, and the F-net 1D model (Kubo et al. 2002). The source models were the optimal solutions of 3D CMT inversion for two earthquakes (Events a and b in Table 3). Figure 6 compares the simulated and observed vertical velocity seismograms. Two F-net stations were selected, and other simulation results were stored at https ://doi. org/10.5281/zenod o.39268 88. We found that simulation results using the JIVSM and the JIVSM without sediments reproduced observed F-net seismograms, with the exception of the N.JIZF seismograms for the JIVSM without sediments. This suggests that the effects of lowvelocity sediments around the Kanto region on CMT inversion using long-period (25-100 s) seismograms are minor. On the other hand, in the Nankai region, thick low-velocity accretionary prism has significant influence on seismograms even for periods of 25-100 s. This difference could be explained by differences in thickness between the Kanto basin (~ 3 km) and accretionary prism (~ 5 km). Because the Kanto Basin and marine sediments exist along the path from event b to N.JIZF, the difference in waveforms of N.JIZF for periods of 10-50 s might appear. Around the Nankai Trough, a thicker (> 5 km) accretionary prism has a significant influence on surface waves even for periods longer than 20 s, and consequently the affects the results of CMT inversions and ground motion simulations (e.g., Nakamura et al. 2015;Takemura et al. 2018aTakemura et al. , b, 2019aTakemura et al. , 2019bTakemura et al. , 2020. On the other hand, the amplitudes of simulation seismograms with a similar source and the F-net 1D model were approximately 35-45% of the observed amplitudes. The effects of the Kanto Basin have a minor influence on ground motion at outcrop rock sites (F-net), and differences in mechanisms and depths compared with F-net solutions that are not significant. This difference could be explained by differences in heterogeneities around the seismic source. The 3D CMT solutions of events a and b were located just beneath the upper surface of the oceanic crust layer 2 and near the boundary between oceanic crust layers 2 and 3 of the Pacific Plate, respectively. In the JIVSM (Table 1), the rigidities of source areas for both events were 20.4-34.3. In contrast, the rigidity at depths between 33 and 100 km was a uniform value (63.7 Fig. 4 Differences between 3D CMT and F-net MT catalogs in the Kanto region. a Cross-correlation of P-wave radiation patterns (e.g., Kuge and Kawakatsu 1993;Helffrich 1997) between 3D CMT and F-net MT catalogs. b Depth differences and c Mw differences from F-net MT solutions Takemura et al. Earth, Planets and Space (2021) Table 1). Locations of seismic sources and stations are illustrated in the upper left panel. Plotted focal mechanisms are the optimal solutions of events a and b GPa; Table 2) in the F-net 1D model. The differences in rigidities around source regions between the JIVSM and the F-net 1D model correspond to differences in seismic moments between the 3D CMT and F-net MT solutions (34-42%). We additionally simulated long-period ground motions for CMT solutions with Mw differences equal to or smaller than − 0.2. Simulated earthquakes are listed in Additional file 1: Table S1. Amplitude differences at each station were calculated by dividing maximum simulated amplitudes in the F-net 1D model by those in the JIVSM. Maximum amplitudes were measured using vertical seismograms for periods of 25-100 s. At almost stations, amplitudes using the JIVSM were larger than those using the F-net 1D model (Fig. 7a). Due to structural differences and radiation pattern of the assumed sources, relationship between Mw and amplitude differences exhibited simple linear trend, but amplitude differences tended to increase with decreasing Mw difference (− 0.2 to − 0.3). We also compared the relationship between Mw differences and rigidity structures (Fig. 7b, c). The rigidity differences were evaluated by differences of the JIVSM rigidities at centroid depths from the F-net 1D model. We confirmed the positive correlation between rigidity and Mw differences. In the Nankai region (Fig. 7c), because centroid depths were changed in the cases of offshore earthquakes (Fig. 8 of Takemura et al. 2020), the values of Mw differences were scattered. As such, it may be concluded that the major cause of differences in seismic moments between the 3D CMT and F-net 1D MT solutions is the difference in rigidity around the source areas.
For the Nankai Trough, both overestimations and underestimations of seismic moments compared to the F-net catalog were observed (Fig. 5). Large Mw differences only appeared in the offshore region, where many intraslab and interplate earthquakes occurred. In particular, intraslab earthquakes along the Nankai Trough occurred within the low-velocity oceanic crust and highvelocity oceanic mantle (see Figs. 5 and 6 of Takemura et al. 2020), not modeled in the F-net 1D model. The difference in Mw values along the Nankai Trough could also be explained by the differences in heterogeneous structures between the 3D and 1D models.
In the F-net routine system, the origin times and epicenters were fixed as those in the JMA unified hypocenter catalog, and time shifts between observed and synthetic seismograms at each station were enabled. Miyoshi et al. (2017) notes that prior to estimating structural properties, the re-evaluation of centroid times for F-net MT solutions should be required to obtain suitable waveform inversion results. In this study, we found that the estimation of seismic moments was affected by the rigidity structure around the source region. The difference in the estimation of seismic moments directly impacts the amplitude of ground motion simulations. The amplitude of ground motion simulation is important to evaluate seismic hazards and estimate structural properties along propagation paths. The 3D CMT catalog allow us to directly conduct 3D ground motion simulations without adjusting source parameters, such as centroid depth, focal mechanism and seismic moment.
The earthquakes with large Mw differences tend to be located within the oceanic crust and mantle (Figs. 4 and 5). The precise fault size evaluation is one of important issues for evaluating seismicity within the slab and expecting seismic hazard due to intraslab earthquakes. The 3D CMT catalog can provide precise values of seismic moments of interaslab earthquakes because the effects of 3D heterogeneous rigidity structures can be evaluated by numerical simulations during preparing Green's function dataset.

Long-period ground motion simulations in the Kanto region
By using our 3D CMT solutions based on the JIVSM, we conducted numerical simulations of long-period ground motions and compared with the observed seismograms. For the SNR of the MeSO-net for periods longer than 5 s, three earthquakes were selected with an Mw equal to or larger than 5.5 for simulations of long-period ground motion in the Kanto Basin. The source parameters of selected events (A-C) are listed in Table 3. Complete files of simulated velocity waveforms and wavefields are available online https ://doi.org/10.5281/zenod o.39268 88. Figure 8 shows an example of simulated vertical velocity wavefields for the simulation of Event A at 40, 60, 80, 100, 120, 140, 160, and 180 s from the earthquake origin (movie file is also available from https ://doi.org/10.5281/ zenod o.39268 88). The seismic waves radiated from the source complicatedly propagate through the Kanto region. In the Kanto, Niigata, and offshore regions, the wavelengths and propagation speeds of the Rayleigh waves became shorter and slower due to low-velocity sediments. The energy of these shorter-wavelength components (i.e., long-period ground motion) was trapped within low-velocity sediments. Thus, the duration of long-period ground motion was elongated in the Kanto, Niigata, and offshore regions (lapse time of 180 s). Peak ground velocities (PGVs) were calculated by the vector sum of three-component filtered seismograms at the F-net and MeSO-net stations; the passband period was 5-30 s. Figure 9 shows the spatial distributions of PGVs for each event. With the exception of Event C, the simulations were able to roughly reproduce the observed Fig. 7 Amplitude and Mw differences between the JIVSM and F-net 1D model. a Amplitude differences of simulation results from the JIVSM and F-net 1D models. Circle, diamond, inverse triangle, square, and star represent amplitude ratio or station location of N.ASIF, N.TSKF, N,ONSF, N.TTOF, and N.JIZF, respectively. Amplitudes ratios were calculated by dividing maximum amplitudes of simulations using the F-net 1D model by those using the JIVSM. Maximum amplitudes were measured using vertical seismograms for periods of 25-100 s. The relationships between Mw differences and rigidity structures in (b) Kanto and (c) Nankai regions. The rigidity differences were evaluated by differences of the JIVSM rigidities at centroid depths from the F-net 1D model. The Mw differences in (c) was made from the data of Takemura et al. (2020) features of PGVs. Large PGVs appeared in regions with bedrock depths greater than 3 km. Figures 10, 11, and 12 compare the filtered seismograms (15-30, 10-20, and 5-16 s) between observations and simulations. We selected two F-net stations and three MeSO-net stations. The selected MeSO-net stations are located at the site with deeper (> 3 km) bedrock depths. Simulations reproduced the observed seismograms at two F-net stations (N.ASIF and N.JIZF), with the exception of the simulation results for Event C for periods of 5-16 s. These results suggest that these 3D CMT solutions have the ability to reproduce observed ground motion with sufficient accuracy for periods longer than 5 s at F-net stations, which are deployed at outcrop rock sites. Although the observed seismograms at MeSOnet stations were reproduced by the simulated seismograms for periods longer than 10 s with the exception of the later phases at E.YROM, the simulation results for 5-16 s periods could not explain the observed seismograms in the Kanto Basin. This period band includes the

Lon. [°]
Lat dominant period (~ 6 s) of long-period ground motions within the Kanto Basin (e.g., Yoshimoto and Takemura 2014). These seismogram discrepancies at the MeSOnet stations could be caused by the JIVSM sedimentary structure. For Event C, because the centroid depth was 8 km, the ground motion for 10-20 and 5-16 s periods were affected by both the Kanto sedimentary basin and oceanic sediments from the epicenter to coastal regions. Thus, these sedimentary structures in the offshore region may decrease waveform fitness for this event. The overestimation of PGVs (Fig. 9c) may also be attributed to the models of the Kanto Basin and the oceanic sediments along propagation paths. We also conducted simulations of long-period ground motions using F-net MT solutions. The correlation coefficients of P-wave radiation patterns for Events A-C are 0.99, 0.91, and 0.91, respectively, and depth differences from F-net solutions are 9, − 2, and − 6 km, respectively. Simulation results of F-net solutions (Additional file 1: Figures S5-7) could not reproduce observed seismograms, and discrepancies from observations could not be corrected by simple adjustment of seismic moments. This result implies that although differences in focal mechanisms and depths between our and F-net 1D catalogs were not so significant, these differences improved the reproducibility of broadband (5-100 s) observed seismograms at stations with/without the Kanto Basin. The 3D CMT catalog allows us to directly conduct the precise forward and inverse modeling of long-period ground motion without adjusting source models, Fig. 9 Peak ground velocity (PGV) for periods of 5-30 s. a Event A occurred on August 4, 2019, b Event B occurred on January 3, 2020, and c Event C occurred on February 6, 2020. Detail source parameters described in Events a-c of Table 3. The top and bottom panels are the synthetic and observed PGVs, respectively. PGVs were evaluated by calculating the vector sum of three-component filtered seismograms at the F-net and MeSO-net stations. Plotted focal mechanisms are the optimal solutions of the 3D CMT inversions for the corresponding events. The JIVSM bedrock depths are illustrated in Figs. 1, 7 and 9, 10 and 11

Conclusions
We conducted CMT inversions of moderate earthquakes in the Kanto region from April 2017 to March 2020. The estimated focal mechanisms and depths using the 3D CMT method were not significantly different from the corresponding F-net MT solutions. However, the Mw values were systematically smaller than those in the F-net catalog. Earthquakes with large Mw differences tended to be located within the subducting plate, i.e., intraslab earthquakes. Using numerical simulations with 3D and 1D velocity models, we concluded that the major cause for the Mw differences is the difference in rigidity between the 1D and 3D velocity models. The 3D subducting oceanic crust and mantle could not be modeled in the 1D CMT system. The 3D CMT method provide precise values of seismic moments including the effects of 3D rigidity structure and thus, the 3D CMT catalog is important for evaluating characteristics of intraslab seismicity. The differences in the estimation of seismic moments directly affect the amplitude of ground motion simulations. The 3D simulation of an intraslab earthquake using the 1D CMT catalog could cause overestimations in the amplitude, even at outcrop rock sites. The 3D CMT solutions should be adopted for precise forward and inverse modeling of long-period ground motion simulations. The 3D CMT solutions can be directly used in ground motion simulations without adjusting source parameters, such as centroid depth, focal mechanism and seismic moment. Simulations using the 3D CMT solutions in this study and the JIVSM were able to reproduce ground motion for periods longer than 5 s at outcrop rock sites. On the other hand, discrepancies between observations and simulations using F-net solutions could not be corrected by simple adjustment of seismic moments. These results mean that the 3D CMT inversion works well in the Kanto region, and 3D CMT solutions are suitable for modeling long-period (> 5 s) ground motion. However, while simulations at stations within the Kanto Basin reproduced observed seismograms for periods longer than 10 s, the reproducibility of these simulations decreased for periods shorter than 10 s.
In the Kanto region, because the predominant period of long-period ground motion is approximately 6 s, a wellconstrained sedimentary model is required to evaluate long-period ground motion for observed and anticipated large earthquakes. Recently, other regional/local velocity structure models of the sedimentary basin and subducting oceanic plate have been released (e.g., Hirose et al. 2008;Takemura et al. 2015; Headquarters for Earthquake Research Promotion 2017; Ito et al. 2019). The reproducibility of long-period ground motion in the Kanto Basin can be improved through the utilization of other published or adjusted models based on ground motion Fig. 11 Comparisons between observed and simulated seismograms for periods of 10-20 s. We selected two F-net (N.ASIF and N.JIZF) and three MeSO-net (E.YYIM, OK.KTGM, and E.YROM) stations. Amplitudes at each station were normalized by the maximum amplitudes of observed seismograms for periods of 10-20 s at each station. a Event A, b Event B, and c Event C. Gray and blue dotted lines are the observed and simulated seismograms, respectively. Left panels show used stations and bedrock depths of the JIVSM. Color for bedrock depths is similar as Fig. 1 simulations using 3D CMT solutions; this should be the primary focus of future research.

Supplementary information
The online version contains supplementary material available at https ://doi. org/10.1186/s4062 3-020-01348 -2. Table S1. Figure S1 shows CMT solutions at the optimal source grid as a function of centroid time. Figures S2 and S3 illustrate examples of CMT solutions of a shallow earthquake and broader period case. Figure S4 shows comparisons of percentages of each component between our 3D and F-net catalogs. Simulated waveforms of Events A-C using F-net MT solutions are illustrated in Figure S5-S7. The list of used earthquakes in Figure 7a is in Tabel S1.

Fig. 12
Comparisons between observed and simulated seismograms for periods of 5-16 s. We selected two F-net (N.ASIF and N.JIZF) and three MeSO-net (E.YYIM, OK.KTGM, and E.YROM) stations. Amplitudes at each station were normalized by the maximum amplitudes of observed seismograms for periods of 5-16 s at each station. a Event A, b Event B, and c Event C. Gray and blue dotted lines are the observed and simulated seismograms, respectively. Left panels show used stations and bedrock depths of the JIVSM. Color for bedrock depths is similar as Fig. 1