علوم الأرض والفضاء الجامعية — السنة الثالثة · 本科三年级
6الحمل في الوشاح والقوى المحرّكة للصفائح
يرتفع ساحل شمال السويد من البحر بنحو سنتيمتر واحد في السنة، بعد عشرة آلاف سنة من ذوبان الغطاء الجليدي الذي كان يغطيه؛ وترتفع سواحل خليج هدسون بالسرعة نفسها. لا يدفعها إلى أعلى سوى الوشاح، يتدفّق ببطء عائدًا أسفل الموضع الذي ضغطه الجليد نحو الأسفل. ويدلّ معدّل ارتداد الأرض على لزوجة الوشاح، وتحدّد اللزوجة نفسها سرعة حمل الوشاح وحركة الصفائح. يتساءل هذا الفصل لماذا يحمل الوشاح الصلب أصلًا، وبأيّ سرعة يجب أن يتدفّق لينقل حرارة الأرض، وما القوى التي تحرّك الصفائح، وماذا يقول ارتداد الأرض عن لزوجة الوشاح.
ما تعرفه من قبل
اللزوجة، أي الجريان الزاحف لمائع لزج (معادلات ستوكس ∇p=η∇2u، ∇⋅u=0)، وقانون ستوكس للكرة، وعدد رينولدز والطفو واردة في المجلّدات الجامعية لكتب الفيزياء. وقد ورد الحمل في الوشاح والارتداد الأيزوستاتي في مجلّد المرحلة الثانوية (الصف 12)، وقوة دفع الظهرة وشدّ اللوح المندس في جزئه الخاص بالصف 11. الطبقة الحدّية الحرارية والغلاف الصخري الحراري والمنحنى الأدياباتي للوشاح موجودة في مجلد السنة الأولى، في فصلها عن ميزانية الحرارة الأرضية؛ والبقع الساخنة، وأعمدة الوشاح والأقاليم البركانية الكبرى في فصلها عن البقع الساخنة. ويرد نموذج تبريد نصف الفضاء في الفصل 5، وانثناء الغلاف الصخري في الفصل 3، والتصوير المقطعي الزلزالي في الفصل 1.
6.1 لماذا يحدث الحمل في الوشاح
الوشاح صلب، إذ ينقل الموجة S. غير أنّه على مدى ملايين السنين يزحف زحفًا كمائع بالغ اللزوجة، وطبقةٌ من مائع تُسخَّن من الأسفل وتُبرَّد من الأعلى قد تبقى ساكنة، وتعبرها الحرارة بالتوصيل، أو تنقلب. كتلةٌ من مائع ساخن حجمها d، أخفّ من محيطها بمقدار ραΔT، ترتفع بسرعة من رتبة ρgαΔTd2/η؛ وتفقد حرارتها بالانتشار خلال زمن من رتبة d2/κ. فإن استطاعت أن ترتفع عبر الطبقة قبل أن تبرد، تنقلب الطبقة.
تعريف 6.1(عدد رايلي)
بالنسبة لطبقة مائعة سُمكها d، وكثافتها ρ، ومعامل التمدّد الحراري α، والانتشارية الحرارية κ واللزوجة η، تحت تأثير الجاذبية g، وبفرق درجة الحرارة ΔT بين قاعدتها وقمّتها، فإنّ عدد رايلي هو
Ra=κηρgαΔTd3,
أي نسبة زمن الانتشار d2/κ إلى زمن الارتفاع بالطفو η/(ρgαΔTd).
تعريف 6.2(عدد رايلي الحرج)
عدد رايلي الحرجRac لطبقة ما هي قيمة عدد رايلي التي فوقها يكبر الاضطراب الصغير في حالة السكون فيتحوّل إلى حمل، وأدنى منها يتلاشى.
قضية 6.3(بداية الحمل)
لطبقة خلا سطحاها العلوي والسفلي من إجهاد القصّ، وثُبِّتت درجتا حرارتهما ثابتتين، ينمو اضطراب ذو عدد موجي أفقي k (بوحدات 1/d) عندما Ra>(π2+k2)3/k2. وتبلغ أدنى قيمة للطرف الأيمن، عند k=π/2، ما يعطي
Rac=427π4≈657.5,
وأولى الخلايا التي تظهر بعرض 22d. وبين حدود صلبة، Rac≈1708.
برهان.نقبله في هذا المستوى.∎
تحليل الاستقرار الخطّي الذي يقوم عليه، وهو اضطراب لمعادلتي الحرارة وستوكس حول الحالة التوصيلية، وهو وارد في معالجة كتب الفيزياء للحمل؛ ويتحقّق التمرين 11 من الحدّ الأدنى.
يُظهر الشكل أدناه الطبقة التي يحدث فيها الحمل. وقيم الوشاح، مع عمق الحدّ بين اللبّ والوشاح، d=2891km، والكثافة المتوسطة للوشاح، 4476kgm−3، وg بين 9.8 و10.7ms−2، وكلّها مأخوذة من نموذج الأرض الزلزالي PREM؛ وα=3.1×10−5K−1 وΔT=1450K من نموذج الصفيحة في الفصل 5؛ وκ=10−6m2s−1؛ وη بين 1021 و1022Pas (انظر القسم 6.5)، فتعطي Ra بين 4.9×106 و4.9×107، نحو سبعة آلاف إلى سبعين ألف مرّة قيمة الحدّ الحرج. فلا يستطيع الوشاح إلّا أن يحمل.
الحمل في طبقة تُسخَّن من الأسفل وتُبرَّد من أعلى. تعبر الحرارة الطبقات الحدّية الرقيقة، ذات السماكة δ، بالتوصيل؛ في ما بينهما يحمل الداخل جيّد الخلط الحرارة بالجريان. تنفصل صفائح باردة عن الطبقة العليا وتغوص؛ وتصعد صفائح حارّة من السفلى.
ويُسخَّن الوشاح كذلك من داخله، بتحلّل العناصر المشعّة، ومن الأسفل، بفعل اللبّ. الطبقة المُسخَّنة من داخلها وحدها لا تملك طبقة حدّية حارّة عند قاعدتها: يُدفع جريانها بصفائح باردة تهبط من الأعلى، بينما يكون الجريان الصاعد عودةً واسعةً سلبية. والألواح المندسة من هذه الصفائح، والأعمدة تحتاج إلى الحرارة القادمة من اللبّ.
6.2 الطبقات الحدّية وقوّة الحمل
تعريف 6.4(عدد نوسلت)
يُعرَّف عدد نوسلت لطبقة يجري فيها الحمل بأنّه نسبة التدفق الحراري عبرها إلى التدفق الذي يحمله التوصيل وحده عبر الطبقة نفسها، Nu=qd/(kΔT).
مبرهنة 6.5(قياس الحمل بالطبقة الحدّية)
في طبقة يجري فيها الحمل عند Ra≫Rac، إذا نمت كلّ طبقة حدّية حرارية حتى يبلغ عدد رايلي الخاص بها، المبنيّ على سماكتها δ وفرق درجة الحرارة ΔT/2 عبرها، يبلغ Rac، فإنّ
حيث u هو مقياس سرعة الجريان، المعرَّف باشتراط أن تنمو الطبقة إلى δ بينما يعبر المائع الخليّة. وهذا هو قياس الحمل بالطبقة الحدّية.
برهان. يكون الداخل جيّد الخلط عند متوسط درجة الحرارة، فتحمل كلّ طبقة حدّية نصف ΔT. وعدد رايلي لها هو ρgα(ΔT/2)δ3/(κη)=Ra(δ/d)3/2؛ وبجعله مساويًا Rac فنحصل على δ. ويُوصَّل التدفق الحراري عبر الطبقة، q=k(ΔT/2)/δ، فيكون Nu=qd/(kΔT)=d/(2δ). تنمو طبقة حدّية بالانتشار حتى السماكة δ خلال زمن δ2/κ؛ ويعبر المائع عند القمة الخليّة، التي عرضها من رتبة d، في زمن d/u. وبمساواة الزمنين، u=κd/δ2=(κ/d)(d/δ)2. ∎
لا يزداد التدفق الحراري إلّا بمقدار الجذر التكعيبي للعدد Ra: فلزوجة أدنى بمائة مرّة تُعطي حرارة أكثر بأقلّ من خمس مرّات. ومع عدد رايلي للوشاح Ra=4.9×107، يُعطي القياس δ≈87km، Nu≈17، وتدفّقًا حراريًا قدره NukΔT/d≈26mWm−2 مع k=3.1Wm−1K−1، وu≈1.2cm/yr (الشكل أدناه). وتساوي الطبقة الحدّية سماكة الصفيحة المحيطية في الفصل 5، 95km؛ أمّا الجذر التربيعي المتوسّط لسرعات الصفائح اليوم، فيتراوح من 1.1cm/yr (أمريكا الجنوبية) إلى 7.5cm/yr (كوكوس)، و4.4cm/yr للغلاف الصخري كلّه. ويبلغ متوسط التدفق الحراري للأرض 92mWm−2، ثلاثة إلى أربعة أضعاف القياس: فالحرارة المتولّدة داخل الوشاح، التي تُغفلها طبقة مُسخَّنة من الأسفل، والثوابت من رتبة الواحد التي يُغفلها القياس كلاهما يسهم في الفرق. وأمّا قانون بمدخل واحد، فرتب المقادير صحيحة: فالصفائح هي الطبقة الحدّية العليا لوشاح ذي Ra∼107.
قياس الحمل. (أ) دون Rac تنقل الطبقة الحرارة بالتوصيل (Nu=1)؛ وأعلى منه بكثير، فيُرسم قانون الطبقة الحدّية، من Ra=16Rac≈1.1×104، حيث يُعطي Nu=1 (وبين الحالتين يبدأ الحمل تدريجيًا). (ب) مقياس السرعة لعمق الوشاح وκ=10−6m2s−1، مع مدى سرعات الصفائح اليوم (بخطّ شرطي-نقطي). وتحدّ الخطوط المتقطّعة عدد رايلي للوشاح عند لزوجتين 1022 و1021Pas.
6.3 القوى المؤثّرة في الصفيحة
تتحرّك الصفيحة بسرعة ثابتة، فتتوازن القوى المؤثّرة فيها. وتدفعها ثلاث قوى أو تقاومها من حوافّها وقاعدتها (الشكل أدناه). والقوّة التي يؤثّر بها اللوح المندس الغائص في صفيحته، أي شدّ اللوح المندس، تنشأ عن الكثافة الزائدة للوح مندس بارد في الوشاح الأدفأ. أمّا قوّة دفع الظهرة، وهي دفع الظهرة المرتفعة الساخنة للصفيحة الأبرد المجاورة لها، فتُستنتج من نموذج تبريد نصف الفضاء في الفصل 5.
القوى المؤثّرة في صفيحة محيطية: دفع الظهرة وشدّ اللوح المندس تدفعها (بالأحمر)، ويقاومها السحبُ عند قاعدتها حيث تتحرّك الصفيحةُ أسرعَ من الوشاح الذي تحتها، وكذلك مقاومةُ الوشاح اللزجةُ لغوص اللوح المندس (بالأزرق).
قضية 6.6(قوة دفع الظهرة)
في نموذج تبريد نصف الفضاء، مع درجة حرارة الوشاح Tm ودرجة حرارة السطح T0 والتعويض الأيزوستاتي، فإنّ القوّة الأفقية لكلّ وحدة طولٍ من الظهرة التي يؤثّر بها العمودُ عند الظهرة في الصفيحة ذات العمر t، وهي قوة دفع الظهرة، وتُعطى بالعلاقة:
برهان. ليكن y العمقَ تحت قاع البحر عند الظهرة، وw الهبوطَ لقاع عمره t، وθ(z)=Tm−T(z) عجزَ درجة الحرارة عند العمق z=y−w في صفيحته، وΘ(z)=∫0zθ. والعمودُ عند الظهرة يحوي وشاحًا كثافتُه ρm ابتداءً من y=0؛ أمّا العمود القديم فيحوي ماءً حتى العمق w، ثم صفيحةً كثافتُها ρm(1+αθ). وتختلف الضغوط الليثوستاتية بمقدار Pr−Pt=(ρm−ρw)gy عندما y<w، وأسفل ذلك،
Pr−Pt=(ρm−ρw)gw−ρmgαΘ(z)=ρmgα(Θ(∞)−Θ(z)),
بالتوازن الأيزوستاتي، (ρm−ρw)w=ρmαΘ(∞). والقوّةُ هي تكاملُ فرق الضغط على العمق؛ وبمكاملة الجزء الثاني بالتجزئة،
Frp=21(ρm−ρw)gw2+ρmgα∫0∞zθ(z)dz.
أمّا نصف الفضاء، فإنّ θ=(Tm−T0)erfc(z/2κt) و∫0∞xerfcxdx=1/4، فيكون التكامل يساوي (Tm−T0)κt؛ وw=2ρmα(Tm−T0)κt/π/(ρm−ρw). وبجمع الحدّين نحصل على الصيغة. ∎
يتزايد دفع الظهرة خطيًّا مع العمر: فهو ليس دفعًا عند المحور، بل تدرّجًا في الضغط موزّعًا على الصفيحة المتبرّدة كلّها. وبالقيم الواردة في الفصل 5، يبلغ الحدّ بين القوسين 1.04 وFrp=3.8×1012Nm−1 عند 100Myr. وفي نموذج الصفيحة يتوقّف النموّ مع العمق، عند نحو 2.3×1012Nm−1 (الشكل أدناه). وشدّ اللوح المندس أكبر: لوحٌ مندس سُمكُه 95km، وهو أبرد في المتوسّط بمقدار (Tm−T0)/2 من الوشاح، ومن ثَمّ أكثف بمقدار ρmα(Tm−T0)/2≈74kgm−3، ويمتدّ حتى 670km، ويزن 4.6×1013Nm−1 أكثر من الوشاح الذي يزيحه، أي اثني عشر ضعفًا لقوّة دفع الظهرة. ويُنفَق جزءٌ كبير من هذا الوزن على الوشاح اللزج الذي يقاوم غوص اللوح المندس، ولا يصل إلى الصفيحة إلّا جزءٌ منه.
قوّة دفع الظهرة بدلالة عمر الصفيحة، استنادًا إلى نموذجَي التبريد في الفصل 5، بالتكامل العددي لفرق الضغط (منحنى نصف الفضاء هو الصيغة المغلقة للقضية).
تعريف 6.7(السحب القاعدي)
السحب القاعدي هو قوّة القصّ لكلّ وحدة مساحة التي يؤثّر بها الوشاح تحت الصفيحة في قاعدتها. وبالنسبة إلى صفيحة تتحرّك بسرعة u فوق طبقة سُمكها h ولزوجتها η وقاعدتها ساكنة، فإنّها تساوي ηu/h، وهي معاكسة للحركة؛ وتدفع الصفيحة حيث يتحرّك الوشاح الذي تحتها أسرع من الصفيحة.
وأيّ القوى هي المهمّة يمكن قراءته من الصفائح نفسها. ويبيّن توازن العزوم على كلّ صفيحة، والمطابَق مع حركات الصفائح، أنّ القوى المؤثّرة في اللوح المندس الغائص هي التي تتحكّم في سرعة الصفائح المحيطية، وهي أقوى بنحو عشرة أضعاف من أيّ قوّة أخرى، وأنّ المقاومة عند القاعدة أكبر تحت القارّات منها تحت المحيطات. وتتّفق السرعات: فأسرع الصفائح، كوكوس بسرعة 75mm/yr، وناسكا بسرعة 70، والمحيط الهادئ بسرعة 66mm/yr، صفائحُ محيطيةٌ تغوص حوافّها في خنادق طويلة؛ وأبطؤها هي التي تحمل قارّات كبيرة: أوراسيا بسرعة 23، وأمريكا الشمالية بسرعة 20، وأمريكا الجنوبية بسرعة 11mm/yr.
6.4 أعمدة الوشاح
تنتظم الطبقات الحارّة الصاعدة من قاعدة الوشاح في أعمدة، لها رأس كبير وذيل رفيع، وهي أعمدة الوشاح المذكورة في مجلد السنة الأولى (الشكل أدناه). يقيس الانتفاخُ فوق البقعة الساخنة، وهو ارتفاعٌ واسع لقاع البحر حول البراكين وفي اتجاه حركة الصفيحة بعدها، مقدارَ الطفو الذي يجلبه عمودُها.
تعريف 6.8(دفق طفو العمود)
دفق طفو العمودB هو العجز في الكتلة، بالنسبة إلى الوشاح المحيط، الذي يحمله العمود إلى أعلى في وحدة الزمن؛ ولعمود ذي دفق حجمي Qv وزيادةٍ في درجة الحرارة ΔT، B=ραΔTQv.
قضية 6.9(دفق الطفو من انتفاخ)
إذا تحرّكت صفيحة بالسرعة u فوق عمود يحمل انتفاخًا مساحةُ مقطعه العرضي A عبر المسار، ويُسنَد أيزوستاتيًّا بمادة العمود التي تحته، فإنّ دفق الطفو للعمود هو
B=(ρm−ρw)Au.
برهان. بالتوازن الأيزوستاتي، تُسنَد شريحةٌ من الانتفاخ طولها وحدة واحدة على طول المسار، ترفع قاع البحر على مساحة المقطع العرضي A وتستبدل الماء بالصخر، بعجزٍ في الكتلة مقداره (ρm−ρw)A في العمود الذي تحتها. تحمل الصفيحة الشريحة بعيدًا بالسرعة u، ولذلك يجب في الحالة المستقرّة أن يوفّر العمود العجز لطول u من الانتفاخ الجديد في وحدة الزمن. ∎
استنادًا إلى المقاطع العرضية للانتفاخات، تتراوح دفوق الطفو للبقع الساخنة ضمن عامل قدره عشرون: هاواي، وهي الأكبر، بمقدار 8.7Mg/s؛ وآيسلندا بمقدار 1.4Mg/s؛ و37 بقعة ساخنة معًا نحو 50Mg/s. وكلّ كيلوغرام من الطفو يعادل cp/α جول من الحرارة الزائدة، فيكون التدفق الحراري cpB/α؛ مع cp=k/(ρmκ)=1170Jkg−1K−1 وα=3.1×10−5K−1، تحمل 37 بقعة ساخنة نحو 1.9TW، أي 3.7mWm−2 على سطح الأرض البالغ 5.1×1014m2، وهو قريب من 4mWm−2 المنشورة مع الدفوق: بضعة في المئة من 47TW التي تفقدها الأرض. الأعمدة شريك ثانوي في ميزانية الحرارة؛ أما الباقي فتحمله الصفائح.
عمود فتيّ (يسارًا) يصعد برأس كبير تغذّيه ذيل رفيع؛ وعندما ينتشر الرأس تحت الغلاف الصخري، يرتفع السطح، وقد تتفجّر البازلتات الفيضية. أما العمود الأقدم (يمينًا) فهو ذيل تحت صفيحة متحرّكة يغذّي انتفاخًا ومسار بقعة ساخنة. ليس بمقياس الرسم.
يمثّل رأس العمود كرةً من وشاح طافٍ صاعدًا عبر الوشاح اللزج: يعطي قانونُ ستوكس المذكور في كتب الفيزياء، v=2Δρga2/(9η)، سرعتَه بالرتبة العددية، وهي كسرٌ من المتر في السنة لرأس يبلغ عرضه بضع مئات الكيلومترات (التمرين 8)، فيعبر الوشاح في ملايين السنين.
6.5 التكيّف الأيزوستاتي الجليدي ولزوجة الوشاح
تعريف 6.10(التكيّف الأيزوستاتي الجليدي، لزوجة الوشاح)
التكيّف الأيزوستاتي الجليدي (GIA) هو التشوّه البطيء للأرض الصلبة، وتغيّرُ مجالها الثقالي، ومستوى سطح البحر، أثناء عودتها نحو التوازن الأيزوستاتي بعد أن ينمو الغطاء الجليدي أو يذوب. أما لزوجة الوشاح فهي اللزوجة الفعّالة η التي تربط الزحف البطيء المستقرّ للوشاح بالإجهادات التي تحرّكه، على مدى آلاف السنين وأكثر.
عندما يذوب الغطاء الجليدي، يكون الانخفاض الذي يتركه خارج التوازن الأيزوستاتي: فيكون ضغط الوشاح تحته أقلّ منه حوله، فيتدفّق الوشاح إليه ليرفع القاع من جديد. اليوم يرتفع ساحل خليج بوثنيا بمقدار يصل إلى 10.3mm/yr قرب أومو، ويرتفع خليج هدسون بنحو 10mm/yr، حيث كان الجليد أسمك ما يكون، في حين تهبط الأرض جنوب البحيرات العظمى بمقدار من 1 إلى 2mm/yr، إذ يهبط الانتفاخ الذي رفعه الجليد حولها. وعلى السواحل تقف خطوط الشاطئ القديمة الآن فوق مستوى سطح البحر، وهي الشواطئ المرتفعة الموضّحة في صندوق الحقل أدناه.
مبرهنة 6.11(استرخاء نصف فضاء لزج)
ليكن الانخفاض w(x,t)=wm(t)cos(2πx/λ) في سطح نصف فضاء ذي لزوجة منتظمة η وكثافة ρ، يُترك دون حمولة، فتتناقص السعة بحسب wm(t)=wm(0)e−t/τ مع
τ=ρgλ4πη.
برهان. نأخذ z موجّهًا نحو الأسفل، وk=2π/λ، ونبحث عن تدفّق ستوكس يتناقص مع العمق، بسرعة رأسية (نحو الأسفل) v=(C+Dkz)e−kzcoskx. ويُعطي عدم الانضغاطية السرعة الأفقية u=−(D−C−Dkz)e−kzsinkx. تحقّق هذه الدوال ∇4=0، كما تفعل كلّ تدفّقات ستوكس في بُعدين. والسطح الحرّ لا يحمل إجهاد قصّ: ∂u/∂z+∂v/∂x=2k(D−C)sinkx=0 عند z=0، فيكون D=C، ومنه
ويُستنتج الضغط من ∂p/∂x=η∇2u=−2ηCk2e−kzsinkx (ويتحقّق ∂p/∂z=η∇2v). وعند السطح ∂v/∂z=0، فإنّ الإجهاد العمودي هناك هو −p. وتحت انخفاضٍ w، يُخفّض وزنُ الصخر المُزاح الضغطَ بمقدار ρgw: p(0)=−ρgw، فيكون 2ηCk=−ρgwm. ويتحرّك السطح مع المائع: dwm/dt=v(0)=C=−ρgwm/(2ηk)، الذي يتناقص حلّه بالزمن τ=2ηk/(ρg)=4πη/(ρgλ). ∎
تعريف 6.12(زمن الارتخاء الأيزوستاتي)
زمن الارتخاء الأيزوستاتيτ لحمولة ذات مقياس أفقي λ هو الزمن الذي يتناقص فيه انحرافها عن التوازن الأيزوستاتي بالعامل e.
يرتخي نصف الفضاء المنتظم الانخفاضات الواسعة أسرع من الضيّقة، والتدفّق تحت حمولة طول موجتها λ يتناقص مع العمق وفق e−kz، على عمق λ/2π: يُحرّك الغطاء الجليدي الكبير الوشاح عميقًا، بينما لا تُحرّك بحيرة صغيرة إلّا قمّته. بالنسبة إلى انخفاض عرضه 3000km وη=1021Pas، مع ρ=3300kgm−3، τ=4.1kyr (الشكل أدناه): بعد عشرة آلاف سنة، يبقى e−2.4=9% من الانخفاض، ولأنّ w˙=−w/τ، فإنّ رفعًا بمقدار 10mm/yr يعني أنّ نحو 40m من الارتفاع ما زال قادمًا.
ارتخاء نصف فضاء لزج منتظم (ρ=3300kgm−3، g=9.8ms−2) لثلاث قيم للزوجة. (أ) نسبة الانخفاض المتبقية مع الزمن. (ب) زمن الارتخاء 4πη/(ρgλ)، أطول للحمولات الضيّقة.
لا يكون الوشاح الحقيقي منتظمًا، ويُقرأ ارتداده بنماذج تُقسّمه إلى طبقات وتُضيف غلافًا صخريًا مرنًا يحمل الحمولات الضيّقة بالانثناء. وقد أعطى التقدير الأول، المستمدّ من ارتداد فينوسكانديا، 1021Pas، وتناقص الرفع عند نهر أنغرمان في السويد وعند خليج ريتشموند على خليج هدسون يُعطي متوسطًا من 0.65 إلى 1.1×1021Pas، ابتداءً من قاعدة الغلاف الصخري حتى عمق نحو 1400km، ويكون الوشاح العلوي أقلّ لزوجة من السفلي بنحو عشرة أضعاف. وتُعطي انعكاسات معدّلات الرفع بنظام GNSS في فينوسكانديا لزوجات الوشاح العلوي من 5×1020 إلى 1021Pas في دراسة، ومن 3.4×1020 إلى 5.0×1020Pas في دراسة أخرى، ولزوجات الوشاح السفلي من 5×1021 إلى 5×1022Pas، وسماكة الغلاف الصخري من 93 إلى 110km. تختلف التقديرات بعوامل من اثنين إلى ثلاثة، لكنها تتفق على رتبة القدر وعلى الزيادة مع العمق.
تعريف 6.13(زمن ماكسويل)
زمن ماكسويل لمادّة لزوجتها η ومعامل القصّ μ هو τM=η/μ: الحمولات المطبّقة لمدّة أقصر بكثير من τM تلقى جسمًا مرنًا، والحمولات المطبّقة لمدّة أطول بكثير منه تلقى مائعًا لزجًا.
بأخذ μ=ρvS2=7.7×1010Pa على عمق 300km في PREM وη=1021Pas، τM=1.3×1010s، أي نحو أربعمئة سنة. أمّا الموجة الزلزالية، التي تدوم ثوانٍ، فتعبر وشاحًا مرنًا؛ والعصر الجليدي، والحمل، يُشوّهان مائعًا. والصخر نفسه يكون الاثنين، بحسب المدّة التي يُدفع فيها.
في الميدان— شواطئ مرتفعة
الشواطئ المرتفعة في بالي هيلين، مالين هيد، شمال أيرلندا: حواجز من الحصى لخطوط شواطئ قديمة رُفعت فوق البحر مع ارتفاع الأرض بعد ذوبان غطائها الجليدي. الصورة: A. F. Borchert, CC BY-SA 4.0.
على سواحل الأراضي التي غطّاها الجليد سابقًا، يجد الجيومورفولوجي عمل البحر القديم مُخلَّفًا في الداخل: حواجز من الحصى المستدير، وتجاويف، ومنصّات مسطّحة نحتتها الأمواج، وتقف الآن فوق مدى أعلى المدّ الحالي، وتتّخذ أحيانًا شكل درج من عدّة مستويات. وكان كلّ مستوى شاطئًا حين كان البحر واقفًا عنده؛ وقد ارتفعت الأرض منذ ذلك الحين. وتُؤرَّخ الأصداف وقطع الخشب المنجرف في الشواطئ القديمة بالكربون المشعّ، وتُعطي ارتفاعات المستويات المؤرَّخة، مُسوّاةً بعلامة مرجعية، الرفعَ بدلالة الزمن. وعند تتبّع شاطئ قديم على امتداد ساحل، يُلاحظ أنّ الشاطئ القديم مائل، وأعلاه حيث كان الجليد أسمك: فالميل يرسم شكل الارتداد.
6.6 الحمل الوشاحي الكامل وبنيته العميقة
هل يتحرّك الوشاح حملًا في طبقة واحدة، أم في طبقتين، فوق الانقطاع عند 670km وتحته؟ يُجيب التصوير المقطعي (انظر الفصل 1) جزئيًا عن ذلك: تتسطّح بعض اللوحات المندسّة فوق الانقطاع، بينما تغوص أخرى عبره إلى الوشاح السفلي، فتعبر الكتلة منه، ولا تكون الطبقتان معزولتين إحداهما عن الأخرى. يجعل ازدياد اللزوجة مع العمق الوشاح السفلي يتدفّق ببطء أكبر، وتتباطأ اللوحات المندسّة وتتثخّن عند دخولها إليه.
عند قاع الوشاح، يُظهر التصوير المقطعي منطقتين واسعتين ذواتي سرعة منخفضة لموجات S تحت أفريقيا وتحت المحيط الهادئ، وهما الأقاليم الكبيرة ذات السرعة المنخفضة للقصّ (LLSVPs). يُرجَّح أنّها متمايزة كيميائيًا عن الوشاح الذي يعلوها وأكثف منه، فهي أكوام لا أعمدة حارّة. أمّا الأقاليم البركانية النارية الكبيرة والكيمبرلايتات في المئات القليلة الأخيرة من ملايين السنين اندفعت في معظمها فوق حوافّها، كأنّ الأعمدة صعدت من أطراف الأكوام، حيث تُجرف الطبقة الحدّية الحارّة عند قاع الوشاح باتجاهها.
6.7 التمارين
تمرين 6.1★
احسب عدد رايلي لطبقة من الوشاح العلوي سمكها 670km مع ρ=3500kgm−3، وg=9.9ms−2، وα=3×10−5K−1، وΔT=1000K، وκ=10−6m2s−1، وη=5×1020Pas. كم مرّة يتجاوز عدد رايلي الحرج، بالنسبة إلى حدود انزلاق حرّ؟
حل
حل التمرين 6.1.
Ra=3500×9.9×3×10−5×1000×(6.7×105)3/(10−6×5×1020)=6.25×105، نحو 950 ضعف عدد رايلي الحرج 657.5: تحمل الطبقة حملًا قويًّا.
تمرين 6.2★
احسب زمن ماكسويل لوشاح لزوجته 1021Pas ومعامل القصّ 7.7×1010Pa، بالسنوات. هل يستجيب الوشاح مرنًا أم لزجًا للمدّ الأرضي نصف اليومي، ولذوبان الغطاء الجليدي، ولغوص لوح مندسّ؟
حل
حل التمرين 6.2.
τM=1021/7.7×1010=1.30×1010s، نحو 412yr. أمّا المدّ الذي يدوم ساعات فيلتقي بوشاح مرن؛ وأمّا ذوبان الغطاء الجليدي على مدى آلاف السنين، وهبوط اللوح المندس على مدى ملايين السنين، فيلتقيان بوشاح لزج.
تمرين 6.3★
مع η=1021Pas، وρ=3300kgm−3، وg=9.8ms−2، احسب زمن الارتخاء الأيزوستاتي لنصف فضاء منتظم لأحمال طول موجتها 1000 و3000km. أيّهما يتعافى أسرع؟
حل
حل التمرين 6.3.
τ=4π×1021/(3300×9.8λ): 3.89×1011s = 12.3kyr من أجل 1000km، و4.1kyr من أجل 3000km. أمّا الحمل الأوسع فيتعافى أسرع.
تمرين 6.4★
مع g=9.8ms−2، وρm=3300kgm−3، وρw=1030kgm−3، وα=3.1×10−5K−1، وTm−T0=1450K، وκ=0.8×10−6m2s−1، و1Myr=3.156×1013s، احسب قوة دفع الظهرة عند 50Myr. كم مرّة تكون أصغر من شدّ لوح مندسّ قدره 4.6×1013Nm−1؟
حل
حل التمرين 6.4.
gρmα(Tm−T0)κt=9.8×3300×3.1×10−5×1450×0.8×10−6×50×3.156×1013=1.835×1012Nm−1؛ والقوس يساوي 1+2×3300×3.1×10−5×1450/(π×2270)=1.0416، لذا Frp=1.91×1012Nm−1، أصغر بنحو 24 مرّة من شدّ اللوح المندس.
تمرين 6.5★★
كتلة من الوشاح عرضها 10km، أسخن من محيطها بمقدار 300K، تصعد بالسرعة ρgαΔTd2/η (خذ ρ=3300kgm−3، وg=9.8ms−2، وα=3×10−5K−1، η=1021Pas، وκ=10−6m2s−1). (أ) احسب زمن صعودها عبر حجمها نفسه وزمن انتشارها d2/κ، وبيّن أنّ نسبتهما هي عدد رايلي لها. (ب) هل تصعد أولًا أم تبرد أولًا؟ (ج) كم ينبغي أن يكون حجمها حتى يبلغ عدد رايلي لها 657.5؟
حل
حل التمرين 6.5.
(أ) v=3300×9.8×3×10−5×300×(104)2/1021=2.91×10−11ms−1؛ زمن الصعود عبر 10km هو d/v=3.44×1014s، 10.9Myr، وزمن الانتشار d2/κ=1014s، 3.2Myr. ونسبتهما هي (d2/κ)/(d/v)=vd/κ=ρgαΔTd3/(κη)=Ra=0.29. (ب) يبرد أوّلًا: الكتل الصغيرة لا يحدث فيها الحمل. (ج) Ra∝d3، لذا d=10×(657.5/0.291)1/3=131km.
تمرين 6.6★★
لطبقة عمقها 2891km عند Ra=107، مع κ=10−6m2s−1، وk=3.1Wm−1K−1، وΔT=1450K، وRac=657.5، احسب سمك الطبقة الحدّية، وعدد نوسلت، ومقياس السرعة بوحدة cm/yr، وتدفّق الحرارة التدفق الحراري.
لدينا انتفاخ افتراضي للبقعة الساخنة مساحة مقطعه العرضي 8×108m2 عبر المسار، والصفيحة فوقه تتحرّك بسرعة 9cm/yr. مع ρm−ρw=2270kgm−3، احسب دفق طفوه، ودفق حرارته مع cp=1170Jkg−1K−1 وα=3.1×10−5K−1. ما الكسر الذي يمثّله من 46.7TW للأرض؟
حل
حل التمرين 6.7.
u=0.09/3.156×107=2.85×10−9ms−1 وB=2270×8×108×2.85×10−9=5.18×103kgs−1=5.2Mg/s، نحو ثلاثة أخماس تدفّق هاواي. ويساوي التدفّق الحراري 1170×5178/3.1×10−5=1.95×1011W، 0.2TW، أو 0.42% من الفقد الحراري للأرض.
تمرين 6.8★★
رأس عمود نصف قطره 400km، أخفّ من الوشاح بمقدار 30kgm−3، يصعد عبر وشاح لزوجته 1021Pas. باستخدام قانون ستوكس، احسب سرعته بوحدة m/yr، والزمن الذي يستغرقه لصعود 2000km.
حل
حل التمرين 6.8.
v=2×30×9.8×(4×105)2/(9×1021)=1.05×10−8ms−1=0.33m/yr؛ ويستغرق قطع 2000km زمنًا 2×106/1.05×10−8=1.9×1014s، نحو 6Myr.
تمرين 6.9★★
انخفضت منطقة بمقدار 300m عندما ذاب جليدها، قبل 10kyr، وتتعافى بثابت زمني τ=4.4kyr. احسب الانخفاض المتبقّي اليوم ومعدّل الرفع الحالي. قارن بالرفع 10.3mm/yr المقيس قرب أوميو.
حل
حل التمرين 6.9.
w=300e−10/4.4=30.9m وw˙=w/τ=30.9/4400=7.0mm/yr، وهو ثلثا المعدّل المقيس قرب أوميو. ونصف الفضاء المتجانس بهذا τ يقع في النطاق الصحيح، لكنّ الارتداد الحقيقي أسرع، أو أنّ انخفاضه أعمق.
تمرين 6.10★★★
(أ) من تعبير Frp في برهان القضية 6.6، بيّن أنّ الحد الأول، 21(ρm−ρw)gw2، هو الحد الثاني داخل القوس. (ب) بيّن أنّ ∫0∞xerfcxdx=1/4، بالتكامل بالتجزئة. (ج) باستخدام قيم التمرين 4، احسب الحدّين عند 100Myr.
(ب) بالتكامل بالتجزئة، ∫0∞xerfcxdx=[2x2erfcx]0∞+∫0∞2x2π2e−x2dx=0+π1⋅4π=41. (ج) عند 100Myr، ρmgα∫zθ=3.67×1012Nm−1 وبوضع w=3705m، 21(ρm−ρw)gw2=1.53×1011Nm−1: حدّ الماء لا يمثّل سوى 4% من القوّة.
تمرين 6.11★★★
بيّن أنّ f(k)=(π2+k2)3/k2 يبلغ أصغر قيمة له عند k=π/2، حيث يساوي 27π4/4. استنتج عرض زوج من خلايا الحمل عند البدء، بوحدات d، وعرض خلية واحدة.
حل
حل التمرين 6.11.
lnf=3ln(π2+k2)−2lnk، ومشتقّتها 6k/(π2+k2)−2/k تنعدم عندما 3k2=π2+k2، k2=π2/2. وعندئذٍ f=(3π2/2)3/(π2/2)=27π4/4. يضمّ طول الموجة 2π/k=22d=2.83d زوجًا من الخلايا المتعاكسة الدوران، عرض كلّ منها 2d=1.41d.
تمرين 6.12★★★
صفيحة طولها 6000km تتحرّك بسرعة 6cm/yr فوق طبقة من الوشاح سمكها 150km، لزوجتها 1020Pas، وقاعدتها ساكنة. (أ) احسب السحب القاعدي بوحدة MPa وقوة السحب لكل وحدة طول من الصفيحة. (ب) قارنها بدفع ظهرة قدره 3.8×1012Nm−1 وشدّ لوح مندسّ قدره 4.6×1013Nm−1. (ج) إذا وصل عُشر شدّ اللوح المندس إلى الصفيحة، فبأيّ سرعة يتعادل السحب مع قوى الدفع؟
حل
حل التمرين 6.12.
(أ) u=1.90×10−9ms−1، ηu/h=1020×1.90×10−9/1.5×105=1.27MPa، وعلى مدى 6000km، 7.6×1012Nm−1. (ب) ضعف قوة دفع الظهرة وسدس شدّ اللوح المندس: يستطيع السحب القاعدي أن يقاوم قوة دفع الظهرة وحدها، لا اللوح المندس. (ج) ηuL/h=3.8×1012+4.6×1012 يعطي u=8.4×1012×1.5×105/(6×106×1020)=2.1×10−9ms−1=6.6cm/yr، وهي سرعة الصفائح المحيطية السريعة.
6.8 مسألة: ارتداد بحيرة أردو
مسألة 6.1
مسألة نهاية الأسبوع — ارتداد حوض بحيرة متخيَّلة، لزوجة الوشاح أسفله، والحمل الذي تُتيحه تلك اللزوجة
بحيرة أردو، وهي بحيرة متخيَّلة، ملأت حوضًا عرضه 250km، حجزها غطاءٌ جليدي، وبقيت بعمق 300m نحو خمسين ألف سنة. وقبل أربعة عشر ألف سنة انهار السدّ الجليدي، وتصرّفت البحيرة في بضعة أسابيع. وأعلى خط شاطئ لها، وكان أفقيًّا حين تكوّن، فهو الآن مائل؛ وخطوط الشاطئ المؤرَّخة للبحيرات الصغيرة المتبقّية في الحوض تعطي رفع المركز منذ التصرّف، بالنسبة إلى الحافة، التي بالكاد تحرّكت:
احسب هبوط القاع عند التوازن الأيزوستاتي، وبيّن لماذا تظهر فيه كثافة الوشاح لا كثافة القشرة.
عامِل البحيرة قرصًا واحسب كتلة مائها.
فسّر لماذا أعلى خط شاطئ مائل الآن، واحسب متوسط ميله من الحافة إلى المركز، بالمتر لكل كيلومتر.
عامِل الهبوط توافقيًّا واحدًا نصفه المنخفض، من −λ/4 إلى λ/4، يغطي البحيرة. أعطِ λ وk=2π/λ.
عند أيّ عمق يتناقص التدفق تحت حمولة كهذه بالمعامل e؟
الجزء الثاني — زمن الارتخاء.
بافتراض أنّ البحيرة بلغت التوازن، احسب الهبوط المتبقّي عند كلّ تاريخ في الجدول.
بيّن أنّه إذا تناقص الهبوط أُسّيًّا، فإنّ لوغاريتمه يتناقص خطيًّا مع الزمن.
احسب τ من التاريخين 1 و14kyr.
تنبّأ بالرفع عند 6kyr باستعمال τ المحسوبة، وقارنه مع الجدول.
احسب الهبوط المتبقّي اليوم ومعدّل الرفع الحالي للمركز، بوحدة mm/yr.
كم من الرفع ما زال باقيًا؟
بعد كم سنة أخرى سينخفض معدّل الرفع دون 0.5mm/yr؟
الجزء الثالث — اللزوجة.
احسب لزوجة نصف الفضاء من τ.
قارنها بلزوجات الوشاح العلوي المحصورة بين 3.4×1020 و1×1021Pas التي وُجدت أسفل فينوسكانديا. أيّ جزء من الوشاح يسبره ارتداد بحيرة أردو؟
بتلك اللزوجة، احسب زمن الارتخاء لحمولة 3000km طول موجتها. لماذا يرتخي غطاء جليدي حقيقي بهذا الحجم أبطأ من هذا؟
ماذا يفعل الغلاف الصخري المرن، الذي أُهمل هنا، في ارتداد حمولة لا يتجاوز عرضها 250km؟
احسب زمن ماكسويل للوشاح باستخدام μ=7.7×1010Pa، وبيّن إن كان النموذج اللزج مبرَّرًا.
تحقّق من الفرضية في السؤال 7: ما كسر هبوط التوازن الذي بلغته خلال بقاء البحيرة؟
بمقدار كم سيتغيّر log10η لو أُخذ λ يساوي 250km؟
الجزء الرابع — الحمل. ليكن طبقة من الوشاح العلوي عمقها 670km وكثافتها ρ=3700kgm−3، وg=9.9ms−2، وα=3.1×10−5K−1، وΔT=1450K، وκ=10−6m2s−1، وk=3.1Wm−1K−1 واللزوجة الواردة في الجزء الثالث.
4. تشكّل خط الشاطئ أفقيًّا، حين كان القاع أشدّ انخفاضًا عند المركز؛ فقد ارتفع المركز أكثر من الحافّة منذ ذلك الحين، لذا يصعد الخطّ الآن نحو المركز، بمقدار 77.2m على امتداد 125km: 0.62m/km.
15. أدنى قليلًا من الطرف الأدنى لمدى فينوسكانديا. يتلاشى الجريان أسفل حمل طول موجته 500km على امتداد 80km: يُعاين الارتداد أضعف الوشاح وأقربه إلى السطح، أسفل ما هو في الأرض الحقيقية جزئيًّا، داخل الغلاف الصخري.
16.τ=4π×3.0×1020/(3300×9.8×3×106)=3.9×1010s=1.2kyr. تحت حمل بهذا الاتّساع يصل الجريان إلى عمق 477km فأكثر، نحو وشاح سفليّ أكثر لزوجةً بعشر مرّات، وهو ما يُبطئه.
17. يحمل جزءًا من حمل ضيّق بالانثناء: الانخفاض عند التوازن أقلّ عمقًا وأوسع من الانخفاض الأيزوستاتي البالغ 90.9m، وجزءٌ من الارتداد مرنٌ وفوريّ. نصفُ فضاءٍ لزِجٍ متجانسٍ يُلائِم الرفعَ ثم يُعطي لزوجةً تمزج التأثيرين معًا.
18.τM=3.0×1020/7.7×1010=3.9×109s، نحو 124yr، أقصر بستّين مرّة من τ: على مدى آلاف السنين يستجيب الوشاح استجابةَ مائعٍ لزِج.
19.1−e−50/7.40=0.9988: بلغ 99.9% من الانخفاض عند التوازن، والافتراض صحيح.
24.u=(10−6/6.7×105)×10.782=1.73×10−10ms−1=0.55cm/yr، أبطأ من الصفائح بما بين مرّتين وأربع عشرة مرّة.
25.q=5.4×3.1×1450/6.7×105=36mWm−2. تُهمل الطبقةُ الحرارةَ المتولّدةَ داخل الوشاح والحرارةَ الآتية من اللبّ والوشاح السفليّ، والقياسُ يُسقط ثوابتَه من رتبة الواحد، بينما تفقد الصفائح، إذ تبرد بوصفها أنصاف فضاءات منطلقةً من الظهرات، حرارةً أكبر بكثير على القاع الفتيّ مما كانت ستفقده طبقة حدّية متوسّطة.