地球与空间科学 · 第 5 册 · 本科三年级

Ilmu Kebumian dan Antariksa Universitas — Tahun 3

书籍暂未推出简体中文版——当前显示 Bahasa Indonesia 版。

Ilmu Kebumian dan Antariksa Universitas — Tahun 3 · 本科三年级

6Konveksi Mantel dan Gaya Penggerak Lempeng

Pesisir utara Swedia naik dari laut sekitar satu sentimeter per tahun, sepuluh ribu tahun setelah lembaran es yang menutupinya mencair; pantai Teluk Hudson naik sama cepatnya. Tidak ada yang mendorong mereka ke atas selain mantel, yang perlahan mengalir kembali di bawah tempat yang semula ditekan es. Laju ketika tanah memantul kembali menakar viskositas mantel, dan viskositas yang sama juga menentukan seberapa cepat mantel berkonveksi dan lempeng bergerak. Bab ini menanyakan mengapa mantel padat berkonveksi sama sekali, seberapa cepat ia harus mengalir untuk membawa panas Bumi, gaya apa yang menggerakkan lempeng, dan apa yang diungkap pantulan tanah tentang viskositas mantel.

Yang sudah kamu ketahui

Viskositas, aliran merayap fluida kental (persamaan Stokes ∇p=η∇2u\nabla p = \eta\nabla^2\mathbf u, ∇⋅u=0\nabla\cdot\mathbf u = 0), hukum Stokes untuk bola, bilangan Reynolds, dan daya apung terdapat dalam volume universitas buku fisika. Konveksi mantel dan pantulan isostatik dibahas dalam volume sekolah menengah atas (kelas 12); dorongan punggung dan tarikan lempengan di bagian kelas 11. Lapisan batas termal, litosfer termal, dan adiabat mantel terdapat dalam volume Tahun 1, pada bab tentang anggaran panas Bumi; titik panas, plume mantel dan provinsi beku besar dalam bab tentang titik panas. Model pendinginan setengah-ruang terdapat pada Bab 5, lenturan litosfer ada pada Bab 3, serta tomografi seismik pada Bab 1.

6.1 Mengapa mantel berkonveksi

Mantel bersifat padat: ia meneruskan gelombang S. Namun, dalam jutaan tahun, ia merayap seperti fluida dengan viskositas sangat besar, dan lapisan fluida yang dipanaskan dari bawah dan didinginkan dari atas dapat tetap diam, dengan panas melintasinya secara konduksi, atau berpusar. Gumpalan fluida panas berukuran dd, yang lebih ringan daripada lingkungannya sebesar ρα ΔT\rho\alpha\,\Delta T, naik dengan laju berorde ρgα ΔTd2/η\rho g\alpha\,\Delta T d^2/\eta; ia kehilangan panas melalui difusi dalam waktu berorde d2/κd^2/\kappa. Jika ia dapat naik melintasi lapisan sebelum mendingin, lapisan itu berpusar.

Definisi 6.1 (Bilangan Rayleigh)

Untuk lapisan fluida dengan ketebalan dd, massa jenis ρ\rho, koefisien muai termal α\alpha, difusivitas termal κ\kappa, dan viskositas η\eta, di bawah gravitasi gg, dengan selisih suhu ΔT\Delta T antara dasarnya dan puncaknya, bilangan Rayleigh adalah

Ra=ρgα ΔTd3κη,\mathrm{Ra} = \frac{\rho g\alpha\,\Delta T d^3}{\kappa\eta},

rasio waktu difusi d2/κd^2/\kappa terhadap waktu naik karena daya apung η/(ρgα ΔTd)\eta/(\rho g\alpha\,\Delta Td).

Definisi 6.2 (Bilangan Rayleigh kritis)

Yang disebut bilangan Rayleigh kritis Rac\mathrm{Ra}_c suatu lapisan adalah nilai bilangan Rayleigh yang bila terlampaui, gangguan kecil pada keadaan diam tumbuh menjadi konveksi, dan di bawahnya gangguan itu mereda.

Proposisi 6.3 (Permulaan konveksi)

Untuk lapisan yang puncak dan dasarnya bebas dari tegangan geser dan dijaga pada suhu tetap, gangguan dengan bilangan gelombang horizontal kk (dalam satuan 1/d1/d) tumbuh bila Ra>(π2+k2)3/k2\mathrm{Ra} > (\pi^2 + k^2)^3/k^2. Nilai minimum ruas kanan itu, pada k=π/2k = \pi/\sqrt2, memberikan

Rac=27π44≈657.5,\mathrm{Ra}_c = \frac{27\pi^4}{4} \approx 657.5,

dan sel pertama yang muncul selebar 22 d2\sqrt2\,d. Di antara batas kaku, Rac≈1708\mathrm{Ra}_c \approx 1708.

Bukti. Diterima tanpa bukti pada tingkat ini. ∎

Analisis stabilitas linear yang mendasarinya, yakni perturbasi persamaan panas dan Stokes terhadap keadaan konduktif, terdapat dalam pembahasan konveksi di buku fisika; latihan 11 memeriksa nilai minimumnya.

Gambar di bawah memperlihatkan lapisan yang berkonveksi. Nilai untuk mantel, beserta kedalaman batas inti–mantel, d=2891 kmd = 2891\,\mathrm{km}, massa jenis rata-rata mantel, 4476 kg m−34476\,\mathrm{kg}\,\mathrm{m}^{-3}, dan gg antara 9.8 dan 10.7 m s−210.7\,\mathrm{m}\,\mathrm{s}^{-2}, semuanya dari model Bumi seismik PREM; α=3.1×10−5 K−1\alpha = 3.1 \times 10^{-5}\,\mathrm{K}^{-1} dan ΔT=1450 K\Delta T = 1450\,\mathrm{K} dari model lempeng pada Bab 5; κ=10−6 m2 s−1\kappa = 10^{-6}\,\mathrm{m}^{2}\,\mathrm{s}^{-1}; dan η\eta antara 102110^{21} dan 1022 Pa s10^{22}\,\mathrm{Pa}\,\mathrm{s} (Bagian 6.5), memberikan Ra\mathrm{Ra} antara 4.9×1064.9 \times 10^{6} dan 4.9×1074.9 \times 10^{7}, sekitar tujuh ribu hingga tujuh puluh ribu kali nilai kritis. Mantel tak mungkin tidak berkonveksi.

Konveksi dalam lapisan yang dipanaskan dari bawah dan didinginkan dari atas. Panas melintasi lapisan batas yang tipis, setebal , secara konduksi; di antaranya, bagian dalam yang tercampur baik membawa panas dengan aliran. Lembaran dingin lepas dari lapisan atas dan tenggelam; yang panas naik dari bawah.
Konveksi dalam lapisan yang dipanaskan dari bawah dan didinginkan dari atas. Panas melintasi lapisan batas yang tipis, setebal δ\delta, secara konduksi; di antaranya, bagian dalam yang tercampur baik membawa panas dengan aliran. Lembaran dingin lepas dari lapisan atas dan tenggelam; yang panas naik dari bawah.

Mantel juga dipanaskan dari dalam, oleh peluruhan unsur-unsur radioaktif, dan dari bawah, oleh inti. Lapisan yang hanya dipanaskan dari dalam tidak memiliki lapisan batas panas di dasarnya: alirannya digerakkan oleh lembaran dingin yang tenggelam dari atas, sedangkan aliran naik berupa arus balik yang lebar dan pasif. Lempengan termasuk lembaran semacam itu, dan plume membutuhkan panas dari inti.

6.2 Lapisan batas dan kekuatan konveksi

Definisi 6.4 (Bilangan Nusselt)

Yang disebut bilangan Nusselt suatu lapisan yang berkonveksi adalah aliran panas yang melewatinya dibagi aliran panas yang dibawa konduksi saja pada lapisan yang sama, Nu=q d/(k ΔT)\mathrm{Nu} = q\,d/(k\,\Delta T).

Teorema 6.5 (Penskalaan lapisan batas konveksi)

Dalam lapisan yang berkonveksi pada Ra≫Rac\mathrm{Ra} \gg \mathrm{Ra}_c, jika setiap lapisan batas termal tumbuh hingga bilangan Rayleigh-nya sendiri, yang dihitung dari tebalnya δ\delta dan beda suhu ΔT/2\Delta T/2 melintasinya, mencapai Rac\mathrm{Ra}_c, maka

δ=d(2RacRa)1/3,Nu=12(Ra2Rac)1/3,u=κd(Ra2Rac)2/3,\delta = d\Bigl(\frac{2\mathrm{Ra}_c}{\mathrm{Ra}}\Bigr)^{1/3},\qquad \mathrm{Nu} = \frac12\Bigl(\frac{\mathrm{Ra}}{2\mathrm{Ra}_c}\Bigr)^{1/3}, \qquad u = \frac\kappa d\Bigl(\frac{\mathrm{Ra}}{2\mathrm{Ra}_c}\Bigr)^{2/3},

Di sini uu adalah skala kecepatan aliran, yang didefinisikan dengan mensyaratkan bahwa lapisan tumbuh hingga δ\delta selagi fluida melintasi sel. Inilah yang disebut penskalaan lapisan batas konveksi.

