1909 年库帕河谷地震:走时拐点怎样画出莫霍面
The 1909 Kupa Valley Earthquake and the Moho
1909 年 10 月 8 日,格林尼治时间约 9 时 59 分,萨格勒布以南库帕河谷发生地震。烟囱倒塌,严重破坏有限。萨格勒布气象台台长安德里亚·莫霍洛维奇开始收集记录。一年前,维歇尔特地震仪刚替换了台里的维琴蒂尼仪器,还有一座里夫勒时钟。震中离萨格勒布约 30 公里。他向欧洲各台索要到时,最终用了约 36 个台的读数。1910 年,克罗地亚文与德文双语长文发表在气象台年报上,题为《1909 年 10 月 8 日的地震》。对象是这批走时,不是后来教科书里光滑的地球分层蛋糕。
先破除一个误解
案例不是“地球有地壳”这句话的插图。当时并不清楚地壳与其下介质之间存在全球可比的速度界面。莫霍洛维奇看到的是:在特定震中距上,P 波和 S 波似乎各有两套到时。直达波与另一支更快的波不能用恒定速度的直线射线同时说圆。他放弃匀速直线,改设速度随深度按指数增加,再试参数,直到理论走时表贴上观测。界面深度落在约 54 公里——那是与资料最合的上介质厚度,不是钻机读出的米数。
地球事件无法被一台仪器完整看见。每台地震仪给出某地的振动轨迹。到时是人从记录上读出的钟面。震中位置还依赖宏观烈度。研究者把这些代理量转换成射线路径和界面深度。每次转换都引入读数误差、台站分布和速度模型结构的不确定性。“观测到莫霍面”与“由走时模型推断出界面”必须分开。
资料是谁生产的,比较对象是什么
萨格勒布台给出近震记录。远一些的欧洲台给出远震相。赫拉克等人 2010 年复查库帕河谷地震时指出:莫霍洛维奇与阿加门农、本多夫、康拉德、盖格、米尔恩等人通信,信件仍在档案里。比较对象不是今天的全球定位震中,而是 1909 年能够对钟的那张台网。台网偏欧洲。射线主要照亮东北方向的地壳。把 54 公里写成全球地壳厚度,是把区域走时表升级成行星常数。大陆莫霍通常约 30–50 公里,大洋约 5–10 公里。莫霍洛维奇测到的是他的射线所穿过的那一块。
他先假设速度 这类随深度增加的形式,射线变成向下凹的曲线。再用试错匹配: 约 5.60 km/s,震源深度一度取 25 公里,$k=3$。界面取使走时拟合最好的上介质厚度,得到约 54 公里。现代重定位把这次地震的震源深度估在大约 14±7 公里,与他为拟合而设的 25 公里不是同一个量。震源深度是事件参数。莫霍深度是界面参数。科普若把两个 50 公里量级的数字合成一个,会把震源和界面叠死。
命名来得更晚。其他人在别的区域看到同类速度跳变之后,界面才被称为莫霍洛维奇不连续面,简称 Moho。名字是荣誉,也是简化。简化之后,学生容易以为 1909 年已经画出全球三维莫霍。画出的是一张区域走时表和一个可迁移的解释:速度在某个深度附近明显加快,快到需要界面,而不是平滑的梯度。
推断怎样从折线走到分层
若地壳与地幔只是温度慢慢升高,波速可以平滑增加,走时曲线不必出现第二支。第二支 P、第二支 S,要求射线在某深度被折射,走一条更快的下路,在足够远的震中距上追上直达波。追上的距离,约束界面有多深、下伏介质有多快。莫霍洛维奇的贡献是把这句射线几何写进可计算的走时方程,而不是只在记录上圈两个波至。
地震 条目把 P、S、面波当作地球的 X 光。本剖析要保留的是 1909 年那一次曝光:近、仪器刚够用、欧洲台网、两套到时。没有这次近且够强的地震,同一个人也许仍能用更早的地中海地震走时做经验表。文献指出他确实用过更早的走时。1909 年事件的特殊处,是近台记录把两套到时差推到必须用界面说话。
地球内部结构 把莫霍写成化学成分分界的习惯画法:地壳岩石与地幔橄榄岩质不同。1909 年的资料直接约束的是弹性波速度,不是手标本。速度跳变与成分跳变,是后来岩石学与实验矿物物理接上的桥。桥可以走,但 1910 年论文并不分析石榴子石。把莫霍直接写成“钻到玄武岩就结束”,会把地球化学后验塞进走时表。
岩石圈与软流圈按力学划分,莫霍按速度(并常被解释成成分)划分。两者不必重合。板块带走的是岩石圈,里面已经包含地壳和地幔盖层。1909 年还没有板块学说。用板块回写莫霍洛维奇,只能写:他的界面后来成为地壳厚度图的骨干,不是他已经在画板块。
反事实与稳健性压力测试
若只有宏观烈度而没有维歇尔特记录,两套到时可能根本读不出来,界面会继续停在猜测。仪器是条件。条件不是命运:别的区域、别的地震,几年内也会被人读到类似折线。1909 年是这条推断的一个清楚实例,不是唯一可能的日期。
若台网全在近处,折射支走不远,界面深度会漂。若台网全在远处,直达地壳波变弱,折线的近端会糊。36 个台的分布决定了 54 这个数的稳度。稳度要用残差和替代速度律来测。单一指数律被后来的分层模型取代。取代之后,莫霍仍在,深度图变细。细,不等于 1910 年的方法作废。方法的核心是:走时的第二支要求速度不连续或近不连续。
若主张“有了全球台网就可以不用界面”,应能指出哪一张平滑速度模型同时拟合近震直达与远震折射。现有资料给不出这张模型。成熟的剖析不是把 54 公里关进定论,而是留下接口:区域深度可修订,界面这类物体在走时里仍必要。
古登堡后来用核–幔边界的影区工作,莱曼用内核。那是同一门 X 光的更深切片。本事件停在地壳底部。切片不能互相替代。把莫霍、古登堡、莱曼收成同一句“他们发现了地球分层”,会把三次不同的走时异常叠成一次英雄叙事。
深一层:速度律、区域厚度与不该用钻机去要求的证明
指数速度律是 1910 年论文里的工作假设,不是今天的地球模型。它的用处是:在还没有高速计算机的时候,给出可微的射线方程,让理论走时能被手算出来并与读数比较。后来的分层模型、走时层析,换了参数化。换参数化之后,界面仍然被需要,因为平滑梯度很难同时解释直达支与折射支的交叉。需要界面,不等于 54 公里这个数被冻结。克罗地亚东北、阿尔卑斯、克拉通,莫霍深度可以差十几公里。把 54 写成全球均值,是把一次拟合的上介质厚度升级成行星常数。
读数本身是人的劳动。模拟记录上的波至,要在噪声和后续振动里挑出第一动。不同分析员会差零点几秒。零点几秒在近震上对应数公里的深度漂。莫霍洛维奇用多台、两套波(P 与 S)互相约束,是为了不让单台误读决定界面。稳健性来自冗余,不是来自一次漂亮的图。今天重做,应报告到时拾取的不确定性,而不是只报告最后那个深度整数。整数是教学方便。方便不能冒充资料精度。
宏观烈度给震中,仪器给到时。两者不一致时,必须声明以谁为准。1909 年的库帕河谷宏观中心与萨格勒布台的几何,让近台记录特别清楚。换一次远在大西洋的地震,萨格勒布台可能根本读不出两套地壳波。对象是这次地震与这张台网的耦合,不是“凡地震都能发现莫霍”。发现依赖几何。几何是方法条件。
钻探不能替代这次推断,也不被这次推断取消。科拉超深钻约 12 公里,仍停在地壳里。大洋钻探可以接近莫霍,因为大洋地壳薄,却仍很少真正穿过。穿过与否,检验的是岩石学解释:界面两侧是否真是玄武质与橄榄岩质。走时检验的是弹性对比。两套检验可以汇合,也可以在局部失配——例如速度梯度带被画成界面,或蛇纹石化把速度压低。失配是研究前沿,不是 1910 年论文的丑闻。1910 年没有义务提供手标本。
灾害解读也要克制。库帕河谷地震对萨格勒布的破坏有限,却是该震源区已知最强事件之一,进入当地危险性评估。危险性评估用的是复发与强度,不是莫霍深度。把“发现地球分层的地震”写成减灾寓言,会把对象从走时换成道德。走时仍然值得单独成篇。减灾有自己的事件剖析,不该借用这次的科学史光环。
教学剖面上的光滑莫霍,是把许多区域界面平均之后的漫画。漫画有用,只要标明它是漫画。真实的莫霍可以起伏、可以成带、可以在造山带下深凹。起伏进入重力与均衡解释,那是后续史。1909 年只要求:至少在他的射线照亮的块体里,速度在某个深度附近明显加快。加快这一句,已经足够改写“地球从地表到地心均匀变热”的十九世纪图景。
S 波不能穿过液体,这一性质在 1909 年的地壳工作里还不是主角。主角是同一套固体里的速度跳变。跳变之后,下伏介质仍传播 S,所以它当时被读成固体,而不是熔融层。熔融的软流圈是更浅、更靠力学衰减的另一回事,不能用莫霍来冒充。冒充会发生在教学简图:把莫霍画成岩石圈底,再把软流圈画成地幔的别名。简图省事,板块会画错。1909 年的论文没有义务预防这种简图。今天的剖析有义务把速度界面和力学界面拆开。拆开之后,走时对象仍然干净:第二支到时要求更快的下路。下路在固体里。固体里的更快,足以定义一个界面,不足以定义板块是否会滑。滑是另一门课。课可以链接,对象不要合并。
1910 年双语文本很长,叙述语气也重,不像后来期刊的短通讯。长度本身说明:他必须同时交代仪器、通信、速度律和拟合。今天的学生若只记住“莫霍 = 50 公里”,就把这四层劳动叠成一个整数。整数可以进填空题。填空题不是剖析。剖析要能回答:哪一张记录上出现了第二支?用了多少台?速度律是什么?54 公里约束的是哪一块地壳?四问都答得上,对象才从名词变回一次推断。推断可以被后人用数字记录重做。重做得不同的深度,并不丢脸。那说明界面深度是模型参数。参数可修订。需要界面这一句,至今仍站得住。站得住的是物体类型,不是某一个整数。
维歇尔特仪器相对维琴蒂尼的进步,是让近震的初动更可读。可读不是“更真”,是带宽与阻尼更适合这类近震。仪器史进入剖析,是为了说明资料生产条件。没有这些条件,第二支可能被当成噪声或仪器故障。故障叙事在科学史里常见:新物相先被怀疑是设备。莫霍洛维奇用多台、P 与 S 成对出现来对抗这种怀疑。成对出现提高了界面解释的优先级。优先级仍是推断,不是照片。照片要等后来的转换波与接收函数。接收函数是另一代资料。资料世代不能回写 1909 年,只能说明同一物体后来被怎样重新照明。
跨域连接
- 地震:P 与 S 的到时是钟面。1909 年的第二支到时,把钟面从“有地震”推进到“有界面”。震级、烈度与烟囱,解释不了这条折线。折线来自射线在更快介质里走了下路,并在足够远的震中距上追上直达波。追上的距离,才是深度的杠杆。
- 地球内部结构:分层图是走时与影区的解,不是钻探日志。莫霍在图上是地壳底。图上的平滑球壳,来自后来把许多区域界面平均到一张教学剖面。平均有用,但 54 公里仍是一次区域拟合,不是行星半径的固定比例。
- 板块构造:漂移的是岩石圈,里面已经包含地壳和地幔盖层。1909 年的界面给地壳厚度,不给板块运动学。运动学要等海底扩张把漂移改写成可检验的洋底年龄之后才有力。用板块回写莫霍洛维奇,只能写他提供了厚度骨干,不能写他已经在画俯冲。
- 莱曼:1936 年用影区里不该出现的波推断固态内核。方法仍是走时异常,对象换到地核,几何完全不同。把莫霍、古登堡与莱曼收成同一句“发现分层”,会把三次不同的异常叠成一次英雄叙事,并让学生无法说出各自到底看见了哪一种波。
- 材料强度:波速依赖弹性模量与密度。界面两侧速度跳变,要求物性跳变。跳变可以是成分,也可以是相变或蛇纹石化。走时先给出弹性对比,实验矿物物理再给出候选。顺序不能倒过来:不能先打开教科书矿物表,再宣称已经“看见”了 1910 年的模拟记录。
参考文献
- Mohorovičić, A. “Potres od 8. X 1909.” Godišnje izvješće Zagrebačkog meteorološkog opservatorija za godinu 1909, 9(4), 1910, 1–56. 德文本 Das Beben vom 8. X. 1909 同期;英译见 Geofizika 9, 1992, 3–55.
- Herak, D. & Herak, M. “The Kupa Valley (Croatia) Earthquake of 8 October 1909—100 Years Later.” Seismological Research Letters, 81(1), 2010. http://www.gfz.hr/PDF/CC20-Herak-SRL%202010.pdf
- Skoko, D. & Mokrović, J. 对莫霍洛维奇 1909 年分析的地球物理回顾,见萨格勒布地球物理研究所公开说明。http://www.gfz.hr/sobe-en/discontinuity.htm
- Tkalčić, H. 英译评述:“Andrija Mohorovičić – an extraordinary scientist who pointed his binoculars down.” 2016.
延伸阅读
- Gutenberg, B. Physics of the Earth’s Interior. Academic Press, 1959. 核–幔边界与走时,须与 1909 年地壳界面分开读。
- Watts, A. B. Isostasy and Flexure of the Lithosphere. Cambridge University Press, 2001. 莫霍深度进入均衡解释的后续史。