DOI: https://doi.org/10.1029/2025ms005349
تاريخ النشر: 2026-02-01
المؤلف: Angel Peinado وآخرون
الموضوع الرئيسي: الظواهر الجوية والمحاكاة
نظرة عامة
يتناول هذا القسم من ورقة البحث استخدام تجارب الكواكب المائية لتقييم التقارب الإحصائي لنموذج العواصف العالمية (GSRM) نموذج أيكوسايدرا غير الهيدروستاتيكي (ICON) من خلال تحسينات متتالية في تباعد الشبكة الأفقية، تتراوح من 160 كم إلى 1.25 كم. يستخدم المؤلفون منهجية الاستقراء ريتشاردسون، مستفيدين من التناظر نصف الكروي للكوكب المائي لتحديد التقارب. يقومون بتحليل إحصائيات المناخ المختلفة، مع التركيز بشكل خاص على المقاييس العالمية المتعلقة بالدورة العامة، بما في ذلك الدورة الدورانية العرضية، وهيكل منطقة التقارب الاستوائية، وميزانيات الطاقة والمياه.
تشير النتائج إلى أن التقارب الإحصائي للدورة الدورانية العرضية والإحصائيات العالمية المتوسطة يتحقق عند تباعد شبكة أفقية قدرها 10 كم. ومن الجدير بالذكر أنه بينما يقترب الهيكل الاستوائي من قيمه الحدية عند دقة على مستوى الكيلومتر، فإن شدته، إلى جانب مقاييس مثل مسار العواصف وخصائص تيار النفاث، تتقارب بشكل أبطأ. بالإضافة إلى ذلك، مع زيادة تباعد الشبكة، يتأثر التقارب في ميزانيات الطاقة والمياه بتقليل السحب والتوزيع العرضي لبخار الماء. بشكل عام، تستنتج الدراسة أن نموذج ICON GSRM يظهر تقاربًا إحصائيًا في مقاييس المناخ المختلفة عند دقة 10 كم، خاصة في الهياكل العمودية الكبيرة والنطاق الاستوائي، بينما تظهر بعض المقاييس ميلاً نحو التقارب عند دقات أدق.
مقدمة
تستعرض المقدمة التحديات والمنهجيات المرتبطة بتقييم التقارب في نماذج المناخ، مع التركيز بشكل خاص على نماذج العواصف العالمية (GSRMs) التي تستخدم معادلات تفاضلية جزئية غير خطية مرتبطة لمحاكاة العمليات الجوية. تؤكد على أنه بينما ينبغي أن يظهر نموذج عددي جيد التحديد تقارب الكميات الإحصائية مع انخفاض تباعد الشبكة، فإن الحاجة إلى التوصيف في نماذج المناخ تعقد هذه العلاقة. أظهرت الدراسات السابقة درجات متفاوتة من التقارب في GSRMs، مع نتائج تشير إلى أن فيزياء النموذج يمكن أن تؤثر بشكل كبير على نتائج التقارب. على سبيل المثال، أفادت الدراسات بالتقارب عند تباعدات شبكة معينة، مثل 1 كم لعدد خلايا الحمل، لكنها أيضًا سلطت الضوء على تباين كبير بين النماذج.
لتحقيق تقارب أكثر موثوقية، يقترح المؤلفون استخدام تجارب الكواكب المائية، التي تبسط التفاعلات الحدودية مع الحفاظ على الميزات الأساسية للدورة العامة للأرض. تهدف هذه الدراسة إلى تحديد ما إذا كان نموذج ICON يظهر تقاربًا إحصائيًا عبر تباعدات شبكة أفقية مختلفة، تحديدًا من 160 إلى 1.25 كم، مع استخدام مجموعة محدودة من التوصيفات. من خلال تجنب توصيف الحمل واستخدام إغلاق بسيط للتوربيد، تستهدف البحث حد محاكاة الدوامة الكبيرة (LES). تشمل المنهجية طريقة الاستقراء ريتشاردسون وتركز على الخصائص المناخية الرئيسية مثل منطقة التقارب الاستوائية (ITCZ) والهيكل العمودي للاستوائيات. ستوفر النتائج رؤى حول الاستقرار الهيكلي لحلول النموذج وتؤسس قاعدة للدراسات المستقبلية.
النتائج
في قسم النتائج، تقيم الدراسة تقارب المحاكاة مع تباعدات شبكة أفقية متغيرة، مع التركيز بشكل خاص على شبكة 1.25 كم وفترة أخذ عينات محدودة لمدة 20 يومًا. يتم تعريف التقارب على أنه الاقتراب المنهجي لكمية نحو تقديرها الحدّي ضمن أشرطة عدم اليقين الخاصة بها، حيث تعتبر شبكة 1.25 كم متقاربة عندما تقع ضمن هذه الأشرطة لمقاييس مأخوذة بشكل كافٍ. تحدد التحليل المقاييس التي تظهر أنماط تقارب منهجية، على الرغم من أن بعضها لم يصل بعد إلى تقديراتها الحدية. بالإضافة إلى ذلك، يتم ملاحظة التقارب التذبذبي، حيث قد تشير أشرطة عدم اليقين إلى تقارب محتمل مع أوقات أخذ عينات أطول.
تصنف الدراسة أيضًا التقديرات على أنها متباينة إذا زادت الانحرافات مع تباعد الشبكة الأكثر دقة، بينما تظهر أنماط تقارب غير واضحة عندما تكون أشرطة عدم اليقين كبيرة والثقة منخفضة (|1 – ζ| > 0.25). يتم استخدام تباعد الشبكة 1.25 كم كمعيار لتقييم تقديرات التقارب. ستتناول الأقسام التالية التقارب للإحصائيات العالمية المتوسطة (القسم 3.1)، والمقاييس المتعلقة بخصائص الدورة العامة (القسم 3.2)، وتلك التي تصف الهيكل العمودي في الاستوائيات (القسم 3.3).
المناقشة
في هذا القسم، يناقش المؤلفون المنهجية وإعداد التجارب لدراستهم باستخدام نموذج أيكوسايدرا غير الهيدروستاتيكي (ICON)، مع التركيز على محاكاة الكواكب المائية عبر تباعدات شبكة أفقية مختلفة. يستخدم النموذج معادلات غير هيدروستاتيكية على شبكة مثلثية، مما يسمح بشبكة عالمية شبه موحدة. تستخدم التجارب تكوين كوكب مائي عالمي مع درجة حرارة سطح البحر (SST) محددة واعتدال دائم، مما يحافظ على حالة مناخية متناظرة. يقدم المؤلفون طريقتين للت initialization: بدء بارد، يتطلب فترة طويلة من التهيئة، وبدء ساخن، يقلل بشكل كبير من وقت التهيئة باستخدام شروط أولية من محاكاة تم تهيئتها مسبقًا. يظهرون أن كلا الطريقتين تتقارب إلى نفس الحالة الثابتة، مما يثبت صحة نهج البدء الساخن.
يتناول المؤلفون أيضًا وقت أخذ العينات المطلوب لتحقيق الاستقرار في الحالة المناخية، مشددين على الحاجة إلى فترات أخذ عينات مختلفة بناءً على المتغير الذي يتم تحليله. يستخدمون الارتباط الذاتي لتحديد أوقات أخذ العينات المناسبة، كاشفين أن مقاييس مثل هطول الأمطار تتطلب فترات أقصر مقارنة بخصائص الدورة الكبيرة. يتم إجراء تحليل التقارب باستخدام طريقة الاستقراء ريتشاردسون، التي تقدر أخطاء التقطيع وتقييم موثوقية أنماط التقارب عبر تباعدات الشبكة المختلفة. تشير النتائج إلى معدلات تقارب متفاوتة لكميات مختلفة، حيث تتطلب بعض المقاييس تباعدات شبكة أدق لتحقيق تقارب موثوق. ومن الجدير بالذكر أن الدراسة تجد أن الانتقال من منطقة تقارب استوائية واحدة إلى منطقتين يحدث حول تباعد شبكة قدره 40 كم، مع آثار على الدورة الجوية العامة وديناميات الطاقة.
DOI: https://doi.org/10.1029/2025ms005349
Publication Date: 2026-02-01
Author(s): Angel Peinado et al.
Primary Topic: Meteorological Phenomena and Simulations
Overview
This section of the research paper discusses the use of aquaplanet experiments to assess the statistical convergence of the Global Storm-Resolving model (GSRM) ICOsahedral Nonhydrostatic (ICON) model through successive refinements in horizontal grid spacing, ranging from 160 km to 1.25 km. The authors employ a Richardson extrapolation methodology, leveraging the hemispherical symmetry of the aquaplanet to quantify convergence. They analyze various climate statistics, particularly focusing on global-mean metrics related to general circulation, including the meridional overturning circulation, the structure of the Inter-Tropical Convergence Zone, and the energy and water budgets.
The findings indicate that statistical convergence for the meridional overturning circulation and global mean statistics is achieved at a horizontal grid spacing of 10 km. Notably, while the tropical structure approaches its asymptotic values at kilometer-scale resolution, its intensity, along with metrics like storm-track and jet-stream characteristics, converges more slowly. Additionally, as grid spacing increases, the convergence in energy and water budgets is influenced by cloud reduction and the zonal distribution of water vapor. Overall, the study concludes that the ICON GSRM exhibits statistical convergence in various climate metrics at a 10 km resolution, particularly in large-scale and tropical vertical structures, while some metrics show a tendency toward convergence at finer resolutions.
Introduction
The introduction outlines the challenges and methodologies associated with assessing convergence in climate models, particularly focusing on Global Storm-Resolving Models (GSRMs) that utilize non-linear coupled partial differential equations (PDEs) to simulate atmospheric processes. It emphasizes that while a well-posed numerical model should exhibit convergence of statistical quantities as grid spacing decreases, the necessity for parameterization in climate models complicates this relationship. Previous studies have shown varying degrees of convergence in GSRMs, with findings indicating that model physics can significantly influence convergence outcomes. For instance, studies have reported convergence at specific grid spacings, such as 1 km for convective cell numbers, but also highlighted substantial inter-model variability.
To investigate convergence more robustly, the authors propose using aquaplanet experiments, which simplify boundary interactions while preserving essential features of Earth’s general circulation. This study aims to determine whether the ICON model demonstrates statistical convergence across various horizontal grid spacings, specifically from 160 to 1.25 km, while employing a minimal set of parameterizations. By avoiding convective parameterization and utilizing a simple turbulence closure, the research targets the Large Eddy Simulation (LES) limit. The methodology includes the Richardson extrapolation method and focuses on key climate characteristics such as the Inter-Tropical Convergence Zone (ITCZ) and the vertical structure of the tropics. The results will provide insights into the structural stability of the model’s solutions and establish a foundation for future studies.
Results
In the Results section, the study evaluates the convergence of simulations with varying horizontal grid spacings, specifically focusing on a 1.25 km grid and a limited 20-day sampling period. Convergence is defined as the systematic approach of a quantity towards its asymptotic estimate within its uncertainty bars, with the 1.25 km grid being considered converged when it falls within these bars for sufficiently sampled metrics. The analysis identifies metrics that exhibit systematic convergence patterns, although some have not yet reached their asymptotic estimates. Additionally, oscillatory convergence is noted, where uncertainty bars may suggest potential convergence with longer sampling times.
The study also categorizes estimates as diverging if deviations increase with finer grid spacing, while unclear convergence patterns arise when uncertainty bars are large and confidence is low (|1 – ζ| > 0.25). The 1.25 km grid spacing is utilized as a benchmark for assessing convergence estimates. The subsequent sections will delve into convergence for global mean statistics (Section 3.1), metrics related to general circulation features (Section 3.2), and those characterizing the vertical structure in the tropics (Section 3.3).
Discussion
In this section, the authors discuss the methodology and experimental setup for their study using the ICOsahedral Nonhydrostatic (ICON) model, focusing on aquaplanet simulations across various horizontal grid spacings. The model employs non-hydrostatic equations on a triangular grid, allowing for a quasi-uniform global mesh. The experiments utilize a global aquaplanet configuration with a prescribed sea surface temperature (SST) and a perpetual equinox, maintaining a symmetric climate state. The authors introduce two initialization methods: a cold start, which requires a lengthy spin-up period, and a hot start, which significantly reduces the spin-up time by using initial conditions from a previously spun-up simulation. They demonstrate that both methods converge to the same stationary state, validating the hot start approach.
The authors also address the sampling time required to achieve stationarity in the climate state, highlighting the need for different sampling intervals based on the variable being analyzed. They employ autocorrelation to determine appropriate sampling times, revealing that metrics like precipitation require shorter intervals compared to large-scale circulation characteristics. The convergence analysis is conducted using the Richardson extrapolation method, which estimates discretization errors and assesses the reliability of convergence patterns across different grid spacings. Results indicate varying convergence rates for different quantities, with some metrics requiring finer grid spacings for reliable convergence. Notably, the study finds that the transition from a single to a double Intertropical Convergence Zone (ITCZ) occurs around a grid spacing of 40 km, with implications for the overall atmospheric circulation and energy dynamics.