High-precision sea-level proxies along the Tyrrhenian coasts provide new constraints on regional isostatic adjustment and future projections by a Bayesian framework
تستخدم هذه الدراسة إطاراً بايزياً وعلامات أثرية عالية الدقة من ساحل تيريني لعزل التعديل التوازني الإقليمي عن حركة الأرض الرأسية المحلية، مما ينتج عنه إعادة بناء محسنة لمستوى سطح البحر خلال عصر الهولوسين تعمل على تحسين معايرة توقعات مستوى سطح البحر المستقبلية بشكل كبير بما يتجاوز المعايير الحالية للهيئة الحكومية الدولية المعنية بتغير المناخ.
لا يرتفع المحيط أو ينخفض ببساطة في خط منتظم؛ بل يتحرك استجابةً لحوار معقد بين وزن الجليد الذائب، والوزن المتغير للمياه نفسها، والتنفس المرن البطيء للأرض الصلبة تحتنا. فعندما ذابت الصفائح الجليدية الضخمة من العصر الجليدي الأخير، خففت الضغط عن اليابسة، مما سمح للقشرة بالارتداد ببطء نحو الأعلى، بينما دفع وزن الماء المضاف مناطق أخرى نحو الأسفل. هذه العملية، المعروفة باسم التعديل التوازني الجليدي (glacial isostatic adjustment)، مستمرة حتى اليوم، وتُشكل ارتفاع سطح البحر بالنسبة للأرض. ومع ذلك، فإن قياس هذه الحركة أمر صعب لأن الأرض نفسها تتحرك غالباً لأسباب أخرى، مثل استقرار رواسب الأنهار الناعمة أو النشاط البركاني المحلي. وللتنبؤ بمقدار غمر السواحل في المستقبل، يحتاج العلماء إلى فصل الإشارة العالمية لتعديل الأرض عن هذه الخصائص المحلية الغريبة. وبدون هذا التمييز، يمكن أن تكون التوقعات لمناطق محددة مضللة، مما يترك المجتمعات غير مستعدة للحجم الحقيقي لارتفاع المياه.
ركز فريق من الباحثين على ساحل تيريني (Tyrrhenian) في إيطاليا، وهو منطقة حيث اليابسة عبارة عن رقعة من المنحدرات الصخرية، والجزر البركانية، والسهول الرسوبية، لفك تشابك هذه الحركات المتداخلة. لقد لجأوا إلى الماضي بحثاً عن أدلة، باستخدام قياسات عالية الدقة لعلامات مستوى البحر القديمة التي خلفتها الأنشطة البشرية. هذه العلامات، الموجودة في البيئات الساحلية الصخرية مثل المرافئ القديمة وأحواض الأسماك، تعمل كلقطات مجمدة لمكان وقوف البحر قبل آلاف السنين. ولأن هذه الهياكل بُنيت على صخور صلبة بدلاً من الرمال أو الطين المتحرك، فهي تقدم سجلاً أوضح لارتفاع البحر الحقيقي، خالياً من التشويه الناتج عن انضغاط الرواسب. جمع الفريق بيانات من هذه المواقع الصخرية وطبقوا طريقة إحصائية متطورة مصمما لإيجاد نمط مشترك بين القياسات مع تحديد المواقع التي تتصرف بشكل مختلف.
بدأ التحليل بالنظر في جميع البيانات معاً لمعرفة ما إذا كانت تروي قصة واحدة. وجد الباحثون أنه بينما تجمعت معظم المواقع الصخرية حول مستوى ثابت، برزت نقطة بيانات واحدة محددة كقيمة شاذة واضحة. نبعت هذه النقطة من جزيرة بونزا (Ponza)، وهي بيئة بركانية حيث أشارت العلامات القديمة إلى مستوى بحر أقل بكثير، ويرجع ذلك على الأرجح إلى حركات الأرض البركانية-التكتونية المحلية. وبعد استبعاد هذا الشذوذ البركاني، فحص الفريق المواقع المتبقية واكتشف نمطاً متميزاً: فبينما تجمعت عدة قطاعات صخرية معاً بإحكام، أظهرت شبه جزيرة سورينتو (Sorrento) باستمرار إشارة مختلفة. أشارت العلامات القديمة هناك إلى مستوى بحر أقل بكيعاً من القطاعات الصخرية المستقرة الأخرى، مما كشف أن الأرض في شبه جزيرة سورينتو كانت تغرق بالنسبة للمتوسط الإقليمي. ومن خلال عزل القطاعات الصخرية المستقرة حقاً وإزالة كل من الشذوذ البركاني وشبه جزيرة سورينتو الغارقة، نجح الفريق في تنقيح صورتهم لمستوى البحر الإقليمي. وقد حددوا أنه قبل حوالي 2000 عام، كان البحر يقع بنحو 0.57 متراً تحت مستواه الحالي في هذه المناطق المستقرة. ومن هنا، حسبوا أن الأرض في هذه القطاعات المستقرة كانت ترتفع ببطء شديد بمعدل 0.29 ملم في السنة، وهي حركة مدفوعة بتعديل الأرض طويل الأمد للجليد الذائب من الماضي البعيد.
هذا الرقم المنقح أكثر دقة بكثير من التقديرات الواسعة المستخدمة حالياً في تقارير المناخ العالمية. وبينما تضع التوقعات الدولية غالباً نطاقاً واسعاً من عدم اليقين لمراعاة الاختلافات المحلية، توفر هذه الدراسة قيمة محددة بدقة لهذا الجزء من البحر الأبيض المتوسط. ثم استخدم الباحثون هذا الخط المرجعي المستقر لقياس سلوك المناطق الساحلية الأخرى. ووجدوا أن شبه جزيرة سورينتو تغرق بالفعل، وكذلك السهول الساحلية المسطحة حيث تتقلص الرواسب الناعمة تحت وزنها الخاص. وفي المقابل، أظهرت جزيرة بونزا البركانية الحركة الأكثر دراماتيكية، حيث كانت تغرق بمعدل يتفق مع طبيعتها الجيولوجية النشطة. كما قارنت الدراسة هذه المعدلات القديمة بالقياسات الحديثة المأخوذة من الأقمار الصناعية التي تتبع سطح الأرض. تطابقت مجموعتا البيانات في نمطهما، حيث أظهرت أن نفس المناطق التي تغرق اليوم كانت تغرق أيضاً خلال الألفي عام الماضيين، وإن كانت الأقمار الصناعية الحديثة قد رصدت حركة أسرع، مما يعني أنها تلتقط التحولات قصيرة المدى التي يعمل المتوسط طويل الأمد على تسويتها.
يُظهر هذا العمل أنه حتى المناطق التي يُفترض تقليدياً أنها مستقرة يمكن أن تتأثر بحركات الأرض المحلية الدقيقة التي تشوه رؤيتنا للبحر. ومن خلال عزل القطاعات الصخرية المستقرة حقاً، قدم الباحثون معياراً أوضح لكيفية تعديل الأرض للمناخ عبر آلاف السنين. وهذا الوضوح يسمح بفصل أكثر دقة لارتفاع مستوى سطح البحر العالمي عن غرق الأرض المحلي. وبالنسبة للمستقبل، يعني هذا أن التنبؤات بفيضانات السواحل في منطقة تيريني يمكن أن تستند إلى فهم أكثر دقة للأرض التي تحتها، بدلاً من الاعتماد على تصحيحات عامة وواسعة النطاق. تسلط الدراسة الضوء على أنه للتنبؤ بمستقبل سواحلنا، يجب علينا أولاً فهم التاريخ المحدد والهادئ للأرض تحت أقدامنا.
ملخص تقني: توفر مؤشرات مستوى سطح البحر عالية الدقة على طول سواحل تيريني ضوابط جديدة للتعديل التوازني الإقليمي والتوقعات المستقبلية عبر إطار عمل بايزي
بيان المشكلة يعد إعادة بناء تغيرات مستوى سطح البحر النسبي (RSL) في الماضي أمراً ضرورياً للتمييز بين التعديل الجليدي-المائي-الايستاتيكي (GIA) وحركة الأرض الرأسية (VLM)، مما يؤدي إلى تحسين توقعات مستوى سطح البحر في المستقبل. ومع ذلك، فإن قواعد بيانات مستوى سطح البحر النسبي الموجودة في حوض المتوسط تعتمد غالباً على السجلات الطبقية من السهول الساحلية، والتي تكون عرضة لانحيازات منهجية مثل انضغاط الرواسب أو الهبوط المحلي. يمكن لهذه الانحيازات أن تحجب الإشارة الإيستاتيكية الإقليمية، لا سيما عند مقياس السنتيمترات، مما يعقد معايرة نماذج GIA والتحقق من صحة المواقع المرجعية المستقرة تكتونياً. علاوة على ذلك، ورغم أن العلامات الأثرية عالية الدقة في البيئات الساحلية الصخرية توفر قيوداً دقيقة للارتفاع، إلا أن اتساقها الداخلي وقدرتها على عزل إشارة إقليمية متماسكة عبر قطاعات متباعدة جغرافياً لم يتم اختبارها بصرامة باستخدام أطر إحصائية قادرة على قياس التباين المكاني.
المنهجية تستخدم الدراسة إطار عمل هرمي بايزي مكون من أربع خطوات مطبق على مجموعة بيانات عالية الدقة من علامات مستوى سطح البحر الجيوأثرية (SLIPs) من أواخر العصر الهولوسيني والمحفوظة في بيئات ساحلية صخرية على طول ساحل تيريني في إيطاليا. تستبعد مجموعة البيانات السهول الرسوبية لتقليل آثار الانضغاط، مع التركيز على العلامات المقاسة عبر نظام GNSS-RTK.
التقدير الإقليمي (الخطوة 1): يقدر نموذج هرمي بايزي مستوى سطح البحر النسبي الإقليمي عند حوالي 2000 سنة قبل الحاضر (yr BP). يعامل النموذج كل ملاحظة كتحقيق ضوضائي لمتوسط إقليمي كامن (μ) مع عدم يقين في القياس (σi) ومعلمة تشتت إضافية (τ) تأخذ في الاعتبار التباين بين المواقع.
اختبار الاستقرار (الخطوة 2): يتم تحليل البواقي (Residuals) بين القيم المرصودة والتقدير الإقليمي الأولي حسب القطاع الجغرافي. تُبنى نماذج بايزية خاصة بكل قطاع لاختبار فرضية الاستقرار التكتوني. تُحسب احتمالات الحالة اللاحقة الزوجية (P(μA>μB)) لتحديد ما إذا كانت القطاعات تشترك في إشارة واحدة أو تظهر انحرافات منهجية.
التقدير المنقح والمقارنة مع GIA (الخطوة 3): يتم استبعاد القطاعات التي تم تحديدها إحصائياً على أنها غير متسقة مع الإشارة الإقليمية. ويتم حساب تقدير منقح لمستوى سطح البحر النسبي الإقليمي باستخدام القطاعات المستقرة فقط. وتتم مقارنة هذه القيمة المنقحة مع أغلفة التنبؤ لثلاثة نماذج GIA رئيسية (ICE-5G، وICE-6G، وANICE-SELEN) لتقييم الاتساق مع افتراضات ريولوجيا الأرض.
التحليل القائم على القطاعات واستدلال VLM (الخطوة 4): باستخدام التقدير الإقليمي المنقح كخط أساس، يتم حساب معدلات ارتفاع مستوى سطح البحر الظاهري لجميع القطاعات. وتُحول الانحرافات عن المعدل الإقليمي إلى تقديرات حركة الأرض الرأسية (VLM) النسبية للقطاع. وتتم مقارنة هذه التقديرات الجيولوجية مع سرعات GNSS الرأسية لعقود من الزمن لتقي وتقييم أنماط التشوه.
النتائج الرئيسية
تحديد القطاعات المستقرة مقابل غير المستقرة: كشف التحليل الأولي عن تباين مكاني كبير (τ≈0.23 م). حدد اختبار الاستقرار شبه جزيرة سورينتو (القطاع رقم 7) كقطاع ينحرف بشكل منهجي عن الإشارة الإقليمية، حيث كانت قيم مستوى سطح البحر النسبي أقل بنحو 0.5 م من القطاعات الصخرية الأخرى. وقد تم استبعاد هذا القطاع من إعادة البناء الإقليمي النهائي.
مستوى سطح البحر النسبي الإقليمي المنقح ومعدل GIA: ينتج النموذج البايزي المنقح، القائم على القطاعات الصخرية المستقرة (لاتيوم الشمالية، القطاع رقم 2؛ تشيلينتو، القطاع رقم 8؛ والقطاع رقم 3)، تقديراً إقليمياً لمستوى سطح البحر النسبي قدره −0.57 م عند ~2000 سنة قبل الحاضر (فاصل ثقة 95%: −0.66 إلى −0.48 م). وهذا يقابل معدل ارتفاع إيستاتيكي منقح قدره 0.29 ملم في السنة (فاصل ثقة 95%: 0.24–0.33 ملم في السنة).
قيود النماذج: يقع التقدير المنقح ضمن أغلفة التنبؤ لجميع نماذج GIA المختبرة، ولكنه يشغل الجزء العلوي المركزي من أغلفة ICE-5G وICE-6G، مما يوفر قيداً أكثر إحكاماً على بارامترات ريولوجيا الأرض المحتملة مقارنة بمجموعة البيانات الأولية.
توصيف حركة الأرض الرأسية:
القطاعات المستقرة: تظهر القطاعات رقم 2 و3 و8 معدلات ارتفاع ظاهرة متسقة مع الإشارة الإيستاتيكية الإقليمية، مما يؤكد استقرارها التكتوني.
القطاعات الهابطة: تظهر شبه جزيرة سورينتو (القطاع رقم 7) انحرافاً موجباً قدره ~0.28–0.33 ملم في السنة، مما يشير إلى هبوط محلي. وتظهر السهول الرسوبية (القطاعات رقم 5 و6) انحرافات أكبر تتفق مع انضغاط الرواسب.
القطاع البركاني: يسجل جزيرة بونزا (القطاع رقم 4) أكبر شذوذ، مدفوعاً بالعمليات البركانية-التكتونية.
المقارنة الجيولوجية والجيوديسية: يوجد ارتباط إيجابي قوي (r=0.93) بين تقديرات VLM الجيولوجية على المقياس الألفي وسرعات GNSS لعقود من الزمن. ومع ذلك، فإن معدلات GNSS أكبر في المقدار بشكل منهجي (عامل قياس ~1.88)، مما يعكس التقاط التشوهات العابرة أو غير المستقرة في السجلات الجيوديسية قصيرة المدى مقارنة بالمتوسط طويل الأمد.
الأهمية والادعاءات يزعم البحث أن سجلات مستوى سطح البحر النسبي من أواخر العصر الهولوسيني المعتمدة من علامات السواحل الصخرية عالية الدقة يمكن أن تقلل بشكل كبير من عدم اليقين في مكون حركة الأرض الرأسية لتوقعات مستوى سطح البحر. ومن خلال عزل إشارة إيستاتيكية إقليمية متماسكة، توفر الدراسة خط أساساً مقيداً جيولوجياً يعد أكثر دقة بكثير من التصحيحات الإقليمية الواسعة المعتمدة حالياً في توقعات التقرير السادس للهيئة الحكومية الدولية المعنية بتغير المناخ (IPCC AR6) (والتي تستخدم نطاقاً يتراوح بين −0.90 و +1.10 ملم في السنة لهذه المنطقة).
ويرى المؤلفون أن نهجهم يقدم منهجية قابلة للنقل لفصل المساهمات المناخية والإيستاتيكية والتكتونية في تغير مستوى سطح البحر النسبي. وتحديداً، توضح الدراسة أن:
القطاعات التي يُفترض تقليدياً أنها مستقرة تكتونياً (مثل سورينتو) قد تحتوي على حركة رأسية محلية طفيفة ولكنها مهمة، مما يستليزم التحقق الكمي قبل استخدامها في معايرة GIA.
يمكن للنمذجة الهرمية البايزية تحديد واستبعاد القطاعات ذات الإشارات المحلية المتميزة بفعالية، مما يؤدي إلى تحسين معدل الارتفاع الإيستاتيكي الإقليمي.
إن استبدال تصحيحات VLM الإقليمية العامة بتقديرات جيولوجية محددة لكل قطاع ومقيدة تجريبياً يمكن أن يحسن المعايرة المحلية لتوقعات مستوى سطح البحر المستقبلية وتقييمات التأثير الساحلي.
وتخلص الورقة إلى أن دمج مثل هذه البيانات الجيولوجية عالية الدقة والمختبرة للاستقرار في أطر التوقعات التشغيلية يوفر أساساً أكثر قوة لاستراتيجيات التكيف الساحلي في المناطق ذات التباين التكتوني مثل منطقة البحر الأبيض المتوسط.