隱沒系統常伴隨著弧後擴張，但並非每一個隱沒系統皆有弧後擴張的發生。弧後擴張發生的時機和影響因素至今仍是一個待解的問題。隨著板塊的隱沒，含水礦物在一定溫壓條件下釋放水分並降低了周圍岩漿的熔點，使其更容易部分熔融產生岩漿。這些液體與熔體大大降低了地函楔(mantle wedge)內的強度，形成一低黏滯性地函楔(low viscosity wedge, LVW)。LVW的存在誘發強烈的朝向海溝方向且向上的地函流。上升的高溫地函物質使上覆地殼被拉薄，最後產生弧後擴張。本研究利用FLAC (Fast Lagrangian Analysis of Continua)的演算方式做數值模擬，並以參數化模擬熔體與液體造成地函楔與其上板塊弱化的效應，以探討影響弧後擴張行為的因素。LVW屈服應力和黏滯度會依據當區體積熔化比例下降，且分別降至f_ys和f_v倍。在此研究中對f_ys和f_v的數值進行系統化測試。
　　模擬結果顯示，弧後張裂活動的時間點和隱沒板塊的總長度有關。在板塊隱沒至地函楔後，大量的流體在地函楔中累積並降低當區強度形成LVW。持續回捲的板塊和LVW促使地函流上升產生並拉薄上覆地殼，最終產生弧後擴張中心。弧後擴張中心會隨時間離海溝越來越遠，且在到達一定距離後轉移到距海溝較近的地方。
　　我們可根據f_ys和f_v的不同將弧後擴張模型結果分成四個類型。(一)無弧後擴張(f_ys≥0.5、f_v≥0.5)：此類型中岩石圈的強度下降幅度小，弧後延伸與擴張等行為不易發生；(二)弧前擴張(f_ys>0.5、0.5>f_v≥0.01)：破裂帶沿著隱沒界面發展並延伸到地面，擴張中心發生在弧前；(三)一次弧後擴張(0.5>f_ys≥0.1)和(四)二次弧後擴張(f_ys<0.1)：由於流體的大幅弱化，上覆岩石圈強度驟減並且容易脆性破裂，擴散中心將在火山島弧區形成。比較四個類型與現實中的弧後擴張行為，表明在實際的弧後擴張階段弧後岩石圈的屈服應力需要急劇下降到50％或以下。這項研究可以幫助我們了解弧後區岩石圈的真實強度，同時也為未來的弧後擴散模型提供了可靠的岩石圈強度範圍。

　　Back-arc spreading is usually accompanied by slab rollback and trench retreat, which make overriding plate easier to extend. However, the occurrence and development of back-arc spreading have not been well studied. With the continually subduction of the slab, the aqueous mineral reaches higher temperature and pressure then releases water. The released water hydrates surrounding mantle and consequently induces partial melting of the surrounding mantle. The high concentrations of water in the mantle wedge could decrease the strength of the mantle, forming a low-viscosity wedge (LVW). The presence of LVW provides an opportunity for back-arc spreading to occur because of weak enough strength and induces a strongly upward mantle flow which also toward the trench. The rising hot mantle material causes the overlying crust to be thinned and eventually causes back-arc spreading.
　　In this study, we parameterize the bark-arc spreading processes using thermo-mechanical models with FLAC (Fast Lagrangian Analysis of Continua) technology in the subduction system, in order to explore the influencing factors affecting the back-arc spreading behavior. The generation of LVW is added to this model parametrically. The viscosity and yield stress is decreased by fv and fys times respectively depending on the regional volume melting ratio. We vary fys and fv systematically.
　　Our results show that the duration of back-arc spreading depends on the total length of subducted slab, regardless of convergence rate. With slab subducted, magma begins to build up and accumulates in mantle wedge, reducing the strength. The presence of the low viscosity wedge and the sinking slab induce strong upwelling in the wedge, which thins the overlying crust and eventually produces back-arc spreading. As time spreads, the spreading center will be farther away from the trench. A new center might form closer to the trench. The models with different fys and fv can roughly divide into 4 types. In type Ⅰ ( fys/ fv > 0.5), lithosphere strength is strong, resulting in no back-arc spreading. In type Ⅱ (fys >0.5, 0.5> fv ≥0.01), overriding plate still has great yield stress so the rupture zone develops along the subducting interface and then extends to the surface, forming a spreading center which is close to the trench. And for the type Ⅲ (0.5>fys≥0.1) and Ⅳ (0.1> fys), overriding plate get really weak due to magma and can rupture easily. The spreading center will form inside the overriding plate. Type 3 and 4 is closer to the reality, suggesting that the yield stress in back-arc need to drop sharply to 50-10-1% in the real back-arc spreading stage. This research can help us to understand the real strength of the back-arc region. It also provides a reliable parameter range for future back-arc spreading models.

