DOI: https://doi.org/10.1007/s10040-025-02994-0
تاريخ النشر: 2026-01-10
المؤلف: Sophie ten Hietbrink وآخرون
الموضوع الرئيسي: أبحاث الجيولوجيا وعلم المناخ القديم
نظرة عامة
تقدم هذه الدراسة أول محاكاة عددية لتخفيف المياه الجوفية في البحر قبالة الهامش القاري النرويجي، كاشفة عن الدور الكبير للجليد في هذه العملية. تشير الأبحاث إلى أنه خلال العصر الجليدي المتأخر والهولوسين، أدى تقدم الغطاء الجليدي الفينوسكاندي إلى إنشاء تدرجات رأس هيدروليكي عمودية سهلت تسرب مياه ذوبان الجليد إلى نظام المياه الجوفية البحرية بمعدل 486 مم في السنة. بحلول 23,000 سنة مضت، توسعت مناطق المياه الجوفية المنعشة بالكامل إلى 21.5 كم²، بينما وصلت المناطق ذات تركيزات الكلوريد أقل من 400 مللي مول إلى 47.4 كم²، وكلاهما أكبر بكثير من التقديرات المستندة فقط إلى تغييرات مستوى سطح البحر الإيستاتي. كما تسلط الدراسة الضوء على أن تراجع الجليد أدى إلى تصريف المياه الجوفية بمعدلات تصل إلى 246 مم في السنة، مما قد يعيق إعادة شحن مياه البحر في المناطق التي كانت مغطاة بالجليد سابقًا.
تتحدى النتائج الافتراضات السابقة المتعلقة بالاتصالات الجانبية بين اليابسة والمحيط كآلية رئيسية لتخفيف المياه البحرية، بدلاً من ذلك تُظهر أن التثبيت العمودي أثناء الجليد هو السائد. تؤكد استمرارية المياه الجوفية المنعشة في العمق خلال الفترات بين الجليدية على التأثير طويل الأمد للجليد على الهيدرولوجيا تحت السطح. تقترح الدراسة أن تدرجات رأس الهيدروليكي الناتجة عن تحميل الجليد تختلف عن تلك الناتجة عن تغييرات مستوى سطح البحر، وأن معدلات تصريف المياه الجوفية البحرية الحالية الملاحظة في موقع لوفوتن-فسترولين لا يمكن تفسيرها فقط من خلال هذه الآليات. تمتد تداعيات هذا البحث عالميًا، حيث قد تحدث عمليات مماثلة في هوامش جليدية أخرى، وقد يؤدي التراجع المستمر للصفائح الجليدية بسبب تغير المناخ إلى تغيير ديناميات المياه الجوفية وتدفقات المواد المذابة إلى المحيطات الساحلية، مما قد يؤدي إلى إطلاق الميثان من الرواسب تحت السطح.
مقدمة
تناقش مقدمة هذه الورقة البحثية الدور الكبير لتصريف المياه الجوفية البحرية (SGD) في تشكيل البيئة الفيزيائية والكيميائية البحرية، خصوصًا في الدورات البيوجيوكيميائية والنظم البيئية الساحلية. تسلط الضوء على كيفية تعزيز الجليد لتدرجات الهيدروليك في المناطق التي تم تجليدها مؤخرًا، مما يمكن أن يزيد من تصريف السوائل البحرية، مما يؤدي إلى تسرب المياه الجوفية من الخزانات الكامنة. هذه العملية حاسمة لنقل الكربون والمواد المغذية وغيرها من المواد المذابة إلى البيئات البحرية، وبالتالي تؤثر على كيمياء قاع البحر وديناميات المواد المغذية. تؤكد الورقة على الحاجة لفهم الضوابط على تدفق المياه الجوفية البحرية وتوزيع المياه العذبة على مدى فترات زمنية جليدية-بين جليدية، خاصة في ضوء التغيرات المعاصرة بسبب تراجع الجليد وارتفاع مستويات سطح البحر.
تشير المقدمة أيضًا إلى أنه بينما يتم التعرف على مبادئ تخفيف المياه البحرية، لا تزال التغيرات المكانية والزمانية لـ SGD غير مفهومة جيدًا. تشير إلى أن الدراسات الحالية قد ركزت بشكل أساسي على الرف القاري الأطلسي قبالة نيو إنجلاند ورف غرينلاند الجنوبي الشرقي، مع وجود فجوة ملحوظة في المحاكاة العددية للهامش النرويجي، على الرغم من وجود أدلة على التخفيف في مواقع متعددة. تدعو الورقة إلى فهم متماسك للآليات التي تحرك التخفيف في الهامش القاري النرويجي، وهو أمر ضروري لفهم أنظمة المياه الجوفية المتأثرة بالديناميات الجليدية وللتطبيقات العملية مثل تخزين الكربون في التكوينات البازلتية.
طرق البحث
في هذه الدراسة، تم تطوير نموذج هيدروجيولوجي باستخدام COMSOL Multiphysics (الإصدار 6.1) لمحاكاة تدفق المياه الجوفية ونقل المواد المذابة في منطقة لوفوتن-فسترولين (LV). يستخدم النموذج قانون دارسي لوصف تدفق المياه الجوفية عبر الوسائط المسامية، ممثلة بالمعادلة \( u = -\nabla p \)، حيث \( u \) هو مجال سرعة دارسي، \( \kappa \) هو النفاذية، \( \mu \) هو اللزوجة الديناميكية، و \( p \) هو ضغط المسام. تتضمن معادلة الاستمرارية قابلية الانضغاط \( S \)، وكثافة السائل \( \rho \)، وعبارة مصدر الكتلة \( Q_m \)، التي تم تعيينها على صفر بسبب عدم وجود مصادر سائلة إضافية. يتم نمذجة نقل المواد المذابة، وبشكل خاص الكلوريد، من خلال ربط مجال سرعة دارسي بمعادلة نقل تأخذ في الاعتبار تراكم المواد المذابة، والامتزاز، والتشتت الهيدروديناميكي، والحمل.
تستند هندسة النموذج إلى التفسيرات الزلزالية وتمتد على مدى أفقي يبلغ 160 كم مع مدى عمودي يبلغ حوالي 6 كم، مقسمة إلى ثمانية وحدات طبقية. تم إجراء تحليلات حساسية لتقييم تأثير تغيير المعلمات مثل المسامية، والنفاذية، وقابلية الانضغاط، والتشتت، مع تحديد نطاقات واقعية لهامش لوفوتن-فسترولين القاري. تتماشى القيم المختارة للمسامية (0.15 لوحدة القاعدة و0.2 للوحدات غير القاعدية) وقابلية الانضغاط (5 × 10⁻⁶ با) مع الأدبيات الموجودة. يفترض النموذج هندسة طبقية ثابتة على مدى زمني قدره 120,000 سنة ويستخدم شبكة مثلثية حرة من أجل الكفاءة الحاسوبية، مع نهج خطوة زمنية تكيفية لحل التغيرات الزمنية بشكل فعال.
النتائج
في هذه الدراسة، تم تحليل سيناريوهين نموذجيين لتقييم تأثير تحميل الجليد على ارتفاع مستوى سطح البحر النسبي. السيناريو الأول، المشار إليه باسم GLAC، يتضمن تأثيرات تحميل الجليد، ممثلة بالمنحنى الأرجواني في الشكل 4b. السيناريو الثاني، المسمى EUS، يتبع منحنى مستوى سطح البحر الإيستاتي دون أي تحميل جليدي، موضحًا بالمنحنى الأخضر في الشكل 4b. من خلال مقارنة هذين السيناريوهين، تسلط الأبحاث الضوء على المساهمة الكبيرة لتحميل الجليد في التغيرات في مستويات سطح البحر النسبي.
المناقشة
في هذا القسم، يناقش المؤلفون المنهجيات والنتائج المتعلقة بشروط حدود رأس الهيدروليك وديناميات المياه الجوفية المتأثرة بتحميل الجليد في منطقة الصفيحة الجليدية الأوراسية. تم اشتقاق بيانات سمك الصفيحة الجليدية من محاكاة الدورات الجليدية، مع التركيز بشكل خاص على الحد الأقصى للجليد الأخير (LGM) وتأثيره على رأس الهيدروليك ومستوى سطح البحر النسبي (RSL). تسلط الدراسة الضوء على التباين الكبير بين RSL ومنحنيات مستوى سطح البحر الإيستاتي خلال الجليد، المدفوعة بالاكتئاب الإيزوستاتي، وتحدد الطرق المستخدمة لاستقراء بيانات RSL إلى ما بعد الـ 22,000 سنة المتاحة. تم حساب رأس الهيدروليك الناتج عن كتلة الجليد، وتم تطبيق شرط حدود ديريشليه لنمذجة ديناميات رأس الهيدروليك، مما يربط بين توسع الصفيحة الجليدية وزيادة رؤوس الهيدروليك على الهامش القاري.
تشير النتائج إلى أنه خلال الجليد، تتطور تدرجات رأس هيدروليكي عمودية قوية، تؤثر بشكل كبير على سرعات تدفق المياه الجوفية وتوزيعات الملوحة. يؤدي سيناريو GLAC، الذي يتميز بتحميل الجليد، إلى تخفيف واسع النطاق في البحر بسبب إعادة شحن المياه الممطرة وتدرجات هيدروليكية بارزة، بينما يعكس سيناريو EUS بشكل أساسي تقلبات مستوى سطح البحر مع تسرب محدود للمياه العذبة. تؤكد الدراسة أيضًا على السلوك غير المتناظر لنظام المياه الجوفية، حيث تختلف معدلات التخفيف خلال الجليد عن معدلات الملوحة خلال ذوبان الجليد، مما يشير إلى أن نظام المياه الجوفية يحتفظ بذاكرة للظروف السابقة. بالإضافة إلى ذلك، يتم فحص تأثير انضغاط الرواسب على النفاذية وتسلل المياه العذبة، مما يكشف أنه بينما قد يعيق الانضغاط تسرب المياه العذبة العميقة، إلا أنه لا يحد بشكل كبير من إعادة شحن المياه العذبة المدفوعة بالجليد تحت الظروف المودلة.
DOI: https://doi.org/10.1007/s10040-025-02994-0
Publication Date: 2026-01-10
Author(s): Sophie ten Hietbrink et al.
Primary Topic: Geology and Paleoclimatology Research
Overview
This study presents the first numerical simulations of offshore freshening of groundwater along the Norwegian continental margin, revealing the significant role of glaciation in this process. The research indicates that during the Late Pleistocene and Holocene, the advance of the Fennoscandian ice sheet onto the continental shelf created vertical hydraulic head gradients that facilitated the infiltration of glacial meltwater into the offshore groundwater system at a rate of 486 mm year$^{-1}$. By 23,000 years ago, areas of fully freshened groundwater expanded to 21.5 km$^{2}$, while regions with chloride concentrations below 400 mM reached 47.4 km$^{2}$, both substantially larger than estimates based solely on eustatic sea-level changes. The study also highlights that glacial retreat led to groundwater discharge at rates up to 246 mm year$^{-1}$, potentially inhibiting seawater recharge in previously glaciated areas.
The findings challenge previous assumptions regarding lateral land-ocean connections as the primary mechanism for offshore freshening, instead demonstrating that vertical emplacement during glaciation is predominant. The persistence of freshened groundwater at depth into interglacial periods underscores the long-term influence of glaciation on subsurface hydrology. The study suggests that the hydraulic head gradients resulting from glacial loading are distinct from those caused by sea-level changes, and current submarine groundwater discharge rates observed at the Lofoten-Vesterålen site cannot be solely explained by these mechanisms. The implications of this research extend globally, as similar processes may occur in other glaciated margins, and the ongoing retreat of ice sheets due to climate change could significantly alter groundwater dynamics and solute fluxes into coastal oceans, potentially releasing methane from subsurface sediments.
Introduction
The introduction of this research paper discusses the significant role of submarine groundwater discharge (SGD) in shaping the offshore physicochemical environment, particularly in coastal biogeochemical cycles and ecosystems. It highlights how glacially enhanced hydraulic gradients in recently glaciated regions can amplify offshore fluid discharge, leading to the exfiltration of groundwater from underlying aquifers. This process is crucial for transporting carbon, nutrients, and other solutes into marine environments, thereby influencing near-seafloor chemistry and nutrient dynamics. The paper emphasizes the need to understand the controls on submarine groundwater flow and freshwater distribution over glacial-interglacial timescales, especially in light of contemporary changes due to glacial retreat and rising sea levels.
The introduction also notes that while the principles of offshore freshening are recognized, the spatial and temporal variability of SGD remains poorly understood. It points out that existing studies have primarily focused on the Atlantic continental shelf off New England and the southeastern Greenland shelf, with a notable gap in numerical simulations for the Norwegian margin, despite evidence of freshening at multiple sites. The paper calls for a coherent understanding of the mechanisms driving freshening in the Norwegian continental margin, which is essential for comprehending groundwater systems influenced by glacial dynamics and for practical applications such as carbon storage in basaltic formations.
Methods
In this study, a hydrogeological model was developed using COMSOL Multiphysics (version 6.1) to simulate groundwater flow and solute transport in the Lofoten-Vesterålen (LV) region. The model employs Darcy’s Law to describe groundwater flow through porous media, represented by the equation \( u = -\nabla p \), where \( u \) is the Darcy velocity field, \( \kappa \) is permeability, \( \mu \) is dynamic viscosity, and \( p \) is pore pressure. The continuity equation incorporates compressibility \( S \), fluid density \( \rho \), and a mass source term \( Q_m \), which was set to zero due to the absence of additional fluid sources. The solute transport, specifically chloride, is modeled by coupling the Darcy velocity field to a transport equation that accounts for solute accumulation, adsorption, hydrodynamic dispersion, and convection.
The model geometry is based on seismic interpretations and spans a horizontal extent of 160 km with a vertical extent of approximately 6 km, divided into eight stratigraphic units. Sensitivity analyses were conducted to assess the impact of varying parameters such as porosity, permeability, compressibility, and dispersivity, with realistic ranges defined for the Lofoten-Vesterålen continental margin. The chosen values for porosity (0.15 for the basement unit and 0.2 for non-basement units) and compressibility (5 × 10⁻⁶ Pa⁻¹) align with existing literature. The model assumes a constant stratigraphic geometry over a timescale of 120,000 years and utilizes a free triangular mesh for computational efficiency, with an adaptive time-stepping approach to resolve temporal changes effectively.
Results
In this study, two modeled scenarios were analyzed to assess the impact of glacial loading on relative sea-level rise. The first scenario, referred to as GLAC, incorporates glacial loading effects, represented by the purple curve in Figure 4b. The second scenario, labeled EUS, follows the eustatic sea-level curve without any glacial loading, depicted by the green curve in Figure 4b. By comparing these two scenarios, the research highlights the significant contribution of glacial loading to changes in relative sea levels.
Discussion
In this section, the authors discuss the methodologies and findings related to hydraulic head boundary conditions and groundwater dynamics influenced by glacial loading in the Eurasian ice sheet region. Ice sheet thickness data were derived from simulations of glacial cycles, specifically focusing on the Last Glacial Maximum (LGM) and its impact on hydraulic head and relative sea level (RSL). The study highlights the significant divergence between RSL and eustatic sea-level curves during glaciation, driven by isostatic depression, and outlines the methods used to extrapolate RSL data beyond the available 22,000 years. The hydraulic head induced by the ice mass was calculated, and a Dirichlet boundary condition was applied to model the hydraulic head dynamics, linking ice sheet expansion to increased hydraulic heads on the continental margin.
The results indicate that during glaciation, strong vertical hydraulic head gradients develop, significantly influencing groundwater flow velocities and salinity distributions. The GLAC scenario, characterized by glacial loading, leads to extensive offshore freshening due to meteoric water recharge and pronounced hydraulic gradients, while the EUS scenario primarily reflects sea-level fluctuations with limited freshwater infiltration. The study also emphasizes the asymmetrical behavior of the groundwater system, where freshening rates during glaciation differ from salinization rates during deglaciation, suggesting that the groundwater system retains a memory of past conditions. Additionally, the influence of sediment compaction on permeability and freshwater infiltration is examined, revealing that while compaction may inhibit deep freshwater infiltration, it does not substantially limit glacially driven freshwater recharge under the conditions modeled.