Bukti. Bagian dalam tercampur baik, pada suhu rata-rata, sehingga setiap lapisan batas membawa setengah dari ΔT\Delta T. Bilangan Rayleigh-nya adalah ρgα(ΔT/2)δ3/(κη)=Ra (δ/d)3/2\rho g\alpha(\Delta T/2) \delta^3/(\kappa\eta) = \mathrm{Ra}\,(\delta/d)^3/2; dengan menyamakannya dengan Rac\mathrm{Ra}_c memberikan δ\delta. Aliran panas dihantarkan melintasi lapisan, q=k(ΔT/2)/δq = k(\Delta T/2)/\delta, sehingga Nu=qd/(kΔT)=d/(2δ)\mathrm{Nu} = qd/(k\Delta T) = d/(2\delta). Lapisan batas tumbuh secara difusi hingga tebal δ\delta dalam waktu δ2/κ\delta^2/\kappa; fluida di bagian atas melintasi sel yang lebarnya berorde dd, dalam waktu d/ud/u. Dengan menyamakan kedua waktu itu, u=κd/δ2=(κ/d)(d/δ)2u = \kappa d/\delta^2 = (\kappa/d)(d/\delta)^2. ∎

Aliran panas hanya tumbuh sebagai akar pangkat tiga dari Ra\mathrm{Ra}: viskositas seratus kali lebih rendah, panas kurang dari lima kali lebih banyak. Untuk mantel dengan Ra=4.9×107\mathrm{Ra} = 4.9 \times 10^{7}, penskalaan memberikan δ≈87 km\delta \approx 87\,\mathrm{km}, Nu≈17\mathrm{Nu} \approx 17, aliran panas sebesar Nu kΔT/d≈26 mW m−2\mathrm{Nu}\,k\Delta T/d \approx 26\,\mathrm{mW}\,\mathrm{m}^{-2} dengan k=3.1 W m−1 K−1k = 3.1\,\mathrm{W}\,\mathrm{m}^{-1}\,\mathrm{K}^{-1}, dan u≈1.2 cm/yru \approx 1.2\,\mathrm{cm}/\mathrm{yr} (gambar di bawah). Tebal lapisan batas sama dengan lempeng samudra pada Bab 5, 95 km95\,\mathrm{km}; kecepatan akar rata-rata kuadrat lempeng masa kini berkisar dari 1.1 cm/yr1.1\,\mathrm{cm}/\mathrm{yr} (Amerika Selatan) hingga 7.5 cm/yr7.5\,\mathrm{cm}/\mathrm{yr} (Cocos), 4.4 cm/yr4.4\,\mathrm{cm}/\mathrm{yr} untuk seluruh litosfer. Aliran panas rata-rata Bumi, 92 mW m−292\,\mathrm{mW}\,\mathrm{m}^{-2}, adalah tiga sampai empat kali hasil penskalaan: panas yang dihasilkan di dalam mantel, yang tidak diperhitungkan lapisan yang dipanaskan dari bawah, serta konstanta berorde satu yang diabaikan penskalaan; keduanya ikut berperan. Untuk hukum dengan satu masukan, orde besarannya tepat: lempeng merupakan lapisan batas atas mantel dengan Ra∼107\mathrm{Ra} \sim 10^7.

Penskalaan konveksi. (a) Di bawah Ra_c lapisan menghantarkan panas secara konduksi (Nu = 1); jauh di atasnya, hukum lapisan batas, yang digambar dari Ra = 16\, Ra_c 1.1 × 104, tempat hukum itu memberikan Nu = 1 (di antaranya, konveksi berangsur mulai). (b) Skala kecepatan untuk kedalaman mantel dan = 10-6\, m2\, s-1, dengan rentang kecepatan lempeng masa kini (garis titik-strip). Garis putus-putus membatasi bilangan Rayleigh mantel untuk viskositas 1022 dan 1021\, Pa\, s.
Penskalaan konveksi. (a) Di bawah Rac\mathrm{Ra}_c lapisan menghantarkan panas secara konduksi (Nu=1\mathrm{Nu} = 1); jauh di atasnya, hukum lapisan batas, yang digambar dari Ra=16 Rac≈1.1×104\mathrm{Ra} = 16\,\mathrm{Ra}_c \approx 1.1 \times 10^{4}, tempat hukum itu memberikan Nu=1\mathrm{Nu} = 1 (di antaranya, konveksi berangsur mulai). (b) Skala kecepatan untuk kedalaman mantel dan κ=10−6 m2 s−1\kappa = 10^{-6}\,\mathrm{m}^{2}\,\mathrm{s}^{-1}, dengan rentang kecepatan lempeng masa kini (garis titik-strip). Garis putus-putus membatasi bilangan Rayleigh mantel untuk viskositas 102210^{22} dan 1021 Pa s10^{21}\,\mathrm{Pa}\,\mathrm{s}.

6.3 Gaya-gaya pada lempeng

Lempeng bergerak dengan kecepatan tetap: gaya-gaya padanya berada dalam keseimbangan. Tiga gaya mendorong atau menahannya dari tepi dan dasarnya (gambar di bawah). Gaya yang diberikan lempengan yang tenggelam pada lempengnya, yakni tarikan lempengan, berasal dari kelebihan massa jenis lempengan dingin di dalam mantel yang lebih panas. Dorongan punggung, yakni dorongan punggungan yang tinggi dan panas terhadap lempeng yang lebih dingin di sebelahnya, diturunkan dari model pendinginan setengah-ruang pada Bab 5.

Gaya-gaya pada lempeng samudra: dorongan punggung dan tarikan lempengan mendorong lempeng (merah); seretan di dasarnya, tempat lempeng bergerak lebih cepat daripada mantel di bawahnya, dan hambatan kental mantel terhadap lempengan yang tenggelam menahannya (biru).
Gaya-gaya pada lempeng samudra: dorongan punggung dan tarikan lempengan mendorong lempeng (merah); seretan di dasarnya, tempat lempeng bergerak lebih cepat daripada mantel di bawahnya, dan hambatan kental mantel terhadap lempengan yang tenggelam menahannya (biru).

Proposisi 6.6 (Gaya dorongan punggung)

Pada model pendinginan setengah-ruang, dengan suhu mantel TmT_m, suhu permukaan T0T_0 dan kompensasi isostatik, gaya horizontal per satuan panjang punggungan yang dikerjakan kolom di punggungan pada lempeng berumur tt, yakni gaya dorongan punggung, adalah

Frp=gρmα(Tm−T0)κt[1+2ρmα(Tm−T0)π(ρm−ρw)].F_{rp} = g\rho_m\alpha(T_m - T_0)\kappa t\Bigl[1 + \frac{2\rho_m\alpha(T_m - T_0)}{\pi(\rho_m - \rho_w)}\Bigr].

Bukti. Misalkan yy kedalaman di bawah dasar laut di punggungan, ww penurunan dasar laut berumur tt, θ(z)=Tm−T(z)\theta(z) = T_m - T(z) defisit suhu pada kedalaman z=y−wz = y - w dalam lempengnya, dan Θ(z)=∫0zθ\Theta(z) = \int_0^z\theta. Kolom punggungan memuat mantel bermassa jenis ρm\rho_m dari y=0y = 0; kolom tua memuat air hingga kedalaman ww, lalu lempeng bermassa jenis ρm(1+αθ)\rho_m(1 + \alpha\theta). Tekanan litostatik berbeda sebesar Pr−Pt=(ρm−ρw)gyP_r - P_t = (\rho_m - \rho_w)gy untuk y<wy < w, dan di bawahnya,

Pr−Pt=(ρm−ρw)gw−ρmgα Θ(z)=ρmgα(Θ(∞)−Θ(z)),P_r - P_t = (\rho_m - \rho_w)gw - \rho_mg\alpha\,\Theta(z) = \rho_mg\alpha\bigl(\Theta(\infty) - \Theta(z)\bigr),

karena isostasi, (ρm−ρw)w=ρmαΘ(∞)(\rho_m - \rho_w)w = \rho_m\alpha\Theta(\infty). Gaya adalah integral selisih tekanan terhadap kedalaman; dengan mengintegralkan bagian kedua secara parsial,

Frp=12(ρm−ρw)gw2+ρmgα∫0∞z θ(z) dz.F_{rp} = \tfrac12(\rho_m - \rho_w)gw^2 + \rho_mg\alpha\int_0^\infty z\,\theta(z)\,\mathrm dz .

Untuk setengah-ruang, θ=(Tm−T0)erfc⁡(z/2κt)\theta = (T_m - T_0)\operatorname{erfc}(z/2\sqrt{\kappa t}) dan ∫0∞xerfc⁡x dx=1/4\int_0^\infty x\operatorname{erfc}x\,\mathrm dx = 1/4, sehingga integralnya bernilai (Tm−T0)κt(T_m - T_0)\kappa t; dan w=2ρmα(Tm−T0)κt/π/(ρm−ρw)w = 2\rho_m\alpha(T_m - T_0)\sqrt{\kappa t/\pi}/(\rho_m - \rho_w). Menjumlahkan kedua suku memberikan rumus. ∎

