العمليات المتنافسة تحدد التأثير طويل الأمد لتحديد معلمات الاحتكاك القاعدي لفقدان الكتلة في القارة القطبية الجنوبية
Competing processes determine the long-term impact of basal friction parameterizations for Antarctic mass loss

المجلة: ˜The œcryosphere، المجلد: 20، العدد: 2
DOI: https://doi.org/10.5194/tc-20-1217-2026
تاريخ النشر: 2026-02-16
المؤلف: Tim van den Akker وآخرون
الموضوع الرئيسي: دراسات وملاحظات الكريوسفير

نظرة عامة

تبحث الدراسة في تأثير قوانين الاحتكاك القاعدي على توقعات مستوى سطح البحر باستخدام نموذج الجليد المجتمعي (CISM)، مشددة على أن هذه الحساسية تعتمد على التكوين الهندسي والظروف الأولية لنموذج الجليد. تدرس الدراسة بشكل خاص خليج أموندسن، حيث يتم محاكاة انهيار الأنهار الجليدية ثويتس وباين آيلاند تحت تأثير المناخ الحالي. تؤدي حالتان مختلفتان من التهيئة، تختلفان في الضغط الفعال للجليد على الصخور الأساسية، إلى مسارات تطورية مختلفة لصفائح الجليد في القارة القطبية الجنوبية. في السيناريو الأول، يؤدي انهيار نهر ثويتس الجليدي إلى تعزيز الدعم الذي يخفف من تأثير تغييرات قوانين الاحتكاك القاعدي، مما ينتج عنه عدم حساسية النموذج لهذه المعلمات. على العكس من ذلك، يكشف السيناريو الثاني، حيث ينهار نهر باين آيلاند الجليدي أولاً، عن حساسية ملحوظة للاحتكاك القاعدي بسبب نقص الدعم القوي.

تؤكد النتائج على أن ترتيب انهيار الأنهار الجليدية، المتأثر بخيارات التهيئة، يؤثر بشكل كبير على حساسية النموذج لتحديد معلمات الاحتكاك القاعدي. تشير كلا المحاكاة إلى أن القوة الحرارية الحالية للمحيط كافية لتحفيز انهيار كلا النهرين الجليديين، ومع ذلك فإن الاختلافات الطفيفة في التهيئة تقدم عدم يقين يبلغ حوالي 500 عام في توقيت الانهيار. تقترح الدراسة أن الاستكشاف الإضافي باستخدام نماذج مختلفة لصفائح الجليد واختبارات موحدة يمكن أن يعزز فهم التباينات في الحساسية لقوانين الاحتكاك القاعدي التي لوحظت في محاكاة مختلفة. بالإضافة إلى ذلك، يمكن أن تساعد البيانات المراقبة المحسنة حول الانزلاق مقابل التدفق المهيمن على التشوه في تحسين عمليات تهيئة النموذج.

مقدمة

تناقش مقدمة ورقة البحث المساهمات المتوقعة لمنطقة القارة القطبية الجنوبية في ارتفاع مستوى سطح البحر العالمي، مقدرة الزيادات بين 0.03-0.27 م تحت سيناريو SSP1.9 و0.03-0.34 م تحت سيناريو SSP8.5 بحلول عام 2100. بعد هذه النقطة، تتصاعد عدم اليقين بسبب العمليات الديناميكية، مما قد يؤدي إلى ارتفاعات متعددة الأمتار بحلول عام 2300، مدفوعة بشكل خاص بإزالة الجليد الذاتية الاستدامة لصفائح الجليد في غرب القارة القطبية الجنوبية (WAIS). يتم تسليط الضوء على خليج أموندسن (ASE) كمنطقة حاسمة للتغيرات الديناميكية في الكتلة، حيث تم تحديد نهر باين آيلاند الجليدي (PIG) ونهر ثويتس الجليدي (TG) كمصادر هامة لعدم اليقين بشأن سلوك صفائح الجليد المستقبلية، خاصة بسبب تعرضها لعدم استقرار صفائح الجليد البحرية (MISI).

تؤكد الورقة على دور تحديدات الاحتكاك القاعدي في نمذجة ديناميات صفائح الجليد، والتي يمكن أن تؤثر بشكل كبير على توقعات فقدان الجليد وارتفاع مستوى سطح البحر. يتم مناقشة نوعين رئيسيين من قوانين الاحتكاك القاعدي: احتكاك قانون القوة، الذي يربط الاحتكاك بسرعة القاعدة، واحتكاك كولوم، الذي يصبح مستقلاً عن السرعة للجليد سريع التدفق. يهدف المؤلفون إلى دراسة حساسية تراجع خط التثبيت وفقدان الكتلة الجليدية في WAIS لهذه القوانين باستخدام نموذج الجليد المجتمعي (CISM)، موسعين المحاكاة على مدى 2000 عام. تتضمن الدراسة معدلات التغير في الكتلة الملاحظة للتهيئة وتستكشف كيف يمكن أن تؤثر قوانين الاحتكاك المختلفة وتأثيرات الدعم على نتائج النموذج، مشددة على أهمية التطور الهندسي في تحديد حساسية ديناميات صفائح الجليد لاختيارات قوانين الاحتكاك.

