岩石圈弯曲的黏弹性松弛及其对确定弹性岩石圈厚度的影响

doi:10.7523/j.issn.2095-6134.2020.05.008

摘要/Abstract

摘要： 确定弹性岩石圈厚度在地球科学中具有重要意义。基于弹性板弯曲解析解结合岩石圈弯曲外缘隆起位置反演弹性岩石圈厚度的传统方法具有深远影响。利用变黏度可压缩黏弹性有限元模型定量研究岩石圈在表面垂向载荷作用下弯曲的黏弹性松弛问题。模拟结果显示，外缘隆起位置具有先向弯曲中心运动而后向弯曲外侧移动的特点，反向运动发生在载荷作用约9 Ma后。与基于松弛时间定义的弹性岩石圈厚度相比，基于弹性板弯曲解由外缘隆起位置约束弹性岩石圈厚度的方法仅在载荷作用后0.2~5.0 Ma的时间段内给出较为合理的结果。

关键词: 岩石圈弯曲, 弹性岩石圈, 黏弹性松弛, 数值模拟

Abstract: It is of significant importance to determine the thickness of elastic lithosphere in geoscience. The traditional method of inverting the thickness of elastic lithosphere constrained by the location of peripheral bulge in the flexural lithosphere based on the analytical solution of elastic plate flexure has far-reaching influence. In this work, the viscoelastic relaxation of the lithosphere flexure under vertical surface loads is studied quantitatively by means of a compressible viscoelastic finite element model with variable viscosities. The simulation results show that the position of the peripheral bulge moves first to the bending center and then to the outside of the bending, and the reverse motion occurs about 9 Ma after the loading. Compared with the elastic lithospheric thickness defined by relaxation time, the method of determining the thickness of the elastic lithosphere by the position of the peripheral bulge based on elastic plate bending solution gives reasonable results only in the period of 0.2-5.0 Ma after loading.

Key words: lithospheric flexure, elastic lithosphere, viscoelastic relaxation, numerical simulation

中图分类号:

邵翠法, 王世民. 岩石圈弯曲的黏弹性松弛及其对确定弹性岩石圈厚度的影响[J]. 中国科学院大学学报, 2020, 37(5): 640-649.

SHAO Cuifa, WANG Shimin. Viscoelastic relaxation of lithospheric flexure and its influence on determining thickness of elastic lithosphere[J]. , 2020, 37(5): 640-649.

参考文献