致謝 i
摘要 iii
Abstract v
第一章 緒論 1
1.1 研究動機 1
1.2 弧後擴張數值模型 7
1.3 弧後擴張實例：馬里亞納隱沒系統 9
第2章 研究方法 14
2.1 數值模擬方法 14
2.1.1運動方程式 14
2.1.2流變行為 14
2.1.3相變、脫水與水合行為 15
2.1.4部分熔融事件 19
2.1.5火山噴發與潛熱作用 21
2.2 初始模型設計 25
2.3 模型使用參數 27
第3章 結果與討論 32
3.1 自發性隱沒 32
3.2 參考模型結果 36
3.3 其餘參數測試 53
3.3.1岩漿分布範圍測試 53
3.3.2 岩漿噴發比例測試 53
3.3.3 潛熱值測試 55
3.4 弧後擴張模型 56
第4章 討論 79
4.1 弧後擴張發生的類型 79
4.2 弧後擴張的發生強度 81
4.3 弧後擴張的島弧—盆地系統 86
4.4 弧後擴張的岩漿來源 88
4.5 熱構造的影響 90
第5章 結論 91
參考文獻 92

Anderson, M. O., Chadwick Jr, W. W., Hannington, M. D., Merle, S. G., Resing, J. A., Baker, E. T., ... & Augustin, N. (2017). Geological interpretation of volcanism and segmentation of the M ariana back‐arc spreading center between 12.7° N and 18.3° N. Geochemistry, Geophysics, Geosystems, 18(6), 2240-2274.
Barazangi, M., & Isacks, B. (1971). Lateral variations of seismic‐wave attenuation in the upper mantle above the inclined earthquake zone of the Tonga Island Arc: Deep anomaly in the upper mantle. Journal of Geophysical Research, 76(35), 8493-8516.
Billen, M. I., & Gurnis, M. (2001). A low viscosity wedge in subduction zones. Earth and Planetary Science Letters, 193(1-2), 227-236.
Brenan, J. M., Shaw, H. F., & Ryerson, F. J. (1995). Experimental evidence for the origin of lead enrichment in convergent-margin magmas. Nature, 378(6552), 54.
Calvert, A. J., Klemperer, S. L., Takahashi, N., & Kerr, B. C. (2008). Three‐dimensional crustal structure of the Mariana island arc from seismic tomography. Journal of Geophysical Research: Solid Earth, 113(B1).
Clark, S. R., Stegman, D., & Müller, R. D. (2008). Episodicity in back-arc tectonic regimes. Physics of the Earth and Planetary Interiors, 171(1-4), 265-279
Cundall, P. A. (1989). Numerical experiments on localization in frictional materials. Ingenieur-archiv, 59(2), 148-159.
Faccenna, C., Davy, P., Brun, J. P., Funiciello, R., Giardini, D., Mattei, M., & Nalpas, T. (1996). The dynamics of back-arc extension: an experimental approach to the opening of the Tyrrhenian Sea. Geophysical Journal International, 126(3), 781-795.
Faccenna, C., Holt, A. F., Becker, T. W., Lallemand, S., & Royden, L. H. (2018). Dynamics of the Ryukyu/Izu-Bonin-Marianas double subduction system. Tectonophysics, 746, 229-238.
Gerya, T. V., & Meilick, F. I. (2011). Geodynamic regimes of subduction under an active margin: effects of rheological weakening by fluids and melts. Journal of Metamorphic Geology, 29(1), 7-31.
Gribble, R. F., Stern, R. J., Bloomer, S. H., Stüben, D., O''Hearn, T., & Newman, S. (1996). MORB mantle and subduction components interact to generate basalts in the southern Mariana Trough back-arc basin. Geochimica et Cosmochimica Acta, 60(12), 2153-2166.
Grove, T. L., Till, C. B., Lev, E., Chatterjee, N., & Médard, E. (2009). Kinematic variables and water transport control the formation and location of arc volcanoes. Nature, 459(7247), 694.
Hacker, B. R. (1996). Eclogite Formation and the Rheology, Buoyancy, Seismicity, and H~ 2O Content of Oceanic Crust. GEOPHYSICAL MONOGRAPH-AMERICAN GEOPHYSICAL UNION, 96, 337-346.
Hall, C. E., Gurnis, M., Sdrolias, M., Lavier, L. L., & Müller, R. D. (2003). Catastrophic initiation of subduction following forced convergence across fracture zones. Earth and Planetary Science Letters, 212(1-2), 15-30.
Hasegawa, A., Zhao, D., Hori, S., Yamamoto, A., & Horiuchi, S. (1991). Deep structure of the northeastern Japan arc and its relationship to seismic and volcanic activity. Nature, 352(6337), 683.
Holt, A. F., & Becker, T. W. (2016). The effect of a power-law mantle viscosity on trench retreat rate. Geophysical Journal International, ggw392.
Jamieson, R. A., Unsworth, M. J., Harris, N. B., Rosenberg, C. L., & Schulmann, K. (2011). Crustal melting and the flow of mountains. Elements, 7(4), 253-260.
Jicha, B. R., & Jagoutz, O. (2015). Magma production rates for intraoceanic arcs. Elements, 11(2), 105-111.
Jolivet, L., Tamaki, K., & Fournier, M. (1994). Japan Sea, opening history and mechanism: A synthesis. Journal of Geophysical Research: Solid Earth, 99(B11), 22237-22259.
Katz, R. F., Spiegelman, M., & Langmuir, C. H. (2003). A new parameterization of hydrous mantle melting. Geochemistry, Geophysics, Geosystems, 4(9).
Katz, R. F., & Weatherley, S. M. (2012). Consequences of mantle heterogeneity for melt extraction at mid-ocean ridges. Earth and Planetary Science Letters, 335, 226-237.
Kincaid, C., & Hall, P. S. (2003). Role of back arc spreading in circulation and melting at subduction zones. Journal of Geophysical Research: Solid Earth, 108(B5).
Lallemand, S., Heuret, A., & Boutelier, D. (2005). On the relationships between slab dip, back‐arc stress, upper plate absolute motion, and crustal nature in subduction zones. Geochemistry, Geophysics, Geosystems, 6(9).
Nakakuki, T., & Mura, E. (2013). Dynamics of slab rollback and induced back-arc basin formation. Earth and Planetary Science Letters, 361, 287-297.
Nichols, G. T., Wyllie, P. J., & Stern, C. R. (1994). Subduction zone melting of pelagic sediments constrained by melting experiments. Nature, 371(6500), 785.
Roth, E. G., Wiens, D. A., & Zhao, D. (2000). An empirical relationship between seismic attenuation and velocity anomalies in the upper mantle. Geophysical Research Letters, 27(5), 601-604.
Schellart, W. P., Freeman, J., Stegman, D. R., Moresi, L., & May, D. (2007). Evolution and diversity of subduction zones controlled by slab width. Nature, 446(7133), 308.
Sdrolias, M., & Müller, R. D. (2006). Controls on back‐arc basin formation. Geochemistry, Geophysics, Geosystems, 7(4).
Stern, R. J. (2002). Subduction zones. Reviews of geophysics, 40(4), 3-1.
Straub, S. M. (2003). The evolution of the Izu Bonin‐Mariana volcanic arcs (NW Pacific) in terms of major element chemistry. Geochemistry, Geophysics, Geosystems, 4(2).
Takahashi, E., & Kushiro, I. (1983). Melting of a dry peridotite at high pressures and basalt magma genesis. American Mineralogist, 68(9-10), 859-879.
Tan, E., Lavier, L. L., Van Avendonk, H. J., & Heuret, A. (2012). The role of frictional strength on plate coupling at the subduction interface. Geochemistry, Geophysics, Geosystems, 13(10).
Tatsumi, Y., Sakuyama, M., Fukuyama, H., & Kushiro, I. (1983). Generation of arc basalt magmas and thermal structure of the mantle wedge in subduction zones. Journal of Geophysical Research: Solid Earth, 88(B7), 5815-5825.
Taylor, B., & Martinez, F. (2003). Back-arc basin basalt systematics. Earth and Planetary Science Letters, 210(3-4), 481-497.
Ulmer, P. (2001). Partial melting in the mantle wedge—the role of H2O in the genesis of mantle-derived ‘arc-related’magmas. Physics of the Earth and Planetary Interiors, 127(1-4), 215-232.
Ulmer, P., & Trommsdorff, V. (1995). Serpentine stability to mantle depths and subduction-related magmatism. Science, 268(5212), 858-861.
Wolf, S. G., & Huismans, R. S. Mountain building or backarc extension in ocean‐continent subduction systems‐a function of backarc lithospheric strength and absolute plate velocities. Journal of Geophysical Research: Solid Earth.
Yaxley, G. M., & Green, D. H. (1994). Experimental demonstration of refractory carbonate-bearing eclogite and siliceous melt in the subduction regime. Earth and Planetary Science Letters, 128(3-4), 313-325.
Zheng, Y. F., & Hermann, J. (2014). Geochemistry of continental subduction-zone fluids. Earth, Planets and Space, 66(1), 93.