We present a 2D thermo-elastic model of a magma body formation in a granitic crust by injection of rhyolitic or basaltic dikes and sills. Elastic analytical solutions enable the computation of rock displacement in response to magma intrusion and evacuation during volcanic eruptions. Phase diagrams for magma and rocks govern melting/crystallization behavior and temperature evolution in 16 million computational cells. Calculated temperature histories are used to predict zircon crystallization/dissolution, their ages, and isotopic ratios within individual batches of magma and rocks. Incremental dike injection naturally generates magma batches of melt that form kilometer-wide clusters appearing in different parts of the domain with diverse shapes, horizontal and vertical interconnectivity. The volume of eruptive melt strongly depends on magma influx rates Q, the width of the injection region W, and eruptions. For example, rhyolitic dike injection with Q = 0.125 m³/s with W = 5 km during 100 ka generates ∼50 km³ of eruptive melt while no significant melt forms if W = 10 km. Injection of basaltic dikes into the granitic crust generates comparable amounts of rhyolitic melt from dike-residual and country rocks. High injection intensities produce magma reservoirs capable of large eruptions, while repetitive eruptions lead to shrinkage of magma bodies. High heat input causes host rock zircons to lose a significant portion of their old cores. Magmatic zircons in the periphery form quickly due to rapid cooling of intruded dikes while in the central part crystals can grow during several hundreds of ka. The model predicts highly heterogeneous O and Hf isotopic ratios recorded in zircons.

中文翻译：

---

岩脉吸积和地壳熔化形成的岩浆室：二维热成分模型，重点是锆石记录的喷发和意义

我们通过注入流纹质或玄武质岩脉和窗台，展示了花岗岩地壳中岩浆体形成的二维热弹性模型。弹性分析解决方案能够计算火山喷发期间岩浆侵入和疏散时响应的岩石位移。岩浆和岩石的相图控制着 1600 万个计算单元的熔化/结晶行为和温度演化。计算出的温度历史用于预测个别批次的岩浆和岩石中的锆石结晶/溶解、年龄和同位素比率。渐增的堤防注入自然会产生岩浆批次的熔体，这些熔体形成千米宽的集群，以不同的形状、水平和垂直的相互联系出现在域的不同部分。W，和喷发。例如， 在 100 ka 期间，Q  = 0.125 m ³ /s 且W = 5 km 的流纹岩脉注入会产生约 50 km ³的喷发熔体，而如果W = 10 公里。将玄武质岩脉注入花岗岩地壳，会从岩脉残岩和围岩中产生相当数量的流纹岩熔体。高注入强度产生能够大规模喷发的岩浆储层，而重复喷发导致岩浆体收缩。高热输入导致主岩锆石失去大部分旧核。由于侵入岩脉的快速冷却，外围的岩浆锆石迅速形成，而在中心部分的晶体可以在数百 ka 期间生长。该模型预测了锆石中记录的高度异质的 O 和 Hf 同位素比率。