DOI: https://doi.org/10.5194/tc-20-1997-2026
تاريخ النشر: 2026-04-10
المؤلف: Dorothée Vallot وآخرون
الموضوع الرئيسي: دراسات وملاحظات الكريوسفير
نظرة عامة
في هذا القسم، يحقق المؤلفون في ذوبان الجليد عند الواجهة بين المحيط وقاعدة الجرف الجليدي، باستخدام نموذج ذوبان يتحكم فيه القص يأخذ في الاعتبار درجة الحرارة، والملوحة، وسرعة الاحتكاك. على عكس الطريقة التقليدية المستخدمة في نموذج المحيط NEMO، التي تعتمد على معامل سحب ثابت وسرعة محيط متوسطة رأسياً داخل طبقة لوش، تتبنى هذه الدراسة نهجاً لوغاريتمياً. تستخدم هذه الطريقة طول خشونة هيدروغرافي ثابت لتحديد معامل السحب، مأخوذ من أول خلية رطبة، بهدف تقليل الاعتماد على الدقة الرأسية والسماح بمعامل سحب متغير يعكس طبوغرافيا قاعدة الجرف الجليدي على مقياس دون الشبكة. تستند الفكرة وراء هذا النهج إلى الملاحظات التي تشير إلى أن الميزات الطبوغرافية الخشنة، مثل الشقوق وقنوات الذوبان، تشهد ذوباناً أكبر مقارنةً بالمناطق الأكثر تسطحاً.
تظهر النتائج من التجارب التي أجريت في بحر أماندسن عند دقة 1/12° أن النهج الجديد أقل حساسية للدقة الرأسية الخشنة، حيث يظهر فقط فرقاً بنسبة 6% في معدلات الذوبان مقارنةً بفارق 22% تم ملاحظته مع طريقة طبقة لوش. بالإضافة إلى ذلك، تلتقط المعايرة الجديدة معدلات ذوبان أعلى بفعالية في المناطق ذات الطبوغرافيا الخشنة بينما تحافظ على معدلات ذوبان الجرف الجليدي الكلية ضمن النطاقات الملاحظة. تستكشف الدراسة أيضاً تداعيات زيادة طول الخشونة الهيدروغرافي، مشيرةً إلى أنه إذا وصلت الخشونة إلى أعلى القيم الملاحظة، فقد يزيد ذوبان الجرف الجليدي الكلي بنسبة 16%. وهذا يشير إلى حلقة تغذية راجعة إيجابية محتملة حيث يؤدي زيادة الذوبان إلى مزيد من تلف الجليد والخشونة، مما يعزز معدلات الذوبان بشكل أكبر.
مقدمة
تسلط المقدمة الضوء على الدور الحاسم للجرف الجليدي العائم حول القارة القطبية الجنوبية في استقرار الصفائح الجليدية القارية وتنظيم تدفق الجليد. تشير الدراسات الحديثة إلى أن هذه الجرف الجليدي تشهد فقداناً كبيراً للكتلة، مما له تداعيات عميقة على المناخ العالمي وارتفاع مستوى سطح البحر. قد يؤدي ضعف الجرف الجليدي، الذي تفاقم بسبب الذوبان تحت الرف الذي تحركه مياه المحيط الأكثر دفئاً، إلى حالة تفكك تتميز بمزيج من الجليد، كما لوحظ في مناطق مثل لسان الجليد الغربي ثويتس.
يؤكد النص على العلاقة بين خشونة الجليد ومعدلات الذوبان، مشيراً إلى أن مناطق الجليد الأكثر خشونة تظهر ذوباناً معززاً مقارنةً بالمناطق الأكثر سلاسة. تؤثر ميزات مختلفة، مثل الشقوق القاعدية والقنوات، على دوران المحيط وأنماط الذوبان، مما يبرز ضرورة وجود بيانات طبوغرافية عالية الدقة لنمذجة هذه التفاعلات بدقة. غالباً ما تتجاهل نماذج المحيط الحالية التمييز بين الجليد السلس والخشن، وهو ما يمكن معالجته من خلال معايرة عمليات السحب دون الشبكة. تناقش المقدمة أيضاً تعقيدات طبقة الحدود بين الجليد والمحيط، حيث يتم عادةً نمذجة التبادلات الديناميكية الحرارية باستخدام معايرة ثلاثية المعادلات، على الرغم من أن اختيار سمك الطبقة يمكن أن يؤثر بشكل كبير على نتائج هذه النماذج.
النتائج
تقدم النتائج المعروضة في الجدول 2 ملخصاً شاملاً لاختلافات معدلات الذوبان الكلية عبر تجارب مختلفة أجريت على جرف جليدي مختلفة ضمن منطقة شبه الجزيرة القطبية الجنوبية (ASE). تبرز البيانات الاختلافات النسبية في معدلات الذوبان، والتي تعتبر حاسمة لفهم ديناميات استقرار الجرف الجليدي ومساهماتها في ارتفاع مستوى سطح البحر.
ستتناول المناقشات اللاحقة التحليل المقارن لهذه التجارب، مقدمةً رؤى حول العوامل التي تؤثر على معدلات الذوبان الملاحظة وتداعياتها على سيناريوهات المناخ المستقبلية.
المناقشة
في هذا القسم، يناقش المؤلفون التعقيدات والشكوك المرتبطة بنمذجة طبقة الحدود بين الجليد والمحيط، مع التركيز بشكل خاص على معايرة معاملات السحب، $C_d$، وتأثيرها على معدلات الذوبان تحت الجرف الجليدي. يبرزون ضرورة وجود نهج أكثر دقة في المعايرة يأخذ في الاعتبار التغيرات المكانية والزمنية في $C_d$، بالإضافة إلى الآليات الفيزيائية المعنية في واجهة الجليد والمحيط. يقترح المؤلفون استخدام “قانون الجدار” لوصف ملف السرعة بالقرب من قاعدة الجرف الجليدي، الممثل بالمعادلة \( U(z) = u_* \kappa \ln \left( \frac{z}{z_0} \right) \)، حيث \( \kappa \) هو ثابت فون كارمان و \( z_0 \) هو طول الخشونة الهيدروغرافي. يقترحون أن خشونة هيدروغرافية متغيرة مكانياً قد تلتقط بشكل أفضل ديناميات عمليات الذوبان، خاصةً في المناطق ذات الطبوغرافيا المعقدة.
تستخدم الدراسة نموذج المحيط NEMO4.2 لمحاكاة الظروف في بحر أماندسن، مقارنةً بين النهج اللوغاريتمي المقترح وطريقة طبقة لوش التقليدية. تشير النتائج إلى أنه بينما تكون معدلات الذوبان الكلية متشابهة عبر كلا النهجين، توجد اختلافات مكانية كبيرة، خاصةً بالنسبة لجرف جليدي محدد مثل ثويتس، حيث يمكن أن تتفاوت معدلات الذوبان بشكل كبير. يجد المؤلفون أن النهج اللوغاريتمي أقل حساسية للتغيرات في الدقة الرأسية مقارنةً بطريقة طبقة لوش، مما يمكن أن يؤدي إلى عدم دقة في توقعات الذوبان تحت الدقة الخشنة. يخلصون إلى أن دمج طول خشونة هيدروغرافي متغير في معايرات الذوبان قد يعزز دقة النماذج التي تتنبأ بمعدلات الذوبان القاعدية المستقبلية، خاصةً في ضوء الزيادات المتوقعة في تلف الجرف الجليدي وتغير ظروف المحيط.
DOI: https://doi.org/10.5194/tc-20-1997-2026
Publication Date: 2026-04-10
Author(s): Dorothée Vallot et al.
Primary Topic: Cryospheric studies and observations
Overview
In this section, the authors investigate the melting of ice at the interface between the ocean and the ice shelf base, employing a shear-controlled melt parameterization that accounts for temperature, salinity, and friction velocity. Unlike the conventional method used in the NEMO ocean model, which relies on a constant drag coefficient and vertically averaged ocean velocity within the Losch layer, this study adopts a logarithmic approach. This method utilizes a constant hydrographic roughness length to determine the drag coefficient, sampled from the first wet cell, aiming to reduce dependency on vertical resolution and to allow for a variable drag coefficient that reflects subgrid-scale ice shelf basal topography. The rationale for this approach stems from observations indicating that rough topographic features, such as crevasses and melt channels, experience greater melting compared to flatter areas.
The findings from experiments conducted in the Amundsen Sea at a resolution of 1/12° reveal that the new approach is less sensitive to coarser vertical resolutions, showing only a 6% difference in melt rates compared to a 22% difference observed with the Losch layer method. Additionally, the new parameterization effectively captures higher melt rates in regions with rough topography while maintaining total ice shelf melt rates within observed ranges. The study also explores the implications of increased hydrographic roughness length, suggesting that if roughness were to reach the highest observed values, overall ice shelf melting could increase by 16%. This indicates a potential positive feedback loop where increased melting leads to greater ice damage and roughness, further enhancing melt rates.
Introduction
The introduction highlights the critical role of floating ice shelves around Antarctica in stabilizing the continental ice sheet and regulating ice flow. Recent studies indicate that these ice shelves are experiencing significant mass loss, which has profound implications for global climate and sea level rise. The weakening of ice shelves, exacerbated by sub-shelf melting driven by warmer ocean waters, could lead to a disintegration state characterized by ice mélange, as observed in regions like the Thwaites Western Ice Tongue.
The text emphasizes the relationship between ice roughness and melting rates, noting that rougher ice regions exhibit enhanced melting compared to smoother areas. Various features, such as basal crevasses and channels, influence ocean circulation and melt patterns, underscoring the necessity for high-resolution topographical data to accurately model these interactions. Current ocean models often overlook the distinction between smooth and rough ice, which can be addressed through the parameterization of subgrid drag processes. The introduction also discusses the complexities of the ice-ocean boundary layer, where thermodynamic exchanges are typically modeled using a three-equation parameterization, although the choice of layer thickness can significantly affect the outcomes of these models.
Results
The results presented in Table 2 provide a comprehensive summary of the total melt rate differences across various experiments conducted on different ice shelves within the Antarctic Peninsula region (ASE). The data highlights the relative variations in melt rates, which are critical for understanding the dynamics of ice shelf stability and their contributions to sea level rise.
Subsequent discussions will elaborate on the comparative analysis of these experiments, offering insights into the factors influencing the observed melt rates and their implications for future climate scenarios.
Discussion
In this section, the authors discuss the complexities and uncertainties associated with modeling the ice-ocean boundary layer, particularly focusing on the parameterization of drag coefficients, $C_d$, and their influence on melt rates beneath ice shelves. They highlight the necessity for a more nuanced approach to parameterization that accounts for spatial and temporal variability in $C_d$, as well as the physical mechanisms at play in the ice-ocean interface. The authors propose using the “law of the wall” to describe the velocity profile near the ice shelf base, represented by the equation \( U(z) = u_* \kappa \ln \left( \frac{z}{z_0} \right) \), where \( \kappa \) is von Kármán’s constant and \( z_0 \) is the hydrographic roughness length. They suggest that a spatially varying hydrographic roughness could better capture the dynamics of melt processes, especially in regions with complex topography.
The study employs the NEMO4.2 ocean model to simulate conditions in the Amundsen Sea, comparing the proposed logarithmic approach to the traditional Losch layer method. Results indicate that while total melt rates are similar across both approaches, significant spatial differences exist, particularly for specific ice shelves like Thwaites, where melt rates can vary dramatically. The authors find that the logarithmic approach is less sensitive to vertical resolution changes compared to the Losch layer method, which can lead to inaccuracies in melt predictions under coarse resolutions. They conclude that incorporating a variable hydrographic roughness length into melt parameterizations could enhance the accuracy of models predicting future basal melt rates, especially in light of anticipated increases in ice shelf damage and changing ocean conditions.