Dorongan punggung tumbuh secara linear terhadap umur: ia bukan dorongan di sumbu, melainkan gradien tekanan yang tersebar di seluruh lempeng yang mendingin. Dengan nilai-nilai dari Bab 5, faktor dalam kurung bernilai 1.04 dan Frp=3.8×1012 N m−1F_{rp} = 3.8 \times 10^{12}\,\mathrm{N}\,\mathrm{m}^{-1} pada 100 Ma100\,\mathrm{Ma}. Pada model lempeng, gaya ini berhenti bertambah seiring kedalaman, sekitar 2.3×1012 N m−12.3 \times 10^{12}\,\mathrm{N}\,\mathrm{m}^{-1} (gambar di bawah). Tarikan lempengan lebih besar: lempengan setebal 95 km95\,\mathrm{km}, rata-rata (Tm−T0)/2(T_m - T_0)/2 lebih dingin daripada mantel, sehingga lebih padat sebesar ρmα(Tm−T0)/2≈74 kg m−3\rho_m\alpha(T_m - T_0)/2 \approx 74\,\mathrm{kg}\,\mathrm{m}^{-3}, menggantung hingga 670 km670\,\mathrm{km}, beratnya 4.6×1013 N m−14.6 \times 10^{13}\,\mathrm{N}\,\mathrm{m}^{-1} lebih besar daripada mantel yang dipindahkannya, dua belas kali dorongan punggung. Sebagian besar berat itu dipakai untuk mengatasi mantel kental yang menahan penurunan lempengan, dan hanya sebagian yang mencapai lempeng.

Gaya dorongan punggung terhadap umur lempeng, dari kedua model pendinginan pada , dengan integrasi numerik selisih tekanan (kurva setengah-ruang adalah bentuk tertutup dari proposisi).
Gaya dorongan punggung terhadap umur lempeng, dari kedua model pendinginan pada Bab 5, dengan integrasi numerik selisih tekanan (kurva setengah-ruang adalah bentuk tertutup dari proposisi).

Definisi 6.7 (Seretan basal)

Yang disebut seretan basal adalah gaya geser per satuan luas yang dikerjakan mantel di bawah lempeng pada dasarnya. Untuk lempeng yang bergerak dengan kecepatan uu di atas lapisan setebal hh dengan viskositas η\eta yang dasarnya diam, gayanya ηu/h\eta u/h, berlawanan dengan gerak; gaya ini mendorong lempeng di tempat mantel di bawahnya bergerak lebih cepat daripada lempeng.

Gaya-gaya mana yang berpengaruh dapat dibaca dari lempeng itu sendiri. Keseimbangan torsi pada setiap lempeng, yang dicocokkan dengan gerak lempeng, menyatakan bahwa gaya-gaya yang bekerja pada lempengan yang tenggelam mengendalikan kecepatan lempeng samudra dan sekitar sepuluh kali lebih kuat daripada gaya lain, dan bahwa hambatan di dasar lebih besar di bawah benua daripada di bawah samudra. Kecepatannya sesuai: lempeng tercepat, Cocos pada 75 mm/yr75\,\mathrm{mm}/\mathrm{yr}, Nazca pada 70 dan Pasifik pada 66 mm/yr66\,\mathrm{mm}/\mathrm{yr}, adalah lempeng samudra yang tepinya tenggelam di palung panjang; yang paling lambat membawa benua besar, Eurasia pada 23, Amerika Utara pada 20, Amerika Selatan pada 11 mm/yr11\,\mathrm{mm}/\mathrm{yr}.

6.4 Plume

Lembaran panas yang naik dari dasar mantel berkumpul menjadi kolom: plume mantel yang dibahas di volume Tahun 1, dengan kepala besar dan ekor tipis (gambar di bawah). Tonjolan titik panas, yaitu naiknya dasar laut yang lebar di sekitar dan di hilir gunung api, mengukur seberapa besar daya apung yang dibawa plume-nya.

Definisi 6.8 (Fluks daya apung plume)

Yang disebut fluks daya apung plume BB adalah defisit massa, relatif terhadap mantel di sekitarnya, yang dibawa plume ke atas per satuan waktu; untuk plume dengan fluks volume QvQ_v dan kelebihan temperatur ΔT\Delta T, B=ρα ΔT QvB = \rho\alpha\,\Delta T\,Q_v.

Proposisi 6.9 (Fluks daya apung dari tonjolan)

Jika lempeng yang bergerak dengan kecepatan uu di atas plume membawa tonjolan dengan luas penampang AA melintasi jejak, yang ditopang secara isostatik oleh material plume di bawahnya, maka fluks daya apung plume adalah

B=(ρm−ρw) A u.B = (\rho_m - \rho_w)\,A\,u .

Bukti. Menurut isostasi, sebuah irisan tonjolan dengan panjang satuan sepanjang jejak, yang mengangkat dasar laut di atas luas penampang AA dan menggantikan air dengan batuan, ditopang oleh defisit massa (ρm−ρw)A(\rho_m - \rho_w)A di kolom di bawahnya. Lempeng membawa irisan itu pergi dengan kecepatan uu, sehingga dalam keadaan tunak plume harus menyediakan defisit itu untuk panjang uu tonjolan baru per satuan waktu. ∎

Dari penampang melintang tonjolan, fluks daya apung titik panas berkisar dalam faktor dua puluh: Hawaii, yang terbesar, 8.7 Mg/s8.7\,\mathrm{Mg}/\mathrm{s}; Islandia 1.4 Mg/s1.4\,\mathrm{Mg}/\mathrm{s}; 37 titik panas bersama-sama, sekitar 50 Mg/s50\,\mathrm{Mg}/\mathrm{s}. Setiap kilogram daya apung setara dengan cp/αc_p/\alpha joule panas berlebih, sehingga fluks panasnya cpB/αc_pB/\alpha; dengan cp=k/(ρmκ)=1170 J kg−1 K−1c_p = k/(\rho_m\kappa) = 1170\,\mathrm{J}\,\mathrm{kg}^{-1}\,\mathrm{K}^{-1} dan α=3.1×10−5 K−1\alpha = 3.1 \times 10^{-5}\,\mathrm{K}^{-1}, ke-37 titik panas membawa sekitar 1.9 TW1.9\,\mathrm{TW}, atau 3.7 mW m−23.7\,\mathrm{mW}\,\mathrm{m}^{-2} bila dibagi ke seluruh permukaan Bumi 5.1×1014 m25.1 \times 10^{14}\,\mathrm{m}^{2}, mendekati 4 mW m−24\,\mathrm{mW}\,\mathrm{m}^{-2} yang dipublikasikan bersama fluks-fluks itu: hanya beberapa persen dari 47 TW47\,\mathrm{TW} yang hilang dari Bumi. Plume adalah mitra kecil dalam anggaran panas; lempenglah yang membawa sisanya.

Plume muda (kiri) naik sebagai kepala besar yang diisi oleh ekor tipis; ketika kepala menyebar di bawah litosfer, permukaan naik dan basal banjir dapat meletus. Plume tua (kanan) adalah ekor di bawah lempeng yang bergerak, memasok tonjolan dan jejak titik panas. Tidak sesuai skala.
Plume muda (kiri) naik sebagai kepala besar yang diisi oleh ekor tipis; ketika kepala menyebar di bawah litosfer, permukaan naik dan basal banjir dapat meletus. Plume tua (kanan) adalah ekor di bawah lempeng yang bergerak, memasok tonjolan dan jejak titik panas. Tidak sesuai skala.

Kepala plume adalah bola mantel berdaya apung yang naik melalui mantel kental: hukum Stokes dari buku fisika, v=2Δρ ga2/(9η)v = 2\Delta\rho\,ga^2/(9\eta), memberi kecepatannya dalam orde besaran, sebagian kecil meter per tahun untuk kepala dengan lebar beberapa ratus kilometer (latihan 8), sehingga ia melintasi mantel dalam jutaan tahun.

6.5 Penyesuaian isostatik glasial dan viskositas mantel

Definisi 6.10 (Penyesuaian isostatik glasial, viskositas mantel)

Yang disebut penyesuaian isostatik glasial (GIA) adalah deformasi lambat Bumi padat, dan perubahan medan gravitasinya serta permukaan laut, ketika ia kembali menuju kesetimbangan isostatik setelah lembaran es tumbuh atau mencair. Viskositas mantel adalah viskositas efektif η\eta yang menghubungkan rayapan mantel yang lambat dan tunak dengan tegangan yang menggerakkannya, dalam rentang waktu ribuan tahun atau lebih.

