পৃথিবী ও মহাকাশ বিজ্ঞান · বই 5 · স্নাতক তৃতীয় বর্ষ
علوم الأرض والفضاء الجامعية — السنة الثالثة
বইগুলো এখনো বাংলা ভাষায় পাওয়া যায় না — العربية সংস্করণ দেখানো হচ্ছে।
علوم الأرض والفضاء الجامعية — السنة الثالثة · স্নাতক তৃতীয় বর্ষ
3الجاذبية وشكل الأرض
يدور قمران صناعيان أحدهما خلف الآخر حول الأرض، على ارتفاع يزيد على 500km ويفصل بينهما نحو 220km، ويقيسان المسافة بينهما باستمرار بدقّة تصل إلى ميكرومتر واحد.
عندما يمرّان فوق مكان تتركّز فيه الكتلة، يُسحب القمر المتقدّم إلى الأمام أوّلًا فيتّسع الفاصل بينهما؛ ثم يُسحب القمر المتأخّر بدوره فيضيق الفاصل. ومن هذه التغيّرات الصغيرة، شهرًا بعد شهر منذ عام 2002، رُسم حقل جاذبية الأرض وهو يتغيّر: فالمياه المخزّنة في الأحواض النهرية ترتفع وتنخفض مع الفصول، وتفقد الأغطية الجليدية الكبرى في غرينلاند والقارة القطبية الجنوبية كتلتها عامًا بعد عام. وتَزِن الجاذبية ما تعجز العين عن رؤيته. يبني هذا الفصل حقل الجاذبية لكوكب دوّار انطلاقًا من جهده، يستخرج منه التفلطح وعزم القصور الذاتي، ويتتبّع شذوذاته حتى المدار، ثم ينتقل إلى أحد ما تكشفه الجاذبية عند السطح: انثناء الغلاف الصخري تحت ثقل بركان.
ما تعرفه من قبل
قانون نيوتن للتجاذب الكوني، وطاقة الكمون الجذبي، والقوة الطاردة المركزية في إطار دوّار هي من كتب الفيزياء؛ وعزوم القصور الذاتي لجسم صلب وكثيرات حدود ليجندر P2(x)=(3x2−1)/2 واردة في المجلّدات الجامعية لكتب الفيزياء والرياضيات، ومعامل يونغ E ونسبة بواسون ν في مجلّد السنة الثالثة لكتب الفيزياء. وقد عُرفت الجيوديسيا والتفلطح f=(a−b)/a للإهليلجي ذي نصفي المحورين a>b في الصف العاشر. وعُرّف التوازن الأيزوستاتي، والتعويض الأيري، والجاذبية المرجعية، والجيود، وارتفاع الجيود، وشذوذات الجاذبية في مجلّد السنة الأولى، في فصله المخصّص للتوازن الأيزوستاتي وشذوذات الجاذبية.
3.1 الجهد الأرضي
على كوكب دوّار يشعر الجسم الساكن بجذب كتلة الكوكب وبالقوة الطاردة المركزية في الإطار الدوّار معه. وكلاهما يشتقّ من جهدٍ واحد.
تعريف 3.1(الجهد الأرضي)
الجهد الأرضي هو طاقة الكمون لكلّ وحدة كتلة، في الإطار الدوّار مع الأرض بالسرعة الزاوية ω، للجذب والقوة الطاردة المركزية معًا:
حيث s هو البعد عن محور الدوران. والجاذبية g=−∇W، والسطوح W=ثابت هي سطوحه المتساوية الجهد: الجاذبية عمودية عليها في كلّ مكان.
يتدلّى الخيط الشاقولي على امتداد −∇W، ويستقرّ الماء الساكن على سطح متساوي الجهد. ومنه يأخذ كوكب مائع دوّار شكل أحد سطوحه المتساوية الجهد: فالكمون الطارد المركزي، الأعمق عند خطّ الاستواء، يُضخّم خطّ الاستواء ويُفلطح القطبين.
سطح متساوي الجهد لكوكب دوّار (تخطيطيًّا، مع تضخيم التفلطح تضخيمًا كبيرًا). عند كلّ نقطة، تكون الجاذبية، أي مجموع الجذب والقوة الطاردة المركزية، عمودية على سطح متساوي الجهد؛ فهي لا تتّجه نحو المركز إلّا عند القطبين وعلى خطّ الاستواء.
خارج الكوكب يحقّق الجهد الجذبي معادلة لابلاس، وحلوله المنتظمة عند اللانهاية هي التوافقيات الكروية الصلبة.
قضية 3.2(نشر الجهد الأرضي بالتوافقيات الكروية)
خارج جسم كتلته M ونصف قطره الاستوائي a، ومبدؤه عند مركز كتلته، ومحوره z على امتداد محور دورانه، يكون الجهد الجذبي عند المسافة r والزاوية القطبية θ وخط الطول λ هو
مع كثيرات حدود ليجندر Pn والدوال المرافقة Pnm، والمعاملات Jn، Cnm، Snm أعداد تُعيَّن بالمطابقة مع الرصد. وبالنسبة إلى جسم متناظر حول محوره، لا يبقى إلّا الحدود المنطقية Jn. ولا يظهر حدّ من الدرجة 1.
برهان. يُعطي فصل المتغيّرات في معادلة لابلاس في الإحداثيات الكروية الحلول المنتظمة عند اللانهاية r−(n+1) مضروبًا في Pnm(cosθ)e±imλ (المجلّدات الجامعية لكتب الرياضيات)؛ وأيّ جهد خارج الكتل هو مجموع هذه الحلول. والحدّ من الدرجة 0 هو −GM/r، أي جهد الكتلة كلّها. والحدود من الدرجة 1 تتناسب مع العزوم الأولى ∫r′dm′، وتنعدم هذه عندما يكون المبدأ عند مركز الكتلة. ويُزيل التناظر المحوري كلّ حدّ يعتمد على λ. ∎
تعريف 3.3(معامل الشكل الديناميكي)
معامل الشكل الديناميكيJ2 لكوكبٍ ما هو معامل الحدّ الزونالي التوافقي من الدرجة 2، أي أكبر انحراف لجهده الخارجي عن جهد الكرة، وذلك بسبب تفلطحه. وبالنسبة إلى الأرض، J2=1.08263×10−3.
مبرهنة 3.4(صيغة ماكولا)
لتكن A، B، C عزوم القصور الذاتي الرئيسية لجسمٍ كتلته M حول محاور تمرّ بمركز كتلته، وI عزم قصوره الذاتي حول المستقيم الواصل من مركز الكتلة إلى نقطة خارجية على بعد r. وحتى الرتبة r−3 يُعطى الجهد الثقالي عند تلك النقطة وفق صيغة ماكولا
V=−rGM−2r3G(A+B+C−3I).
وبالنسبة إلى جسم متماثل حول المحور z (A=B)، V=−rGM+r3G(C−A)P2(cosθ)، ومنه
J2=Ma2C−A.
برهان. وليكن r^ متجهة الوحدة الموجَّهة نحو النقطة، وr′≪r،
∣r−r′∣1=r1+r2r^⋅r′+2r33(r^⋅r′)2−r′2+O(r−4).
بالتكامل على الكتلة، يُعطي الحدّ الأوّل M/r، وينعدم الحدّ الثاني (مركز الكتلة). وفي الحدّ الثالث، ∫r′2dm′=(A+B+C)/2، لأنّ A+B+C=∫2r′2dm′، و∫(r^⋅r′)2dm′=∫r′2dm′−I، ومربّع المسافة إلى المستقيم هو r′2−(r^⋅r′)2. ويساوي تكامل القوس 3(A+B+C)/2−3I−(A+B+C)/2=A+B+C−3I. وعند A=B، يكون العزم حول مستقيم عند الزاوية المتمّمة لخط العرض θ هو I=Asin2θ+Ccos2θ، و2A+C−3I=(C−A)(1−3cos2θ)=−2(C−A)P2(cosθ). وبمقارنة ذلك مع القضية 3.2 نجد J2Ma2=C−A. ∎
ومن ثمّ يقيس الحدّ من الدرجة 2 في الجهد الفرقَ بين عزمَي القصور الذاتي، لا العزوم نفسها. أمّا تقدّم محور الأرض، الذي تدفعه عزوم الشمس والقمر المؤثّرة في الانتفاخ الاستوائي، يقيس الإهليلجية الديناميكية H=(C−A)/C؛ وبالجمع بين J2 وH نحصل على C/(Ma2)، وهو عامل عزم القصور الذاتي، ويساوي 0.3308 للأرض،
أدنى بكثير من 0.4 الخاصّة بالكرة المتجانسة: إذ تتركّز الكتلة نحو المركز.
3.2 الإهليلجي الدوّار
تعريف 3.5(الإهليلجي المرجعي)
الإهليلجي المرجعي هو إهليلجي دوراني، ذا نصف المحور الاستوائي a والتفلطح f، ويُختار بحيث يكون هو نفسه سطح تساوي جهدٍ نموذجيّ: جهدُ جسمٍ دوّار كتلته M بالسرعة الزاوية ωومعامل الشكل الديناميكيJ2 للأرض. وتحدّد هذا الإهليلجي أربعةُ ثوابت؛ ويأخذ نظام المرجع الجيوديسي 1980 القيم الآتية:
ومنها يُشتقّ التفلطح f=0.00335281068 (1/f=298.257222) ونصف المحور الصغير b=6356752.3141m.
يُعرف جهدُ النموذج، أي الجهد المرجعي، على صورة مغلقة، وكذلك جاذبيتُه على الإهليلجي، أي الجاذبية المرجعية في مجلد السنة الأولى.
قضية 3.6(صيغة سوميليانا)
على إهليلجي مرجعي ذي نصفي المحورين a وb، والجاذبية المرجعية γe عند خط الاستواء وγp عند القطبين تُعطى عند خط العرض الجيوديسي φ بالعلاقة صيغة سوميليانا
γ(φ)=a2cos2φ+b2sin2φaγecos2φ+bγpsin2φ.
بالنسبة إلى نظام المرجع الجيوديسي 1980، γe=9.7803267715ms−2 وγp=9.8321863685ms−2، والسلسلة γ=9.780327(1+0.0053024sin2φ−0.0000058sin22φ)ms−2 تُطابق الصيغة إلى حدود 0.1mGal.
برهان.نقبله في هذا المستوى.∎
ما مقدار التفلطح الذي ينبغي أن يبلغه كوكبٌ دوّار؟ إذا كان مائعًا في التوازن الهيدروستاتيكي، فإنّ تفلطحه يتحدّد بدورانه وبكيفية توزّع كتلته في العمق.
قضية 3.7(علاقتا كلّيرو ودارون–رادو)
لكوكب في توازن هيدروستاتيكي، ولتكن m=ω2a2b/(GM) هي نسبة التسارع النابذ إلى التسارع الثقالي عند خط الاستواء، وحتى الرتبة الأولى في f وm:
الكوكب المتجانس (C/Ma2=0.4) يتحقّق فيه f=5m/4؛ أمّا كوكب تتركّز كتلته نحو المركز فيكون أقلّ تفلطحًا.
برهان.نقبله في هذا المستوى.∎
تُحوّل علاقةُ كلّيرو، مع m=0.00344979 الخاصّ بالأرض، قيمةَ J2 المقاسة إلى f=0.0033488، أو 1/298.6: ضمن 0.12% من تفلطح الإهليلجي. وأمّا علاقة دارون–رادو مع عامل عزم القصور الذاتي، تُعطي التفلطحَ الذي كانت ستبلغه الأرض لو كانت مائعًا في توازن: 1/299.9، مقابل 1/298.26 المقاسة. والأرض أكثر تفلطحًا قليلًا من كوكب مائع بدورانه وتوزيع كتلته؛ فانتفاخها الاستوائي يعلو بنحو 116m فوق ما يتنبّأ به التوازن الهيدروستاتيكي. والأرض ليست تمامًا مائعًا في توازن، والفرقُ الصغير بحدّ ذاته قياسٌ للإجهادات التي يتحمّلها باطنها.
تفلطح كوكب مائع بدوران الأرض، مقابل عامل عزم القصور الذاتي له، من علاقة دارون–رادو (محسوب من Geodetic Reference System 1980). وكلّما ازداد تركّز الكتلة نحو المركز صار الكوكب أكثر استدارة. وعند قيمة الأرض 0.3308 يكون التفلطح الهيدروستاتيكي 1/299.9، أقلّ قليلًا من 1/298.26 المقاسة.
3.3 الجيود والجهد المُضطرِب
الجيود هو السطح المتساوي الجهد للجهد الأرضي الحقيقي الذي ينطبق على متوسط مستوى سطح البحر؛ أمّا الإهليلجي المرجعي فهو سطح متساوي الجهد الطبيعي U ذي القيمة نفسها، U0=W0. أمّا فصلهما، أي ارتفاع الجيود N، فيتحدّد بفرق الجهدين.
قضية 3.8(صيغة برونس)
ليكن T=W−U الجهد المُضطرِب. يُعطى ارتفاع الجيود عند نقطة وفق صيغة برونس:
N=−γT,
حيث يُحسب T على الجيود، وγ هي الجاذبية المرجعية هناك. أمّا فائض الكتلة الذي يُعمّق الجهد (T<0)، فيرفع الجيود.
برهان. ليكن P على الجيود وQ النقطةَ الواقعةَ على الإهليلجي تحتها، على بُعد N. لأنّ U يتزايد نحو الخارج بالمعدّل γ، U(P)=U(Q)+γN=U0+γN. ثمّ W0=W(P)=U(P)+T(P)=U0+γN+T، وW0=U0 يُعطي N=−T/γ. ∎
يقابل ارتفاعَ الجيود 50m فرقُ جهدٍ قدرُه فقط 490Jkg−1، أي جزءًا واحدًا من 105 من الجهد الأرضي. وتتراوح ارتفاعات الجيود للأرض من نحو +85m في آيسلندا إلى −106m في جنوب الهند،
أي مدًى دون 200m على كوكب يختلف نصفا قطريه الاستوائي والقطبي بمقدار 21km.
3.4 قياس الجاذبية بالأقمار الصناعية
تعريف 3.9(قياس الجاذبية بالأقمار الصناعية)
قياس الجاذبية بالأقمار الصناعية يحدّد حقل الجاذبية الأرضي من حركة الأقمار الصناعية: من الاضطرابات طويلة الأمد لمداراتها، ومن تغيّر المسافة بين قمرين صناعيين على المدار نفسه، أو من مقاييس التدرّج التي تقيس فروق الجاذبية عبر المركبة الفضائية.
أكبر اضطراب لمدار قريب من الأرض يأتي من J2. ويمارس الانتفاخ الاستوائي عزمًا على المدار المائل، فيدور ببطء حول محور الأرض.
قضية 3.10(التقدّم العقدي)
إلى الرتبة الأولى في J2، يدور خطّ العقد لمدار دائري نصف قطره r وميله i بالمعدّل
Ω˙=−23nJ2(ra)2cosi,n=GM/r3,
يتجه غربًا في المدار المباشر (i<90∘)، ويتجه شرقًا في المدار الرجعي.
برهان. عبر دورات كثيرة يتصرّف القمر الصناعي كحلقة كتلتها ms موزّعة على مداره. يُكتب الجزء ذو الدرجة 2 من الجهد، (GMJ2a2/r3)P2(cosθ) لكلّ وحدة كتلة، متوسّطًا على دائرة مائلة بميل i، يتضمّن ⟨cos2θ⟩=21sin2i، فإنّ ⟨P2⟩=43sin2i−21، وطاقة الحلقة هي E(i)=ms(GMJ2a2/r3)(43sin2i−21). ويكون العزم حول خطّ العقد هو τ=−dE/di=−23ms(GMJ2a2/r3)sinicosi. وهو عمودي على كمية الحركة الزاوية للحلقة L=msnr2، لذلك يدير L حول المحور القطبي دون أن يغيّر طولها: LsiniΩ˙=τ. ومع GM=n2r3، فإنّ Ω˙=−23nJ2(a/r)2cosi. ∎
يبلغ المعدّل عدّة درجات في اليوم لمدار منخفض قرب خطّ الاستواء؛ وعند ميل قليلًا بعد 90∘ يتطابق مع الحركة السنوية للشمس، ويحتفظ المدار بالتوجّه نفسه نحو الشمس طوال السنة (انظر التمرين 12). وعند القياس على أقمار صناعية كثيرة بميول مختلفة، تُعطي هذه التقدّمات J2 والتوافقيات المنطقية الأعلى.
يؤدّي زوج من الأقمار الصناعية أفضل: إذ يقيس تغيّر المسافة بينهما فرق الجاذبية التي يشعران بها على طول المدار. طار زوج GRACE من 2002 إلى 2017، وخلَفه GRACE Follow-On، منذ إطلاقه في 22 مايو 2018، وتقيس وصلته الميكروويفية المسافة بدقّة في حدود الميكرومتر، ووصلة ليزرية تختبر دقّةً تصل إلى عشرين مرّة أفضل.
يتغيّر الحقل من شهر إلى آخر لأنّ الماء يتحرّك: الثلج ورطوبة التربة، والمياه الجوفية، والبحيرات، والمحيطات، والغطاء الجليدي.
مبدأ زوج الأقمار الصناعية (تخطيطي، ليس بمقياس الرسم). عند الاقتراب من كتلة زائدة، يشعر القمر المتقدّم بجذبها أولًا ويُعجَّل نحو الأمام، فيزداد المدى بينهما؛ وعندما يمرّ القمر المتأخّر فوق الكتلة فيتناقص المدى من جديد. ومعدّل تغيّر المدى المقاس باستمرار يُعكس لاستنتاج حقل الجاذبية.
القمران الصناعيان GRACE Follow-On، تصوّر فنّي (المسافة بينهما مختصرة كثيرًا). يرمز الشعاع الأحمر إلى وصلة القياس بالليزر. NASA/JPL-Caltech، ملك عام.
تعريف 3.11(ارتفاع الماء المكافئ)
يمكن تمثيل تغيّر حقل الجاذبية الذي يتغيّر ببطء على مسافات أكبر بكثير من سمك الطبقة التي تحرّكت فيها الكتلة بتغيّر في الكثافة السطحية σ (kgm−2) على قشرة رقيقة عند سطح الأرض. يُسمّى ارتفاع الماء المكافئ هو سمك الماء العذب الذي له الكثافة السطحية نفسها:
hw=ρwσ.
متر واحد من الماء فوق حوض مساحته 1×106km2 يساوي 1×1015kg، أو 1000Gt. مستوى سطح البحر هو الطرف الآخر من الميزانية: 360Gt من الجليد يُضاف إلى المحيط يرفعه بنحو ملّيمتر واحد.
قياس الجاذبية بالأقمار الصناعية، وقياس الارتفاع بالأقمار الصناعية، وحصيلة تساقط الثلج مقابل تصريف الجليد، تقيس كلٌّ منها الحصيلة الكتلية للغطاء الجليدي؛ وبعد التوفيق بينها، تُعطي السلسلة الآتية.
التغيّر التراكمي للكتلة في الغطاءين الجليديين لغرينلاند وأنتاركتيكا منذ 1992، مع نطاقات بانحراف معياري واحد (تقديرات IMBIE 2021 المصالَحة من قياس الارتفاع والجاذبية وطريقة المدخلات–المخرجات). بحلول نهاية عام 2020 فقدت غرينلاند 4890±460Gt وأنتاركتيكا 2670±530Gt؛ تعطي الملاءمات الخطّية −176Gtyr−1 و−91Gtyr−1 خلال الفترة 1992–2020، و−246Gtyr−1 و−127Gtyr−1 خلال سنوات GRACE.
شذوذ الجاذبية المقاس عند السطح هو مجموع توافقيات؛ وعند ارتفاع قمر صناعي، يضعف كلٌّ منها بعامل يعتمد على طول موجته.
قضية 3.12(الاستمرار إلى الأعلى)
فوق الكتل المسبّبة له، يكون شذوذ الجاذبية الذي يتغيّر على طول السطح بالصيغة Δg0coskx، حيث k=2π/λ، عند الارتفاع z (أي الاستمرار إلى الأعلى)
Δg(x,z)=Δg0e−kzcoskx.
برهان. فوق الكتل، يحقّق الجهد الشاذّ، ومعه الشذوذ، معادلة لابلاس ∂x2Δg+∂z2Δg=0. بفصل المتغيّرات، Δg=Z(z)coskx مع Z′′=k2Z، وحلّها المحدود عندما z→∞ هو Z=Δg0e−kz. ∎
عند ارتفاع 400km يحتفظ شذوذ طوله 1000km بنسبة 8% من سعته، والذي طوله 100km يحتفظ بنحو 1×10−11 منها: ترى الأقمار الصناعية الأطوال الموجية الطويلة، وترصد المسوح عند السطح القصيرة.
3.5 انثناء الغلاف الصخري
عوّض مجلد السنة الأولى حملًا محليًا، إذ يطفو كلّ عمود على حدة. لا يُعوَّض بركانٌ واقعٌ على صفيحة محيطية تحته وحده: تنحني الصفيحة تحته على عرض مئات الكيلومترات، ويهبط قاع البحر حول البركان في خندقٍ محيطيّ، وما وراءه يرتفع في قوس منخفض.
تعريف 3.13(انثناء الغلاف الصخري، صلابة الانثناء، السمك المرن الفعّال)
انثناء الغلاف الصخري هو انحناء الغلاف الصخري تحت حمل، ويُنمذَج بوصفه انحناء صفيحة مرنة رقيقة طافية على وشاح مائع. ومقاومة الصفيحة للانحناء هي صلابة الانثناء
D=12(1−ν2)ETe3,
وTe هو سمك الصفيحة المرنة المتجانسة التي تنحني كما ينحني الغلاف الصخري الحقيقي، وهو السمك المرن الفعّال.
مبرهنة 3.14(انثناء صفيحة رقيقة)
صفيحة مرنة رقيقة صلابتها D، تنحني في اتجاه واحد بفعل حمل q(x) لكلّ وحدة مساحة، تنحرف نحو الأسفل بمقدار w(x)، وتطفو على وشاح كثافته ρm، والتجويف مملوء بمادة كثافتها ρi، فتحقّق معادلة الانثناء
Ddx4d4w+(ρm−ρi)gw=q(x).
برهان. في الصفيحة المنحنية، تتقلّص الألياف الواقعة فوق المستوى المنتصف، وتتمدّد الواقعة تحته، بما يتناسب مع الانحناء w′′؛ وبتكامل إجهاد الانحناء على السمك نحصل على عزم الانحناء M=−Dw′′، مع الصلابة الواردة في التعريف (لا تتمدّد الصفيحة جانبيًا، ومن ثمّ المعامل 1−ν2؛ انظر مجلد السنة الثالثة من كتب الفيزياء). شريحةٌ dx من الصفيحة تحمل قوة القصّ V على وجهيها والعزم M. التوازن الشاقولي: الحمل qdx نحو الأسفل، وقوة الطفو للوشاح المُزاح (ρm−ρi)gwdx نحو الأعلى، فيعطي dV/dx=−q+(ρm−ρi)gw. وتوازن العزوم يعطي dM/dx=V. ومنه −Dw′′′′=M′′=V′=−q+(ρm−ρi)gw. ∎
تعريف 3.15(معامل الانثناء، الانتفاخ الأمامي)
معامل الانثناء لصفيحة ما هو الطول
α=((ρm−ρi)g4D)1/4,
الذي يحدّد المقياس الأفقي لانثنائها. الانتفاخ الأمامي هو القوس المنخفض من الانثناء نحو الأعلى الذي ترفعه الصفيحة المنحنية خارج الخندق المحيطي أو الحوض المجاور لحملها.
قضية 3.16(حمل خطّي على صفيحة متّصلة)
يُحدث حمل خطّي V0 لكلّ وحدة طول على امتداد x=0 انحرافًا للصفيحة مقداره
w(x)=w0e−∣x∣/α(cosα∣x∣+sinα∣x∣),w0=8DV0α3.
ينعدم الانحراف عند ∣x∣=3πα/4، أي حافّة الخندق المحيطي، وتقع قمّة الانتفاخ الأمامي عند ∣x∣=πα، وارتفاعها w0e−π=0.0432w0. وتوازن قوةُ الطفو لكامل الانحراف الحملَ: ∫(ρm−ρi)gwdx=V0.
برهان. من أجل x>0 تكون المعادلة w′′′′+4w/α4=0؛ وجذور r4=−4/α4 هي (±1±i)/α، والحلول التي تنعدم بعيدًا هي e−x/α(Acos(x/α)+Bsin(x/α)). وتفرض التناظر w′(0)=0، وهذا يعطي A=B=w0. وبتكامل التوازن الشاقولي عبر الحمل، تقفز قوة القصّ V=−Dw′′′ بمقدار −V0؛ وبالتناظر w′′′(0−)=−w′′′(0+)، وw′′′(0+)=4w0/α3، فيكون 8Dw0/α3=V0. وينعدم الانحراف حيث tan(x/α)=−1، أوّلًا عند x=3πα/4؛ ومشتقّته، −(2w0/α)e−x/αsin(x/α)، تنعدم بعد ذلك عند x=πα. وأخيرًا ∫0∞e−u(cosu+sinu)du=1، إذن ∫wdx=2w0α=V0α4/(4D)=V0/((ρm−ρi)g). ∎
قضية 3.17(صفيحة مكسورة)
صفيحةٌ مكسورة عند x=0 ويُدفع طرفها الحرّ نحو الأسفل، كما هو الحال عند خندقٍ بحريّ عميق، تنحرف وفق w=w0e−x/αcos(x/α) من أجل x≥0: ينعدم الانحراف عند x=πα/2، وتقع قمّة الانتفاخ الأمامي، أي الارتفاع الخارجي، عند x=3πα/4، وارتفاعها w0e−3π/4/2=0.067w0.
برهان. الحلول المتناقصة نفسها، لكن مع طرف حرّ لا يحمل عزمًا، w′′(0)=0، وهذا يعطي B=0. ويقع الصفر عند cos(x/α)=0؛ والمشتقّة −(w0/α)e−x/α(cos+sin)(x/α) تنعدم عند x=3πα/4، حيث w=w0e−3π/4cos(3π/4). ∎
ومن ثمّ فالمسافات هي مفتاح قوّة الصفيحة: يُقاس موضع انتهاء الخندق المحيطي أو موضع قمّة الانتفاخ، فيُستنتج α، ثم D، ثم Te. ولمّا كان Te يدخل بوصفه α4/3، فإنّ تقدير المسافة تقديرًا معقولًا يعطي تقديرًا جيدًا لسمك الصفيحة.
انثناء صفيحة مرنة (محسوب؛ E=100GPa، ν=0.25، الوشاح 3300kgm−3، الماء في الحفرة المحيطة؛ مقياس رأسي مبالَغ فيه كثيرًا). (أ) تحت حمل ظهرة بركانية (4.17×1012Nm−1)، تغوص صفيحة رقيقة عميقًا وضيقًا، وتغوص صفيحة سميكة ضحلًا وعريضًا؛ ينتقل الانتفاخ الأمامي من 112km إلى 316km. (ب) صفيحة مكسورة عند خندق ومدفوعة إلى أسفل بمقدار 1km: يصل الارتفاع الخارجي إلى 67m، على بُعد 167km من الخندق.
تعريف 3.18(التعويض الأيزوستاتي الإقليمي)
التعويض الأيزوستاتي الإقليمي هو إسناد حمل بطفو انحراف ممتدّ على منطقة أوسع من الحمل، بواسطة صلابة الانثناء في الغلاف الصخري؛ أمّا في تعويض إيري أو برات، فيُعوَّض كل عمود تحت ذاته على العكس.
يوازن التعويضان المحلي والإقليمي الكتلة الكلية نفسها، لكنّهما يوزّعانها بطريقة مختلفة: تحت حمل ضيق على صفيحة جاسئة يكون الجذر ضحلًا وعريضًا، ويكون الشذوذ الجاذبي فوق الحمل موجبًا بقوة، يكتنفه من جانبيه شذوذات سالبة فوق الحفرة المحيطة. وحول جزر ظهرة هاواي تعطي المقاطع الجاذبية والزلزالية سمكًا مرنًا فعّالًا من 17 إلى 37km، و10 إلى 20km على طول الجبال البحرية الأقدم في سلسلة الإمبراطور، وهو فرق يُتوقَّع من أعمار البراكين وأعمار الصفيحة التي حمّلتها؛ وتتبع القيم عمق خط متساوي الحرارة قرب 350∘C في الصفيحة المحيطية المتبرّدة.
يفتح الانثناء نفسه، تحت وزن حزام جبلي، أحواض المقدّمة في مجلد السنة الثانية، في فصله عن الأوروجين الألبي.
جزيرة بركانية على صفيحة محيطية (مقطع عرضي تخطيطي، مقياس رأسي مبالَغ فيه). تنحني الصفيحة تحت البركان؛ الفجوة المحيطة به، أي الحفرة المحيطة، تمتلئ بالفتات البركاني والرواسب؛ أمّا ما وراءها فترتفع الصفيحة في انتفاخ أمامي منخفض.
في الميدان— مسح جاذبي نسبي
مقياس جاذبية ميداني على قاعدة التسوية الخاصة به، مع هوائي قمر صناعي في الخلفية (رسم توضيحي مولَّد بالذكاء الاصطناعي).
يُسقط مقياس الجاذبية المطلق مرآةً مكعبة الزوايا في غرفة مفرغة مرارًا وتكرارًا، ويقيس زمن سقوطها بمقياس تداخل ليزري وساعة ذرية؛ فيعطي الجاذبية عند محطة واحدة، لكنّه ثقيل وبطيء الإعداد. ويُنجَز المسح نفسه بمقياس الجاذبية النسبي، وهو ميزان زنبركي صغير يقيس فروق الجاذبية. يُسوّيه الجيوفيزيائي على لوحة قاعدة عند كل محطة على امتداد خطّ، فيقرأه ويدوّن الزمن والارتفاع وينتقل إلى التالية؛ وكل بضع ساعات يعود المسار إلى محطة أساس، فيُقاس الانجراف البطيء للزنبرك ومدّ الأرض الصلبة ويُزال. ويعطي مستقبِل القمر الصناعي على عمود قريب الموضعَ والارتفاع، لأنّ الجاذبية تتغيّر بنحو 0.3mGal مع كل متر من الارتفاع.
3.6 التمارين
تمرين 3.1★
انطلاقًا من قيم نظام المرجع الجيوديسي 1980 a=6378137m وf=0.00335281068، احسب نصف المحور الأصغر b وa−b. مع J2=1.08263×10−3 وm=0.00344979، احسب التفلطح الذي تتنبّأ به علاقة كليرو، وقارِن.
حل
حل التمرين 3.1.
b=a(1−f)=6378137×0.99664719=6356752.3m، وa−b=21384.7m، أي نحو 21.4km. كليرو: f=(3J2+m)/2=(0.00324789+0.00344979)/2=0.0033488، أو 1/298.61، مقابل 1/298.26: أصغر بمقدار 0.12%، أي بقدر حدود الرتبة الثانية التي تُهمِلها العلاقة.
تمرين 3.2★
في سنةٍ جافة، يفقد حوض نهري مساحته 5.0×105km2 كتلةً قدرها 60Gt من الماء. احسب التغيّر في ارتفاع الماء المكافئ. بكم سيرفع الكتلة نفسها مستوى سطح البحر، علمًا بأنّ 360Gt لكل مليمتر؟
حل
حل التمرين 3.2.
σ=6.0×1013kg/5.0×1011m2=120kgm−2، إذن hw=120/1000=0.12m: 12cm من الماء. وعند تصريفه في المحيط، 60/360=0.17mm من مستوى سطح البحر.
تمرين 3.3★
باستخدام المتسلسلة γ=9.780327(1+0.0053024sin2φ−0.0000058sin22φ)ms−2، احسب الجاذبية المرجعية عند خطّ الاستواء، وعند 45∘ وعند القطبين. عبّر عن الفرق بين القطب والاستواء بالميليغال.
حل
حل التمرين 3.3.
عند خط الاستواء γ=9.780327ms−2. وعند 45∘، sin2φ=0.5 وsin22φ=1: γ=9.780327×1.0026454=9.806200ms−2. وعند القطبين، 9.780327×1.0053024=9.832186ms−2. الفرق هو 0.051859ms−2، أو 5186mGal.
تمرين 3.4★
احسب صلابة الانثناءومعامل الانثناء للصفائح المحيطية حيث Te=10km وTe=30km (E=100GPa، ν=0.25، ρm−ρw=2270kgm−3، g=9.81ms−2). أين تقع حافة الحفرة المحيطة والانتفاخ الأمامي تحت حمل خطّي؟
حل
حل التمرين 3.4.
D=ETe3/(12(1−ν2))=1011Te3/11.25: 8.89×1021Nm من أجل 10km، 2.40×1023Nm من أجل 30km. مع (ρm−ρw)g=22269Nm−3، α=(4D/22269)1/4: 35.5km و81.0km. تنتهي الحفرة المحيطة عند 3πα/4، 84km و191km؛ تبلغ قمم الانتفاخ الأمامي عند πα، 112km و255km.
تمرين 3.5★★
(أ) ارتفاع الجيود +50m: احسب جهد الاضطراب هناك، مستخدمًا γ=9.80ms−2. (ب) على امتداد 100km يرتفع الجيود بمقدار 10m. بأيّ زاوية يختلف خيط الشاقول هناك عن العمودي على الإهليلجي؟ أعطِ الجواب بالثواني القوسية.
حل
حل التمرين 3.5.
(أ) T=−γN=−9.80×50=−490Jkg−1. (ب) يميل الجيود بمقدار 10/100000=1.0×10−4 راديان بالنسبة إلى الإهليلجي، ويميل الشاقول وهو العمودي على الجيود، بالزاوية نفسها بالنسبة إلى عمودي الإهليلجي: 1.0×10−4×206265=20.6′′.
H=J2/(C/Ma2)=1.08263×10−3/0.3308=0.0032728، و1/H=305.55. J2 هو نفسه بالفرض؛ ولو كانت الأرض متجانسة لكان C/(Ma2)=0.4 ومنه H=J2/0.4=0.0027066: لكانت المبادرة أبطأ. إنّ قياس J2 وH معًا هو ما يكشف تركّز الكتلة نحو المركز.
تمرين 3.7★★
يدور قمر صناعي في مدار دائري على ارتفاع 500km فوق نصف القطر الاستوائي (a=6378.137km، GM=3.986005×1014m3s−2)، مائلًا بزاوية 89∘. احسب دورَه ومعدّل تقدّم العقدة بالدرجات في اليوم. لماذا تختار بعثات الجاذبية مدارات شبه قطبية؟
حل
حل التمرين 3.7.
r=6878.137km؛ n=GM/r3=1.10678×10−3rads−1، أي دور قدره 2π/n=94.6min. Ω˙=−1.5×1.10678×10−3×1.08263×10−3×(6378.137/6878.137)2×cos89∘=−2.70×10−8rads−1، أو −0.134∘ في اليوم، غربًا. يمرّ مدار شبه قطبي فوق كل خط عرض، فتُرسم الأرض كلها؛ كما يتقدّم المدار ببطء، فتملأ مسارات السطح الفراغات بانتظام.
تمرين 3.8★★
لدينا شذوذان قيمة كلٍّ منهما 10mGal عند السطح، وطولا موجتيهما 100km و1000km. احسب سعتيهما عند ارتفاع 500km.
حل
حل التمرين 3.8.
e−2π×500/100=e−31.4=2.3×10−14: 2×10−13mGal، أي لا شيء. e−2π×500/1000=e−3.14=0.0432: 0.43mGal.
تمرين 3.9★★
نحو البحر من الخندق، تقع قمّة الانتفاخ الخارجي على بعد 120km من محور الخندق. باعتماد نموذج الصفيحة المنكسرة مع وجود الماء فوقها، احسب α، D وTe، والمسافة التي تقطع عندها الصفيحة مستواها غير المنثني.
حل
حل التمرين 3.9.
3πα/4=120km يعطي α=50.9km. D=(ρm−ρw)gα4/4=22269×(50930)4/4=3.75×1022Nm، وTe=(11.25D/E)1/3=16.2km. تعبر الصفيحة مستواها غير المنثني عند بُعد πα/2=80km من الخندق.
تمرين 3.10★★★
(أ) تحقّق من أنّ w=w0e−u(cosu+sinu)، حيث u=x/α، يحقّق Dw′′′′+(ρm−ρi)gw=0 من أجل x>0، مع α كما هو معرَّف. (ب) بيّن أنّ عزم الانحناء −Dw′′ يبلغ أكبر قيمة له تحت الحمل ويغيّر إشارته عند x=πα/4. (ج) أين يُرجَّح أكثر أن تتشقّق الصفيحة عند سطحها العلوي، ولماذا؟
حل
حل التمرين 3.10.
(أ) بوضع u=x/α، d4/dx4=α−4d4/du4. بكتابة cosu+sinu=2cos(u−π/4)، فتكون الدالة هي 2w0Ree(−1+i)u−iπ/4، و(−1+i)4=−4: إذن w′′′′=−4w/α4. ثم Dw′′′′+(ρm−ρi)gw=(−4D/α4+(ρm−ρi)g)w=0، لأنّ α4=4D/((ρm−ρi)g). (ب) w′′=−(2w0/α2)e−u(cosu−sinu)، إذن M=−Dw′′=(2Dw0/α2)e−u(cosu−sinu): أكبر ما يكون في المقدار المطلق عند u=0، وينعدم حيث tanu=1، أي x=πα/4. (ج) تحت الحمل تكون الصفيحة مقعّرة نحو الأعلى، وسطحها العلوي في انضغاط؛ وبعد πα/4 تصبح محدّبة نحو الأعلى، وسطحها العلوي في تمدّد. وهناك M′∝−e−ucosu ينعدم عند u=π/2: يبلغ التمدّد أقصاه عند x=πα/2، على المنحدر الخارجي للحفرة المحيطة نحو الانتفاخ الأمامي، حيث ∣M∣ هو e−π/2=21% من قيمته تحت الحمل. وهذا هو الموضع الذي تفتح فيه الصدوع العادية والشقوق في أعلى الصفيحة.
تمرين 3.11★★★
باعتماد m=0.00344979 وC/(Ma2)=0.3308، اعكس علاقة داروين–راديو لإيجاد التفلطح الهيدروستاتيكي. قارن مع f=0.00335281: احسب الفرق والارتفاع الزائد لخطّ الاستواء فوق الشكل الهيدروستاتيكي. ماذا يدلّ ذلك عن الوشاح؟
حل
حل التمرين 3.11.
1+η=2.5(1−1.5×0.3308)=1.2595، إذن 1+η=1.5863 و25m/f=2.5863: fh=2.5×0.00344979/2.5863=0.0033346، 1/fh=299.9. والفرق f−fh=1.82×10−5؛ وبضربه في a، يرتفع خط الاستواء بمقدار 116m فوق الشكل الهيدروستاتيكي. والأرض ليست تمامًا مائعًا في حالة توازن: إذ يتحمّل باطنها إجهادات، ناتجة عن متانة الغلاف الصخري أو عن فروق التدفّق والكثافة في الوشاح، تحفظ الشكل مبتعدًا عنه قليلًا.
تمرين 3.12★★★
يُراد لقمر صناعي على مدار دائري ارتفاعه 700km أن يحافظ على التوجّه نفسه بالنسبة إلى الشمس طوال العام: يجب أن تتقدّم عقدته بمقدار 360∘ في سنة مدارية مدّتها 365.242 يومًا. احسب ميل المدار المطلوب.
حل
حل التمرين 3.12.
يجب أن تتقدّم العقدة بمقدار 2π/(365.242×86400)=1.991×10−7rads−1. عند r=7078.137km، n=1.0602×10−3rads−1 و23nJ2(a/r)2=1.398×10−6، إذن cosi=−1.991×10−7/1.398×10−6=−0.1424 وi=98.2∘: مدار قليل الارتداد، متزامن مع الشمس.
3.7 مسألة: الحفرة المحيطة بظهرة أورو
مسألة 3.1
مسألة نهاية الأسبوع — شكل الأرض، وحمل بركاني، حفرته الانثنائية وما تراه الأقمار الصناعية منها
أورو ظهرة بركانية متخيَّلة في محيط: سلسلة مستقيمة من براكين متلاحمة طولها 400km، مبنية فوق الغلاف الصخري المحيطي القديم. ويُعتبر مقطعها مثلثًا عرضه 80km عند قاعدته وارتفاعه 6km، وكثافته 2800kgm−3، واقفًا على قاع بحر عمقه 5km. كثافة الوشاح 3300kgm−3، وكثافة مياه البحر 1030kgm−3، E=100GPa، ν=0.25، g=9.81ms−2، G=6.674×10−11m3kg−1s−2.
الجزء الأول — شكل الأرض. استخدم نظام Geodetic Reference System 1980: a=6378137m، b=6356752.3m، GM=3.986005×1014m3s−2، J2=1.08263×10−3، ω=7.292115×10−5rads−1، f=0.00335281؛ وC/(Ma2)=0.3308.
احسب m=ω2a2b/(GM).
احسب التفلطح الذي تعطيه علاقة كليرو انطلاقًا من J2 وm، ومقلوبه.
احسب التفلطح الهيدروستاتيكي من علاقة داروين–راديو، ومقلوبه.
قارن كليهما بالتفلطح المرصود. بكم مترًا يعلو خطّ الاستواء فوق الشكل الهيدروستاتيكي؟
احسب الانفراج الديناميكي H=(C−A)/C.
الجزء الثاني — الحمل.
احسب مساحة المقطع العرضي للظهرة.
احسب الحمل الخطّي V0، معتبرًا الظهرة كلّها مغمورة، بحيث يكون وزنها (ρv−ρw)g لكلّ وحدة حجم.
لو عُوّضت الظهرة محليًّا على طريقة إيري بجذرٍ من صخرها الخاص في الوشاح، فما عمق الجذر تحت القمّة؟
لو كان التعويض موضعيًّا، لانتهى الجذر عند سفح الظهرة، على بعد 40km من محورها. تُظهر البيانات الباثيمترية حفرةً تمتدّ حتى 165km. ماذا يُستنتج من ذلك؟
الجزء الثالث — الحفرة الانثنائية. الحفرة مملوءة بالماء (ρi=ρw). يجد مقطع زلزالي عبر الظهرة أعلى القشرة المحيطية عائدًا إلى عمقه الإقليمي على بعد 165km من المحور، وقمّة الانتفاخ الأمامي على بعد 220km من المحور.
تحقّق من أنّ الطفو لكامل الانخفاض، 2w0α(ρm−ρw)g، يوازن V0.
قارن Te بالقيم المرصودة على طول ظهرة هاواي، من 17 إلى 37km. قارن الانخفاض تحت المحور مع جذر إيري في السؤال 8.
الجزء الرابع — الجاذبية فوق الحفرة المحيطة. خذ توافقيًّا واحدًا من الانثناء، طول موجته 2πα وسعته W=500m، مُنحنيًا كلًّا من قاع البحر (القشرة 2900kgm−3 تحت الماء، على عمق 5km تحت مستوى سطح البحر) والموهو (الوشاح تحت القشرة، على عمق 12km تحت مستوى سطح البحر). صفيحة رقيقة ذات كثافة سطحية σ تُنتج الشذوذ 2πGσ مباشرةً فوقها.
احسب طول الموجة والعدد الموجي k للتوافقي.
احسب سعة الشذوذ عند مستوى سطح البحر الناتج عن قاع البحر المنحني بعد الاستمرار إلى الأعلى من عمق 5km.
احسب سعة الموهو المنحني، والمجموع الكلّي.
احسب المجموع عند ارتفاع قمر صناعي، 400km فوق مستوى سطح البحر.
بأيّ عامل يُضعَف توافقي طول موجته 80km، وهو عرض الظهرة، على ذلك الارتفاع؟ وماذا يمكن للقمر الصناعي أن يميّزه؟
إلى أيّ درجة توافقية كروية n≈2πR/λ يناظر التوافقي في السؤال 18 (R=6371km)؟
4. قيمة كليرو، 1/298.61، تقع ضمن 0.12% من القيمة المرصودة 1/298.26؛ أمّا القيمة الهيدروستاتيكية، 1/299.88، فهي أصغر بمقدار Δf=1.82×10−5، وaΔf=116m: يعلو خطّ الاستواء 116m فوق الشكل الهيدروستاتيكي.
5.H=1.08263×10−3/0.3308=0.0032728 (1/305.6).
6.21×80km×6km=240km2=2.40×108m2.
7.V0=(2800−1030)×9.81×2.40×108=4.17×1012Nm−1.
8. إيري: (ρm−ρv)r=(ρv−ρw)h، ومنه r=6×1770/500=21.2km تحت القمّة.
9. ينحني قاع البحر انحناءً كبيرًا إلى ما بعد قاعدة الظهرة: الحمل لا يُعوَّض عمودًا فعمودًا، بل بانثناء الصفيحة عبر منطقة أوسع أربع مرات. للصفيحة صلابة: التعويض الأيزوستاتي الإقليمي.
10.α=165/(3π/4)=70.0km.
11.α=220/π=70.0km: القيمة نفسها، كما تقتضيه صفيحة مرنة واحدة.
16.2×1334×7.00×104×2270×9.81=4.16×1012Nm−1، يساوي V0 ضمن حدود تقريب D.
17.24.7km تقع ضمن نطاق ظهرة هاواي. الصفيحة تغوص 1.3km تحت المحور بدلًا من أن ينمو لها جذر بسمك 21km: إذ يتوزّع الطفو نفسه على 330km من الخندق المحيط.
18.λ=2π×70.0km=440km, k=1/α=1.43×10−5m−1.
19.σ=(2900−1030)×500=9.35×105kgm−2، ومنه 2πGσ=3.92×10−4ms−2=39.2mGal عند قاع البحر؛ مضروبًا في e−k×5km=0.931، 36.5mGal عند مستوى سطح البحر (سالبة حيث تُدفع الصفيحة إلى أسفل).
20.σ=(3300−2900)×500=2.0×105kgm−2: 8.39mGal، مضروبًا في e−k×12km=0.843، 7.07mGal. ويكون المجموع 43.6mGal.
21.e−k×400km=e−5.71=0.0033: 0.14mGal.
22.e−2π×400/80=e−31.4=2.3×10−14: القمر الصناعي لا يرى من الظهرة ذاتها شيئًا، بل الانثناء الواسع فحسب، وحتى هذا يضعف بمقدار ثلاثمئة مرّة؛ ويجب قياس الشذوذ الخاص بالظهرة في البحر.
23.n≈2π×6371/440=91.
24. من المسائل 10 إلى 13، السمك المرن الفعّال للصفيحة تحت ظهرة أورو هو 24.7km.