DOI: https://doi.org/10.1038/s44304-026-00175-x
تاريخ النشر: 2026-02-27
المؤلف: Coleman Blakely وآخرون
الموضوع الرئيسي: أبحاث الأعاصير الاستوائية وغير الاستوائية
نظرة عامة
تدرس هذه الدراسة ضعف خليج البنغال (BoB) تجاه أحداث المد العاصف الناتجة عن الأعاصير الاستوائية (TCs)، نظرًا لارتفاعه المنخفض وكثافة السكان العالية على طول الساحل. باستخدام نموذج هيدروديناميكي ونموذج موجي مترابط عالي الدقة، يقوم الباحثون بمحاكاة عدة آلاف من الأعاصير الاستوائية الاصطناعية تحت سيناريوهات المناخ الحالية والمستقبلية. تركز التحليل على التفاعلات غير الخطية بين المد العاصف، والمد، وتصريف الأنهار، وارتفاع مستوى سطح البحر، مع تحديد ارتفاع الموجات، وتفاعلات المد-العاصف، وارتباط مستويات البحر المتوسطة مع المد العاصف كأهم العمليات المؤثرة على ديناميات المد العاصف.
تشير النتائج إلى تباين مكاني ملحوظ في تأثيرات المد العاصف، مع إنشاء مجموعة بيانات شاملة لمستويات المد العاصف في مواقع البنية التحتية الحيوية. تكشف التوقعات طويلة الأجل (أكثر من 100 عام) عن انخفاض في خطر الفيضانات عبر دلتا غانغس-براهمابوترا-ميغنا الشرقية؛ ومع ذلك، تواجه الساحل الشرقي للهند زيادة في المخاطر. على وجه التحديد، من المتوقع أن يشهد مشروع كوفادا للطاقة النووية زيادة تصل إلى 78% في مستويات المد العاصف لمدة 5000 عام بحلول منتصف القرن، مما يبرز الحاجة الملحة لتدابير التكيف في تخطيط البنية التحتية الساحلية.
الطرق
في هذا القسم، يصف المؤلفون المنهجية المستخدمة لتحليل القيم القصوى، مع التركيز على تقييم مستويات المد العاصف القصوى على مدى فترة 1000 عام عبر نماذج مناخية مختلفة (pc، cmcc، cnrm، ecearth، وhadgem). تم أخذ عينات من كل نموذج 200 مرة مع الاستبدال، واختيار اضطرابات العواصف لتمثيل ظروف المد وتدفق الأنهار المختلفة. تم تحديد مستويات المياه القصوى من أعلى القيم ضمن دائرة نصف قطرها 5 كم من كل موقع من المواقع المعنية. تم استخدام توزيع باريتو العام (GPD) لنمذجة الذيل الثقيل للعينات الناتجة، مع التركيز بشكل خاص على القيم التي تتجاوز عتبة محددة باستخدام نهج القيم القصوى فوق العتبة (POT).
لضمان ملاءمة دقيقة لـ GPD لذيل التوزيع، تم اعتماد منهجية آلية لاختيار العتبة، باستخدام اختبار أندرسون-دارلينغ المعدل (A²R) وتقليل قيمة p. تم تقدير معلمات GPD—الشكل ($k$)، والموقع ($\mu$)، والمقياس ($\sigma$)—باستخدام لحظات L، مما وفر وسيلة فعالة حسابيًا للتعامل مع العدد الكبير من العينات. تضمنت العملية فرز مستويات المياه القصوى، وتحديد العتبات المحتملة، وتناسب GPD لكل عتبة، واختيار العتبة التي تقلل من متمم قيمة p. أخيرًا، تم حساب فترات العودة المتوسطة وفترات الثقة من ملاءمات GPD باستخدام طرق إعادة أخذ العينات غير المعلمية، مما يسمح بتقدير قوي للأحداث القصوى للمد العاصف.
النتائج
يستعرض قسم النتائج تحليلين رئيسيين تم إجراؤهما لتقييم خطر الفيضانات الساحلية الناتجة عن الأعاصير الاستوائية في خليج البنغال، مع التركيز بشكل خاص على مستويات المياه القصوى للمد العاصف. يستكشف التحليل الأول أهمية التفاعلات غير الخطية بين مختلف محركات المد العاصف من خلال مقارنة قائمة على الانحدار مع التركيب الخطي. يستخدم التحليل الثاني تحليل القيم القصوى لتقدير مستويات العودة للمد العاصف لفترات طويلة تحت سيناريوهات المناخ الحالية وبمنتصف القرن. تعرف الدراسة “المد العاصف” بأنه أقصى مستوى مياه إجمالي محاكى في مواقع محددة، مع تضمين الارتفاعات الجوية، والمد الفلكي، وضغط إشعاع الموجات، وتصريف الأنهار، وانحراف مستوى سطح البحر المتوسط.
تُستخدم العواصف الاصطناعية في التحليل باستخدام نموذج توليد الأعاصير الاستوائية الاصطناعية (STORM)، مع توافق الظروف الحالية مع البيانات التاريخية من الأرشيف الدولي لأفضل مسار لرعاية المناخ (IBTrACS) وسيناريوهات المستقبل المستمدة من نماذج المناخ في منتصف القرن من مشروع مقارنة نماذج المناخ الإصدار 6 (CMIP6). تُجرى المحاكاة باستخدام نموذج ADCIRC+SWAN المترابط، الذي يأخذ في الاعتبار تأثيرات الموجات. يتم استكشاف التفاعلات غير الخطية من خلال محاكاة متعددة لكل عاصفة تحت مراحل المد المختلفة وحالات تصريف الأنهار. بالإضافة إلى ذلك، يتم إجراء محاكاة مساعدة لعزل مساهمات آليات القوة الفردية، والتي تتم مقارنتها مع نتائج المجموعة الكاملة لتقييم أهمية التفاعلات غير الخطية. يتم تقدير مستويات العودة طويلة الأجل باستخدام توزيع باريتو العام (GPD) مع نهج القيم القصوى فوق العتبة (POT)، مع تحديد عدم اليقين من خلال إعادة أخذ العينات. يتم تقديم منهجيات مفصلة في قسم “الطرق” من الورقة.
المناقشة
في هذه الدراسة، أجرينا تحليل حساسية لتقييم التفاعلات غير الخطية بين المد العاصف، والمد، وتصريف الأنهار، ومستوى سطح البحر المتوسط على مستويات العودة للمد العاصف في خليج البنغال (BoB). من خلال مقارنة مستويات المياه الإجمالية من مجموعة كاملة من العواصف مع تلك المستمدة من التركيب الخطي، حددنا اختلافات كبيرة تُعزى إلى العمليات غير الخطية مثل تفاعلات المد-العاصف وإعداد الموجات. تم استخدام إطار انحدار لتحديد هذه التأثيرات، مما يكشف أن نموذج “الملاءمة الكاملة”، الذي يتضمن تفاعلات غير خطية مختلفة، يتفوق بشكل كبير على النماذج الخطية الأبسط. من الجدير بالذكر أن التحليل أشار إلى أن إعداد الموجات له تأثير بارز على مستويات المد العاصف، خاصة على الساحل الشرقي للهند، بينما تلعب تفاعلات المد-العاصف أيضًا دورًا حاسمًا في تقييمات مخاطر الفيضانات.
أظهر تحليل القيم القصوى لدينا أن مستويات العودة للمد العاصف تختلف بشكل كبير عبر مواقع مختلفة، حيث تشهد دلتا GBM عمومًا انخفاضًا في مستويات العودة طويلة الأجل تحت سيناريوهات المناخ المستقبلية، بينما تواجه الساحل الشرقي للهند مخاطر متزايدة. تشير النتائج إلى أن الاعتماد على التركيب الخطي يمكن أن يقلل أو يبالغ في تقدير مستويات العودة بنسبة تصل إلى 40%. تؤكد الدراسة على ضرورة تقييم مخاطر الفيضانات المحددة حسب المنطقة ودمج أطر النمذجة عالية الدقة في تصميم البنية التحتية واستراتيجيات الاستعداد للكوارث، خاصة مع تغير الظروف المناخية التي تغير مخاطر الفيضانات في المناطق الساحلية. تشمل اتجاهات البحث المستقبلية تحسين دقة النموذج من خلال دمج تأثيرات هبوط الأرض واستخدام تقنيات التعلم الآلي لتحسين كفاءة محاكاة العواصف.
DOI: https://doi.org/10.1038/s44304-026-00175-x
Publication Date: 2026-02-27
Author(s): Coleman Blakely et al.
Primary Topic: Tropical and Extratropical Cyclones Research
Overview
This study investigates the vulnerability of the Bay of Bengal (BoB) to storm-tide events driven by tropical cyclones (TCs), given its low elevation and high population density along the coast. Utilizing a high-resolution, coupled hydrodynamic and wave model, the researchers simulate several millennia of synthetic TCs under both present and future climate scenarios. The analysis focuses on the nonlinear interactions among storm surge, tides, river discharge, and sea-level rise, identifying wave runup, tide-surge interactions, and the coupling of mean sea levels with surge as the most significant processes influencing storm-tide dynamics.
The findings indicate notable spatial variability in storm-tide impacts, with the generation of a comprehensive dataset of storm-tide levels at critical infrastructure sites. Long-term projections (over 100 years) reveal a decrease in flood risk across the eastern Ganges-Brahmaputra-Meghna delta; however, the eastern coast of India faces increased risk. Specifically, the Kovvada Atomic Power Project is projected to experience up to a 78% increase in 5000-year storm-tide levels by mid-century, highlighting the urgent need for adaptive measures in coastal infrastructure planning.
Methods
In this section, the authors describe the methodology employed for extreme value analysis, focusing on the assessment of peak storm-tide levels over a 1000-year period across various climate models (pc, cmcc, cnrm, ecearth, and hadgem). Each model was sampled 200 times with replacement, selecting storm perturbations to represent different tidal and riverine flow conditions. The peak water levels were determined from the highest values within a 5 km radius of each site of interest. The Generalized Pareto Distribution (GPD) was utilized to model the heavy tails of the resulting samples, specifically focusing on values exceeding a defined threshold using a peaks over threshold (POT) approach.
To ensure accurate fitting of the GPD to the tail of the distribution, an automated methodology was adopted for threshold selection, employing the modified Anderson-Darling (A²R) test and p-value minimization. The parameters of the GPD—shape ($k$), location ($\mu$), and scale ($\sigma$)—were estimated using L-Moments, which provided a computationally efficient means to handle the large number of samples. The process involved sorting maximum water levels, identifying potential thresholds, fitting the GPD for each threshold, and selecting the threshold that minimized the complement of the p-value. Finally, mean return periods and confidence intervals were calculated from the GPD fits using non-parametric bootstrapping methods, allowing for robust estimation of extreme storm-tide events.
Results
The results section outlines two primary analyses conducted to assess the coastal flood risk from tropical cyclones in the Bay of Bengal, specifically focusing on extreme storm-tide water levels. The first analysis investigates the significance of nonlinear interactions among various storm-tide drivers through a regression-based comparison with linear superposition. The second analysis employs extreme value analysis to estimate long-period storm-tide return levels under current and mid-century climate scenarios. The study defines “storm-tide” as the maximum simulated total water level at specific sites, incorporating meteorological surge, astronomical tides, wave radiation stress, riverine discharge, and mean sea level offset.
The synthetic storms utilized in the analysis are generated using the Synthetic Tropical cyclOne geneRation Model (STORM), with present-day conditions aligned with historical data from the International Best Track Archive for Climate Stewardship (IBTrACS) and future scenarios derived from mid-century climate models from the Climate Model Intercomparison Project version 6 (CMIP6). The simulations are conducted using a coupled ADCIRC+SWAN model, which accounts for wave effects. Nonlinear interactions are explored through multiple simulations of each storm under varying tidal phases and riverine discharge states. Additionally, auxiliary simulations are performed to isolate the contributions of individual forcing mechanisms, which are compared to the full ensemble results to assess the importance of nonlinear interactions. Long-term return levels are estimated using a Generalized Pareto Distribution (GPD) with a Peaks Over Threshold (POT) approach, with uncertainty quantified through bootstrap resampling. Detailed methodologies are provided in the “Methods” section of the paper.
Discussion
In this study, we conducted a sensitivity analysis to evaluate the nonlinear interactions among storm surge, tides, riverine discharge, and mean sea level on storm-tide return levels in the Bay of Bengal (BoB). By comparing total water levels from a full storm ensemble with those derived from linear superposition, we identified significant discrepancies attributed to nonlinear processes such as tide-surge interactions and wave setup. A regression framework was employed to quantify these effects, revealing that the “Full Fit” model, which incorporates various nonlinear interactions, significantly outperforms simpler linear models. Notably, the analysis indicated that wave setup has a pronounced impact on storm-tide levels, particularly along the east coast of India, while the tide-surge interaction also plays a critical role in flood risk assessments.
Our extreme value analysis demonstrated that storm-tide return levels vary considerably across different sites, with the GBM Delta generally experiencing decreased long-term return levels under future climate scenarios, while the eastern coast of India faces increased risks. The findings suggest that reliance on linear superposition can underestimate or overestimate return levels by up to 40%. The study emphasizes the necessity of region-specific flood risk assessments and the integration of high-resolution modeling frameworks in infrastructure design and disaster preparedness strategies, particularly as changing climatic conditions alter flood risks in coastal zones. Future research directions include enhancing model accuracy by incorporating land subsidence effects and employing machine learning techniques to improve storm simulation efficiency.