Ketika lembaran es mencair, cekungan yang ditinggalkannya berada di luar kesetimbangan isostatik: di bawahnya, tekanan mantel lebih rendah daripada di sekitarnya, dan mantel mengalir masuk untuk mengangkat kembali dasarnya. Hari ini pesisir Teluk Bothnia naik hingga 10.3 mm/yr10.3\,\mathrm{mm}/\mathrm{yr}, di dekat Umeå; Teluk Hudson sekitar 10 mm/yr10\,\mathrm{mm}/\mathrm{yr}, tempat es paling tebal, sedangkan daratan di selatan Great Lakes turun sebesar 1 hingga 2 mm/yr2\,\mathrm{mm}/\mathrm{yr}, ketika tonjolan yang didorong naik oleh es di sekitarnya menyusut. Di pesisir, garis pantai lama kini berada di atas laut, yaitu pantai terangkat dalam kotak lapangan di bawah.

Teorema 6.11 (Relaksasi setengah-ruang kental)

Cekungan w(x,t)=wm(t)cos⁡(2πx/λ)w(x, t) = w_m(t)\cos(2\pi x/\lambda) pada permukaan suatu setengah-ruang dengan viskositas seragam η\eta dan massa jenis ρ\rho, yang dibiarkan tanpa beban, meluruh sebagai wm(t)=wm(0)e−t/τw_m(t) = w_m(0)e^{-t/\tau} dengan

τ=4πηρgλ.\tau = \frac{4\pi\eta}{\rho g\lambda} .

Bukti. Ambil zz ke arah bawah, k=2π/λk = 2\pi/\lambda, dan cari aliran Stokes yang meluruh terhadap kedalaman, dengan kecepatan vertikal (ke bawah) v=(C+Dkz)e−kzcos⁡kxv = (C + Dkz)e^{-kz}\cos kx. Kondisi tak-mampumampat memberikan kecepatan horizontal u=−(D−C−Dkz)e−kzsin⁡kxu = -(D - C - Dkz)e^{-kz}\sin kx. Fungsi-fungsi ini memenuhi ∇4=0\nabla^4 = 0, seperti halnya setiap aliran Stokes dalam dua dimensi. Permukaan bebas tidak memikul tegangan geser: ∂u/∂z+∂v/∂x=2k(D−C)sin⁡kx=0\partial u/\partial z + \partial v/\partial x = 2k(D - C)\sin kx = 0 pada z=0z = 0, sehingga D=CD = C, dan

v=C(1+kz)e−kzcos⁡kx,u=Ckze−kzsin⁡kx,p=2ηCke−kzcos⁡kx,v = C(1 + kz)e^{-kz}\cos kx,\quad u = Ckze^{-kz}\sin kx,\quad p = 2\eta Cke^{-kz}\cos kx,

tekanan diperoleh dari ∂p/∂x=η∇2u=−2ηCk2e−kzsin⁡kx\partial p/\partial x = \eta\nabla^2u = -2\eta Ck^2e^{-kz}\sin kx (dan ∂p/∂z=η∇2v\partial p/\partial z = \eta\nabla^2v terverifikasi). Pada permukaan ∂v/∂z=0\partial v/\partial z = 0, sehingga tegangan normal di sana adalah −p-p. Di bawah cekungan ww, berat batuan yang hilang menurunkan tekanan sebesar ρgw\rho gw: p(0)=−ρgwp(0) = -\rho gw, sehingga 2ηCk=−ρgwm2\eta Ck = -\rho gw_m. Permukaan bergerak bersama fluida: dwm/dt=v(0)=C=−ρgwm/(2ηk)\mathrm dw_m/\mathrm dt = v(0) = C = -\rho gw_m/(2\eta k), yang solusinya meluruh dengan τ=2ηk/(ρg)=4πη/(ρgλ)\tau = 2\eta k/(\rho g) = 4\pi\eta/(\rho g\lambda). ∎

Definisi 6.12 (Waktu relaksasi isostatik)

Yang disebut waktu relaksasi isostatik τ\tau dari suatu beban dengan skala horizontal λ\lambda adalah waktu ketika simpangannya dari kesetimbangan isostatik meluruh dengan faktor ee.

Setengah-ruang seragam merelaksasi cekungan lebar lebih cepat daripada cekungan sempit, dan aliran di bawah beban dengan panjang gelombang λ\lambda meluruh terhadap kedalaman sebagai e−kze^{-kz}, hingga kedalaman λ/2π\lambda/2\pi: lembaran es besar mengaduk mantel jauh ke dalam, danau kecil hanya bagian atasnya. Untuk cekungan selebar 3000 km3000\,\mathrm{km} dan η=1021 Pa s\eta = 10^{21}\,\mathrm{Pa}\,\mathrm{s}, dengan ρ=3300 kg m−3\rho = 3300\,\mathrm{kg}\,\mathrm{m}^{-3}, τ=4.1 ka\tau = 4.1\,\mathrm{ka} (gambar di bawah): setelah sepuluh ribu tahun, e−2.4=9 %e^{-2.4} = 9\,\% dari cekungan yang tersisa, dan karena w˙=−w/τ\dot w = -w/\tau, pengangkatan sebesar 10 mm/yr10\,\mathrm{mm}/\mathrm{yr} berarti masih ada sekitar 40 m40\,\mathrm{m} kenaikan yang akan datang.

Relaksasi setengah-ruang kental seragam (= 3300\, kg\, m-3, g = 9.8\, m\, s-2) untuk tiga viskositas. (a) Fraksi cekungan yang tersisa seiring waktu. (b) Waktu relaksasi 4π /( g ), lebih panjang untuk beban sempit.
Relaksasi setengah-ruang kental seragam (ρ=3300 kg m−3\rho = 3300\,\mathrm{kg}\,\mathrm{m}^{-3}, g=9.8 m s−2g = 9.8\,\mathrm{m}\,\mathrm{s}^{-2}) untuk tiga viskositas. (a) Fraksi cekungan yang tersisa seiring waktu. (b) Waktu relaksasi 4πη/(ρgλ)4\pi\eta/(\rho g\lambda), lebih panjang untuk beban sempit.

Mantel sungguhan tidak seragam, dan pantulannya dibaca dengan model yang melapisinya dan menambahkan litosfer elastik, yang menahan beban sempit melalui lenturan. Taksiran pertama, dari pantulan Fennoscandia, memberikan 1021 Pa s10^{21}\,\mathrm{Pa}\,\mathrm{s}, dan peluruhan pengangkatan di Sungai Angerman di Swedia dan di Teluk Richmond di Teluk Hudson memberikan rata-rata 0.65 hingga 1.1×1021 Pa s1.1 \times 10^{21}\,\mathrm{Pa}\,\mathrm{s} dari dasar litosfer hingga sekitar 1400 km1400\,\mathrm{km}, dengan mantel atas sekitar sepuluh kali lebih rendah viskositasnya daripada mantel bawah. Inversi laju pengangkatan GNSS Fennoscandia memberikan viskositas mantel atas sebesar 5×10205 \times 10^{20} hingga 1021 Pa s10^{21}\,\mathrm{Pa}\,\mathrm{s} dalam satu studi dan 3.4×10203.4 \times 10^{20} hingga 5.0×1020 Pa s5.0 \times 10^{20}\,\mathrm{Pa}\,\mathrm{s} dalam studi lain, serta viskositas mantel bawah sebesar 5×10215 \times 10^{21} hingga 5×1022 Pa s5 \times 10^{22}\,\mathrm{Pa}\,\mathrm{s}, dan litosfer setebal 93 hingga 110 km110\,\mathrm{km}. Taksiran-taksiran berbeda dengan faktor dua atau tiga, tetapi keduanya sepakat pada orde magnitudo dan pada peningkatan seiring kedalaman.

Definisi 6.13 (Waktu Maxwell)

Yang disebut waktu Maxwell suatu bahan dengan viskositas η\eta dan modulus geser μ\mu adalah τM=η/μ\tau_M = \eta/\mu: beban yang diterapkan dalam waktu jauh lebih singkat daripada τM\tau_M menemui padatan elastik; beban yang diterapkan jauh lebih lama menemui fluida kental.

Dengan μ=ρvS2=7.7×1010 Pa\mu = \rho v_S^2 = 7.7 \times 10^{10}\,\mathrm{Pa} pada 300 km300\,\mathrm{km} dalam PREM dan η=1021 Pa s\eta = 10^{21}\,\mathrm{Pa}\,\mathrm{s}, τM=1.3×1010 s\tau_M = 1.3 \times 10^{10}\,\mathrm{s}, sekitar empat ratus tahun. Gelombang seismik, yang berlangsung beberapa detik, melintasi mantel elastik; zaman es, dan konveksi, mendeformasi fluida. Batuan yang sama adalah keduanya, bergantung pada seberapa lama ia didorong.

Di lapangan — Gisik terangkat

Gisik terangkat di Ballyhillin, Malin Head, di utara Irlandia: punggungan kerikil garis pantai lama terangkat di atas laut seiring naiknya daratan setelah lembaran esnya mencair. Foto A. F. Borchert, CC BY-SA 4.0.