طرق

يستعرض قسم “الطرق” الأساليب التجريبية والتحليلية المستخدمة في الدراسة. يوضح تصميم التجارب، بما في ذلك اختيار الموضوعات، والمواد المستخدمة، والإجراءات المحددة المتبعة لضمان إمكانية التكرار. تم إجراء تحليلات إحصائية لتقييم أهمية النتائج، باستخدام تقنيات مثل ANOVA وتحليل الانحدار لتفسير البيانات بدقة.

بالإضافة إلى ذلك، يصف القسم النماذج الرياضية المطبقة لتحليل العلاقات بين المتغيرات، مع ضمان تلبية افتراضات هذه النماذج. تم اختبار المنهجيات بدقة من حيث الموثوقية والصلاحية، مع تنفيذ ضوابط مناسبة لتخفيف التحيزات المحتملة. بشكل عام، توفر الطرق المستخدمة إطارًا قويًا لفهم النتائج المقدمة في الدراسة.

نتائج

في هذا القسم، يقوم المؤلفون بتحليل نموذج صفائح الجليد الحالي من خلال طريقتين مختلفتين للتهيئة، كما هو موضح في القسم 2.4. يؤكدون على الفروق الحاسمة بين هذه الطرق وتأثيرات خيارات التهيئة على نتائج النموذج. بعد ذلك، يقدم المؤلفون محاكاة مستقبلية غير مجبرة، تفحص كيف تؤثر التغيرات في الاحتكاك القاعدي على استجابة صفائح الجليد. يوضحون أيضًا العلاقة بين هذه الاستجابات والهندسة الأساسية لصفائح الجليد، بالإضافة إلى دور الأرفف الجليدية المدعومة في تعديل هذه الديناميات.

مناقشة

في هذا القسم، يقوم المؤلفون بتقييم أربعة قوانين للاحتكاك القاعدي لفهم تأثيرها على ديناميات صفائح الجليد، مع التركيز بشكل خاص على قانون زويت-إيفرسون (ZI) ومقارنته بثلاثة نماذج احتكاك أخرى: قانون القوة، قانون شوف، وقانون اللدونة الزائفة. يتضمن قانون ZI معلمة ضبط متغيرة مكانيًا، \(C_c\)، التي تعدل النموذج ليتناسب مع سمك الجليد الملاحظ. تسلط الدراسة الضوء على أن قانون القوة وقانون اللدونة الزائفة يظهران زيادة في الاحتكاك القاعدي مع سرعات أعلى، مما يقدم آلية تغذية راجعة سلبية تخفف الاضطرابات بشكل أكثر فعالية من قوانين ZI وشوف، التي تقترب من قيمة احتكاك ثابتة عند السرعات العالية. هذه التمييزات حاسمة لمحاكاة ديناميات الجليد، خاصة بالقرب من خطوط التثبيت حيث يؤثر الضغط الفعال على الاحتكاك القاعدي.

كما يقوم المؤلفون بتحديد قدرة الدعم للأرفف الجليدية باستخدام طريقتين: نهج توازن الإجهاد وتجارب إزالة الأرفف. يشير رقم الدعم، المشتق من مقارنات الإجهاد، إلى درجة الدعم المقدم من الرف الجليدي للجليد المثبت. تقيم الطريقة الثانية تسارع الجليد المثبت بعد إزالة الرف، مما يوفر لمحة زمنية عن تأثيرات الدعم. يتم التحقق من كلا الطريقتين مقابل الحالات النظرية والحالة الحالية لصفائح الجليد في القارة القطبية الجنوبية. يتم استخدام نموذج الجليد المجتمعي (CISM) للمحاكاة، مع استخدام تقريب من الدرجة الأعلى لتوازن الزخم ودمج ضبط المعلمات للاحتكاك القاعدي ومعدلات الذوبان بناءً على البيانات الملاحظة. تم تصميم التهيئات (P05 وP1) لضمان الاستقرار وتقليل الانجراف في مخرجات النموذج، مما يمهد الطريق لتجارب الاستمرار التي تستكشف حساسية تطور صفائح الجليد لقوانين الاحتكاك القاعدي المختلفة.

Journal: ˜The œcryosphere, Volume: 20, Issue: 2
DOI: https://doi.org/10.5194/tc-20-1217-2026
Publication Date: 2026-02-16
Author(s): Tim van den Akker et al.
Primary Topic: Cryospheric studies and observations

Overview