[1] Turcotte D L, Schubert G. Geodynamics[M]. 3rd ed. Cambridge, UK:Cambridge University Press, 2014.
[2] Diao F, Wang R, Wang Y, et al. Fault behavior and lower crustal rheology inferred from the first seven years of postseismic GPS data after the 2008 Wenchuan earthquake[J]. Earth and Planetary Science Letters, 2018, 495:202-212.
[3] Freed A M, Hashima A, Becker T W, et al. Resolving depth-dependent subduction zone viscosity and afterslip from postseismic displacements following the 2011 Tohoku-oki, Japan earthquake[J]. Earth and Planetary Science Letters, 2017, 459:279-290.
[4] Crawford O, Al-Attar D, Tromp J, et al. Forward and inverse modelling of post-seismic deformation[J]. Geophysical Journal International, 2017, 208(2):845-876.
[5] Hu Y, Buergmann R, Banerjee P, et al. Asthenosphere rheology inferred from observations of the 2012 Indian Ocean earthquake[J]. Nature, 2016, 538(7625):368-372.
[6] Hoechner A, Sobolev S V, Einarsson I, et al. Investigation on afterslip and steady state and transient rheology based on postseismic deformation and geoid change caused by the Sumatra 2004 earthquake[J]. Geochemistry Geophysics Geosystems, 2011, 12(7):Q07010, doi:10.1029/2010GC003450.
[7] Klein F W. Lithospheric flexure under the Hawaiian volcanic load:internal stresses and a broken plate revealed by earthquakes[J]. J Geophys Res-Solid Earth, 2016, 121(4):2400-2428.
[8] Caron L, Metivier L, Greff-Lefftz M, et al. Inverting glacial isostatic adjustment signal using Bayesian framework and two linearly relaxing rheologies[J]. Geophysical Journal International, 2017, 209(2):1126-1147.
[9] Tanaka Y, Klemann V, Martinec Z, et al. Spectral-finite element approach to viscoelastic relaxation in a spherical compressible Earth:application to GIA modelling[J]. Geophysical Journal International, 2011, 184(1):220-234.
[10] Zhong S J, Paulson A, Wahr J. Three-dimensional finite-element modelling of Earth's viscoelastic deformation:effects of lateral variations in lithospheric thickness[J]. Geophysical Journal International, 2003, 155(2):679-695.
[11] Capitanio F A, Morra G, Goes S. Dynamic models of downgoing plate-buoyancy driven subduction:subduction motions and energy dissipation[J]. Earth and Planetary Science Letters, 2007, 262(1/2):284-297.
[12] Capitanio F A, Morra G, Goes S. Dynamics of plate bending at the trench and slab-plate coupling[J]. Geochemistry Geophysics Geosystems, 2009, 10(4):Q04002, doi:10.1029/2008gc002348.
[13] Guest A, Schubert G, Gable C W. Stress field in the subducting lithosphere and comparison with deep earthquakes in Tonga[J]. J Geophys Res-Solid Earth, 2003, 108(B6):2288, doi:10.1029/2002jb002161.
[14] Gurnis M, Eloy C, Zhong S J. Free-surface formulation of mantle convection.2. Implication for subduction-zone observables[J]. Geophysical Journal International, 1996, 127(3):719-727.
[15] Hunter J, Watts A B. Gravity anomalies, flexure and mantle rheology seaward of circum-Pacific trenches[J]. Geophysical Journal International, 2016, 207(1):288-316.
[16] Watts A B, Zhong S J, Hunter J. The behavior of the lithosphere on seismic to geologic timescales[J]. Annual Review of Earth and Planetary Sciences, 2013, 41(1):443-468.
[17] Zhong S, Watts A B. Lithospheric deformation induced by loading of the Hawaiian Islands and its implications for mantle rheology[J]. J Geophys Res-Solid Earth, 2013, 118(11):6025-6048.
[18] Watts A B, Zhong S. Observations of flexure and the rheology of oceanic lithosphere[J]. Geophysical Journal International, 2000, 142(3):855-875.
[19] Zhong S J. Dynamics of crustal compensation and its influences on crustal isostasy[J]. J Geophys Res-Solid Earth, 1997, 102(B7):15287-15299.
[20] Gunn R. A quantitative evaluation of the influence of the lithosphere on the anomalies of gravity[J]. Journal of the Franklin Institute, 1943, 236:47-65.
[21] Barrell J. The strength of the earth's crust[J]. Journal of Geology, 1914, 22(8):729-741.
[22] Kim S S, Wessel P. Flexure modelling at seamounts with dense cores[J]. Geophysical Journal International, 2010, 182(2):583-598.
[23] Forsyth D W. Subsurface loading and estimates of the flexural rigidity of continental lithosphere[J]. Journal of Geophysical Research-Solid Earth and Planets, 1985, 90(B14):2623-2632.
[24] Wessel P. A reexamination of the flexural deformation beneath the Hawaiian-islands[J]. J Geophys Res-Solid Earth, 1993, 98(B7):12177-12190.
[25] Contreras-Reyes E, Garay J. Flexural modeling of the elastic lithosphere at an ocean trench:a parameter sensitivity analysis using analytical solutions[J]. Journal of Geodynamics, 2018, 113:1-12.
[26] Contreras-Reyes E, Osses A. Lithospheric flexure modelling seaward of the Chile trench:implications for oceanic plate weakening in the Trench Outer Rise region[J]. Geophysical Journal International, 2010, 182(1):97-112.
[27] Manriquez P, Contreras-Reyes E, Osses A. Lithospheric 3-D flexure modelling of the oceanic plate seaward of the trench using variable elastic thickness[J]. Geophysical Journal International, 2014, 196(2):681-693.
[28] 杨永章, 李金岭, 松本晃治, 等. 利用全月重力/地形导纳估计月球弹性层厚度[J]. 中国科学院大学学报, 2018, 35(4):463-467.
[29] Watts A B. Analysis of isostasy in the world's oceans. 1. Hawaii-Emperor seamount chain[J]. Journal of Geophysical Research, 1978, 83(NB12):5989-6004.
[30] Watts A B, Cochran J R. Gravity anomalies and flexure of lithosphere along Hawaiian-Emperor seamount chain[J]. Geophysical Journal of the Royal Astronomical Society, 1974, 38(1):119-141.
[31] Caldwell J G, Turcotte D L. Dependence of the thickness of the elastic oceanic lithosphere on age[J]. Journal of Geophysical Research, 1979, 84(NB13):7572-7576.
[32] Watts A B, Sandwell D T, Smith W H F, et al. Global gravity, bathymetry, and the distribution of submarine volcanism through space and time[J]. J Geophys Res-Solid Earth, 2006, 111(B8):B08408, doi:10.1029/2005jb004083.
[33] Walcott R I. Flexural rigidity, thickness, and viscosity of lithosphere[J]. Journal of Geophysical Research, 1970, 75(20):3941-3954.
[34] Sleep N H, Snell N S. Thermal contraction and flexure of mid-continent and Atlantic marginal basins[J]. Geophysical Journal of the Royal Astronomical Society, 1976, 45(1):125-154.
[35] Quinlan G M, Beaumont C. Appalachian thrusting, lithospheric flexure, and the paleozoic stratigraphy of the eastern interior of north-America[J]. Canadian Journal of Earth Sciences, 1984, 21(9):973-996.
[36] Mcnutt M K, Parker R L. Isostasy in Australia and evolution of compensation mechanism[J]. Science, 1978, 199(4330):773-775.
[37] Lambeck K, Nakiboglu S M. Seamount loading and stress in the ocean lithosphere:2. Viscoelastic and elastic-viscoelastic models[J]. Journal of Geophysical Research:Solid Earth (1978-2012), 1981, 86(B8):6961-6984.
[38] Courtney R C, Beaumont C. Thermally-activated creep and flexure of the oceanic lithosphere[J]. Nature, 1983, 305(5931):201-204.
[39] Watts A B. Isostasy and flexure of the lithosphere[M]. Cambridge:Cambridge University Press, 2001.
[40] Karato S, Wu P. Rheology of the upper mantle:a synthesis[J]. Science, 1993, 260(5109):771-778.
[41] Parsons B, Sclater J G. Analysis of variation of ocean-floor bathymetry and heat-flow with age[J]. Journal of Geophysical Research, 1977, 82(5):803-827.
[42] Comsol. COMSOL multiphysics programming reference manual[M]. COMSOL INC, 2016.