Di pesisir daratan yang dahulu tertutup es, seorang geomorfolog menemukan pekerjaan lama laut yang terdampar di pedalaman: punggungan kerikil membulat, takik, dan platform datar yang dipahat ombak, kini berada di atas jangkauan pasang surut tertinggi masa kini, kadang berupa tangga beberapa tingkat. Setiap tingkat adalah garis pantai ketika laut berada di sana; daratan sejak itu terangkat. Cangkang dan kayu hanyut di gisik lama diberi penanggalan dengan radiokarbon, dan ketinggian aras yang diberi penanggalan, diratakan ke titik acuan, memberikan pengangkatan terhadap waktu. Ditelusuri sepanjang pesisir, garis pantai lama ternyata miring, paling tinggi di tempat es paling tebal: kemiringannya memetakan bentuk pantulan itu.

6.6 Konveksi seluruh mantel dan struktur dalamnya

Apakah mantel berkonveksi sebagai satu lapisan, atau dua, di atas dan di bawah 670 km670\,\mathrm{km} diskontinuitas? Tomografi (Bab 1) menjawab sebagian: sebagian lempengan mendatar di atas diskontinuitas, sebagian lainnya menembusnya ke mantel bawah, sehingga massa melintasinya dan kedua lapisan tidak tertutup satu sama lain. Peningkatan viskositas terhadap kedalaman membuat mantel bawah mengalir lebih lamban, dan lempengan melambat serta menebal saat memasukinya.

Di dasar mantel, tomografi memperlihatkan dua wilayah luas berkecepatan gelombang S rendah, di bawah Afrika dan di bawah Pasifik, yakni provinsi kecepatan geser rendah berskala besar (LLSVPs). Keduanya kemungkinan besar berbeda secara kimiawi dari mantel di atasnya dan lebih padat, berupa timbunan, bukan plume panas. Provinsi beku besar dan kimberlit dari beberapa ratus juta tahun terakhir meletus terutama di atas tepinya, seakan-akan plume naik dari tepi timbunan, tempat lapisan batas panas di dasar mantel tersapu menuju timbunan itu.

6.7 Latihan

Latihan 6.1 ★

Hitung bilangan Rayleigh dari lapisan mantel atas setebal 670 km670\,\mathrm{km} dengan ρ=3500 kg m−3\rho = 3500\,\mathrm{kg}\,\mathrm{m}^{-3}, g=9.9 m s−2g = 9.9\,\mathrm{m}\,\mathrm{s}^{-2}, α=3×10−5 K−1\alpha = 3 \times 10^{-5}\,\mathrm{K}^{-1}, ΔT=1000 K\Delta T = 1000\,\mathrm{K}, κ=10−6 m2 s−1\kappa = 10^{-6}\,\mathrm{m}^{2}\,\mathrm{s}^{-1} dan η=5×1020 Pa s\eta = 5 \times 10^{20}\,\mathrm{Pa}\,\mathrm{s}. Berapa kali lipat nilai kritisnya untuk batas bebas-geser?

Solusi

Solusi Latihan 6.1.

Ra=3500×9.9×3×10−5×1000×(6.7×105)3/(10−6×5×1020)=6.25×105\mathrm{Ra} = 3500\times9.9\times3 \times 10^{-5}\times1000\times(6.7 \times 10^{5})^3/ (10^{-6}\times5 \times 10^{20}) = 6.25 \times 10^{5}, sekitar 950 kali bilangan Rayleigh kritis 657.5: lapisan berkonveksi kuat.

Latihan 6.2 ★

Hitung waktu Maxwell mantel dengan viskositas 1021 Pa s10^{21}\,\mathrm{Pa}\,\mathrm{s} dan modulus geser 7.7×1010 Pa7.7 \times 10^{10}\,\mathrm{Pa}, dalam tahun. Apakah mantel merespons secara elastis atau secara viskos terhadap pasang surut Bumi dua kali sehari, terhadap mencairnya lembaran es, dan terhadap tenggelamnya lempengan?

Solusi

Solusi Latihan 6.2.

τM=1021/7.7×1010=1.30×1010 s\tau_M = 10^{21}/7.7 \times 10^{10} = 1.30 \times 10^{10}\,\mathrm{s}, sekitar 412 yr412\,\mathrm{yr}. Pasang surut, yang berlangsung berjam-jam, bertemu mantel elastis; pencairan lembaran es, selama ribuan tahun, dan tenggelamnya lempengan selama jutaan tahun, bertemu mantel viskos.

Latihan 6.3 ★

Dengan η=1021 Pa s\eta = 10^{21}\,\mathrm{Pa}\,\mathrm{s}, ρ=3300 kg m−3\rho = 3300\,\mathrm{kg}\,\mathrm{m}^{-3} dan g=9.8 m s−2g = 9.8\,\mathrm{m}\,\mathrm{s}^{-2}, hitung waktu relaksasi setengah-ruang seragam untuk beban dengan panjang gelombang 1000 dan 3000 km3000\,\mathrm{km}. Mana yang pulih lebih cepat?

Solusi

Solusi Latihan 6.3.

τ=4π×1021/(3300×9.8λ)\tau = 4\pi\times10^{21}/(3300\times9.8\lambda): 3.89×1011 s3.89 \times 10^{11}\,\mathrm{s} = 12.3 ka12.3\,\mathrm{ka} untuk 1000 km1000\,\mathrm{km}, dan 4.1 ka4.1\,\mathrm{ka} untuk 3000 km3000\,\mathrm{km}. Beban yang lebih luas pulih lebih cepat.

Latihan 6.4 ★

Dengan g=9.8 m s−2g = 9.8\,\mathrm{m}\,\mathrm{s}^{-2}, ρm=3300 kg m−3\rho_m = 3300\,\mathrm{kg}\,\mathrm{m}^{-3}, ρw=1030 kg m−3\rho_w = 1030\,\mathrm{kg}\,\mathrm{m}^{-3}, α=3.1×10−5 K−1\alpha = 3.1 \times 10^{-5}\,\mathrm{K}^{-1}, Tm−T0=1450 KT_m - T_0 = 1450\,\mathrm{K}, κ=0.8×10−6 m2 s−1\kappa = 0.8 \times 10^{-6}\,\mathrm{m}^{2}\,\mathrm{s}^{-1} dan 1 Ma1\,\mathrm{Ma} =3.156×1013 s= 3.156 \times 10^{13}\,\mathrm{s}, hitung gaya dorongan punggung pada 50 Ma50\,\mathrm{Ma}. Berapa kali lebih kecil dibanding tarikan lempengan sebesar 4.6×1013 N m−14.6 \times 10^{13}\,\mathrm{N}\,\mathrm{m}^{-1}?

Solusi

Solusi Latihan 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×1012 N m−1g\rho_m\alpha(T_m - T_0)\kappa t = 9.8\times3300\times3.1 \times 10^{-5}\times1450 \times0.8 \times 10^{-6}\times50\times3.156 \times 10^{13} = 1.835 \times 10^{12}\,\mathrm{N}\,\mathrm{m}^{-1}; faktor dalam kurung bernilai 1+2×3300×3.1×10−5×1450/(π×2270)=1.04161 + 2\times3300\times3.1 \times 10^{-5}\times1450/(\pi\times2270) = 1.0416, sehingga Frp=1.91×1012 N m−1F_{rp} = 1.91 \times 10^{12}\,\mathrm{N}\,\mathrm{m}^{-1}, 24 kali lebih kecil daripada tarikan lempengan.

Latihan 6.5 ★★

Sebuah gumpalan mantel selebar 10 km10\,\mathrm{km}, 300 K300\,\mathrm{K} lebih panas daripada sekelilingnya, naik dengan laju ρgαΔTd2/η\rho g\alpha\Delta Td^2/\eta (ambil ρ=3300 kg m−3\rho = 3300\,\mathrm{kg}\,\mathrm{m}^{-3}, g=9.8 m s−2g = 9.8\,\mathrm{m}\,\mathrm{s}^{-2}, α=3×10−5 K−1\alpha = 3 \times 10^{-5}\,\mathrm{K}^{-1}, η=1021 Pa s\eta = 10^{21}\,\mathrm{Pa}\,\mathrm{s}, κ=10−6 m2 s−1\kappa = 10^{-6}\,\mathrm{m}^{2}\,\mathrm{s}^{-1}). (a) Hitung waktu naik melintasi ukurannya sendiri dan waktu difusinya d2/κd^2/\kappa, lalu tunjukkan bahwa rasionya sama dengan bilangan Rayleigh-nya. (b) Apakah ia naik atau mendingin lebih dulu? (c) Seberapa besar ukurannya agar bilangan Rayleigh-nya mencapai 657.5?

Solusi

Solusi Latihan 6.5.

(a) v=3300×9.8×3×10−5×300×(104)2/1021=2.91×10−11 m s−1v = 3300\times9.8\times3 \times 10^{-5}\times300\times(10^{4})^2/10^{21} = 2.91 \times 10^{-11}\,\mathrm{m}\,\mathrm{s}^{-1}; waktu naik sepanjang 10 km10\,\mathrm{km} adalah d/v=3.44×1014 sd/v = 3.44 \times 10^{14}\,\mathrm{s}, 10.9 Ma10.9\,\mathrm{Ma}, waktu difusi d2/κ=1014 sd^2/\kappa = 10^{14}\,\mathrm{s}, 3.2 Ma3.2\,\mathrm{Ma}. Rasio keduanya adalah (d2/κ)/(d/v)=vd/κ=ρgαΔTd3/(κη)=Ra=0.29(d^2/\kappa)/(d/v) = vd/\kappa = \rho g\alpha\Delta Td^3/(\kappa\eta) = \mathrm{Ra} = 0.29. (b) Ia mendingin lebih dulu: gumpalan kecil tidak berkonveksi. (c) Ra∝d3\mathrm{Ra} \propto d^3, sehingga d=10×(657.5/0.291)1/3=131 kmd = 10\times(657.5/0.291)^{1/3} = 131\,\mathrm{km}.

