<p>The Median Tectonic Line (MTL) is one of the largest active fault systems in Japan. Whether the active MTL fault is north-dipping or near-vertical is still under debate. In this study, we estimated the MTL fault geometry in eastern Shikoku based on interseismic displacements obtained via Interferometric Synthetic Aperture Radar (InSAR) and Global Navigation Satellite System (GNSS). Detecting interseismic displacements in eastern Shikoku with InSAR is challenging because dense vegetation and steep topography increase InSAR observation noises. Therefore, we used Advanced Land Observing Satellite 2 data, which are less affected by vegetation. For InSAR processing, we applied state-of-the-art atmospheric correction methods and conducted small baseline subset -based InSAR time series analysis for two frames in eastern Shikoku. Obtained InSAR displacement velocities were consistent with GNSS velocities, with a root mean square error of 1.49 mm/year in line-of-sight direction. GNSS and InSAR velocity fields primarily reflected the plate coupling effect on the Nankai subduction zone. The displacement velocity discontinuity was not clear along the MTL, suggesting that the shallow part of the MTL was almost fully locked. We then conducted a parameter study for the MTL to estimate deep slip rate <InlineEquation ID="IEq1"> <EquationSource Format="TEX">${S}_{\text{MTL}}$</EquationSource> <EquationSource Format="MATHML"><math> <msub> <mi>S</mi> <mtext>MTL</mtext> </msub> </math></EquationSource> </InlineEquation>, locking depth <InlineEquation ID="IEq2"> <EquationSource Format="TEX">${D}_{\text{MTL}}$</EquationSource> <EquationSource Format="MATHML"><math> <msub> <mi>D</mi> <mtext>MTL</mtext> </msub> </math></EquationSource> </InlineEquation>, and dip angle <InlineEquation ID="IEq3"> <EquationSource Format="TEX">${\delta }_{\text{MTL}}$</EquationSource> <EquationSource Format="MATHML"><math> <msub> <mi>δ</mi> <mtext>MTL</mtext> </msub> </math></EquationSource> </InlineEquation> using obtained velocity fields. Combinations of best-fit parameters (<InlineEquation ID="IEq4"> <EquationSource Format="TEX">${S}_{\text{MTL}}$</EquationSource> <EquationSource Format="MATHML"><math> <msub> <mi>S</mi> <mtext>MTL</mtext> </msub> </math></EquationSource> </InlineEquation>, <InlineEquation ID="IEq5"> <EquationSource Format="TEX">${D}_{\text{MTL}}$</EquationSource> <EquationSource Format="MATHML"><math> <msub> <mi>D</mi> <mtext>MTL</mtext> </msub> </math></EquationSource> </InlineEquation>, <InlineEquation ID="IEq6"> <EquationSource Format="TEX">${\delta }_{\text{MTL}}$</EquationSource> <EquationSource Format="MATHML"><math> <msub> <mi>δ</mi> <mtext>MTL</mtext> </msub> </math></EquationSource> </InlineEquation>) which explained GNSS and InSAR observations with corresponding PH subduction models were (9.0 mm/year, 13 km, 30°) and (4.0 mm/year, 8 km, 35°), and corresponding root mean square values of GNSS and InSAR residual velocities were 1.20 and 0.57 mm/year, respectively. There was a trade-off between <InlineEquation ID="IEq7"> <EquationSource Format="TEX">${S}_{\text{MTL}}$</EquationSource> <EquationSource Format="MATHML"><math> <msub> <mi>S</mi> <mtext>MTL</mtext> </msub> </math></EquationSource> </InlineEquation> and margin-parallel slip on the Nankai subduction zone, which prevented precise parameter estimation. Different combinations of the <InlineEquation ID="IEq8"> <EquationSource Format="TEX">${S}_{\text{MTL}}$</EquationSource> <EquationSource Format="MATHML"><math> <msub> <mi>S</mi> <mtext>MTL</mtext> </msub> </math></EquationSource> </InlineEquation>, <InlineEquation ID="IEq9"> <EquationSource Format="TEX">${D}_{\text{MTL}}$</EquationSource> <EquationSource Format="MATHML"><math> <msub> <mi>D</mi> <mtext>MTL</mtext> </msub> </math></EquationSource> </InlineEquation> and <InlineEquation ID="IEq10"> <EquationSource Format="TEX">${\delta }_{\text{MTL}}$</EquationSource> <EquationSource Format="MATHML"><math> <msub> <mi>δ</mi> <mtext>MTL</mtext> </msub> </math></EquationSource> </InlineEquation> produced spatially similar interseismic displacement fields due to the MTL, which made it difficult to find an optimal parameter combination. Nevertheless, our results generally supported the idea that the MTL is north-dipping rather than near-vertical, and shallow part of the MTL is almost fully locked in eastern Shikoku, which is consistent with a result from a previous study based on campaign GNSS observations.</p> Graphical Abstract <p></p>

