DOI: https://doi.org/10.3934/geosci.2026015
تاريخ النشر: 2026-01-01
المؤلف: Alexander G. Simakin وآخرون
الموضوع الرئيسي: التحليل الجيولوجي والجيوكيميائي
نظرة عامة
تبحث هذه الدراسة في التفاعل بين الصهارة الريوليتية المائية والسوائل الكربونية، مع التركيز بشكل خاص على دور محتويات H2O-CO2 في الشوائب المحتجزة في الكوارتز (MIs) كمؤشرات على مثل هذه التفاعلات. تكشف الدراسة أن السوائل الغنية بـ CO2 يمكن أن تغسل الصهارة الريوليتية، مما يؤدي إلى تكوين تركيبات ذائبة تظهر أنماطًا مميزة على مخططات H2O-CO2. تشير النمذجة العددية على المقياس الصغير إلى أن ذوبانًا مجففًا بشكل موحد يمكن أن يتطور بسرعة، مع آثار كبيرة على فهم الانفجارات الانفجارية، خاصة في سياق ثوران صهارة هكلبيري ريدج (HRT) في كالديرا يلوستون.
تشمل النتائج الرئيسية تحديد حلقتين متميزتين من الغسل المرتبطتين بالثوران المبكر لـ HRT، حيث تتميز بارتباطات متناقضة بين تركيزات الليثيوم (Li) و CO2 في MIs. تضمنت الحلقة الأولى سائلًا غنيًا بـ CO2 و Li، بينما كانت الحلقة الثانية مميزة بسائل فقير في Li، مما أدى إلى سلوكيات ثورانية مختلفة. بالإضافة إلى ذلك، تقوم الدراسة بت quantifying آثار إزالة الغازات داخل الغرفة وفك الضغط على محتوى المتطايرات في MIs، مما يشير إلى أن ديناميات انتشار الماء و CO2 تلعب دورًا حاسمًا في تطور تركيبات الصهارة أثناء الثوران. تؤكد النتائج على تعقيد العمليات الماجماتية وأهمية تفاعلات السوائل في دفع النشاط البركاني الانفجاري.
مقدمة
تناقش مقدمة ورقة البحث التأثير الكبير لتكوين الفقاعات في الصهارة على خصائصها الفيزيائية وديناميات الثوران. تبرز أن وجود المكونات المتطايرة، وخاصة CO2، يمكن أن يعزز حجم ومدة الثورانات من خلال زيادة كثافة الصهارة وقابليتها للضغط. تشير الدراسة إلى نتائج تجريبية تشير إلى أن ارتفاع فقاعات الغاز عبر طبقات الصهارة يمكن أن يحفز عدم استقرار رايلي-تايلور، مما يؤدي إلى أعمدة حمل تؤثر بشكل أكبر على ضغط غرفة الصهارة وقد تحفز الثورانات.
يؤكد المؤلفون على أهمية نمذجة ضغط الصهارة بدقة، مع الأخذ في الاعتبار عوامل مثل قابلية ضغط الصهارة وقابلية ذوبان المتطايرات تحت ضغوط متغيرة. كما يستكشفون التفاعل بين CO2 و الذوبان المائي السيليكي، مشيرين إلى أن الجفاف السريع يحدث مع الحفاظ على محتوى CO2 أثناء كل من فك الضغط وتدفق CO2. تهدف الورقة إلى تحسين النماذج العددية لتفاعلات CO2-الذوبان، مع التركيز بشكل خاص على فترات زمنية قصيرة وآثار حجم الفقاعات ونسبة حجم السائل على الجفاف الأولي. يتم وضع النتائج في سياق أوسع لأنظمة بركانية، مع إشارات محددة إلى ثورانات يلوستون وتوبا توف، مما يبرز أهمية هذه العمليات في فهم النشاط البركاني.
طرق
في هذا القسم، يستكشف المؤلفون التحكم الانتشاري لـ CO2 و H2O أثناء إزالة الغازات من الصهارة، مشيرين إلى أن CO2 ينتشر بشكل أبطأ من الماء. يؤدي هذا الانتشار التفاضلي إلى فقدان الماء بشكل متسلسل قبل CO2 أثناء إزالة الغازات بفعل فك الضغط، مما ينتج عنه مسارات مميزة على مخططات محتوى المتطايرات. تؤكد الملاحظات التجريبية أن كل من إزالة الغازات غير المتوازنة وغسل CO2 تنتج توزيعًا شبه عمودي لتركيبات الذوبان في مخطط H2O-CO2، مع آثار على جفاف الذوبان وزيادة درجات حرارة السائل للمعادن الماجماتية. يشير المؤلفون إلى أن الشوائب المحتجزة (MIs) في الكوارتز يمكن أن تسجل محتويات المتطايرات، وأثناء الغسل، يمكن أن تظهر مجموعة عمودية من تركيبات MI، تعكس درجات متفاوتة من تشبع CO2.
تستكشف الدراسة أيضًا فقدان الماء من MIs في الكوارتز، مشددة على أنه بينما يعتبر الكوارتز مضيفًا أكثر استقرارًا لـ MIs، يمكن أن ينتشر الماء خارجًا على مدى أيام إلى أسابيع عند درجات حرارة مرتفعة. تظهر التجارب انخفاضًا كبيرًا في محتوى الماء في البداية، يُعزى إلى انتشار الهيدروجين عبر الكوارتز، الذي يحركه تحلل الماء وأكسدة الحديد الثنائي في الذوبان. يستخدم المؤلفون نموذج قشرة كروية لتقدير معدل فقدان الماء، مستخلصين حلاً أسّيًا يتماشى جيدًا مع البيانات التجريبية. تشير نتائجهم إلى أن معاملات انتشار الماء في الكوارتز تتأثر بكثافة تشوه البلورة، مع ملاحظة معدلات انتشار أقل في البلورات الطبيعية المثالية مقارنة بتلك ذات كثافات التشوه الأعلى. بشكل عام، توفر الدراسة رؤى قيمة حول آليات نقل المتطايرات في الأنظمة الماجماتية والعوامل التي تؤثر على انتشار الماء في الكوارتز.
نتائج
تشير نتائج النمذجة العددية إلى وجود علاقة بين الوقت اللازم لتحقيق توزيع موحد للماء ($\tau_u$) ونسبة حجم الفقاعات ونصف قطر الفقاعة الأولي ($R_0$)، كما هو موضح في المعادلة:
\[
\log(\tau_u) = -0.231 + 1.782 \log\left(\frac{R_1}{R_0}\right) = -0.231 + 1.782 \frac{\log\left(\frac{1}{\epsilon_{fl}}\right)}{3}
\]
مع انحراف معياري قدره 0.056. تُظهر التحليلات أنه بالنسبة للفقاعات الصغيرة (مع قيم $R_0$ بين 20-50 ميكرومتر)، يكون $\tau_u$ عادةً بضع دقائق، مما يتماشى مع النتائج التجريبية، مثل تلك التي قدمها يوشيمورا، الذي لاحظ تدرجات كبيرة من CO2 وتوزيع موحد لـ H2O في الذوبانات الريوليتية خلال تجارب استمرت ساعة واحدة.
عند توسيع النموذج ليشمل فقاعات أكبر ونسب حجم أقل، كما هو موجود في غرف الصهارة، تتراوح تقديرات $\tau_u$ من حوالي ساعة إلى أسبوع. بالإضافة إلى ذلك، تقدم الدراسة مسارات لعمليات الغسل غير المتوازنة وإزالة الغازات داخل نظام الذوبان-CO2-H2O، موضحة كيف يؤدي فك الضغط إلى مسارات مميزة لتفكك CO2 و H2O. يتم وضع هذه النتائج في سياق بيانات الشوائب المحتجزة من ثوران المناخ LCT في يلوستون، مما يوفر إطارًا مقارنًا لفهم ديناميات سلوك الغاز في الأنظمة الماجماتية.
مناقشة
في هذا القسم، يناقش المؤلفون الديناميات التي تتحكم فيها الانتشار لـ H2O و CO2 في الصهارة الريوليتية أثناء غسل السوائل الكربونية. يبرزون قيود نموذجهم السابق، الذي افترض معاملات انتشار ثابتة لـ H2O و CO2، ويقدمون نهجًا عدديًا أكثر تعقيدًا يأخذ في الاعتبار معاملات انتشار متغيرة تتأثر بالتركيز. يستخدم النموذج طريقة حجم التحكم لمحاكاة التفاعل بين فقاعات CO2 وذوبان مائي، كاشفًا أن معدلات الانتشار لكل من المكونين تتأثر بلزوجة الذوبان ومحتوى الماء. ومن الجدير بالذكر أن معامل انتشار H2O أعلى بكثير من CO2، خاصة عند تركيزات الماء المنخفضة، مما يؤدي إلى انفصال المتطايرات أثناء صعود الصهارة.
كما يقدم المؤلفون نتائج عددية تشير إلى أن الوقت المطلوب لوصول تركيز الماء في الذوبان إلى التوحيد، المشار إليه بـ $\tau_u$، هو أمر حاسم لفهم ديناميات تدفق الصهارة أثناء الثورانات. يجدون أن $\tau_u$ يقارن بوقت تدفق الصهارة في القنوات، مما يشير إلى أن آثار الانتشار تصبح مهمة عندما تقترب سرعة التدفق من 25-50 م/ث. علاوة على ذلك، تؤكد الدراسة على أهمية فهم سلوك الانتشار للمتطايرات في الشوائب المحتجزة والمجاري، حيث يمكن أن توفر هذه رؤى حول الظروف التي تؤدي إلى الثورانات البركانية الانفجارية. تؤكد النتائج على تعقيدات نقل المتطايرات في الأنظمة الماجماتية وضرورة التفسير الدقيق لبيانات الشوائب المحتجزة في سياق ديناميات الثوران.
DOI: https://doi.org/10.3934/geosci.2026015
Publication Date: 2026-01-01
Author(s): Alexander G. Simakin et al.
Primary Topic: Geological and Geochemical Analysis
Overview
This research investigates the interaction between hydrous rhyolitic magma and carbonic fluids, particularly focusing on the role of H2O-CO2 contents in quartz-hosted melt inclusions (MIs) as indicators of such interactions. The study reveals that CO2-rich fluids can flush rhyolitic magma, leading to the formation of melt compositions that display distinct patterns on H2O-CO2 diagrams. Numerical modeling at the micro-scale indicates that a uniformly dehydrated melt can develop rapidly, with significant implications for the understanding of explosive eruptions, particularly in the context of the Yellowstone caldera’s Huckleberry Ridge Tuff (HRT) eruptions.
Key findings include the identification of two distinct flushing episodes associated with the early HRT eruption, characterized by contrasting correlations between lithium (Li) and CO2 concentrations in MIs. The first episode involved a CO2 and Li-rich fluid, while the second was marked by a Li-poor fluid, leading to different eruptive behaviors. Additionally, the study quantifies the effects of intra-chamber degassing and decompression on the volatile content of MIs, suggesting that the dynamics of water and CO2 diffusion play a critical role in the evolution of magma compositions during eruptions. The results underscore the complexity of magmatic processes and the importance of fluid interactions in driving explosive volcanic activity.
Introduction
The introduction of the research paper discusses the significant impact of magma vesiculation on its physical properties and eruption dynamics. It highlights that the presence of volatile components, particularly CO2, can enhance the volume and duration of eruptions by increasing the density and compressibility of the magma. The study references experimental findings indicating that the rise of gas bubbles through magma layers can induce Rayleigh-Taylor instability, leading to convective plumes that further influence magma chamber pressure and potentially trigger eruptions.
The authors emphasize the importance of accurately modeling magma pressure, taking into account factors such as magma compressibility and the solubility of volatiles under varying pressures. They also explore the interaction between CO2 and hydrous silicic melts, noting that rapid dehydration occurs while maintaining CO2 content during both decompression and CO2 fluxing. The paper aims to refine numerical models of CO2-melt interactions, particularly focusing on short time scales and the effects of bubble size and fluid volume fraction on initial dehydration. The findings are contextualized within the broader framework of volcanic systems, with specific references to the Yellowstone and Toba Tuff eruptions, underscoring the relevance of these processes in understanding volcanic activity.
Methods
In this section, the authors investigate the diffusional control of CO2 and H2O during the degassing of magma, highlighting that CO2 diffuses more slowly than water. This differential diffusion leads to the sequential loss of water before CO2 during decompression degassing, resulting in distinct trajectories on volatile content diagrams. Experimental observations confirm that both nonequilibrium degassing and CO2 flushing produce a subvertical distribution of melt compositions in the H2O-CO2 diagram, with implications for melt dehydration and increased liquidus temperatures of magmatic minerals. The authors note that melt inclusions (MIs) in quartz can record volatile contents, and during flushing, a vertical array of MI compositions can emerge, reflecting varying degrees of CO2 saturation.
The study further explores the water loss from MIs in quartz, emphasizing that while quartz is a more stable host for MIs, water can still diffuse out over days to weeks at elevated temperatures. Experiments demonstrate a significant initial drop in water content, attributed to hydrogen diffusion through quartz, which is driven by water decomposition and the oxidation of ferrous iron in the melt. The authors employ a spherical shell model to estimate the rate of water loss, deriving an exponential solution that aligns well with experimental data. Their findings suggest that the diffusion coefficients of water in quartz are influenced by the crystal’s dislocation density, with lower diffusion rates observed in ideal natural crystals compared to those with higher dislocation densities. Overall, the research provides valuable insights into the mechanisms of volatile transport in magmatic systems and the factors affecting water diffusion in quartz.
Results
The results of the numerical modeling indicate a relationship between the time to achieve uniform water distribution ($\tau_u$) and the bubble volume fraction and initial bubble radius ($R_0$), as described by the equation:
\[
\log(\tau_u) = -0.231 + 1.782 \log\left(\frac{R_1}{R_0}\right) = -0.231 + 1.782 \frac{\log\left(\frac{1}{\epsilon_{fl}}\right)}{3}
\]
with a standard deviation of 0.056. The analysis shows that for small bubbles (with $R_0$ values between 20-50 µm), $\tau_u$ is typically a few minutes, aligning with experimental findings, such as those by Yoshimura, who observed significant CO2 gradients and uniform H2O distribution in rhyolitic melts over one-hour experiments.
When scaling the model to larger bubbles and lower volumetric fractions, as found in magma chambers, $\tau_u$ estimates range from approximately one hour to one week. Additionally, the study presents trajectories of disequilibrium flushing and degassing processes within the melt-CO2-H2O system, illustrating how decompression degassing leads to distinct trajectories for CO2 and H2O exsolution. These findings are contextualized with melt inclusion data from the Yellowstone LCT climatic eruption, providing a comparative framework for understanding the dynamics of gas behavior in magmatic systems.
Discussion
In this section, the authors discuss the diffusion-controlled dynamics of H2O and CO2 in rhyolitic magma during carbonic fluid flushing. They highlight the limitations of their previous model, which assumed constant diffusion coefficients for H2O and CO2, and introduce a more sophisticated numerical approach that accounts for variable diffusion coefficients influenced by concentration. The model employs a control volume method to simulate the interaction between CO2 bubbles and a hydrous melt, revealing that the diffusion rates of both components are affected by the melt’s viscosity and water content. Notably, the diffusion coefficient of H2O is significantly higher than that of CO2, particularly at low water concentrations, which leads to a decoupling of the two volatiles during magma ascent.
The authors also present numerical results indicating that the time required for the water concentration in the melt to reach uniformity, denoted as $\tau_u$, is critical for understanding the dynamics of magma flow during eruptions. They find that $\tau_u$ is comparable to the flow time of magma in conduits, suggesting that diffusion effects become significant when the flow velocity approaches 25-50 m/s. Furthermore, the study emphasizes the importance of understanding the diffusion behavior of volatiles in melt inclusions and embayments, as these can provide insights into the conditions leading to explosive volcanic eruptions. The findings underscore the complexities of volatile transport in magmatic systems and the necessity for careful interpretation of melt inclusion data in the context of eruption dynamics.