Latihan 6.6 ★★

Untuk lapisan berkedalaman 2891 km2891\,\mathrm{km} pada Ra=107\mathrm{Ra} = 10^{7}, dengan κ=10−6 m2 s−1\kappa = 10^{-6}\,\mathrm{m}^{2}\,\mathrm{s}^{-1}, k=3.1 W m−1 K−1k = 3.1\,\mathrm{W}\,\mathrm{m}^{-1}\,\mathrm{K}^{-1}, ΔT=1450 K\Delta T = 1450\,\mathrm{K} dan Rac=657.5\mathrm{Ra}_c = 657.5, hitung ketebalan lapisan batas, bilangan Nusselt, skala kecepatan dalam cm/yr\mathrm{cm}/\mathrm{yr}, dan aliran panas.

Solusi

Solusi Latihan 6.6.

2Rac=13152\mathrm{Ra}_c = 1315 dan Ra/2Rac=7604\mathrm{Ra}/2\mathrm{Ra}_c = 7604, yang akar pangkat tiganya adalah 19.66. δ=2891/19.66=147 km\delta = 2891/19.66 = 147\,\mathrm{km}, Nu=19.66/2=9.8\mathrm{Nu} = 19.66/2 = 9.8, u=(10−6/2.891×106)×19.662=1.34×10−10 m s−1=0.42 cm/yru = (10^{-6}/2.891 \times 10^{6})\times19.66^2 = 1.34 \times 10^{-10}\,\mathrm{m}\,\mathrm{s}^{-1} = 0.42\,\mathrm{cm}/\mathrm{yr}, dan q=9.8×3.1×1450/2.891×106=15 mW m−2q = 9.8\times3.1\times1450/2.891 \times 10^{6} = 15\,\mathrm{mW}\,\mathrm{m}^{-2}.

Latihan 6.7 ★★

Sebuah tonjolan titik panas rekaan memiliki luas penampang 8×108 m28 \times 10^{8}\,\mathrm{m}^{2} melintang jejak, dan lempeng di atasnya bergerak dengan 9 cm/yr9\,\mathrm{cm}/\mathrm{yr}. Dengan ρm−ρw=2270 kg m−3\rho_m - \rho_w = 2270\,\mathrm{kg}\,\mathrm{m}^{-3}, hitung fluks daya apungnya, dan fluks panasnya dengan cp=1170 J kg−1 K−1c_p = 1170\,\mathrm{J}\,\mathrm{kg}^{-1}\,\mathrm{K}^{-1} dan α=3.1×10−5 K−1\alpha = 3.1 \times 10^{-5}\,\mathrm{K}^{-1}. Berapa pecahan dari aliran panas Bumi sebesar 46.7 TW46.7\,\mathrm{TW} yang dihasilkannya?

Solusi

Solusi Latihan 6.7.

u=0.09/3.156×107=2.85×10−9 m s−1u = 0.09/3.156 \times 10^{7} = 2.85 \times 10^{-9}\,\mathrm{m}\,\mathrm{s}^{-1} dan B=2270×8×108×2.85×10−9=5.18×103 kg s−1=5.2 Mg/sB = 2270\times8 \times 10^{8} \times2.85 \times 10^{-9} = 5.18 \times 10^{3}\,\mathrm{kg}\,\mathrm{s}^{-1} = 5.2\,\mathrm{Mg}/\mathrm{s}, sekitar tiga perlima milik Hawaii. Fluks panas adalah 1170×5178/3.1×10−5=1.95×1011 W1170\times5178/3.1 \times 10^{-5} = 1.95 \times 10^{11}\,\mathrm{W}, 0.2 TW0.2\,\mathrm{TW}, atau 0.42 %0.42\,\% dari kehilangan panas Bumi.

Latihan 6.8 ★★

Sebuah kepala plume berjari-jari 400 km400\,\mathrm{km}, lebih ringan daripada mantel sebesar 30 kg m−330\,\mathrm{kg}\,\mathrm{m}^{-3}, naik melalui mantel dengan viskositas 1021 Pa s10^{21}\,\mathrm{Pa}\,\mathrm{s}. Dengan hukum Stokes, hitung lajunya dalam m/yr\mathrm{m}/\mathrm{yr} dan waktu yang dibutuhkan untuk naik sejauh 2000 km2000\,\mathrm{km}.

Solusi

Solusi Latihan 6.8.

v=2×30×9.8×(4×105)2/(9×1021)=1.05×10−8 m s−1=0.33 m/yrv = 2\times30\times9.8\times(4 \times 10^{5})^2/(9\times10^{21}) = 1.05 \times 10^{-8}\,\mathrm{m}\,\mathrm{s}^{-1} = 0.33\,\mathrm{m}/\mathrm{yr}; 2000 km2000\,\mathrm{km} butuh 2×106/1.05×10−8=1.9×1014 s2 \times 10^{6}/1.05 \times 10^{-8} = 1.9 \times 10^{14}\,\mathrm{s}, sekitar 6 Ma6\,\mathrm{Ma}.

Latihan 6.9 ★★

Sebuah wilayah tertekan sedalam 300 m300\,\mathrm{m} ketika lembaran esnya mencair, 10 ka10\,\mathrm{ka} lalu, dan berelaksasi dengan τ=4.4 ka\tau = 4.4\,\mathrm{ka}. Hitung penurunan yang tersisa hari ini dan laju pengangkatan saat ini. Bandingkan dengan 10.3 mm/yr10.3\,\mathrm{mm}/\mathrm{yr} yang diukur dekat Umeå.

Solusi

Solusi Latihan 6.9.

w=300e−10/4.4=30.9 mw = 300e^{-10/4.4} = 30.9\,\mathrm{m} dan w˙=w/τ=30.9/4400=7.0 mm/yr\dot w = w/\tau = 30.9/4400 = 7.0\,\mathrm{mm}/\mathrm{yr}, dua pertiga laju yang diukur dekat Umeå: setengah-ruang seragam dengan τ\tau ini berada pada rentang yang tepat, tetapi pantulan sebenarnya lebih cepat atau cekungannya lebih dalam.

Latihan 6.10 ★★★

(a) Dari ekspresi FrpF_{rp} dalam pembuktian Proposisi 6.6, tunjukkan bahwa suku pertama, 12(ρm−ρw)gw2\tfrac12(\rho_m - \rho_w)gw^2, adalah suku kedua dalam kurung. (b) Tunjukkan bahwa ∫0∞xerfc⁡x dx=1/4\int_0^\infty x\operatorname{erfc}x\,\mathrm dx = 1/4, dengan integrasi parsial. (c) Dengan nilai-nilai dari latihan 4, hitung kedua suku pada 100 Ma100\,\mathrm{Ma}.

Solusi

Solusi Latihan 6.10.

(a) Dengan w=2ρmαΔTκt/π/(ρm−ρw)w = 2\rho_m\alpha\Delta T\sqrt{\kappa t/\pi}/(\rho_m - \rho_w),

12(ρm−ρw)gw2=2gρm2α2ΔT2κtπ(ρm−ρw)=gρmαΔTκt×2ρmαΔTπ(ρm−ρw).\tfrac12(\rho_m - \rho_w)gw^2 = \frac{2g\rho_m^2\alpha^2\Delta T^2\kappa t} {\pi(\rho_m - \rho_w)} = g\rho_m\alpha\Delta T\kappa t\times \frac{2\rho_m\alpha\Delta T}{\pi(\rho_m - \rho_w)} .

(b) Dengan integrasi parsial, ∫0∞xerfc⁡x dx=[x22erfc⁡x]0∞+∫0∞x22 2πe−x2 dx=0+1π⋅π4=14\int_0^\infty x\operatorname{erfc}x\,\mathrm dx = [\tfrac{x^2}2\operatorname{erfc} x]_0^\infty + \int_0^\infty\tfrac{x^2}2\,\tfrac2{\sqrt\pi}e^{-x^2}\,\mathrm dx = 0 + \tfrac1{\sqrt\pi}\cdot\tfrac{\sqrt\pi}4 = \tfrac14. (c) Pada 100 Ma100\,\mathrm{Ma}, ρmgα∫zθ=3.67×1012 N m−1\rho_mg\alpha\int z\theta = 3.67 \times 10^{12}\,\mathrm{N}\,\mathrm{m}^{-1} dan, dengan w=3705 mw = 3705\,\mathrm{m}, 12(ρm−ρw)gw2=1.53×1011 N m−1\tfrac12(\rho_m - \rho_w)gw^2 = 1.53 \times 10^{11}\,\mathrm{N}\,\mathrm{m}^{-1}: suku air hanya 4 %4\,\% dari gaya total.

