利用石英边界的分形分析作为应变速率替代指标以追踪地球应力历史

在普通矿物晶粒中读取地球的过去

山脉会记忆。即便造山的力场已消退，深处的岩石仍然保存着地壳如何被压缩与拉伸的记录。这项研究表明，常见于花岗岩和砂岩等岩石中的普通矿物石英，可以作为记录应力历史的微小档案。通过测量石英晶粒边界变得多么曲折和不规则，作者们提出了一种估算岩石过去变形速率的方法——为活跃山带的隐秘史料开启了一扇新窗。

地下深处的繁忙碰撞带

本研究聚焦于伊朗西南部的查扎尔逆冲带，它属于长期碰撞的扎格罗斯山系，在此两块大陆板块数千万年间相互碰撞。在该地区，古老的火山岩和沉积岩被掩埋、加热并挤压成称为片麻岩的条带状岩石，位于地表几十千米深处。在约420–600 °C和高压条件下，矿物能够缓慢改变形状而不是破裂。由于石英在这些岩石中占比较大并在岩体中连通，其内部织构对碰撞期间地壳流动的记录尤其敏感。

石英晶粒如何响应应力

在高温高压下，石英并非保持刚性。其晶粒会长出新晶体、弯曲并重新排列其内部结构。早期研究显示，不同的变形方式往往在不同温度下出现：在相对较低温度下沿晶界膨胀，在中等条件下形成并旋转亚晶粒，以及在较高温度下晶界的大规模迁移。但更新的研究揭示，这些织构并不仅受温度控制。它们同样强烈响应岩石变形速率、水的含量以及应力的分布。这种复杂性使得将晶粒形状直接换算为精确温度或应力水平变得困难，但也暗示晶粒形状编码了关于整体变形环境的丰富信息。

将不规则晶界变为数字

为了解读这些信息，作者们采用了来自粗糙形状研究的数学工具：分形分析。他们拍摄了八个片麻岩样品中石英的高质量显微镜图像，并手工描绘每个样品至少45个晶粒的外边界。随后在每个轮廓上叠加逐渐变小的方格网，统计与晶界相交的方格数量。将这些计数在对数坐标上对方格尺寸作图，可以揭示边界在不同尺度上的复杂性。该图线的斜率即为“分形维数”，一个介于1到2之间的数值，随着边界越发锯齿状和复杂而增大。利用一条由实验推导出的方程，将分形维数与变形温度和应变速率关联起来，研究组将边界粗糙度转换为岩石在形成这些织构时变形速率的估算值。

这些数字对隐秘变形说明了什么

查扎尔片麻岩中的石英展示了从轻微膨胀到高度锯齿状、呈叶状的完整特征套系，表明岩石经历了若干重叠的变形阶段。分形维数从略高于1.01到约1.21不等，意味着变形强度分布广泛。结合从整体矿物组合和石英织构推断的温度范围，这些数值换算出的大致应变速率约在10⁻¹⁰.⁹到10⁻⁶.⁸每秒之间。这些值高于许多教科书中对大尺度、长期地壳流动的估计，但它们符合这样一种图景：变形并非平稳均匀，而是可以集中在狭窄的带内或以短暂爆发的形式发生，即使在总体缓慢变形的地壳中也能产生局部的高应变速率。

这对理解造山有何意义

通过表明石英晶界的粗糙度可作为半定量的应变速率指标，这项研究为地质学家提供了一个强有力的新证据链。该方法并不宣称能给出温度或应力的完美、唯一解，作者也强调最好将其与传统显微观察和区域地质背景相结合使用。尽管如此，它展示了常见矿物内部微小不规则缝隙能够显示中地壳何时何地发生了更强烈的变形。若将此方法应用于其他山带，可能有助于澄清地壳如何以及何时局部化应变、容纳大陆碰撞，并最终塑造我们在地表所见的地貌。

引用: Abdolzadeh, M., Hosseini, S.R., Rasa, I. et al. Fractal analysis of quartz boundaries as a strain rate proxy for tracing Earth’s stress history. Sci Rep 16, 9759 (2026). https://doi.org/10.1038/s41598-026-40639-8

关键词: 石英变形, 分形分析, 应变速率, 山带, 构造应力