পৃথিবী ও মহাকাশ বিজ্ঞান · বই 5 · স্নাতক তৃতীয় বর্ষ
Ilmu Kebumian dan Antariksa Universitas — Tahun 3
বইগুলো এখনো বাংলা ভাষায় পাওয়া যায় না — Bahasa Indonesia সংস্করণ দেখানো হচ্ছে।
Ilmu Kebumian dan Antariksa Universitas — Tahun 3 · স্নাতক তৃতীয় বর্ষ
19Dinamika Fluida Geofisika
Aduk secangkir teh dan pusarannya hilang dalam semenit; daun teh berkumpul di tengah dasar cangkir, dan tidak ada yang terjadi lagi. Lautan dan atmosfer juga teraduk oleh Matahari dan angin, tetapi di planet yang berotasi dan dalam lapisan dengan massa jenis berbeda. Rotasi dan stratifikasi membuat keduanya tunduk pada aturan yang tak pernah tampak di cangkir teh: angin bertiup sejajar isobar, bukan melintasinya, kolom fluida menolak untuk miring, dan angin di troposfer atas ditentukan oleh suhu di bawahnya. Bab ini menyusun persamaan untuk aliran besar, lambat, dan tipis ini, serta beberapa bilangan yang menentukan suku mana yang penting.
Yang sudah kamu ketahui
Gaya Coriolis dan gaya sentrifugal dalam kerangka berputar, persamaan Navier–Stokes dan persamaan kontinuitas, viskositas, bilangan Reynolds, dan relasi hidrostatik dibahas di buku fisika. Parameter Coriolis f=2Ωsinφ, gaya gradien tekanan, keseimbangan geostrofik, angin geostrofik dan angin gradien, serta arus jet dibahas di volume Tahun 1 pada bab tentang sirkulasi umum atmosfer; lapisan Ekman pada bab tentang sirkulasi samudra; suhu potensial, stabilitas statis, dan osilasi parsel udara yang terangkat pada bab tentang struktur atmosfer dan samudra.
19.1 Persamaan dalam kerangka berputar
Dalam kerangka yang berputar bersama Bumi dengan kecepatan sudut Ω, sebuah parsel fluida berkecepatan u mematuhi
DtDu+2Ω×u=−ρ1∇p+g+ν∇2u,∂t∂ρ+∇⋅(ρu)=0,
dengan D/Dt=∂/∂t+u⋅∇ adalah turunan yang mengikuti parsel dan ν viskositas kinematik. Gaya sentrifugal telah digabungkan ke dalam gravitasi: jumlahnya g adalah gravitasi efektif, tegak lurus terhadap permukaan geopotensial tetap, dan vertikal didefinisikan olehnya. Hanya suku Coriolis, −2Ω×u, yang tersisa dari rotasi, dan suku ini tidak melakukan usaha: ia memutar parsel tanpa mengubah lajunya. Bumi berputar dengan Ω=7.292115×10−5rad/s, sekali per hari sideris sebesar 23.9345h.
Untuk gerak pada jarak yang jauh lebih kecil daripada jari-jari Bumi a, sebuah kerangka Kartesius lokal sudah memadai: x ke timur, y ke utara, z ke atas, pada lintang φ, yakni bidang singgung pada gambar di bawah. Di sana Ω=Ω(0,cosφ,sinφ), dan ketiga komponen persamaan momentum adalah
dengan F adalah gesekan. Untuk aliran tipis yang hampir mendatar di samudra dan atmosfer, suku-suku yang memuat f~ diabaikan: f~w kecil dibandingkan fv karena w≪v, dan f~u bernilai sepersepuluh ribu dari g (lihat Latihan 19.5). Gaya Coriolis kemudian hanya bekerja pada arah horizontal, melalui komponen vertikal lokal rotasi, Ωsinφ.
Definisi 19.1(bidang-f dan bidang-β)
Di dekat lintang acuan φ0, pendekatan bidang-f menganggap parameter Coriolis konstan, f=f0=2Ωsinφ0. Pendekatan bidang-β mempertahankan perubahan orde pertamanya terhadap jarak ke utara y: f=f0+βy, dengan β=df/dy=2Ωcosφ0/a.
Pada 45∘, dengan jari-jari rata-rata a=6371km, f0=1.031×10−4s−1 dan β=1.62×10−11m−1s−1: dalam 1000km sepanjang lintang, f berubah sebesar 16%. Bidang-f cocok untuk gerak yang jauh lebih kecil daripada itu; bidang-β diperlukan untuk arus berskala cekungan dan gelombang planet.
Kiri, kerangka lokal pada lintang φ: z sepanjang vertikal lokal, y ke utara pada bidang singgung, x ke timur masuk ke dalam halaman. Vektor rotasi terurai menjadi Ωsinφ terhadap vertikal, yang memutar gerak horizontal, dan Ωcosφ terhadap sumbu ke utara, diabaikan untuk aliran tipis. Kanan, parameter Coriolis dan bidang-β singgungnya pada 45∘: dekat untuk pita lintang tertentu, jauh darinya di dekat khatulistiwa dan kutub.
Dua pendekatan lagi menyederhanakan persamaan tanpa kehilangan hal yang penting. Yang pertama menyangkut massa jenis.
Definisi 19.2(Pendekatan Boussinesq)
Yang disebut pendekatan Boussinesq menuliskan massa jenis sebagai ρ=ρ0+ρ′, dengan ρ0 suatu konstanta dan ∣ρ′∣≪ρ0, dan mempertahankan ρ′ hanya jika dikalikan dengan g. Persamaan kontinuitas menjadi ∇⋅u=0, dan persamaan momentum
DtDu+fk×u=−ρ01∇p′+bk+ν∇2u,b=−gρ0ρ′,
dengan p′ adalah selisih tekanan terhadap tekanan hidrostatik massa jenis ρ0, k vektor satuan ke atas, dan b daya apung.
Di samudra ρ′ hanya beberapa kg/m3 dari 1025, sebagian kecil persen, dan pendekatannya sangat baik. Di atmosfer massa jenis turun dengan faktor e dalam sekitar 8.5km; pendekatan ini berlaku untuk lapisan yang jauh lebih tipis dari itu, dan aliran yang lebih dalam mengganti ρ0 dengan profil rata-rata ρ0(z), dengan fisika yang sama.
Yang kedua menyangkut arah vertikal. Fluida yang diam berada dalam kesetimbangan hidrostatik; fluida yang bergerak tetap memenuhinya dengan pendekatan sangat baik bila gerakannya jauh lebih lebar daripada dalamnya.
Definisi 19.3(Pendekatan hidrostatik)
Yang disebut pendekatan hidrostatik menggantikan persamaan momentum vertikal fluida yang bergerak dengan relasi hidrostatik, ∂p/∂z=−ρg.
Untuk gerak dengan skala horizontal L, kedalaman H, dan laju U, dengan bilangan RossbyU/(fL) paling besar berorde satu, percepatan vertikal lebih kecil daripada gradien vertikal perturbasi tekanan dengan faktor berorde (U/fL)(H/L)2.
Bukti. Ketika gaya Coriolis setidaknya sama besar dengan percepatan, maka persamaan momentum horizontal menyeimbangkan gradien tekanan horizontal terhadapnya: p′/(ρ0L)∼fU, sehingga p′/ρ0∼fUL, dan gradien vertikalnya adalah p′/(ρ0H)∼fUL/H. Kontinuitas, ∂w/∂z=−(∂u/∂x+∂v/∂y), membatasi kecepatan vertikal oleh W∼UH/L, dan percepatan vertikal oleh Dw/Dt∼UW/L=U2H/L2. Rasionya adalah (U2H/L2)/(fUL/H)=(U/fL)(H/L)2. ∎
Sebuah sistem cuaca selebar 1000km di troposfer setebal 10km memiliki (H/L)2=10−4 dan bilangan Rossby sekitar 0.1: galat dari pendekatan hidrostatik berorde 10−5. Ia hanya gagal pada gerak yang sempit dan dalam: arus naik badai petir, gumpalan konvektif di samudra.
19.2 Skala aliran berotasi
Suku mana dari persamaan momentum yang penting ditentukan oleh dua rasio.
Definisi 19.5(Bilangan Rossby dan Ekman)
Untuk aliran dengan laju U, skala horizontal L, dan kedalaman H pada planet dengan parameter Coriolis f, yang disebut bilangan Rossby adalah Ro=U/(fL), rasio percepatan, yang berorde U2/L, terhadap gaya Coriolis, berorde fU. Yang disebut bilangan Ekman adalah Ek=ν/(fH2), rasio gaya viskos, berorde νU/H2, terhadap gaya Coriolis; dalam aliran turbulen ν digantikan oleh viskositas eddy.
Air bak mandi yang dikuras pada laju 0.1m/s sepanjang 0.3m memiliki bilangan Rossby sebesar beberapa ribu: rotasi Bumi tidak menentukan arah pusarannya. Depresi selebar 1000km dengan angin 15m/s memiliki Ro≈0.14, dan gyre samudra Ro∼10−4. Viskositas bahkan lebih lemah lagi: viskositas molekuler air adalah 10−6m2/s, dan viskositas eddy yang mewakili turbulensi samudra, 10−5 sampai 10−3m2/s dalam arah vertikal, memberikan Ek∼10−5 pada kedalaman 1000m. Kedua bilangan itu berhubungan dengan bilangan Reynolds menurut Re=UL/ν=(Ro/Ek)(L/H)2.
Ketika kedua bilangan kecil, persamaan momentum horizontal tunak hanya mempertahankan gaya Coriolis dan gradien tekanan: aliran berada dalam keseimbangan geostrofik, fk×ug=−(1/ρ)∇hp, yang sudah dipakai di volume Tahun 1 untuk membaca peta cuaca. Geostrofi adalah batas Ro,Ek→0; koreksinya berorde Ro, dan gesekan hanya bekerja pada lapisan tipis dengan ketebalan relatif Ek, yaitu lapisan Ekman di dasar atmosfer dan bagian atas dan bawah samudra. Pada bidang-f, kecepatan geostrofik tidak memiliki divergensi horizontal: ∂ug/∂x+∂vg/∂y=−(1/ρf)(∂2p/∂x∂y−∂2p/∂y∂x)=0 untuk ρ seragam.
Batas sebaliknya, yaitu rotasi tanpa gradien tekanan, adalah gerak paling sederhana yang dapat dilakukan fluida berotasi.
Definisi 19.6(Osilasi inersia)
Yang disebut osilasi inersia adalah gerak horizontal fluida yang tidak dikenai gaya horizontal apa pun, kecuali gaya Coriolis.
Teorema 19.7(Lingkaran inersia)
Pada bidang-f, parsel yang bergerak dengan laju U tanpa gaya horizontal selain gaya Coriolis mempertahankan lajunya dan membentuk lingkaran dengan jari-jari U/∣f∣, searah jarum jam di belahan bumi utara dan berlawanan arah jarum jam di belahan bumi selatan, dengan periode inersia 2π/∣f∣, yaitu setengah hari sideris dibagi sinφ.
Bukti. Persamaannya adalah du/dt=fv, dv/dt=−fu. Dengan w=u+iv persamaan itu menjadi dw/dt=−ifw, sehingga w=Ue−ift untuk awal bergerak ke timur: laju ∣w∣=U konstan dan kecepatan berputar dengan laju f, searah jarum jam bila f>0. Dengan mengintegralkan, posisi x+iy=(iU/f)(e−ift−1) berputar mengelilingi titik −iU/f pada jarak U/∣f∣. Periodenya adalah 2π/∣f∣=2π/(2Ω∣sinφ∣), dan 2π/Ω adalah hari sideris. ∎
Pada 45∘, periode inersia adalah 16.9h, dan arus dengan laju 0.2m/s berputar pada lingkaran berjari-jari 1.94km; pada 30∘ periodenya sama dengan satu hari sideris penuh, dan di khatulistiwa tidak ada lingkaran inersia. Setelah badai menggerakkan lapisan permukaan samudra lalu mereda, air melakukan persis hal ini. Pelampung hanyut ber-drogue pada 15m di Pasifik Utara membentuk lintasan melingkar selama beberapa hari setelah badai pada Oktober 1987, terbawa oleh arus rata-rata. Lintasan ketiga pada gambar di bawah, setelah lingkaran dan spiral.
Lintasan selama empat hari dari parsel yang bergerak ke timur dengan laju 0.2m/s pada 45∘ LU (titik), tanpa gaya selain gaya Coriolis: lingkaran berjari-jari 1.94km yang ditempuh searah jarum jam setiap 16.9h (abu-abu); dengan gesekan yang mereduksi laju menjadi setengahnya dalam 2.8 hari, spiral ke dalam (biru); hal yang sama terbawa oleh arus ke timur yang tetap dengan laju 0.05m/s, berupa deretan lingkaran yang makin kecil (merah).
Di antara geostrofi dan lingkaran inersia terdapat aliran di sekitar pusat tekanan tinggi atau rendah, tempat isobar melengkung dan percepatan menuju pusat tidak lagi dapat diabaikan.
Definisi 19.8(Angin gradien)
Yang disebut angin gradien adalah aliran tetap tanpa gesekan sepanjang isobar melengkung, tempat gaya gradien tekanan, gaya Coriolis, dan percepatan sentripetal berada dalam keseimbangan.
Proposisi 19.9(Koreksi kelengkungan)
Misalkan isobar berupa lingkaran berjari-jari r di belahan bumi utara, dan Vg laju geostrofik yang akan diberikan oleh gradien tekanannya, dengan Rog=Vg/(fr). Di sekitar pusat rendah, angin gradien adalah V=2Vg/(1+1+4Rog), lebih lambat daripada Vg; di sekitar pusat tinggi, V=2Vg/(1+1−4Rog), lebih cepat, dan ia hanya ada bila Rog≤1/4. Hingga orde pertama, V=Vg(1∓Rog).
Bukti. Di sekitar pusat rendah angin berputar berlawanan arah jarum jam; gaya gradien tekanan fVg menuju pusat, gaya Coriolis fV menuju luar, dan selisih keduanya memberikan percepatan sentripetal: fVg−fV=V2/r. Akar positif dari V2+frV−frVg=0 adalah V=(fr/2)(1+4Rog−1); dengan mengalikan (1+4Rog+1)/(1+4Rog+1) diperoleh rumus tersebut. Di sekitar pusat tinggi, angin berputar searah jarum jam dan gaya Coriolis menuju ke dalam: fV−fVg=V2/r, sehingga V2−frV+frVg=0, yang hanya memiliki akar real bila f2r2≥4frVg; akar yang menuju Vg ketika Rog→0 adalah V=(fr/2)(1−1−4Rog). Dengan menguraikan akar-akarnya, V=Vg(1∓Rog+…). ∎
Bilangan Rossby isobar mengukur koreksinya: sekitar sepersepuluh untuk depresi berukuran sedang. Syarat Rog≤1/4 menjelaskan mengapa pusat antisiklon memiliki gradien lemah dan angin ringan, sedangkan yang di pusat depresi dapat memiliki gradien curam dan angin kencang.
19.3 Rotasi tanpa stratifikasi
Fluida berotasi dengan massa jenis seragam, dalam gerak lambat, berperilaku seolah-olah tersusun dari kolom kaku yang sejajar dengan sumbu rotasi.
Teorema 19.10(Taylor–Proudman)
Pada aliran tetap fluida tak-kental dengan massa jenis seragam, yang berotasi dengan Ω, dengan bilangan Rossby kecil, kecepatan tidak berubah sepanjang arah sumbu rotasi: (Ω⋅∇)u=0.
Bukti. Untuk Ro→0, tetap dan tak-kental, persamaan momentum tereduksi menjadi 2Ω×u=−∇(p/ρ0+Φ), dengan Φ adalah geopotensial. Kurl suatu gradien bernilai nol, sehingga ∇×(Ω×u)=0. Untuk Ω konstan, kurl ini adalah Ω(∇⋅u)−(Ω⋅∇)u, dan ∇⋅u=0 untuk massa jenis seragam. ∎
Dengan Ω searah z, ∂u/∂z=0: ketiga komponen kecepatan sama di setiap ketinggian. Jika dasar kaku melarang w di sana, w=0 di mana-mana, dan aliran horizontal sama di setiap tingkat. Oleh karena itu, rintangan di dasar tangki yang berputar terasa sampai ke permukaan: fluida di atasnya tidak dapat naik melampauinya, karena hal itu akan meregangkan atau memampatkan kolom, dan mengelilinginya pada setiap ketinggian. Kolom fluida yang tetap berada di atas rintangan disebut kolom Taylor.
Kolom Taylor (skematis). (a) Dalam tangki berisi air bermassa jenis seragam yang berputar cepat, aliran lambat menjumpai rintangan rendah di dasar; fluida di atas rintangan (abu-abu) tertahan di tempatnya sampai ke permukaan, dan aliran pada setiap ketinggian mengelilingi kolom seolah-olah rintangan itu mencapai puncak. (b) Aliran yang terlihat dari atas sama di setiap tingkat.
Atmosfer dan samudra terstratifikasi, dan stratifikasi melonggarkan kendala itu: fluida terstratifikasi dapat bervariasi terhadap ketinggian, sebagaimana ditentukan oleh teorema angin termal.
19.4 Stratifikasi dan angin termal
Volume Tahun 1 mengangkat parsel udara dalam atmosfer stabil dan mendapati bahwa parsel itu berosilasi di sekitar tingkatnya. Hal yang sama berlaku pada fluida terstratifikasi mana pun, dan frekuensi osilasi tersebut adalah ukuran stratifikasi.
Definisi 19.11(Frekuensi daya apung)
Yang disebut frekuensi daya apungN suatu fluida dalam keadaan diam didefinisikan oleh
N2=θgdzdθ(udara),N2=−ρ0gdzdρθ(air laut),
dengan θ adalah suhu potensial udara dan ρθ massa jenis potensial air laut, yakni massa jenis yang akan dimilikinya jika dibawa secara adiabatik ke tekanan acuan.
Teorema 19.12(Osilasi parsel yang tergeser)
Sebuah parsel yang digeser secara vertikal sejauh jarak kecil ζ dari tingkatnya dalam fluida diam, secara adiabatik dan tanpa pencampuran, memenuhi ζ¨=−N2ζ. Jika N2>0 parsel berosilasi pada frekuensi sudut N; jika N2<0 simpangannya membesar dan fluida tidak stabil.
Bukti. Dipindahkan secara adiabatik dan tanpa pencampuran, parsel mempertahankan massa jenis potensialnya ρθ(z0). Pada tingkat z0+ζ ia berada pada tekanan lingkungannya, yang massa jenis potensialnya adalah ρθ(z0)+ζdρθ/dz; dua parsel pada tekanan yang sama berbeda massa jenisnya sebesar selisih massa jenis potensialnya, hingga orde pertama. Gaya daya apung per satuan massa adalah −g(ρp−ρe)/ρ0=gζ(dρθ/dz)/ρ0=−N2ζ. Untuk udara pada tekanan tertentu massa jenis berubah sebagai 1/θ, sehingga −dρθ/ρ0=dθ/θ dan kedua bentuk N2 sama. ∎
Di troposfer atmosfer standar, yang mendingin sebesar 6.49K/km dibandingkan dengan 9.77K/km pada adiabat kering (g/cp, dengan cp=γR/(γ−1), γ=1.4 dan R=286.9Jkg−1K−1), N rata-rata bernilai 0.0113s−1, yang berarti periode 9.3min; di stratosfer bawah yang isotermal, nilainya 0.0210s−1. Di samudra, kolom air di Atlantik Utara pada 40.5∘ LU, 30.5∘ BB, dalam rata-rata tahunan, kolom air memiliki stratifikasi terkuat di piknoklin musiman, N=0.0106s−1 pada sekitar 30m, enam siklus per jam, dan nilai yang jauh lebih lemah di bawahnya: beberapa siklus per jam paling banyak, sebagaimana teramati di seluruh samudra bagian atas.
Frekuensi daya apung. Kiri, atmosfer standar: hampir seragam di troposfer, dua kali lebih besar di stratosfer bawah yang isotermal, dan lebih besar lagi di atas 25km, tempat suhu meningkat. Kanan, kolom rata-rata tahunan Atlantik Utara (suhu dan salinitas dari World Ocean Atlas 2023, massa jenis potensial dari EOS-80): maksimum tajam di piknoklin musiman, stratifikasi lemah di bawahnya.
Stratifikasi memungkinkan aliran berubah menurut ketinggian, dan cara perubahannya bergantung pada apakah permukaan massa jenis konstan mengikuti permukaan tekanan konstan.
Definisi 19.13(Fluida barotropik dan baroklinik)
Suatu fluida disebut barotropik jika massa jenisnya hanya fungsi tekanan, sehingga permukaan massa jenis konstan berimpit dengan permukaan tekanan konstan; fluida disebut baroklinik jika kedua permukaan itu saling berpotongan.
Di atmosfer, gradien suhu horizontal sepanjang permukaan tekanan membuat udara baroklinik. Keseimbangan hidrostatik kemudian membuat permukaan tekanan semakin curam terhadap ketinggian, dan keseimbangan geostrofik mengubah kemiringan yang membesar menjadi angin yang menguat.
Geometri angin termal (skematik, penampang meridional dengan sisi hangat di selatan, p1>p2>p3>p4). Udara hangat memuai, sehingga lapisan di antara dua permukaan tekanan lebih tebal di sisi hangat daripada di sisi dingin: setiap permukaan miring turun ke arah sisi dingin lebih curam daripada permukaan di bawahnya. Angin geostrofik sebanding dengan kemiringan, dan menguat seiring ketinggian (lingkaran, makin besar untuk angin yang lebih kuat), ke timur, keluar dari bidang gambar.
Teorema 19.14(Angin termal)
Dalam atmosfer yang berada dalam keseimbangan hidrostatik dan geostrofik, angin geostrofik berubah terhadap tekanan sebagai
∂lnp∂ug=fR(∂y∂T)p,∂lnp∂vg=−fR(∂x∂T)p.
Di antara dua aras tekanan pb>pt, angin termal, yakni selisih angin geostrofik, adalah
ug(pt)−ug(pb)=fRlnptpbk×∇pTˉ,
dengan Tˉ adalah suhu rata-rata lapisan: angin bertiup sepanjang isoterm lapisan, dengan udara dingin di sebelah kirinya dalam belahan bumi utara.
Bukti. Pada permukaan tekanan konstan, dp=(∂p/∂y)dy+(∂p/∂z)dz=0, dan dengan ∂p/∂z=−ρg, ketinggian zp permukaan tersebut memenuhi g∂zp/∂y=(1/ρ)∂p/∂y. Dalam bentuk geopotensialΦ=gzp permukaan itu, keseimbangan geostrofik dinyatakan sebagai fug=−∂Φ/∂y dan fvg=∂Φ/∂x, dengan turunan diambil sepanjang permukaan itu. Hubungan hidrostatik dan hukum gas ideal memberikan ∂Φ/∂p=−1/ρ=−RT/p, yaitu ∂Φ/∂lnp=−RT. Dengan menurunkan fug=−∂Φ/∂y terhadap lnp menghasilkan f∂ug/∂lnp=R∂T/∂y, dan dengan cara serupa untuk vg. Dengan mengintegralkan dari pb ke pt, lnp berkurang sebesar ln(pb/pt), sehingga ug(pt)−ug(pb)=−(R/f)ln(pb/pt)∂Tˉ/∂y, dan vg(pt)−vg(pb)=(R/f)ln(pb/pt)∂Tˉ/∂x: keduanya adalah komponen bentuk vektor, karena k×∇Tˉ=(−∂Tˉ/∂y,∂Tˉ/∂x). ∎
Akibat 19.15(Angin termal dalam ketinggian)
Dalam koordinat ketinggian, untuk atmosfer, ∂ug/∂z=−(g/fT)∂T/∂y+(ug/T)∂T/∂z; untuk samudra dengan pendekatan Boussinesq, ∂ug/∂z=(g/fρ0)∂ρ/∂y dan ∂vg/∂z=−(g/fρ0)∂ρ/∂x.
Bukti. Di udara, fug=−(1/ρ)∂p/∂y=−RT∂lnp/∂y dan ∂lnp/∂z=−g/(RT). Lalu
f∂z∂(ug/T)=−R∂y∂∂z∂lnp=g∂y∂(1/T)=−T2g∂y∂T,
Dengan menjabarkan ruas kiri diperoleh rumusnya. Di samudra, fug=−(1/ρ0)∂p/∂y dan ∂p/∂z=−ρg, sehingga f∂ug/∂z=(g/ρ0)∂ρ/∂y; demikian pula untuk vg. ∎
Suku kedua bentuk atmosfer, yang berasal dari penurunan suhu terhadap ketinggian, lebih kecil daripada suku pertama dan lenyap dalam lapisan isotermal; diintegralkan melintasi lapisan dengan suhu rata-ratanya, suku pertama saja hampir memberikan bentuk tekanan, karena Δz≈(RTˉ/g)ln(pb/pt) (Soal, Bagian IV). Di samudra, teorema yang sama mengubah penampang massa jenis melintasi arus menjadi geseran vertikalnya: dasar metode dinamis para oseanograf.
Arus jet pada volume Tahun 1 adalah angin termal. Pada Januari rata-rata zonal 1991–2020, angin barat di 30∘ LU bertambah dari 4.3m/s pada 850hPa hingga 40.4m/s pada 250hPa; dengan f di 30∘ dan R, geseran ini adalah ciri lapisan yang suhu rata-ratanya turun ke utara sebesar 7.5K per 1000km. Pada 45∘ LU geseran dari 6.0 sampai 23.2m/s memberi 5.0K per 1000km. Gambar di bawah menunjukkan profil dan, untuk setiap lintang, penurunan suhu ke arah kutub yang dibaca dari geseran: terbesar di subtropik belahan bumi musim dingin, tempat arus jet subtropik berada, 43.9m/s pada 200hPa dan 30∘ LU pada Januari, dan 41.0m/s pada 30∘ LS pada Juli. Di atas 200hPa pada 30∘ LU angin melemah lagi: teorema angin termal yang dibaca terbalik menyatakan bahwa di sana, di stratosfer bawah, suhu naik ke arah kutub.
Kiri, angin zonal rata-rata pada Januari (NCEP–NCAR Reanalysis 1, rata-rata 1991–2020) pada tiga lintang: angin bertambah seiring ketinggian di seluruh troposfer. Kanan, penurunan suhu rata-rata lapisan ke arah kutub antara 850hPa dan 250hPa dibaca dari geseran angin dengan teorema angin termal (di bawah 25∘ ke arah khatulistiwa, tempat geostrofi gagal, dikeluarkan): lebih besar di belahan bumi musim dingin.
19.5 Jari-jari deformasi Rossby
Rotasi dan stratifikasi masing-masing menetapkan suatu skala; bersama-sama keduanya menetapkan suatu panjang.
Definisi 19.16(Jari-jari deformasi Rossby)
Yang disebut jari-jari deformasi Rossby dari lapisan dengan kedalaman H dan frekuensi daya apungN adalah LR=NH/∣f∣.
Proposisi 19.17(Energi potensial dan kinetik aliran seimbang)
Dalam aliran hidrostatik dan geostrofik berskala horizontal L di dalam lapisan berkedalaman H dan frekuensi daya apungN, rasio energi potensial yang tersimpan dalam permukaan massa jenis yang tergeser terhadap energi kinetik berorde (L/LR)2.
Bukti. Geostrofi memberikan gangguan tekanan p′/ρ0∼fUL; hubungan hidrostatik memberikan daya apung b∼p′/(ρ0H)∼fUL/H. Daya apung b dalam fluida berfrekuensi daya apung N berarti permukaan massa jenis tergeser secara vertikal sejauh η=b/N2, yang menyimpan energi potensial per satuan massa berorde N2η2=b2/N2∼f2U2L2/(N2H2). Dibagi dengan energi kinetik yang berorde U2, hasilnya (fL/NH)2=(L/LR)2. ∎
Gerak yang jauh lebih lebar daripada LR menyimpan sebagian besar energinya dalam permukaan massa jenis yang miring, gerak yang jauh lebih sempit dalam kecepatannya; pada LR keduanya sama, dan ketidakstabilan aliran baroklinik, yang memanfaatkan energi potensial permukaan yang miring, tumbuh paling cepat di dekat skala itu. Untuk troposfer, dengan N rata-rata atmosfer standar, H=11km dan f di 45∘, LR=1210km: skala depresi dan antisiklon pada peta cuaca. Untuk kilometer teratas kolom Atlantik Utara di 40.5∘ LU, LR=38km, tiga puluh kali lebih kecil: samudra memiliki cuacanya sendiri, eddy selebar puluhan kilometer, dan model yang mampu menyelesaikannya memerlukan kisi yang jauh lebih halus. Pada skala dekat LR, f berubah sebesar βLR/f, 19% untuk atmosfer dan kurang dari satu persen untuk samudra: eddy samudra hidup pada bidang-f, sistem cuaca merasakan bidang-β.
Di laboratorium— Planet di atas meja putar
Bak yang berputar: pewarna berpilin mengelilingi silinder, direkam oleh kamera yang ikut berputar bersama meja (ilustrasi buatan AI).
Di laboratorium dinamika fluida geofisika, sebuah bak air berdiri di atas meja putar; kamera ikut berputar bersamanya dan melihat air seperti pengamat di planet. Sebelum setiap percobaan, air dibiarkan mengalami spin-up hingga berputar sebagai benda tegar, yang memakan waktu beberapa menit, banyak putaran meja. Lalu para peneliti menyuntikkan zat pewarna, atau menggerakkan penghalang perlahan di sepanjang dasar: pewarna naik membentuk kolom vertikal di atas penghalang dan aliran terbelah mengelilinginya di setiap level. Dengan lapisan air tawar di atas air asin, atau dinding luar yang dipanaskan dan pusat yang didinginkan, bak menjadi baroklinik, dan arus jet, eddy dan front berkelok yang muncul, ukurannya diprediksi oleh jari-jari deformasi Rossby bak itu.
19.6 Latihan
Latihan 19.1★
Ambil Ω=7.292115×10−5rad/s. (a) Hitung parameter Coriolis dan periode inersia pada 30∘, 45∘ dan 60∘ LU. (b) Hitung jari-jari lingkaran inersia arus sebesar 0.20m/s pada lintang yang sama. (c) Bandingkan periode inersia pada 30∘ dengan hari sideris, 23.9345h.
Solusi
Solusi Latihan 19.1.
(a) f=2Ωsinφ: 7.292×10−5, 1.031×10−4, dan 1.263×10−4s−1; 2π/f: 23.93, 16.92, dan 13.82h. (b) U/f: 2.74, 1.94, dan 1.58km. (c) Pada 30∘, sinφ=1/2 dan 2π/f=2π/Ω: periode inersia adalah hari sideral, 23.93h.
Latihan 19.2★
Pada 45∘ LU, hitung bilangan Rossby untuk pusaran debu (10m/s, 10m), angin laut (5m/s, 20km), depresi (15m/s, 1000km) dan gyre samudra (0.05m/s, 2000km). Manakah di antara mereka yang mendekati keseimbangan geostrofik?
Solusi
Solusi Latihan 19.2.
Dengan f=1.031×10−4s−1: pusaran debu 9.7×103, angin laut 2.4, depresi 0.15, gyre 2.4×10−4. Depresi dan gyre mendekati keseimbangan geostrofik; pusaran debu berputar ke kedua arah, tidak menyadari rotasi Bumi, dan angin laut merasakannya hanya selama sehari.
Latihan 19.3★
Pada 45∘ LU suhu rata-rata troposfer, 250K, turun ke arah utara sebesar 5K per 1000km. Dengan menggunakan suku pertama angin termal terhadap ketinggian, hitung geser vertikal angin barat dan pertambahan angin sepanjang 10km.
Solusi
Solusi Latihan 19.3.
∂ug/∂z=(g/fT)∣∂T/∂y∣=9.81×5×10−6/(1.031×10−4×250)=1.90×10−3s−1; sepanjang 10km angin barat meningkat sebesar 19m/s.
Latihan 19.4★
Stabilitas E=−(1/ρ)dρθ/dz pada kilometer teratas samudra berada antara 50×10−8 dan 1000×10−8m−1; di palung dalam nilainya 10−8m−1. Hitung N=gE, periode osilasi daya apung dan N dalam siklus per jam untuk ketiga nilai tersebut.
Solusi
Solusi Latihan 19.4.
E=50×10−8: N=2.21×10−3s−1, periode 47min, 1.27 siklus per jam. E=1000×10−8: N=9.90×10−3s−1, 10.6min, 5.7 siklus per jam. E=10−8: N=3.13×10−4s−1, 5.6h, 0.18 siklus per jam.
Latihan 19.5★★
Suatu sistem cuaca memiliki U=10m/s, L=1000km, H=10km, pada 45∘ LU. (a) Perkirakan kecepatan vertikal dan percepatan vertikalnya, dan bandingkan yang terakhir dengan g dan dengan gradien vertikal gangguan tekanan, fUL/H. (b) Hitung suku Coriolis yang diabaikan f~u dan bandingkan dengan g.
Solusi
Solusi Latihan 19.5.
(a) W∼UH/L=0.1m/s dan Dw/Dt∼UW/L=10−6m/s2, 1.0×10−7 kali g. Suku tekanan, fUL/H=0.103m/s2, lebih besar 1×105 kali: rasionya, 9.7×10−6, adalah Ro(H/L)2 dengan Ro=0.097. (b) f~u=2Ωcos45∘×10=1.03×10−3m/s2, 1.05×10−4 kali g.
Latihan 19.6★★
Ambil f=10−4s−1 dan isobar melingkar berjari-jari 300km. (a) Gradien tekanan di sekitar tekanan rendah memberikan Vg=20m/s: hitung angin gradien. (b) Di sekitar tekanan tinggi, Vg=5m/s: hitung angin gradien. (c) Berapa Vg terbesar yang mungkin di sekitar tekanan tinggi, dan berapa anginnya saat itu?
Solusi
Solusi Latihan 19.6.
fr=30m/s. (a) Rog=0.667, V=40/(1+3.667)=13.7m/s, sepertiga lebih lambat daripada geostrofi. (b) Rog=0.167, V=10/(1+0.333)=6.34m/s. (c) Vg≤fr/4=7.5m/s, dan kemudian V=fr/2=15m/s, dua kali kecepatan geostrofik.
Latihan 19.7★★
Sebuah bak berkedalaman 20cm berputar sekali setiap 6s; airnya memiliki ν=10−6m2/s. (a) Hitung f=2Ω, lalu bilangan Rossby untuk aliran 1cm/s yang melewati silinder berdiameter 5cm, serta bilangan Ekman bak itu. Apakah kolom Taylor akan terbentuk? (b) Teori memberikan waktu spin-up H/νΩ: hitung dalam menit dan dalam putaran meja.
Solusi
Solusi Latihan 19.7.
(a) Ω=2π/6=1.047rad/s, f=2.094s−1. Ro=0.01/(2.094×0.05)=0.095 dan Ek=10−6/(2.094×0.04)=1.19×10−5: keduanya kecil, sehingga aliran mematuhi teorema Taylor–Proudman dan kolom Taylor terbentuk. (b) 0.2/10−6×1.047=195s, 3.3min, 33 putaran meja.
Latihan 19.8★★
Pada rata-rata zonal Januari 1991–2020, angin barat sebesar 4.27m/s pada 850hPa dan 40.39m/s pada 250hPa di 30∘ LU, 5.99m/s dan 23.17m/s di 45∘ LU. Dengan R=286.9Jkg−1K−1, hitung gradien meridional suhu rata-rata lapisan pada kedua lintang, dan jelaskan mengapa arus jet lebih kuat pada 30∘ LU.
Solusi
Solusi Latihan 19.8.
∂Tˉ/∂y=−fΔu/(Rln(850/250)), dengan ln(850/250)=1.224. Pada 30∘ LU, −7.29×10−5×36.12/(286.9×1.224)=−7.50K per 1000km; pada 45∘ LU, −1.031×10−4×17.18/351.1=−5.05K per 1000km. Geser sebanding dengan gradien dibagi f: gradiennya 1.5 kali lebih besar pada 30∘ dan f 1.4 kali lebih kecil, sehingga geser 2.1 kali lebih besar.
Latihan 19.9★★
(a) Hitung jari-jari deformasi Rossby troposfer (N=0.0113s−1, H=11km, 45∘ LU) dan samudra bagian atas (N=0.00357s−1, H=1000m, 40.5∘ LU). (b) Dengan a=6371km, hitung β pada kedua lintang dan perubahan relatif parameter Coriolis, βLR/f, melintasi setiap jari-jari.
Solusi
Solusi Latihan 19.9.
(a) LR=0.0113×11000/1.031×10−4=1206km; 0.00357×1000/9.47×10−5=37.7km, 32 kali lebih kecil. (b) β=2Ωcosφ/a: 1.62×10−11 dan 1.74×10−11m−1s−1. βLR/f=0.19 untuk atmosfer, 0.007 untuk samudra.
Latihan 19.10★★★
Sebuah parsel samudra permukaan yang bergerak dengan kecepatan U diperlambat oleh gesekan linear: du/dt=fv−ru, dv/dt=−fu−rv. (a) Selesaikan dengan w=u+iv dan jelaskan lintasannya. (b) Parsel juga terbawa arus tetap u0 ke timur: berikan posisinya. (c) Pada 45∘ LU dengan U=0.2m/s, r−1=4 hari dan u0=0.05m/s, berikan kecepatannya relatif terhadap arus setelah empat hari dan jarak yang telah dihanyutkannya.
Solusi
Solusi Latihan 19.10.
(a) dw/dt=−(r+if)w, sehingga w=Ue−rte−ift: kecepatan berputar dengan laju f sementara magnitudonya meluruh sebagai e−rt. Posisi x+iy=U(1−e−(r+if)t)/(r+if) berspiral ke titik U/(r+if). (b) x+iy=U(1−e−(r+if)t)/(r+if)+u0t: spiral itu tertarik memanjang menjadi gelung. (c) 0.2e−1=0.074m/s; pusat telah bergeser u0t=0.05×4×86400=17.3km ke timur.
Latihan 19.11★★★
Melintasi front samudra pada 40∘ LU, massa jenis air meningkat ke utara sebesar 1kg/m3 sepanjang 100km, pada lapisan atas setebal 1000m; di bawahnya, air seragam dan diam. Dengan ρ0=1025kg/m3, hitung geser vertikal arus dan arus di permukaan. Ke arah mana arus itu mengalir?
Solusi
Solusi Latihan 19.11.
f=9.37×10−5s−1; ∂u/∂z=(g/fρ0)∂ρ/∂y=9.81×10−5/(9.37×10−5×1025)=1.02×10−3s−1. Dengan u=0 pada 1000m arus permukaan adalah 1.02m/s, ke timur: di belahan bumi utara arus mengalir dengan air yang lebih ringan di sebelah kanannya.
Latihan 19.12★★★
Kolom rata-rata tahunan pada 40.5∘ LU, 30.5∘ BB memiliki massa jenis potensial σθ=ρθ−1000kg/m3 sebesar 26.029 di permukaan, 26.697 pada 100m, 27.079 pada 500m dan 27.628 pada 1000m. (a) Hitung N pada ketiga lapisan, dalam s−1 dan dalam siklus per jam. (b) Lapisan mana yang paling stabil, dan mengapa? (c) Bandingkan nilai-nilainya dengan rentang lapisan atas samudra pada Latihan 19.4.
Solusi
Solusi Latihan 19.12.
(a) N2=(g/ρ)Δσθ/Δz. 0–100m: 9.81×0.668/(1026.4×100), N=7.99×10−3s−1, 4.6 siklus per jam; 100–500m: 3.02×10−3s−1, 1.7 siklus per jam; 500–1000m: 3.24×10−3s−1, 1.9 siklus per jam. (b) 100m teratas: lapisan ini mengandung air permukaan yang hangat di atas piknoklin musiman. Antara 100 dan 500m air lebih seragam daripada di bawahnya. (c) E=N2/g adalah 650×10−8, 93×10−8, dan 107×10−8m−1, semuanya berada dalam kisaran samudra bagian atas.
19.7 Soal: Angin di atas Stepa Orrin
Soal 19.1
Soal akhir pekan — rotasi, stratifikasi dan angin di dataran rekaan pada 48∘ LU
Stepa Orrin adalah dataran datar pada 48∘ LU dengan sebuah danau, Danau Orrin, lebar 12km, dan sebuah laboratorium kecil dengan tangki berputar. Gunakan Ω=7.292115×10−5rad/s, g=9.81m/s2, R=286.9Jkg−1K−1 dan cp=1004Jkg−1K−1. Semua nilai stepa, danaunya, dan anginnya adalah rekaan.
Bagian I — Apakah rotasi penting?
Hitung parameter Coriolis dan periode inersia di Orrin.
Hitung bilangan Rossby dari lima aliran: air yang mengalir keluar dari bak mandi (0.1m/s sepanjang 0.3m); tangki laboratorium yang berputar sekali setiap 10s (dengan f=2Ωtangki; 5mm/s sepanjang 0.2m); angin sore dari danau (5m/s, 20km); sebuah depresi (15m/s, 1000km); gyre dari sebuah cekungan samudra pada lintang yang sama (0.05m/s, 2000km).
Manakah dari aliran ini yang merasakan rotasi yang menggerakkannya, dan mana yang mendekati geostrofik?
Hitung bilangan Ekman untuk tangki (kedalaman 0.2m, ν=10−6m2/s), depresi (viskositas eddy 10m2/s, kedalaman 10km), dan gyre (viskositas eddy 10−3m2/s, kedalaman 1000m).
Lapisan Ekman memiliki ketebalan δ=2ν/f. Tunjukkan bahwa δ/H=2Ek dan hitung δ untuk tangki dan untuk udara di stepa.
Apa yang dikatakan bilangan-bilangan ini tentang tempat gesekan berpengaruh?
Bagian II — Angin di permukaan tanah.
Di permukaan tanah, tekanannya 1012.9hPa dan suhunya 288.1K. Hitung massa jenis udara.
Tekanan permukaan laut turun ke utara sebesar 4hPa per 300km. Hitung angin geostrofik dan berikan arahnya.
Gradien yang sama mengelilingi sebuah depresi yang isobarnya memiliki jari-jari 600km. Hitung bilangan Rossby dari isobar dan angin gradien.
Gradien yang sama mengelilingi sebuah antisiklon dengan jari-jari yang sama. Hitung angin gradien dan gradien tekanan terbesar, dalam hPa per 300km, yang masih dapat dipertahankan antisiklon ini.
Bagian III — Stratifikasi.
Suhu troposfer turun sebesar 6.49K/km. Hitung laju penurunan adiabatik kering g/cp, dan N serta periode osilasi daya apung pada 5km, dengan suhu 255.7K.
Di atas 11km udara bersifat isotermal pada 216.6K. Tunjukkan bahwa N2=g2/(cpT) di sana, dan hitung N serta rasionya terhadap nilai pada pertanyaan 11.
Hitung jari-jari deformasi Rossby troposfer Orrin, dengan N dari pertanyaan 11 dan H=10km, lalu bandingkan dengan ukuran depresi pada Bagian I.
Pada musim panas, 20m lapisan atas Danau Orrin lebih ringan daripada air dalam sebesar 1.5kg/m3 (ρ0=1000kg/m3). Hitung gravitasi tereduksi g′=gΔρ/ρ0, kecepatan g′h dari gelombang internal lambat pada lapisan atas, dan jari-jari deformasi Rossbyg′h/f.
Apakah rotasi Bumi akan membentuk gerakan internal danau?
Bagian IV — Angin termal.
Dua stasiun radiosonde berjarak 600km pada meridian yang sama. Suhu rata-rata lapisan dari permukaan tanah, 1013hPa hingga 265hPa, sekitar 10km, adalah 253.0K di stasiun selatan dan 249.4K di stasiun utara. Hitung gradien suhu meridional dalam K per 1000km.
Ketebalan lapisan di antara dua tekanan adalah Δz=(RTˉ/g)ln(pb/pt). Hitung di kedua stasiun, dan kemiringan permukaan 265hPa relatif terhadap permukaan tanah.
Turunkan kenaikan angin geostrofik dari permukaan tanah hingga 265hPa, dari kemiringan itu.
Hitung kenaikan yang sama dengan teorema angin termal dalam koordinat tekanan.
Hitung kenaikan dari suku pertama angin termal dalam ketinggian, dengan suhu rata-rata troposfer standar antara 0 dan 10km, 255.7K, dan bandingkan.
Hitung angin pada 541hPa, sekitar 5km, dengan menganggap bahwa gradien suhu sama di setiap aras.
Bilangan Richardson N2/(∂u/∂z)2 harus turun di bawah 0.25 agar geser memicu turbulensi. Hitung untuk troposfer Orrin.
Pada Januari, angin barat rata-rata zonal pada 250hPa adalah 23.2m/s di 45∘ LU dan 40.4m/s di 30∘ LU. Bandingkan dengan angin di Orrin, dan jelaskan mengapa angin lokal dapat melampaui rata-rata zonal.
Berikan angin geostrofik pada 10km di atas Stepa Orrin.
2. Bak mandi 3.1×103; tangki (f=4π/10=1.257s−1) 0.020; angin laut 2.3; depresi 0.14; gyre 2.3×10−4.
3. Bak mandi sama sekali tidak merasakan rotasi Bumi; angin laut merasakannya hanya setelah beberapa jam. Tangki, depresi, dan gyre memiliki bilangan Rossby kecil dan mendekati geostrofik.
6. Gesekan bekerja pada lapisan yang tebalnya beberapa perseribu kedalaman di tangki dan samudra, beberapa perseratus di atmosfer; di atasnya aliran hampir tanpa gesekan.
7.ρ=p/(RT)=101290/(286.9×288.1)=1.225kg/m3.
8.ug=(1/ρf)∣∂p/∂y∣=(400/300000)/(1.225×1.084×10−4)=10.0m/s, dari barat, dengan tekanan rendah berada di utara, di sebelah kiri.
9.Rog=10.04/(1.084×10−4×600000)=0.154; V=2×10.04/(1+1.618)=8.84m/s, 12% lebih lambat.
10.V=2×10.04/(1+0.382)=12.4m/s. Nilai Vg terbesar adalah fr/4=16.3m/s, gradien ρfVg=2.16×10−3Pa/m, 6.5hPa per 300km.
11.g/cp=9.77K/km. N2=(9.81/255.7)×(9.77−6.49)×10−3=1.26×10−4s−2, N=0.0112s−1, dengan periode 9.3min.
12. Dengan Γ=0, N2=(g/T)(g/cp)=g2/(cpT)=4.43×10−4s−2, N=0.0210s−1, 1.9 kali nilai troposfer.
13.LR=0.0112×104/1.084×10−4=1035km, kira-kira sama dengan ukuran depresi.
15. Ya: danau selebar 12km, lebih dari dua kali jari-jari deformasi Rossby-nya, sehingga gelombang dan arus internal di dalamnya dibelokkan oleh gaya Coriolis.
16.(249.4−253.0)/600=−6.0K per 1000km.
17.ln(1013/265)=1.3409. Selatan: 286.9×253.0×1.3409/9.81=9922m; utara 9781m. Permukaan 265hPa lebih rendah 141m di utara, dengan kemiringan 2.35×10−4 relatif terhadap permukaan bawah.
19.(R/f)ln(pb/pt)∣∂Tˉ/∂y∣=(286.9/1.084×10−4)×1.3409×6.0×10−6=21.3m/s, hasilnya sama: kedua perhitungan itu sebenarnya satu.
20.Baroklinik: suhu, dan karena itu massa jenis, berubah sepanjang permukaan tekanan. Barotropik, angin tidak akan berubah terhadap ketinggian: 10.0m/s.
21.(g/fTˉ)∣∂T/∂y∣Δz=9.81×6.0×10−6×104/(1.084×10−4×255.7)=21.2m/s, dalam 0.3%.
23.∂u/∂z=21.3/104=2.13×10−3s−1; Ri=1.26×10−4/4.54×10−6=28, jauh di atas 0.25: geser tidak membuat troposfer bebas menjadi turbulen.
24. Angin Orrin, 31m/s, lebih kuat daripada rata-rata zonal pada 45∘ LU dan lebih lemah daripada pada 30∘ LU. Arus jet front kutub bergeser ke utara dan selatan; pada sembarang saat gradien dan anginnya terkonsentrasi di sebuah pita, dan rata-rata sepanjang lingkaran lintang menyebarkannya.
25.10.0+21.3: angin geostrofik pada 10km di atas stepa adalah 31m/s, dari barat (31.3m/s menurut perhitungan).