DOI: https://doi.org/10.5194/gmd-19-261-2026
تاريخ النشر: 2026-01-08
المؤلف: Nicolas Mokus وآخرون
الموضوع الرئيسي: ديناميات الجليد في القطب الشمالي والقطب الجنوبي
نظرة عامة
تدرس الدراسة الانفصال الناتج عن الموجات في الجليد البحري، وهي عملية حاسمة لفهم تشكيل منطقة الجليد الهامشية في المحيطات القطبية. إن النمذجة الدقيقة لديناميات الجليد البحري ضرورية لكل من الأبحاث التشغيلية والمناخية؛ ومع ذلك، هناك نقص في الإجماع حول المعايير المناسبة لتكسير الجليد تحت تأثير الموجات. تقليديًا، تم نمذجة التكسير باستخدام إطار هيدروإيلاستيك مع معيار قائم على الحد الأقصى للإجهاد. تقدم هذه الأبحاث نهجًا بديلًا قائمًا على الطاقة من خلال نموذج أحادي البعد جديد يسمى SWIIFT (أثر الموجات السطحية على الجليد البحري – مجموعة أدوات التكسير)، والذي يستخدم نظرية الألواح الخطية لمحاكاة التكسير الناتج عن الموجات في مجال الزمن.
تم التحقق من صحة نموذج SWIIFT مقابل التجارب المعملية الحالية التي تتضمن تكسير الجليد الناتج عن الموجات لمادة مماثلة، مما يوضح قدرته على إعادة إنتاج أنماط التكسير المرصودة. ومن الجدير بالذكر أن النتائج تكشف أنه تحت ظروف موجية معينة، تتميز بعدد موجي غير بعدي، يفشل النموذج في التنبؤ بالتكسير عند انحناء ثابت. تثير هذه النتيجة تساؤلات حول صحة استخدام معيار الإجهاد الأقصى لتوصيف تكسير الجليد البحري، مما يشير إلى أن نهجًا قائمًا على الطاقة قد يوفر تمثيلًا أكثر دقة للعمليات الفيزيائية الأساسية.
مقدمة
تناقش المقدمة تأثير تغير المناخ على الجليد البحري في القطب الشمالي، مع تسليط الضوء على ظهور مناطق المحيط المفتوح التي تعرض الجليد البحري لنشاط موجي متزايد. يؤدي ذلك إلى جليد أصغر سناً وأرق وأكثر هشاشة، مما يجعله عرضة للتجزئة، مما يؤدي إلى حلقة تغذية راجعة إيجابية حيث يؤدي الانفصال الناتج عن الموجات إلى إنشاء كتل أصغر في منطقة الجليد الهامشية (MIZ). تختلف ديناميات الجليد المجزأ بشكل كبير عن الجليد المتماسك، مما يجعله أكثر حركة وحساسية للذوبان، مما قد يزيد من فقدان الجليد البحري خلال العواصف. تفشل النماذج المناخية الحالية في تمثيل عمليات التجزئة بدقة، مما يستلزم تحسين نمذجة منطقة الجليد الهامشية لتعزيز توقعات تطور الجليد البحري.
تهدف الورقة إلى وضع معيار للتكسير للجليد البحري تحت تأثير الموجات، مع معالجة نقص الإجماع حول هذه القضية داخل المجتمع العلمي. تعتمد النماذج الحالية بشكل أساسي على إجهاد الانحناء المحلي أو عتبات الإجهاد، لكن المؤلفين يقترحون معيار تكسير قائم على الطاقة مستمد من ميكانيكا التكسير التي تأخذ في الاعتبار التأثيرات غير المحلية. يسمح هذا النموذج بمحاكاة انتشار التكسير بمرور الوقت ويتم التحقق من صحته مقابل التجارب المعملية على المواد المماثلة للجليد. تؤكد الدراسة على الحاجة إلى فهم شامل لتشوه الجليد الناتج عن الموجات والتكسير لتحسين دقة نماذج الجليد البحري.
طرق
في هذا القسم، يقدم المؤلفون مقارنة بين نتائجهم العددية والبيانات التجريبية من Auvity et al. (2025)، كما هو موضح في الشكل 10. تشير النتائج إلى أن ملفات السعة الحرجة متسقة عبر خصائص ميكانيكية متنوعة للمادة المحاكية، مع ثابت مضاعف يزداد مع صلابة الانحناء أو الطول. يظهر ترتيب حجم السعة الحرجة واتجاهها المتناقص مع زيادة \( kL_D \) توافقًا نوعيًا بين المحاكاة والتجارب، على الرغم من ظهور اختلافات بسبب الطبيعة غير الخطية للموجات التجريبية، التي تنتج سعات حرجة تختلف كـ \( k^{-1} \) مقارنة بـ \( k^{-2} \) المرصودة في المحاكاة العددية.
كما يبرز المؤلفون الاختلافات في تعريف طول الاسترخاء \( L_\kappa \) المستخدم في محاكاة الموجات الخطية الخاصة بهم مقابل العرض الكامل عند نصف الحد الأقصى (FWHM) الذي استخدمه Auvity et al. (2025) في تجاربهم غير الخطية. لمعالجة ذلك، يمثلون السعات الحرجة التجريبية كدالة لـ \( L_{D_{FWHM}} \) (المربعات الوردية في الشكل 10)، مما يعدل البيانات التجريبية لتتوافق بشكل أفضل مع نتائجهم العددية. يهدف هذا التعديل إلى تسوية الاختلافات الناشئة عن سلوك الموجات المتميز، لا سيما في المناطق التي تشير فيها كلا النموذجين إلى أن الانحناء الحرج يتأثر بطول موجة القوة، مما يتحدى فكرة وجود عتبة إجهاد ثابتة.
نتائج
في هذا القسم، يقدم المؤلفون نتائج حول اكتشاف عتبة التكسير باستخدام صيغة طاقة، مقارنة نتائجهم مع التجارب المعملية التي أجراها Auvity et al. (2025). باستخدام معلمات المادة من Auvity et al. (2025) كما هو موضح في الجدول 1، يمتنع المؤلفون عن ضبط معلمات النموذج للحفاظ على نزاهة تحليلهم. تركز الدراسة على تحديد السعة الحرجة، المشار إليها بـ $a_{cr}$، من خلال طريقة بحث ثنائية النظام، مع التأكيد على فحص الحالات شبه الثابتة دون دمج مصطلح التعب.
تنظم النتائج في قسمين فرعيين: القسم 4.1 يستكشف تأثيرات تغيير عدد الموجات، بينما يوسع القسم 4.2 هذا التحليل من خلال تغيير المعلمات الميكانيكية للوحة أيضًا. يسمح هذا النهج المنظم بفهم شامل للعوامل التي تؤثر على عتبة التكسير، مما يساهم في تقديم رؤى قيمة في مجال علوم المواد وميكانيكا التكسير.
مناقشة
في هذا القسم، يقدم المؤلفون نموذجًا شاملاً لتشوه وتكسير كتل الجليد العائمة تحت تأثير الموجات. يعتمد النموذج على نظرية الألواح الرقيقة لكيرشهوف-لوف، مع اعتبار الجليد مرنًا ومتجانسًا، مع الطفو والجاذبية كالقوى الرئيسية المؤثرة على الكتل. يتم اشتقاق المعادلات الحاكمة لوصف انحراف كتل الجليد استجابةً للاختلالات في سطح السائل، مما يؤدي إلى معادلة تفاضلية تلتقط سلوك الانحناء للجليد. يقدم المؤلفون معيار تكسير قائم على مبادئ تقليل الطاقة من ميكانيكا التكسير، ويقارنونه مع معيار إجهاد تقليدي. يفترض معيار الطاقة أن كتلة ستنكسر إذا تجاوزت الطاقة المرنة المفرج عنها من الانفصال تكلفة الطاقة لإنشاء أسطح جديدة، بينما يسمح معيار الإجهاد بالتكسير عندما يتجاوز إجهاد الانحناء المحلي عتبة حرجة.
يستعرض المؤلفون التنفيذ العددي لنموذجهم، بما في ذلك حساب القوى الناتجة عن الموجات وعملية تكرارية لمحاكاة تشوه الكتل والتكسير المحتمل. يؤكدون على أهمية تقدير معلمات التكسير بدقة، مثل معدل تحرير الطاقة والإجهاد الحرج، التي تتأثر بالظروف البيئية وخصائص الجليد. يختتم القسم بمثال محاكاة يوضح قدرة النموذج على التنبؤ بأحداث التكسير بمرور الوقت، مع تسليط الضوء على الاختلافات بين النهجين القائمين على الطاقة والإجهاد. تشير النتائج إلى أن معيار الطاقة قد يوفر إطارًا أكثر ارتباطًا بالفيزياء لفهم انفصال كتل الجليد تحت ظروف الموجات الديناميكية.
القيود
في قسم القيود، يناقش المؤلفون القيود المفروضة من خطية نموذجهم العددي بالنسبة للنتائج التجريبية لـ Auvity et al. (2025). تتطلب التجارب المعملية وجود موجات غير خطية (تتميز بانحدار موجة $a_k \approx 0.14$) لحدوث التكسير، وهو ما لا يتم تمثيله في نموذج المؤلفين. لتسهيل ملاحظة التكسير ضمن إطارهم، يخففون الشرط على انحدار الموجة، مع تحديد حد أعلى قدره $a_k \leq 0.5$. يضع هذا التعديل تحليلهم خارج النطاق الخطي، الذي هو أساس نموذجهم.
على الرغم من هذه القيود، يجادل المؤلفون بأنها غير ذات أهمية لتقييمهم النوعي لسلوك النموذج. يشيرون إلى أنه بالنسبة للمعلمات النموذجية للجليد البحري – تحديدًا السماكة ومعامل يونغ – يحدث التكسير ضمن نطاق خطي، كما هو موضح في نتائجهم حيث $a_k = 0.015$. بالإضافة إلى ذلك، يوضحون أن انحدار الموجة يتم تعريفه بالنسبة للموجات الهيدروإيلاستيك التي تنتشر تحت اللوحة المرنة، بدلاً من الموجات السطحية الحرة، التي لها علاقة انتشار مختلفة. يبرز هذا التمييز انحدار الموجات الهيدروإيلاستيك مقارنة بموجات الجاذبية السطحية الحرة، مما يضع القيود على نموذجهم في سياق أكبر.
DOI: https://doi.org/10.5194/gmd-19-261-2026
Publication Date: 2026-01-08
Author(s): Nicolas Mokus et al.
Primary Topic: Arctic and Antarctic ice dynamics
Overview
The study investigates the wave-induced breakup of sea ice, a critical process for understanding the formation of the marginal ice zone in polar oceans. Accurate modeling of sea ice dynamics is essential for both operational and climate research; however, there is a lack of consensus on the appropriate criteria for ice fracturing under wave action. Traditionally, fracture has been modeled using a hydroelastic framework with a maximum strain-based criterion. This research introduces an alternative energy-based approach through a new one-dimensional model called SWIIFT (Surface Wave Impact on sea Ice – Fracture Toolkit), which utilizes linear plate theory to simulate wave-induced fracture in the time domain.
The SWIIFT model is validated against existing laboratory experiments involving wave-induced fracture of an analogue material, demonstrating its capability to reproduce observed fracture patterns. Notably, the findings reveal that under certain wave conditions, characterized by a dimensionless wavenumber, the model fails to predict fracture at constant curvature. This outcome raises questions about the validity of using a maximum strain criterion for parametrizing sea ice fracture, suggesting that an energy-based approach may provide a more accurate representation of the underlying physical processes.
Introduction
The introduction discusses the impact of climate change on Arctic sea ice, highlighting the emergence of open ocean areas that expose sea ice to intensified wave activity. This results in younger, thinner, and more fragile ice that is prone to fragmentation, leading to a positive feedback loop where wave-induced breakup creates smaller floes in the marginal ice zone (MIZ). The dynamics of fragmented ice differ significantly from consolidated ice, making it more mobile and sensitive to melting, which can exacerbate sea ice loss during storms. Current climate models fail to accurately represent the fragmentation processes, necessitating improved modeling of the MIZ to enhance predictions of sea ice evolution.
The paper aims to establish a fracture criterion for sea ice under wave forcing, addressing the lack of consensus on this issue within the scientific community. Existing models primarily rely on local flexural stress or strain thresholds, but the authors propose an energy-based fracture criterion derived from fracture mechanics that accounts for non-local effects. This model allows for the simulation of fracture propagation over time and is validated against laboratory experiments on ice analogues. The study emphasizes the need for a comprehensive understanding of wave-induced ice deformation and fracture to improve the accuracy of sea ice models.
Methods
In this section, the authors present a comparison of their numerical results with experimental data from Auvity et al. (2025), as illustrated in Fig. 10. The findings indicate that the critical amplitude profiles are consistent across varying mechanical properties of the simulated material, with a multiplicative constant that increases with flexural rigidity or length. The order of magnitude of the critical amplitude and its decreasing trend with increasing \( kL_D \) show qualitative agreement between simulations and experiments, although discrepancies arise due to the nonlinear nature of the experimental waves, which yield critical amplitudes that vary as \( k^{-1} \) compared to the \( k^{-2} \) observed in the numerical simulations.
The authors also highlight differences in the definition of the relaxation length \( L_\kappa \) used in their linear wave simulations versus the full width at half-maximum (FWHM) employed by Auvity et al. (2025) in their nonlinear experiments. To address this, they represent the experimental critical amplitudes as a function of \( L_{D_{FWHM}} \) (pink squares in Fig. 10), which adjusts the experimental data to better align with their numerical results. This adjustment aims to reconcile the differences arising from the distinct wave behaviors, particularly in regions where both models indicate that critical curvature is influenced by the forcing wavelength, thus challenging the notion of a constant strain threshold.
Results
In this section, the authors present findings on the detection of a fracture threshold using an energy formalism, comparing their results with laboratory experiments conducted by Auvity et al. (2025). Utilizing material parameters from Auvity et al. (2025) as outlined in Table 1, the authors refrain from tuning model parameters to maintain the integrity of their analysis. The study focuses on identifying the critical amplitude, denoted as $a_{cr}$, through a systematic dichotomic search method, emphasizing the examination of quasi-static states without incorporating a fatigue term.
The results are organized into two subsections: Section 4.1 investigates the effects of varying the wavenumber, while Section 4.2 extends this analysis by also altering the mechanical parameters of the plate. This structured approach allows for a comprehensive understanding of the factors influencing the fracture threshold, contributing valuable insights to the field of material science and fracture mechanics.
Discussion
In this section, the authors present a comprehensive model for the deformation and fracture of floating ice floes under wave action. The model is grounded in Kirchhoff-Love thin-plate theory, treating the ice as elastic and isotropic, with buoyancy and gravity as the primary forces acting on the floes. The governing equations are derived to describe the deflection of the ice floes in response to perturbations in the fluid surface, leading to a differential equation that captures the bending behavior of the ice. The authors introduce a fracture criterion based on energy minimization principles from fracture mechanics, contrasting it with a conventional strain criterion. The energy criterion posits that a floe will fracture if the elastic energy released by the breakup exceeds the energy cost of creating new surfaces, while the strain criterion allows for fracture when local bending strain surpasses a critical threshold.
The authors detail the numerical implementation of their model, including the calculation of wave-induced forces and the iterative process for simulating floe deformation and potential fracture. They emphasize the importance of accurately estimating fracture parameters, such as the energy release rate and critical strain, which are influenced by environmental conditions and ice properties. The section concludes with a simulation example demonstrating the model’s capability to predict fracture events over time, highlighting the differences between the energy and strain-based approaches. The results indicate that the energy-based criterion may provide a more physically grounded framework for understanding ice floe breakup under dynamic wave conditions.
Limitations
In the section on limitations, the authors discuss the constraints imposed by the linearity of their numerical model in relation to the experimental findings of Auvity et al. (2025). The laboratory experiments necessitate the presence of nonlinear waves (characterized by a wave slope $a_k \approx 0.14$) for fracture to occur, which is not represented in the authors’ model. To facilitate the observation of fracture within their framework, they relax the condition on the wave slope, setting an upper limit of $a_k \leq 0.5$. This adjustment places their analysis outside the linear regime, which is the basis of their model.
Despite this limitation, the authors argue that it is inconsequential for their qualitative assessment of the model’s behavior. They note that, for parameters typical of sea ice—specifically thickness and Young’s modulus—fracture occurs within a linear regime, as demonstrated in their results where $a_k = 0.015$. Additionally, they clarify that the wave slope is defined concerning the hydroelastic waves propagating beneath the elastic plate, rather than the free surface waves, which have a different dispersion relation. This distinction highlights the relative steepness of hydroelastic waves compared to free surface gravity waves, further contextualizing the limitations of their model.