Estimation of the fault geometry of the Median Tectonic Line (MTL) in Eastern Shikoku, Southwestern Japan, via GNSS and InSAR time series analysis

期刊论文
出版模式：

开放获取
发表日期：

2026年1月2日

查看全文

Kohei Shimotsuma,
Yohei Kinoshita,
Suguru Yabe,
Tadafumi Ochi

摘要

The Median Tectonic Line (MTL) is one of the largest active fault systems in Japan. Whether the active MTL fault is north-dipping or near-vertical is still under debate. In this study, we estimated the MTL fault geometry in eastern Shikoku based on interseismic displacements obtained via Interferometric Synthetic Aperture Radar (InSAR) and Global Navigation Satellite System (GNSS). Detecting interseismic displacements in eastern Shikoku with InSAR is challenging because dense vegetation and steep topography increase InSAR observation noises. Therefore, we used Advanced Land Observing Satellite 2 data, which are less affected by vegetation. For InSAR processing, we applied state-of-the-art atmospheric correction methods and conducted small baseline subset -based InSAR time series analysis for two frames in eastern Shikoku. Obtained InSAR displacement velocities were consistent with GNSS velocities, with a root mean square error of 1.49 mm/year in line-of-sight direction. GNSS and InSAR velocity fields primarily reflected the plate coupling effect on the Nankai subduction zone. The displacement velocity discontinuity was not clear along the MTL, suggesting that the shallow part of the MTL was almost fully locked. We then conducted a parameter study for the MTL to estimate deep slip rate ${S}_{\text{MTL}}$ $S_{MTL}$ , locking depth ${D}_{\text{MTL}}$ $D_{MTL}$ , and dip angle ${\delta }_{\text{MTL}}$ $δ_{MTL}$ using obtained velocity fields. Combinations of best-fit parameters ( ${S}_{\text{MTL}}$ $S_{MTL}$ , ${D}_{\text{MTL}}$ $D_{MTL}$ , ${\delta }_{\text{MTL}}$ $δ_{MTL}$ ) which explained GNSS and InSAR observations with corresponding PH subduction models were (9.0 mm/year, 13 km, 30°) and (4.0 mm/year, 8 km, 35°), and corresponding root mean square values of GNSS and InSAR residual velocities were 1.20 and 0.57 mm/year, respectively. There was a trade-off between ${S}_{\text{MTL}}$ $S_{MTL}$ and margin-parallel slip on the Nankai subduction zone, which prevented precise parameter estimation. Different combinations of the ${S}_{\text{MTL}}$ $S_{MTL}$ , ${D}_{\text{MTL}}$ $D_{MTL}$ and ${\delta }_{\text{MTL}}$ $δ_{MTL}$ produced spatially similar interseismic displacement fields due to the MTL, which made it difficult to find an optimal parameter combination. Nevertheless, our results generally supported the idea that the MTL is north-dipping rather than near-vertical, and shallow part of the MTL is almost fully locked in eastern Shikoku, which is consistent with a result from a previous study based on campaign GNSS observations.

Graphical Abstract