The research investigates the influence of basal friction laws on sea-level projections using the Community Ice Sheet Model (CISM), highlighting that this sensitivity is contingent upon the geometric configuration and initial conditions of the ice sheet model. The study specifically examines the Amundsen Sea Embayment, where the collapse of Thwaites and Pine Island glaciers is simulated under present-day climate forcing. Two distinct initializations, differing in the effective pressure of ice on bedrock, yield different evolutionary paths for the Antarctic Ice Sheet. In the first scenario, the collapse of Thwaites Glacier leads to enhanced buttressing that mitigates the impact of varying basal friction laws, resulting in a model insensitivity to these parameters. Conversely, the second scenario, where Pine Island Glacier collapses first, reveals a pronounced sensitivity to basal friction due to the lack of strong buttressing.

The findings underscore that the order of glacier collapse, influenced by initialization choices, significantly affects the model’s sensitivity to basal friction parameterization. Both simulations indicate that current ocean thermal forcing is adequate to trigger the collapse of both glaciers, yet minor differences in initialization introduce an uncertainty of approximately 500 years in collapse timing. The study suggests that further exploration with various ice sheet models and standardized tests could enhance understanding of the discrepancies in sensitivity to basal friction laws observed in different simulations. Additionally, improved observational data on sliding versus deformation-dominated flow could refine model initialization processes.

Introduction

The introduction of the research paper discusses the projected contributions of the Antarctic region to global mean sea level rise, estimating increases between 0.03-0.27 m under SSP1.9 and 0.03-0.34 m under SSP8.5 by the year 2100. Beyond this point, uncertainties escalate due to dynamic processes, potentially leading to multi-meter rises by 2300, particularly driven by self-sustaining deglaciation of the West Antarctic Ice Sheet (WAIS). The Amundsen Sea Embayment (ASE) is highlighted as a critical area for dynamic mass changes, with the Pine Island Glacier (PIG) and Thwaites Glacier (TG) identified as significant sources of uncertainty regarding future ice sheet behavior, particularly due to their susceptibility to Marine Ice Sheet Instability (MISI).

The paper emphasizes the role of basal friction parameterizations in modeling ice sheet dynamics, which can significantly influence predictions of ice loss and sea level rise. Two primary types of basal friction laws are discussed: power law friction, which relates friction to basal velocity, and Coulomb friction, which becomes independent of velocity for fast-flowing ice. The authors aim to investigate the sensitivity of grounding line retreat and ice mass loss in the WAIS to these friction laws using the Community Ice Sheet Model (CISM), extending simulations over 2000 years. The study incorporates observed mass change rates for initialization and explores how different friction laws and buttressing effects can impact model outcomes, emphasizing the importance of geometric evolution in determining the sensitivity of ice sheet dynamics to friction law choices.

Methods

The “Methods” section outlines the experimental and analytical approaches employed in the study. It details the design of the experiments, including the selection of subjects, materials used, and the specific procedures followed to ensure reproducibility. Statistical analyses were conducted to evaluate the significance of the results, employing techniques such as ANOVA and regression analysis to interpret the data accurately.

Additionally, the section describes the mathematical models applied to analyze the relationships between variables, ensuring that the assumptions of these models were met. The methodologies were rigorously tested for reliability and validity, with appropriate controls implemented to mitigate potential biases. Overall, the methods employed provide a robust framework for understanding the findings presented in the study.

Results

In this section, the authors analyze the present-day ice sheet modeled through two distinct initialization methods, as outlined in Section 2.4. They emphasize the critical differences between these methods and the implications of the initialization choices on the model’s outcomes. Following this, the authors present unforced future simulations, examining how variations in basal friction influence the ice sheet’s response. They further elucidate the relationship between these responses and the underlying ice sheet geometry, as well as the role of buttressed ice shelves in modulating these dynamics.

Discussion

In this section, the authors evaluate four basal friction laws to understand their impact on ice sheet dynamics, particularly focusing on the Zoet-Iverson (ZI) law and its comparison with three other friction models: the power law, Schoof law, and pseudoplastic law. The ZI law incorporates a spatially varying tuning parameter, \(C_c\), which adjusts the model to match observed ice thickness. The study highlights that the power law and pseudoplastic law exhibit increasing basal friction with higher velocities, introducing a negative feedback mechanism that dampens perturbations more effectively than the ZI and Schoof laws, which asymptote to a constant friction value at high velocities. This distinction is crucial for simulating ice dynamics, especially near grounding lines where effective pressure influences basal friction.

The authors also quantify the buttressing capacity of ice shelves using two methods: a stress balance approach and shelf-removal experiments. The buttressing number, derived from stress comparisons, indicates the degree of support provided by the ice shelf to the grounded ice. The second method assesses the acceleration of grounded ice following the removal of the shelf, providing a temporal snapshot of buttressing effects. Both methods are validated against theoretical cases and the current state of the Antarctic Ice Sheet. The Community Ice Sheet Model (CISM) is employed for simulations, utilizing a higher-order approximation to momentum balance and incorporating parameter tuning for basal friction and melt rates based on observed data. The initializations (P05 and P1) are designed to ensure stability and minimize drift in model outputs, setting the stage for continuation experiments that explore the sensitivity of ice sheet evolution to different basal friction laws.