2. 中国地震局地球物理研究所,北京 100081;
3. 中国地质大学(北京)地球物理与信息技术学院,北京 100083
2. Institute of Geophysics, China Earthquake Administration, Beijing 100081, China;
3. School of Geophysics and Information Technology, China University of Geosciences, Beijing 100083, China
青藏高原被称为“世界屋脊",位于中国大陆西南,面积达240km2,是由不同地质历史时期的微大陆拼合起来的,是亚洲乃至全球变形规模最大、形成时间最晚的陆陆碰撞造山带,吸引了全世界地球科学家的极大关注[1-2].目前,国内外众多学者对青藏高原是由印度-欧亚两大板块于50~60 Ma碰撞形成的这一观点已达成共识,但其隆升变形机制仍是长期争论的重要科学问题[3].青藏高原东北缘作为印度-欧亚板块碰撞作用由近南北向向北东、东向转换的重要场所,成为研究高原隆升变形机制问题的关键地区之一.
研究区位于青藏高原东北缘地区,跨越青海、甘肃、宁夏、内蒙等省份部分地区,经纬度范围为东经102°-110°,北纬33°-39°,在我国大地构造上处于华北-塔里木块体、扬子块体和青藏块体等块体的交汇部位,如图 1所示.该区晚新生代以来构造变形十分强裂,广泛发育第三-第四系褶皱、逆冲和走滑断裂,强震频繁发生.
20世纪60年代以来,国内外学者对青藏高原东北缘地区的地壳速度结构、壳内低速层、电性结构和莫霍面形态等问题进行了广泛的研究,并取得许多重要认识与成果[4-10].人工反射地震、天然地震和大地电磁具有较大的探测深度和较高的垂向分辨率,但受施工成本高等因素影响,现阶段该区主要以二维剖面探测为主,较少开展大面积的三维面积性测量,因此无法全面整体地反映出研究区深部结构与构造.
重力测量相对简单,数据采集成本较低,具有覆盖面积广、快速高效的特点,历来是区域地质构造与动力学研究的重要方法之一.经过我国几十年的重力测量积累,到目前为止青藏高原东北缘地区已完全覆盖1∶100 万的重力测量,大部分地区已覆盖1∶50万甚至1∶20 万的重力测量,丰富且高精度的重力资料为研究该区深部结构和构造提供了极其重要的资料.梁桂培[11]、郭守年[12]、李清河[13]和江为为[14]等通过对该区重力资料的研究,得到了关于均衡重力异常特征、莫霍面深度分布和地壳密度结构等方面的研究成果.本文在系统收集该区重力资料和前人研究成果基础上,研究青藏高原东北缘地区不同尺度的重力异常特征,分析其横向构造格架,为探讨青藏高原隆升变形机制问题提供一定的重力学参考.
2 布格重力异常特征图 2显示了研究区布格重力异常,数据比例尺为1∶100万,网度为5km.由图 2可以清楚的看出,研究区内布格重力异常全部为负异常,从东到西呈逐渐减小趋势,幅值为-400mGal~ -100mGal.大致以东经106°线为界,有一条醒目的重力异常梯级带,其东西两侧异常特征在形态和走向上截然不同.
东侧布格重力异常变化平缓,主要呈北东和近东西走向;西侧异常值变化剧烈,走向以北西或北西西为主.该重力异常梯级带即是将中国大陆划分为东、西两大不同地貌和构造部分的贺兰山-六盘山-川滇南北构造带的北段[15].根据布格重力异常特征及随后讨论方便,将研究区划分为Ⅰ区(青藏高原东北缘重力异常区)、Ⅱ区(阿拉善东南缘重力异常区)、Ⅲ区(鄂尔多斯西南缘重力异常区)、Ⅳ 区(古浪、海原、银川边界带重力异常区)和v 区(扬子西北缘重力异常区)五个部分,下面分别介绍各部分布格重力异常特征.
青藏高原为全国布格重力异常值最低的负异常区,高原东部和北部被一条巨大的重力梯级带所围限[16].青藏高原东北缘(Ⅰ区)永登-兰州-武都一线附近的重力梯级带,就是上述巨大梯级带中的一部分,它是青藏高原东北缘在重力场上的反映.其西南侧表现为一系列次一级的重力梯级带,且沿北西向有扭曲现象,显示出高原隆起过程中其东北部经历了推覆、挤压和走滑等作用,物质自西向东“流动"[17].
阿拉善东南缘(Ⅱ区)表现为宽缓平静的负重力异常,幅值较低,与围绕青藏高原的巨型重力异常梯级带相连.由于研究区涉及阿拉善地块有限,本次未进行深入研究.
鄂尔多斯西南缘(Ⅲ区)以北东走向宽缓的高重力异常带为主,从块体西南边界陇县起,斜穿整个鄂尔多斯地块.块体内部定边以北地区,此重力异常高带走向转为近南北向,显示出鄂尔多斯块体南北两部分布格重力异常在延伸方向和走向等方面的差异.
古浪、海原、银川边界带(Ⅳ 区)内海原-六盘山断裂带,西起景泰西,途经海原、固原,向东南延伸至六盘山,呈弧形展布,弧顶在海原附近.在布格重力异常图上该断裂带反映为延伸很长的断续展布的重力异常梯级带,其中局部异常可能是火山岩沿断裂上涌的结果.青铜峡-固原断裂带在布格重力异常图上表现为清晰连续展布的近南北向重力梯级带,其东西两侧重力场特征明显不同.东侧重力异常平缓简单,走向主要为北东或北北东,西侧重力异常复杂多变,有扭曲现象,呈弧形展布,走向以北西、北北西为主,推测此断裂带为深大断裂并可能为鄂尔多斯块体西南边界.银川裂谷位于北东和北西走向的低异常区,布格重力异常特征基本反映其构造走向.
扬子西北缘(v 区)为秦岭北缘造山带,重力异常特征表现为一系列近东西向的高重力异常带,幅值较大,说明扬子板块与华北板块的碰撞、拼合对秦岭造山带进行了强烈的改造作用[18].
3 重力异常多尺度横向构造分析 3.1 异常多尺度分离和梯度分析重力异常是地表到地下深部所有密度不均匀体的综合反映,定性地讲,小尺度的重力异常与规模小、埋藏浅的地质体(比如地质构造或岩矿体)分布有关,而大尺度的重力异常则相反,与规模大、埋藏深的地质体分布有关.要根据重力异常研究某个特定地质体,必须首先从叠加异常中分离出单纯由这个地质体引起的异常,然后用这个异常进行反演解释.
传统向上延拓法是一种常用的异常分离方法,但它存在的主要问题是,在向上延拓时,浅源短波长信息和深源长波长信息同时延拓,即在浅源信息被压制的同时,深源信息也受到压制,使得异常分离不彻底.针对这个问题,意大利学者Pawlowski[19]根据维纳滤波和格林等效层原理提出了优选延拓法,应用该方法进行向上延拓时,可以在压制浅源短波长信息的同时,尽量保持深源长波长信息不衰减,使得异常分离相对彻底.然而,优选延拓法需要已知延拓高度,而实际应用中延拓高度往往是未知的.针对此问题,郭良辉等[20]对优选延拓法进行改进而提出了优化滤波法,它既保留优选延拓法分离异常相对彻底的特性,同时无需已知延拓高度,能实现对某一指定频段或尺度的局部异常有选择地分离.
本文采用优化滤波法对研究区布格重力异常(图 2)进行了多尺度分离,方法原理详见文献[20].根据研究区布格重力异常的径向对数功率谱形状特征(图 3),这里用0~0.0023cycles/km (频段1)、0.0023~0.0047cycles/km (频段2)、0.0047~0.0145cycles/km (频段3)、0.0145 ~ 0.0297cycles/km(频段4)、0.0297~0.0599cycles/km(频段5)和0.0599~0.1cycles/km(频段6)六个等效层分段拟合.其中,频段1~2对应大尺度重力异常,频段3对应中尺度重力异常,频段4~5对应小尺度重力异常,而频段6主要对应高频噪声.图 4a和图 4b分别显示了指定频段1~2 的优化滤波算子和分离出的大尺度重力异常,图 5a和图 5b分别绘制了指定频段3的优化滤波算子和分离出的中尺度重力异常,图 6a和图 6b则分别显示了指定频段4~5的优化滤波算子和分离出的小尺度重力异常.
根据场源视深度估计分析,频段1~2的大尺度重力异常主要反映深度41km 以下的莫霍面起伏和上地幔密度不均匀分布,频段3 的中尺度重力异常主要反映深度12~41km 的中、下地壳密度不均匀分布,频段4~5的小尺度重力异常主要反映深度12km 以上的上地壳及沉积盖层密度不均匀分布.
板块边界、不同地质体边界和断裂带等构造边界往往具有一定密度差异,造成重力异常变化较大,根据这一特点可利用重力方法研究横向线性构造.重力异常线性构造分析主要以重力异常导数(或梯度)为基础,利用异常方向导数、各阶垂直导数或总水平梯度的极值点、零点和其他特征点来识别构造边界.其中,重力异常总水平梯度法是一种常用的方法,其极大值在划分规模较大的线性构造时具有明显优势,但规模较小的线性构造往往会被规模较大的线性构造所淹没.针对此问题,王万银等[21]对总水平梯度法进行改进而提出了归一化总水平梯度垂直导数法,此方法对不同规模大小的线性构造都能较好反映,刻画的细节也较常规方法的丰富.
本文采用归一化总水平梯度垂直导数法对上述优化滤波法分离出的大尺度、中尺度和小尺度重力异常分别进行线性构造分析,结果分别见图 4c、图 5c和图 6c,图中同时显示了本文依据不同尺度重力异常及其线性构造特征综合推断出的五条断裂带,分别为:F1-秦岭地轴北缘断裂带,F2-海原-六盘山断裂带,F3-香山-天景山断裂带,F4-烟筒山断裂带和F5-青铜峡-固原断裂带.
3.2 大尺度异常横向构造特征由图 4b可见,分离出的大尺度重力异常体现了布格重力异常的区域特征,主要反映莫霍面起伏和上地幔深部场源的信息.大尺度重力异常以大致东经106°线为界,东西两侧异常特征在形态和走向上明显不同,说明研究区东西两侧深部密度结构具有显著差异.鄂尔多斯块体内部定边以北,重力异常高带走向由北东向转为近南北向,推测定边附近深部存在一个密度或构造界面,其两侧物质组成和构造特征具有差异.
图 4c所示的大尺度重力异常归一化总水平梯度垂直导数主要反映切穿莫霍面的深大断裂或构造边界.图中可以看出,本文推断的五条断裂带(F1~F5)除了部分与重力异常线性构造(高导数值所指示的)吻合外,大部分分布在低导数值区,推测这五条断裂带大部分没有切穿莫霍面,主要为壳内断裂.
3.3 中尺度异常横向构造特征由图 5b可以清楚的看出,分离出的中尺度重力异常主要反映深度12~41km 的中、下地壳不均匀密度结构.研究区鄂尔多斯地块内部分布近乎平行的四条重力低带和三条重力高带,重力低带、高带相间,异常单调,走向北东,且地块构造西南边界为清晰连续的重力异常梯级带,显示出刚性块体特征.与图 4b不同的是,图 5b中鄂尔多斯地块内部南北两部分重力异常特征相同,表明异常特征的差异不应该由中、下地壳引起的,而主要是由上地幔深部结构引起的.鄂尔多斯地块西南缘存在大范围的重力低带,永登-秦安一线走向北西,向南至宝鸡以南转为近东西向,并在汉中以北突转为北东向延伸,整体上看,负异常呈v 字形状,存在向秦岭造山带插入的趋势,与研究区磁异常特征基本一致[22],推测由于印度-欧亚板块碰撞及随后印度板块持续向北的挤压作用,造成青藏高原东北缘中、下地壳物质在巨大的北东向推挤力和鄂尔多斯刚性块体阻挡的共同作用下,沿着相对软弱的秦岭造山带方向蠕动.
图 5c所示的中尺度重力异常归一化总水平梯度垂直导数主要反映切穿中、下地壳的深大断裂或构造边界.从图可见,海原弧形构造区内海原-六盘山断裂带(F2)和青铜峡-固原断裂带(F5)分布在重力异常线性构造上及其附近,香山-天景山断裂带(F3)西侧分布在重力异常线性构造上,而烟筒山断裂带(F4)分布在低导数值区,因此推测海原-六盘山断裂带和青铜峡-固原断裂带深达下地壳,香山-天景山断裂带属次一级的断裂带,而烟筒山断裂带并不是中、下地壳断裂带.
3.4 小尺度异常横向构造特征由图 6b可见,分离出的小尺度重力异常主要反映深度12km 以上的上地壳及沉积盖层不均匀密度结构.可以看出,海原弧形构造区内断续展布的重力异常梯级带与本文推断的五条断裂带位置有较好的对应关系.图 6b与图 4b不同尺度重力异常对比可见,青藏高原东北部各块体深部边界位置与地表构造分布差异较大,反映出该区构造复杂,演化多期次、多样化,深浅构造不同.
图 6c所示的小尺度重力异常归一化总水平梯度垂直导数主要反映切穿上地壳或基底的断裂带或构造边界.由图可见,祁连-秦岭两大褶皱系的分界-秦岭地轴北缘断裂带(F1)[23]分布在同仁-陇西-西安一线走向北西西向的重力异常线性构造上.海原弧形构造区内由南西到北东的海原-六盘山断裂带(F2)、香山-天景山断裂带(F3)、烟筒山断裂带(F4)和青铜峡-固原断裂带(F5)分别分布在不同的重力异常线性构造带,与地表露头位置相吻合[24].秦岭地轴北缘断裂带(F1)在图 5c中、下地壳重力异常线性构造上显示的位置比图 6c上地壳重力异常线性构造中显示的位置偏向西南,因此推测此断裂带倾向西南,且倾角较缓.海原-六盘山断裂带(F2)在图 5c中、下地壳重力异常线性构造上显示的位置与图 6c上地壳重力异常线性构造中反映的位置基本一致,推断该断裂带倾角较大,近似直立.青铜峡-固原断裂带(F5)在图 5c中、下地壳重力异常线性构造上显示的位置比在图 6c上地壳重力异常线性构造中反映的位置向西移动,因此推测该断裂带倾向西,浅部较陡,且综合中小尺度重力异常特征,推测青铜峡-固原断裂带为研究区内鄂尔多斯地块西边界.
4 结论与讨论研究区大致以东经106°线为界,有一条醒目的重力异常梯级带,即贺兰山-六盘山-川滇南北构造带的北段,其东西两侧布格重力异常特征在形态和走向上截然不同.东侧异常变化平缓,主要呈北东和近东西走向;西侧异常变化剧烈,走向以北西或北西西为主.
通过多尺度重力异常对比分析,青藏高原东北缘各块体深部边界位置与地表构造分布差异较大,反映出该区构造复杂,演化多期次、多样化,深浅构造不同.由于印度-欧亚两大板块碰撞及随后印度板块持续向北的挤压作用,造成青藏高原东北缘中、下地壳物质在巨大的北东向推挤力和鄂尔多斯刚性块体阻挡的共同作用下,沿着相对软弱的秦岭造山带方向蠕动,整体上看,形成v字型负重力异常.
鄂尔多斯地块内部重力异常平缓简单,走向基本一致,与研究区内其他块体异常特征明显不同,且鄂尔多斯地块西南边界为清晰连续的重力异常梯级带,显示出刚性块体特征.块体内部定边向北,重力异常高带走向由北东向转为近南北向,推测定边附近存在一个密度或构造界面,其两侧物质组成和构造特征具有显著差异.对比大尺度重力异常和中尺度重力异常,表明异常特征的差异主要是由上地幔深部结构引起的,但上地幔深部结构差异是由什么引起,还需要近一步更深入的研究.
通过重力异常多尺度横向构造分析,综合推断出五条深断裂带,其中秦岭地轴北缘断裂带深达下地壳,倾向西南,且倾角较缓;海原-六盘山断裂带为深大断裂,倾角较大,近似直立;香山-天景山断裂带只有西侧分布在中尺度重力异常线性构造上,属于次一级的断裂带;烟筒山断裂带仅分布在小尺度重力异常线性构造上,为上地壳内断裂;青铜峡-固原断裂带为下地壳断裂,倾向西,浅部较陡,推测为研究区内鄂尔多斯地块西边界.
[1] | Yin A, Harrison T M. Geologic evolution of the Himalayan-Tibetan orogen. Annu.Rev.Earth Sci. , 2000, 28: 211-280. DOI:10.1146/annurev.earth.28.1.211 |
[2] | 许志琴, 杨经绥, 李海兵, 等. 青藏高原与大陆动力学—地体拼合、碰撞造山及高原隆升的深部驱动力. 中国地质 , 2006, 33(2): 221–238. Xu Z Q, Yang J S, Li H B, et al. The Qinghai-Tibet plateau and continental dynamics: A review on terrain tectonics, collisional orogenesis, and processes and mechanisms for the rise of the plateau. Geology in China (in Chinese) , 2006, 33(2): 221-238. |
[3] | 高锐, 王海燕, 王成善, 等. 青藏高原东北缘岩石圈缩短变形-深地震反射剖面再处理提供的证据. 地球学报 , 2011, 32(5): 513–520. Gao R, Wang H Y, Wang C S, et al. Lithospheric deformation shortening of the northeastern Tibetan Plateau: Evidence from reprocessing of deep seismic reflection data. Acta Geoscientica Sinica (in Chinese) , 2011, 32(5): 513-520. |
[4] | 曾融生, 腾吉文, 阚荣举, 等. 我国西北地区地壳中的高速夹层. 地球物理学报 , 1965, 14(2): 94–106. Tseng J S, Teng C W, Kan Y C, et al. The high velocity imbedded layer in the crust of northwestern China. Acta Geophysica Sinica (Chinese J.Geophys|) (in Chinese) , 1965, 14(2): 94-106. |
[5] | 李松林, 张先康, 张成科, 等. 玛沁—兰州—靖边地震测深剖面地壳速度结构的初步研究. 地球物理学报 , 2002, 45(2): 210–217. Li S L, Zhang X K, Zhang C K, et al. A Preliminary study on the crustal velocity structure of Maqin-Lanzhou-Jingbian by means of deep seismic sounding profile. Chinese Journal of Geophysics (in Chinese) , 2002, 45(2): 210-217. |
[6] | 高锐, 王海燕, 马永生, 等. 松潘地块若尔盖盆地与西秦岭造山带岩石圈尺度的构造关系-深地震反射剖面探测成果. 地球学报 , 2006, 27(5): 411–418. Gao R, Wang H Y, Ma Y S, et al. Tectonic relationships between the Zoigê basin of the Song-Pan block and the West Qinling orogen at lithosphere scale: results of deep seismic reflection profiling. Acta Geoscientica Sinica (in Chinese) , 2006, 27(5): 411-418. |
[7] | Wang C S, Gao R, Yin A, et al. A mid-crustal strain-transfer model for continental deformation: a new perspective from high-resolution deep seismic-reflection profiling across NE Tibet. Earth Planet.Sci.Lett. , 2011, 306(3-4): 279-288. |
[8] | 陈九辉, 刘启元, 李顺承, 等. 青藏高原东北缘—鄂尔多斯地块地壳上地幔S波速度结构. 地球物理学报 , 2005, 48(2): 333–342. Chen J H, Liu Q Y, Li S C, et al. Crust and upper mantle S-wave velocity structure across Northeastern Tibetan Plateau and Ordos block. Chinese Journal of Geophysics (in Chinese) , 2005, 48(2): 333-342. |
[9] | Wei W B, Unsworth M, Jones A, et al. Detection of widespread fluids in the Tibetan crust by magnetotelluric studies. Science , 2001, 292(5517): 716-719. DOI:10.1126/science.1010580 |
[10] | 赵国泽, 汤吉, 詹艳, 等. 青藏高原东北缘地壳电性结构和地块变形关系的研究. 中国科学(D辑:地球科学) , 2004, 43(10): 908–918. Zhao G Z, Tang J, Zhan Y, et al. The relationship of crustal electrical structure and block deformation in the northeastern margin of Qinghai-Tibet plateau. Science in China (Ser.D) Earth Sciences (in Chinese) , 2004, 43(10): 908-918. |
[11] | 梁桂培, 李渭娟.南北地震带北段深部构造与均衡异常.// 大陆多震层研究.北京: 地震出版社, 1992: 183-190. Liang G P, Li W J.The study of Deep structures, isostatic anomaly and continental seismogenic layer in the northern section of the north-south seismic belt .Beijing: Seismological Press, 1992: 183-190.(请核对本条文献信息) |
[12] | 郭守年, 徐敬文, 孟万辉, 等. 天水地区重力剖面测量结果与地壳结构. 西北地震学报 , 1995(1): 36–42. Guo S N, Xu J W, Meng W H, et al. Gravity profile measurement and crustal structure of Tianshui region. Northwestern Seismological Journal (in Chinese) , 1995(1): 36-42. |
[13] | 李清河, 郭守年, 吕德徽. 鄂尔多斯西缘与西南缘深部结构与构造. 北京: 地震出版社, 1999 : 71 -86. Li Q H, Guo S N, Lv D H. The Deep Structure and Tectonics of the Western and Southwestern Margin of Ordos Block (in Chinese). Beijing: Seismological Press, 1999 : 71 -86. |
[14] | 江为为, 郝天珧, 宋海斌. 鄂尔多斯盆地地质地球物理场特征与地壳结构. 地球物理学进展 , 2000, 15(3): 45–52. Jiang W W, Hao T Y, Song H B. Crustal structure and geological and geophysical features of Ordos basin. Progress in Geophysics (in Chinese) , 2000, 15(3): 45-52. |
[15] | Li S L, Mooney W D, Fan J C. Crustal structure of mainland China from deep seismic sounding data. Tectonophysics\ , 2006, 420(1-2): 239-252. DOI:10.1016/j.tecto.2006.01.026 |
[16] | Zhang J S, Gao R, Zeng L S, et al. Relationship between characteristics of gravity and magnetic anomalies and the earthquakes in the Longmenshan range and adjacent areas. Tectonophysics , 2010, 491(1-4): 218-229. DOI:10.1016/j.tecto.2009.12.004 |
[17] | 吴刚, 林伟凡. 青藏高原东北缘深部构造研究. 高原地震 , 1992, 4(2): 1–9. Wu G, Lin W F. A research on the deep structure in the northeast margin of the Qinghai-Xizang plateau. Plateau Earthquake Research (in Chinese) , 1992, 4(2): 1-9. |
[18] | 李曙光, 孙卫东, 张国伟, 等. 南秦岭勉略构造带黑沟峡变质火山岩的年代学和地球化学——古生代洋盆及其闭合时代的证据. 中国科学(D辑:地球科学) , 1996, 26(3): 223–230. Li S G, Sun W D, Zhang G W, et al. The Geochronology and Geochemistry of South Qinling Mianlue tectonic zone with Heigouxia metavolcanic rocks-evidence of Paleozoic oceanic basin and its close time. Science in China (Ser.D) (in Chinese) , 1996, 26(3): 223-230. |
[19] | Pwalowski R S. Preferential continuation for potential-field anomaly enhancement. Geophysics , 1995, 60(2): 390-398. DOI:10.1190/1.1443775 |
[20] | 郭良辉, 孟小红, 石磊, 等. 优化滤波方法及其在中国大陆布格重力异常数据处理中的应用. 地球物理学报 , 2012, 55(12): 4078–4088. Guo L H, Meng X H, Shi L, et al. Preferential filtering method and its application to Bouguer gravity anomaly of Chinese continent. Chinese Journal of Geophysics (in Chinese) , 2012, 55(12): 4078-4088. |
[21] | Wang W Y, Pan Y, Qiu Z Y. A new edge recognition technology based on the normalized vertical derivative of the total horizontal derivative for potential field data. Applied Geophysics , 2009, 6(3): 226-233. DOI:10.1007/s11770-009-0026-x |
[22] | 王涛, 徐鸣洁, 王良书, 等. 鄂尔多斯及邻区航磁异常特征及其大地构造意义. 地球物理学报 , 2007, 50(1): 163–170. Wang T, Xu M J, Wang L S, et al. Aeromagnetic anomaly analysis of Ordos and adjacent regions and its tectonic implications. Chinese Journal of Geophysics (in Chinese) , 2007, 50(1): 163-170. |
[23] | 周民都, 吕太乙, 张元生, 等. 青藏高原东北缘地质构造背景及地壳结构研究. 地震学报 , 2000, 22(6): 645–653. Zhou M D, Lü T J, Zhang Y S, et al. The geological structure background and the crustal structure in the northeastern margin of Qinghai-Tibetan plateau. Acta Seismologica Sinica (in Chinese) , 2000, 22(6): 645-653. |
[24] | Lin X B, Chen H L, Wyrwoll K H, et al. The uplift history of the Haiyuan-Liupan Shan region northeast of the present Tibetan plateau: integrated constraint from Stratigraphy and Thermochronology. Geology , 2011, 119(4): 372-393. DOI:10.1086/660190 |