Latihan 6.11 ★★★

Tunjukkan bahwa f(k)=(π2+k2)3/k2f(k) = (\pi^2 + k^2)^3/k^2 paling kecil pada k=π/2k = \pi/\sqrt2, dengan nilainya 27π4/427\pi^4/4. Turunkan lebar sepasang sel konveksi pada awal konveksi, dalam satuan dd, dan lebar satu sel.

Solusi

Solusi Latihan 6.11.

ln⁡f=3ln⁡(π2+k2)−2ln⁡k\ln f = 3\ln(\pi^2 + k^2) - 2\ln k, yang turunannya 6k/(π2+k2)−2/k6k/(\pi^2 + k^2) - 2/k bernilai nol ketika 3k2=π2+k23k^2 = \pi^2 + k^2, k2=π2/2k^2 = \pi^2/2. Di sana f=(3π2/2)3/(π2/2)=27π4/4f = (3\pi^2/2)^3/ (\pi^2/2) = 27\pi^4/4. Panjang gelombang 2π/k=22 d=2.83 d2\pi/k = 2\sqrt2\,d = 2.83\,d memuat sepasang sel yang berputar berlawanan arah, masing-masing selebar 2 d=1.41 d\sqrt2\,d = 1.41\,d.

Latihan 6.12 ★★★

Sebuah lempeng sepanjang 6000 km6000\,\mathrm{km} bergerak dengan 6 cm/yr6\,\mathrm{cm}/\mathrm{yr} di atas lapisan mantel setebal 150 km150\,\mathrm{km}, dengan viskositas 1020 Pa s10^{20}\,\mathrm{Pa}\,\mathrm{s}, yang dasarnya diam. (a) Hitung seretan basal dalam MPa\mathrm{MPa} dan gaya seret per satuan panjang lempeng. (b) Bandingkan dengan dorongan punggung sebesar 3.8×1012 N m−13.8 \times 10^{12}\,\mathrm{N}\,\mathrm{m}^{-1} dan tarikan lempengan sebesar 4.6×1013 N m−14.6 \times 10^{13}\,\mathrm{N}\,\mathrm{m}^{-1}. (c) Jika sepersepuluh tarikan lempengan mencapai lempeng, pada laju berapa seretan menyeimbangkan gaya penggeraknya?

Solusi

Solusi Latihan 6.12.

(a) u=1.90×10−9 m s−1u = 1.90 \times 10^{-9}\,\mathrm{m}\,\mathrm{s}^{-1}, ηu/h=1020×1.90×10−9/1.5×105=1.27 MPa\eta u/h = 10^{20}\times1.90 \times 10^{-9}/ 1.5 \times 10^{5} = 1.27\,\mathrm{MPa}, dan sepanjang 6000 km6000\,\mathrm{km}, 7.6×1012 N m−17.6 \times 10^{12}\,\mathrm{N}\,\mathrm{m}^{-1}. (b) Dua kali dorongan punggung dan seperenam tarikan lempengan: seretan dapat menahan dorongan punggung saja, tetapi tidak tarikan lempengan. (c) ηuL/h=3.8×1012+4.6×1012\eta uL/h = 3.8 \times 10^{12} + 4.6 \times 10^{12} memberikan u=8.4×1012×1.5×105/(6×106×1020)=2.1×10−9 m s−1=6.6 cm/yru = 8.4 \times 10^{12}\times1.5 \times 10^{5}/(6 \times 10^{6}\times10^{20}) = 2.1 \times 10^{-9}\,\mathrm{m}\,\mathrm{s}^{-1} = 6.6\,\mathrm{cm}/\mathrm{yr}, laju lempeng samudra yang cepat.

6.8 Soal: Pantulan Danau Ardo

Soal 6.1

Soal akhir pekan — pantulan cekungan danau rekaan, viskositas mantel di bawahnya, dan konveksi yang dimungkinkan oleh viskositas itu

Danau Ardo, danau rekaan, mengisi cekungan selebar 250 km250\,\mathrm{km}, dibendung oleh lembaran es, dengan kedalaman 300 m300\,\mathrm{m} selama sekitar lima puluh ribu tahun. Empat belas ribu tahun lalu bendungan es jebol dan danau surut dalam beberapa minggu. Garis pantai tertinggi, yang dulu datar saat terbentuk, kini miring; garis pantai bertanggal dari danau-danau kecil yang tersisa di cekungan memberikan pengangkatan pusat sejak pengeringan, relatif terhadap tepi yang hampir tidak bergerak:

waktu sejak pengeringan (ka)1246914
pengangkatan pusat (m)11.521.638.050.564.077.2

Ambil ρw=1000 kg m−3\rho_w = 1000\,\mathrm{kg}\,\mathrm{m}^{-3}, ρm=3300 kg m−3\rho_m = 3300\,\mathrm{kg}\,\mathrm{m}^{-3}, g=9.8 m s−2g = 9.8\,\mathrm{m}\,\mathrm{s}^{-2}, 1 yr1\,\mathrm{yr} =3.156×107 s= 3.156 \times 10^{7}\,\mathrm{s}, dan perlakukan mantel sebagai ruang setengah kental seragam.

Bagian I — Beban.

  1. Hitung tekanan danau pada dasarnya.
  2. Hitung penekanan dasar pada kesetimbangan isostatik, dan jelaskan mengapa massa jenis mantel, bukan kerak, muncul di dalamnya.
  3. Anggap danau sebagai cakram, lalu hitung massa airnya.
  4. Jelaskan mengapa garis pantai tertinggi kini miring, dan hitung kemiringan rata-ratanya dari tepi ke pusat, dalam meter per kilometer.
  5. Penekanan diperlakukan sebagai satu harmonik yang setengah tertekannya, dari −λ/4-\lambda/4 sampai λ/4\lambda/4, membentang di seluruh danau. Berikan λ\lambda dan k=2π/λk = 2\pi/\lambda.
  6. Pada kedalaman berapa aliran di bawah beban seperti itu meluruh dengan faktor ee?

Bagian II — Waktu relaksasi.

  1. Dengan anggapan danau telah mencapai kesetimbangan, hitung penekanan yang tersisa pada setiap tanggal dalam tabel.
  2. Tunjukkan bahwa, jika penekanan meluruh secara eksponensial, logaritmanya turun secara linear terhadap waktu.
  3. Hitung τ\tau dari tanggal 1 dan 14 ka14\,\mathrm{ka}.
  4. Prediksikan pengangkatan pada 6 ka6\,\mathrm{ka} dengan τ\tau tersebut, dan bandingkan dengan tabel.
  5. Hitung penekanan yang tersisa hari ini dan laju pengangkatan pusat saat ini, dalam mm/yr\mathrm{mm}/\mathrm{yr}.
  6. Berapa banyak pengangkatan yang masih akan terjadi?
  7. Dalam berapa tahun lagi laju pengangkatan akan turun di bawah 0.5 mm/yr0.5\,\mathrm{mm}/\mathrm{yr}?

Bagian III — Viskositas.

  1. Hitung viskositas ruang setengah dari τ\tau.
  2. Bandingkan dengan viskositas mantel atas sebesar 3.4×10203.4 \times 10^{20} sampai 1×1021 Pa s1 \times 10^{21}\,\mathrm{Pa}\,\mathrm{s} yang ditemukan di bawah Fennoskandia. Bagian mantel mana yang tersampel oleh pantulan Danau Ardo?
  3. Dengan viskositas itu, hitung waktu relaksasi beban dengan panjang gelombang 3000 km3000\,\mathrm{km}. Mengapa lembaran es sungguhan sebesar itu akan meluruh lebih lambat daripada ini?
  4. Apa yang dilakukan litosfer elastik, yang diabaikan di sini, terhadap pantulan beban yang hanya selebar 250 km250\,\mathrm{km}?
  5. Hitung waktu Maxwell mantel dengan μ=7.7×1010 Pa\mu = 7.7 \times 10^{10}\,\mathrm{Pa}, lalu nyatakan apakah model kental itu dapat dibenarkan.
  6. Periksa asumsi pertanyaan 7: berapa bagian dari penekanan kesetimbangan yang tercapai selama danau berdiri?
  7. Berapa besar log⁡10η\log_{10}\eta akan berubah jika λ\lambda diambil sebagai 250 km250\,\mathrm{km}?

