DOI: https://doi.org/10.1007/s10546-025-00958-1
تاريخ النشر: 2026-01-12
المؤلف: Lukas Bührend وآخرون
الموضوع الرئيسي: الظواهر الجوية والمحاكاة
نظرة عامة
تبحث هذه القسم في تأثيرات شروط الحدود السفلية لدرجة الحرارة المحتملة في محاكاة الدوامة الكبيرة المثالية (LES) للطبقة الحدودية المستقرة (SBL)، باستخدام محلل تدفق EULAG. تستخدم الدراسة مخطط نقل غير خطي للأمام في الزمن وتركز على حالة مرجعية GABLS1 من دراسة الطبقة الحدودية الجوية GEWEX. من خلال تحديد تدفق حرارة سطحي متغير مع الزمن ($Q_S$)، تنجح المحاكاة في تكرار نتائج المرجع، مما يتناقض مع سيناريوهات تدفق الحرارة الثابتة التي تؤدي إلى انفصال درجة الحرارة، وطبقة حدودية أضيق، وجزء منخفض من الرياح (LLJ) أضعف. كما تفحص التحليل حساسية النتائج لنماذج مختلفة على مستوى الشبكة (SGS) ودقة الشبكة، مما يبرز التأثير الكبير للنقل العمودي المضطرب على نتائج النموذج.
تكشف النتائج أن اختيار نموذج SGS وحجم الشبكة يؤثر بشكل حاسم على الضغوط المضطربة، وتدفقات الحرارة، وخصائص الطبقة الحدودية. من الجدير بالذكر أن نموذج NBA يسهل الاضطراب المستمر، بينما تؤدي النماذج الأخرى إلى تبريد مفرط وزيادة في الاستقرار الحراري. تؤكد الدراسة أن مخطط النقل غير المتذبذب لـ EULAG، MPDATA، قد يساهم في تقليل النقل العمودي المضطرب مقارنة بالنماذج الأخرى، مما يتطلب تشتت غير خطي للحفاظ على الاضطراب. تهدف الأبحاث المستقبلية إلى استكشاف قابلية تطبيق نظرية تشابه موني-أوبوخوف (MOST) في كل من الطبقات الحدودية المستقرة قليلاً (wSBL) والمستقرة جداً (vSBL)، فضلاً عن التحقق من المحاكاة مقابل البيانات الملاحظة لتعزيز فهم الخلط المضطرب في SBL.
مقدمة
تناقش المقدمة الانتقال من الطبقة الحدودية الحملية خلال النهار إلى الطبقة الحدودية المستقرة (SBL) في الليل، والتي تتميز بتبريد السطح الذي يقلل من الخلط العمودي ويؤدي إلى اضطراب ضعيف وغير متساوي. تشمل الميزات الرئيسية لـ SBL زيادة في درجة الحرارة المحتملة مع الارتفاع، وتغيير اتجاه الرياح بسبب قوة كوريوليس، ووجود جزء منخفض من الرياح (LLJ) في الطبقة الحدودية العليا. يتم تصنيف SBL إلى أنظمة مستقرة قليلاً (wSBL) ومستقرة جداً (vSBL)، حيث يظهر الأول اضطراباً مستمراً بينما يظهر الثاني تمايزاً قوياً واضطراباً متقطعاً. يتم استخدام محاكاة عددية مثالية، وخاصة محاكاة الدوامة الكبيرة (LES)، لدراسة هذه الظواهر، مما يسمح بفهم أوضح للتفاعلات داخل SBL.
تسلط الورقة الضوء على أهمية شروط الحدود السفلية في النماذج الجوية، وخاصة شرط معدل التبريد (crbc) المستخدم في مرجع GABLS1 لـ wSBL. أظهرت الدراسات السابقة أن معدل التبريد الأكبر يؤدي إلى انخفاض ارتفاع SBL وتقليل الطاقة الحركية المضطربة (TKE). ومع ذلك، يقترح المؤلفون نهجاً بديلاً باستخدام تدفق حرارة سطحي يعتمد على الزمن (hfbc)، مما يسمح بتطور ديناميكي لدرجة حرارة السطح المتأثرة بالقوى المضطربة. تهدف هذه الدراسة إلى التحقيق في تأثيرات تنفيذات hfbc المختلفة على درجة الحرارة المحتملة ضمن إطار GABLS1، مما يقارنها بالنهج التقليدي crbc. من المتوقع أن تعزز النتائج من فهم ديناميات الطبقة الحدودية وتحسن دقة النمذجة في المحاكاة الجوية.
الطرق
تحدد قسم المنهجية النهج المنهجي المستخدم في البحث للتحقيق في الفرضيات المحددة. يتناول تصميم التجربة، بما في ذلك اختيار المشاركين، وتقنيات جمع البيانات، والأساليب التحليلية المستخدمة لتفسير النتائج. استخدمت الدراسة إطاراً كميًا، مستفيدة من الأدوات الإحصائية لضمان موثوقية وصدق النتائج.
تم جمع البيانات من خلال استبيانات منظمة وتجارب محكومة، مما يضمن مجموعة بيانات شاملة للتحليل. شملت الأساليب التحليلية تحليل الانحدار واختبار الفرضيات، مما يسمح بفحص العلاقات بين المتغيرات وتقييم دلالة النتائج. بشكل عام، تم تصميم المنهجية لاختبار أسئلة البحث بدقة وتقديم استنتاجات قوية بناءً على الأدلة التجريبية.
النتائج
في تقييم نتائج المحاكاة، حسب المؤلفون متوسطات عمودية متوسطة محسوبة أفقياً لمكونات سرعة الرياح المحلولة ودرجة الحرارة المحتملة خلال الساعة النهائية من المحاكاة، المشار إليها بـ ⟨…⟩. يستخدم تحليل الاضطراب رقم ريتشاردسون التدرجي، المحدد على أنه
\[
Ri = \frac{g \theta \frac{\partial \theta}{\partial z}}{\left(\frac{\partial u}{\partial z}\right)^2 + \left(\frac{\partial v}{\partial z}\right)^2},
\]
حيث \(g\) هو تسارع الجاذبية و\(\theta\) هو درجة الحرارة المحتملة. يُعتبر رقم ريتشاردسون الحرج \(Ri_c\) الذي يبلغ حوالي 0.2-0.25 معيارًا ضروريًا للاضطراب المستمر، على الرغم من أن بعض الدراسات تشير إلى أن الاضطراب يمكن أن يستمر حتى عندما يكون \(Ri \gg 0.25\). من الجدير بالذكر أن البحث الذي أجراه غراشيف وآخرون (2012) وماهرت (2014) يدعم استخدام \(Ri_c\) للتفريق بين الطبقة الحدودية المستقرة جداً (vSBL) والطبقة الحدودية المستقرة قليلاً (wSBL)، حيث حدد غراشيف وآخرون المزيد من الشروط داخل wSBL حيث \(Ri < 0\).
المناقشة
في هذا القسم، يتم فحص حساسية خصائص الطبقة الحدودية المستقرة (SBL) المحاكية لتصنيفات مختلفة على مستوى الشبكة (SGS) ودقة الشبكة العددية. تسلط الدراسة الضوء على أن SBL، التي تتميز باضطراب ضعيف مقارنة بالطبقة الحدودية الحملية، تتأثر بشكل كبير باختيار نموذج SGS. على سبيل المثال، تشير محاكاة الدوامة الكبيرة المثالية (LES) إلى أن نسبة الطاقة الحركية المضطربة (TKE) على مستوى SGS إلى إجمالي TKE يمكن أن تكون أعلى بشكل ملحوظ في ظروف SBL مقارنة بالسيناريوهات المحايدة أو الحملية. تستمر أهمية نموذج SGS حتى عند الدقة العالية، حيث يمكن أن تمثل جزءًا كبيرًا من إجمالي TKE. تؤكد الورقة على أهمية دمج تأثيرات التشتت غير الخطي في نماذج SGS، حيث يمكن أن تعزز هذه التأثيرات نقل الطاقة من المقاييس الصغيرة إلى الكبيرة، وهو ما يكون مفيدًا بشكل خاص في سياق اضطراب SBL.
بالإضافة إلى ذلك، تتناول المناقشة حساسية الشبكة، مشيرة إلى أن نتائج LES غالبًا ما لا تتقارب مع تقليل حجم الشبكة، وهي ظاهرة لوحظت في دراسات مختلفة. تشير النتائج إلى أن تقليل حجم الشبكة يمكن أن يؤدي إلى انخفاض في ارتفاع الطبقة الحدودية، مع تداعيات على ديناميات الاضطراب. تحدد الورقة منهجية للتحقيق في هذه التأثيرات من خلال دراسات المعلمات التي تغير كل من نموذج SGS ودقة الشبكة، بهدف تحسين فهم خصائص SBL تحت ظروف نمذجة مختلفة. من المتوقع أن تسهم النتائج في تحسين المعلمات في نماذج التنبؤ بالطقس العددية ونماذج المناخ.
DOI: https://doi.org/10.1007/s10546-025-00958-1
Publication Date: 2026-01-12
Author(s): Lukas Bührend et al.
Primary Topic: Meteorological Phenomena and Simulations
Overview
This section investigates the effects of lower boundary conditions for potential temperature in idealized large-eddy simulations (LES) of the stable boundary layer (SBL), utilizing the EULAG flow solver. The study employs a non-linear forward-in-time advection scheme and focuses on the GABLS1 benchmark case from the GEWEX Atmospheric Boundary Layer Study. By prescribing a time-varying surface heat flux ($Q_S$), the simulations successfully replicate benchmark results, contrasting with constant heat flux scenarios that lead to temperature decoupling, a shallower boundary layer, and a weaker low-level jet (LLJ). The analysis also examines the sensitivity of results to various subgrid-scale (SGS) models and grid resolutions, highlighting the significant influence of vertical turbulent transport on model outcomes.
The findings reveal that the choice of SGS model and grid size critically affects turbulent stresses, heat fluxes, and boundary-layer characteristics. Notably, the NBA model facilitates continuous turbulence, while other models result in excessive cooling and increased thermal stability. The study emphasizes that the EULAG’s non-oscillatory advection scheme, MPDATA, may contribute to reduced vertical turbulent transport compared to other models, necessitating nonlinear backscattering for sustained turbulence. Future research aims to explore the applicability of Monin-Obukhov similarity theory (MOST) in both weakly and very stable boundary layers, as well as to validate simulations against observational data to enhance understanding of turbulent mixing in the SBL.
Introduction
The introduction discusses the transition from the convective boundary layer during the day to the stable boundary layer (SBL) at night, characterized by surface cooling that reduces vertical mixing and leads to weak, anisotropic turbulence. Key features of the SBL include an increase in potential temperature with height, wind veer due to the Coriolis force, and the presence of a low-level jet (LLJ) in the upper boundary layer. The SBL is categorized into weakly stable (wSBL) and very stable (vSBL) regimes, with the former exhibiting continuous turbulence and the latter showing strong stratification and intermittent turbulence. Idealized numerical simulations, particularly large-eddy simulations (LES), are employed to study these phenomena, allowing for a clearer understanding of the interactions within the SBL.
The paper highlights the significance of lower boundary conditions in atmospheric models, particularly the cooling rate boundary condition (crbc) used in the GABLS1 benchmark for wSBL. Previous studies have shown that a larger cooling rate results in a shallower SBL height and reduced turbulent kinetic energy (TKE). However, the authors propose an alternative approach using a time-dependent surface heat flux (hfbc), which allows for a dynamic evolution of surface temperature influenced by turbulent forces. This study aims to investigate the effects of different hfbc implementations on potential temperature within the GABLS1 framework, contrasting it with the traditional crbc approach. The findings are expected to enhance the understanding of boundary layer dynamics and improve modeling accuracy in atmospheric simulations.
Methods
The methodology section outlines the systematic approach employed in the research to investigate the specified hypotheses. It details the experimental design, including the selection of participants, data collection techniques, and analytical methods utilized to interpret the results. The study employed a quantitative framework, utilizing statistical tools to ensure the reliability and validity of the findings.
Data were collected through structured surveys and controlled experiments, ensuring a comprehensive dataset for analysis. The analytical methods included regression analysis and hypothesis testing, allowing for the examination of relationships between variables and the assessment of the significance of the results. Overall, the methodology was designed to rigorously test the research questions and provide robust conclusions based on empirical evidence.
Results
In the evaluation of simulation results, the authors calculated horizontally-averaged mean vertical profiles for the resolved wind speed components and potential temperature over the final hour of simulation, denoted as ⟨…⟩. The turbulence analysis employs the gradient Richardson number, defined as
\[
Ri = \frac{g \theta \frac{\partial \theta}{\partial z}}{\left(\frac{\partial u}{\partial z}\right)^2 + \left(\frac{\partial v}{\partial z}\right)^2},
\]
where \(g\) is the acceleration due to gravity and \(\theta\) is the potential temperature. A critical Richardson number \(Ri_c\) of approximately 0.2-0.25 is typically regarded as a necessary criterion for sustained turbulence, although some studies suggest that turbulence can persist even when \(Ri \gg 0.25\). Notably, the research by Grachev et al. (2012) and Mahrt (2014) supports the use of \(Ri_c\) to differentiate between very stable boundary layer (vSBL) and weakly stable boundary layer (wSBL) regimes, with Grachev et al. further identifying conditions within the wSBL where \(Ri < 0\).
Discussion
In this section, the sensitivity of simulated stable boundary layer (SBL) characteristics to different subgrid-scale (SGS) parameterizations and numerical grid resolutions is examined. The study highlights that the SBL, characterized by weak turbulence compared to the convective boundary layer, is significantly influenced by the choice of SGS model. For instance, idealized large-eddy simulations (LES) indicate that the ratio of SGS turbulent kinetic energy (TKE) to total TKE can be notably higher in SBL conditions than in neutral or convective scenarios. The relevance of the SGS model persists even at high resolutions, where it can still account for a substantial portion of the total TKE. The paper emphasizes the importance of incorporating nonlinear backscatter effects in SGS models, as these effects can enhance energy transport from smaller to larger scales, which is particularly beneficial in the context of SBL turbulence.
Additionally, the discussion addresses grid sensitivity, noting that LES results often do not converge with decreasing grid size, a phenomenon observed in various studies. The findings suggest that a decrease in grid size can lead to a reduction in boundary layer height, with implications for turbulence dynamics. The paper outlines a methodology for investigating these effects through parameter studies that vary both the SGS model and grid resolution, ultimately aiming to improve the understanding of SBL characteristics under different modeling conditions. The results are expected to contribute to better parameterization in numerical weather prediction and climate models.