Bagian IV — Konveksi. Ambil lapisan mantel atas dengan kedalaman 670 km670\,\mathrm{km}, dengan ρ=3700 kg m−3\rho = 3700\,\mathrm{kg}\,\mathrm{m}^{-3}, g=9.9 m s−2g = 9.9\,\mathrm{m}\,\mathrm{s}^{-2}, α=3.1×10−5 K−1\alpha = 3.1 \times 10^{-5}\,\mathrm{K}^{-1}, ΔT=1450 K\Delta T = 1450\,\mathrm{K}, κ=10−6 m2 s−1\kappa = 10^{-6}\,\mathrm{m}^{2}\,\mathrm{s}^{-1}, k=3.1 W m−1 K−1k = 3.1\,\mathrm{W}\,\mathrm{m}^{-1}\,\mathrm{K}^{-1} dan viskositas Bagian III.

  1. Hitung bilangan Rayleigh lapisan itu.
  2. Berapa kali lipat nilai kritisnya (Rac=657.5\mathrm{Ra}_c = 657.5)?
  3. Hitung tebal lapisan batas dan bilangan Nusselt.
  4. Hitung skala kecepatan dalam cm/yr\mathrm{cm}/\mathrm{yr}, dan bandingkan dengan laju lempeng sebesar 1 sampai 7.5 cm/yr7.5\,\mathrm{cm}/\mathrm{yr}.
  5. Hitung aliran panas, dan berikan dua alasan mengapa nilai itu berada di bawah rata-rata Bumi, 92 mW m−292\,\mathrm{mW}\,\mathrm{m}^{-2}.
  6. Berikan viskositas mantel atas di bawah Danau Ardo sebagai log⁡10\log_{10} dari η\eta dalam Pa s\mathrm{Pa}\,\mathrm{s}.
Solusi

Solusi Soal 6.1.

1. p=1000×9.8×300=2.94 MPap = 1000\times9.8\times300 = 2.94\,\mathrm{MPa}.

2. w=ρwh/ρm=300×1000/3300=90.9 mw = \rho_wh/\rho_m = 300\times1000/3300 = 90.9\,\mathrm{m}. Kerak bergerak turun sebagai satu blok; yang mengalir keluar dari bawah beban, lalu kembali, adalah mantel, sehingga kolom diseimbangkan oleh mantel yang terdesak di kedalaman.

3. π×(1.25×105)2=4.91×1010 m2\pi\times(1.25 \times 10^{5})^2 = 4.91 \times 10^{10}\,\mathrm{m}^{2}, dan air menimbang 4.91×1010×300×1000=1.47×1016 kg4.91 \times 10^{10}\times300\times1000 = 1.47 \times 10^{16}\,\mathrm{kg}.

4. Garis pantai terbentuk datar, ketika dasar paling turun di pusat; sejak itu pusat naik lebih banyak daripada tepi, sehingga garis kini naik menuju pusat, sebesar 77.2 m77.2\,\mathrm{m} sepanjang 125 km125\,\mathrm{km}: 0.62 m/km0.62\,\mathrm{m}/\mathrm{km}.

5. λ/2=250 km\lambda/2 = 250\,\mathrm{km}, sehingga λ=500 km\lambda = 500\,\mathrm{km} dan k=2π/5×105=1.257×10−5 m−1k = 2\pi/5 \times 10^{5} = 1.257 \times 10^{-5}\,\mathrm{m}^{-1}.

6. 1/k=λ/2π=79.6 km1/k = \lambda/2\pi = 79.6\,\mathrm{km}.

7. w=90.9−pengangkatanw = 90.9 - \text{pengangkatan}: 79.4, 69.3, 52.9, 40.4, 26.9 dan 13.7 m13.7\,\mathrm{m}.

8. w=w0e−t/τw = w_0e^{-t/\tau} memberikan ln⁡w=ln⁡w0−t/τ\ln w = \ln w_0 - t/\tau, sebuah garis lurus dengan kemiringan −1/τ-1/\tau.

9. τ=(14−1)/ln⁡(79.4/13.7)=13/1.7571=7.40 ka\tau = (14 - 1)/\ln(79.4/13.7) = 13/1.7571 = 7.40\,\mathrm{ka}.

10. w(6)=79.4 e−5/7.40=40.4 mw(6) = 79.4\,e^{-5/7.40} = 40.4\,\mathrm{m}, pengangkatan sebesar 50.5 m50.5\,\mathrm{m}, sesuai dengan tabel.

11. 13.7 m13.7\,\mathrm{m} tersisa; w˙=13.7/7400=1.85 mm/yr\dot w = 13.7/7400 = 1.85\,\mathrm{mm}/\mathrm{yr}.

12. Seluruh 13.7 m13.7\,\mathrm{m}, didekati makin lama makin lambat.

13. Laju turun mengikuti e−t/τe^{-t/\tau}: t=7.40ln⁡(1.85/0.5)=7.40×1.309=9.7 kat = 7.40\ln(1.85/0.5) = 7.40\times1.309 = 9.7\,\mathrm{ka}, sekitar 9700 tahun lagi.

14. τ=7399×3.156×107=2.335×1011 s\tau = 7399\times3.156 \times 10^{7} = 2.335 \times 10^{11}\,\mathrm{s}, dan η=τρmgλ/(4π)=2.335×1011×3300×9.8×5×105/12.566=3.0×1020 Pa s\eta = \tau\rho_mg\lambda/(4\pi) = 2.335 \times 10^{11}\times3300\times9.8\times5 \times 10^{5}/ 12.566 = 3.0 \times 10^{20}\,\mathrm{Pa}\,\mathrm{s}.

15. Tepat di bawah ujung bawah rentang Fennoscandia. Aliran di bawah beban dengan panjang gelombang 500 km500\,\mathrm{km} meluruh dalam 80 km80\,\mathrm{km}: pantulan mengambil sampel mantel terdangkal dan terlemah, di bawah dan sebagian di dalam yang di Bumi nyata adalah litosfer.

16. τ=4π×3.0×1020/(3300×9.8×3×106)=3.9×1010 s=1.2 ka\tau = 4\pi\times3.0 \times 10^{20}/(3300\times9.8\times3 \times 10^{6}) = 3.9 \times 10^{10}\,\mathrm{s} = 1.2\,\mathrm{ka}. Di bawah beban selebar itu aliran mencapai hingga 477 km477\,\mathrm{km} lebih dalam, menuju mantel bawah yang sepuluh kali lebih viskos, yang memperlambatnya.

17. Litosfer menahan sebagian beban sempit melalui lenturan: depresi pada kesetimbangan lebih dangkal dan lebih lebar daripada 90.9 m90.9\,\mathrm{m} isostatik, dan sebagian pantulan itu bersifat elastik dan seketika. Setengah-ruang kental seragam yang dicocokkan dengan pengangkatan kemudian memberikan viskositas yang mencampur kedua efek.

18. τM=3.0×1020/7.7×1010=3.9×109 s\tau_M = 3.0 \times 10^{20}/7.7 \times 10^{10} = 3.9 \times 10^{9}\,\mathrm{s}, sekitar 124 yr124\,\mathrm{yr}, enam puluh kali lebih pendek daripada τ\tau: selama ribuan tahun mantel bereaksi seperti fluida kental.

19. 1−e−50/7.40=0.99881 - e^{-50/7.40} = 0.9988: 99.9 %99.9\,\% dari depresi kesetimbangan telah tercapai, dan asumsi ini berlaku.

20. η∝λ\eta \propto \lambda: memangkas λ\lambda menjadi separuh memangkas η\eta menjadi separuh, dan log⁡10η\log_{10}\eta turun sebesar log⁡102=0.30\log_{10}2 = 0.30.

21. Ra=3700×9.9×3.1×10−5×1450×(6.7×105)3/(10−6×3.0×1020)=1.65×106\mathrm{Ra} = 3700\times9.9\times3.1 \times 10^{-5}\times1450\times (6.7 \times 10^{5})^3/(10^{-6}\times3.0 \times 10^{20}) = 1.65 \times 10^{6}.

22. 1.65×106/657.5≈25001.65 \times 10^{6}/657.5 \approx 2500 kali kritis.

23. (Ra/2Rac)1/3=(1254)1/3=10.78(\mathrm{Ra}/2\mathrm{Ra}_c)^{1/3} = (1254)^{1/3} = 10.78: δ=670/10.78=62 km\delta = 670/10.78 = 62\,\mathrm{km} dan Nu=5.4\mathrm{Nu} = 5.4.

24. u=(10−6/6.7×105)×10.782=1.73×10−10 m s−1=0.55 cm/yru = (10^{-6}/6.7 \times 10^{5})\times10.78^2 = 1.73 \times 10^{-10}\,\mathrm{m}\,\mathrm{s}^{-1} = 0.55\,\mathrm{cm}/\mathrm{yr}, dua sampai empat belas kali lebih lambat daripada lempeng.

25. q=5.4×3.1×1450/6.7×105=36 mW m−2q = 5.4\times3.1\times1450/6.7 \times 10^{5} = 36\,\mathrm{mW}\,\mathrm{m}^{-2}. Lapisan ini mengabaikan panas yang dihasilkan dalam mantel dan panas dari inti serta mantel bawah; dan penskalaan menghilangkan konstanta berorde satu, sedangkan lempeng, yang mendingin sebagai setengah-ruang dari punggungan, kehilangan jauh lebih banyak panas pada dasar samudra muda daripada lapisan batas rata-rata.

26. η=3.0×1020 Pa s\eta = 3.0 \times 10^{20}\,\mathrm{Pa}\,\mathrm{s}: log⁡10η=20.48\log_{10}\eta = \textbf{$20.48$}.

本章定义的术语

在术语表中查看全部 957 个术语