Płyty tectoniczne: fizyka ruchliwej Ziemi

Planeta z miękkim wnętrzem

Zewnątrz Ziemia wydaje się solida, ale jej wnętrze jest warstwowane jak cebula, a jedna z tych warstw jest wystarczająco miękka, aby płynąć. Krusta i górna, twarda część mantelu tworzą lithosferę — warstwę około 70-100 km grubą, podzieloną na około siedemnaście głównych fragmentów, które nazywamy płytami tectonicznymi. Pod nią znajduje się asthenosfera, zasiany w górnej części mantelu obszar połomlony, który jest mekanicznie słaby do tego stopnia, że gęste płyty powyżej mogą na nim przesuwać się. To właśnie na tym cienkim, wolnym warstwie cała zewnętrzna skóra planety spokojnie przemieszcza się.

Dalej w dół mantel zachowuje się jak prawie niezwykle lepkie ciecz — około 10²¹ Pa·s, czyli około triliarda trilionów razy lepsza niż woda — jednak podana ilość lat pozwala mu konwective. Goryczka blisko granicy jądra jest zanurzona i podnosi się; chłodna goryczka blisko podstawy lithosfery opada. Ta kółko przepływu jest motorem, mimo że nie bezpośrednim, za każdym sejsmom, szczytami gór i dolinami oceanów na mapie.

Kiedy zaczyna się konwekcja? Liczba Rayleigh

Czy warstwa cieczy convectuje czy tylko spokojnie przekształca ciepło, jest zarządzana przez pojedyncze bezwymiarowe wyrażenie — liczba Rayleigh:

Ra = (ρ · g · α · ΔT · d³) / (η · κ) ρ gęstość α współczynnik rozszerzalności termicznej g przyspieszenie ziemskie ΔT różnica temperatury w warstwie d grubość warstwy η lepkość κ dyfuzyja termiczna Konwekcja zaczyna się, gdy Ra przekroczy około 1000. Dla mantelu Ziemi, Ra ≈ 10⁷ — intensywnie convectujący. Podstawowa fizyka za tą liczbą to równania Stokes połączone z transportem ciepła: na skalę mantelu moment bezwładności jest niewielki, więc gradienty ciśnienia i siły lepkie dokładnie balansują przyspieszenie ziemskie, podczas gdy temperatura porusza się wraz z przepływem i wolno się rozsuwana jest przez dyfuzję. Rozwiązanie tego systemu w pełni trójwymiarowym jest zadaniem dla superkomputerów — co powoduje, że każda symulacja w czasie rzeczywistym, w tym ta na tej stronie, pracuje z drastycznie uproszczonym podstawowym obrazem.

Ra = (ρ · g · α · ΔT · d³) / (η · κ)

ρ  density          α  thermal expansion coefficient
g  gravity           ΔT temperature difference across the layer
d  layer thickness    η  viscosity
κ  thermal diffusivity

Convection begins once Ra exceeds roughly 1,000.
For Earth's mantle, Ra ≈ 10⁷ — vigorously convecting.

Trzy sposoby spotkania płyt

Gdziekolwiek dwa płyty mają wspólną granicę, ich wzajemne ruch należy do jednej z trzech kategorii, każda ze swoim charakterystycznym śladem geologicznym.

Granice odchylające się oddzielają się nawzajem. Magmę podnoszą na miejsce i solidaryzują w nową krustę oceaniczną — rozszerzanie się dna morskiego. Wyspa Środkowego Atlantyku, 16 000 km niziny podwodnej, jest rezultatem tego procesu. Granice konvergentne kolidują: cięższa płyta oceaniczna jest zmuszona do wędrówki pod lekkojszą kontynentalną, co prowadzi do subdukcji i wywołuje głębokie trzęsienia ziemi oraz eksplozywne vulkanizm arcji — okrąg pacyficzny "Ogniska", oznaczony jest tych obszarów. Gdzie dwa kontynentalne płyty spotykają się, żadna z nich nie chce podchodzić, więc kruesta krumpli się w górę — kolizja ta stworzyła Himalaje, gdy Indie uderzyły w Azję około 40-50 milionów lat temu. Granice transformacyjne przesuwają się wzajemnie poziomo, zamykając i przysypiając w pasach, które wywołują trzęsienia ziemi — Przecinek San Andreas w Kalifornii jest przykładem typowym, gdzie płyta Pacyficzna przesuwa się na zachód względem Ameryki Północnej o około 5 cm rocznie.

Co rzeczywiście tarczy

Jest podziwiać, że ruch płyt jest przypisywany przepływom konwekcyjnym, które je niesą jak półki transportowe, ale dominująca siła jest znacznie bardziej prostym procesem. Tarcie płyt (slab pull) — stara, zimna i gęsta lithosfera oceaniczna opadająca pod własnym ciężarem w obszarach subdukcji — odpowiada za około 90% całkowitej siły tarcia, ciągnąca resztę płyty za sobą jak obracanie się przykrycia na stole. Tarcie górne (ridge push) – spadkowy ruch nowego, podwyższonego odcinka crusty od centralnej doliny oceanicznej – dodaje mniejszą część. Płyn tarczący pod płytą może pomagać lub przeciwdziałać w zależności od geometrii komórki konwekcyjnej pod nią, a lokalne zasuctionowanie podstawowe ze względu na podnoszenie się płomyków lithosfery może rozrywać płyty i wywoływać nowe rzeźby. Nic z tego nie dzieje się szybko – typowe prędkości płyt wynoszą 1-10 cm rocznie, co jest porównywalne do szybkości wzrostu paznokcia, choć GPS może teraz mierzyć to bezpośrednio: Atlantyk rozszerza się o około 2,5 cm rocznie.

Jak symulacja w przeglądarce kłamie przekonująco

Nie rozwiązujemy równań Stoksa w tablii przeglądarki z rozdzielczością planetarną, więc Płyty tectoniczne używa trzech połączeń, które reprezentują zachowanie jakościowe w sposób ekonomiczny. Po pierwsze, powierzchnia globu jest podzielona na diagram Voronoi — każda płyta to obszar siatki najbliżej poruszającego się punktu źródłowego. Po drugie, podlega 2D siatkę o niskiej rozdzielczości przechowywany jest uproszczony pola temperatury i prędkości: na każdym klatce przesuwa temperaturę wzdłuż ruchu, rozpraszając ją za pomocą małego stencilu Laplasa, a następnie wylicza prędkość zanurzeniową z wynikowego gradientu temperatury. To pole prędkości jest to, co przesuwa seed punkty Voronoi, co sprawia, że płyty się poruszają.

function updateBoundaryMesh(plateA, plateB, velocity) { const relSpeed = velocity.dot(boundary.normal); // + konwergentne, - divergentne if (relSpeed > 0) { // konwergencja: podnoszenie terenu wzdłuż granicy, a następnie uśmierzanie boundary.vertices.forEach(v => v.y = Math.min(v.y + relSpeed * dt * MOUNTAIN_SCALE, MAX_HEIGHT)); smoothTerrain(plateA, 3); smoothTerrain(plateB, 3); } else { // divergenca: wstawianie nowych wierzchołków jako młodego, ciepłego pokrycia insertVertices(plateA, interpolateBoundary(boundary, -relSpeed * dt), HOT_CRUST_COLOR); } } Trzecia warstwa zajmuje się samymi granicami: na każdym wspólnym krawędziowym obszarze decyduje o tym, czy region otwiera się, kolide lub zdeformuje prędkość ruchu między sąsiednimi płytami. Krawędzie divergentne wstawiają nowe trójkąty farbowane jako świeże, ciepłe pokrycie; krawędzie konwergentne deforma powierzchnię w górę proporcjonalnie do prędkości zamknięcia, a następnie stosują uśmierzanie Gaussa, aby góry podnoszą się stopniowo zamiast nagle. Krawędzie transformacyjne akumulują naprężenie przemieszczen jako wartość "energii trzęsień ziemi", która rzuca się w partycje, gdy próg jest przekroczony. Przybliżone co sześćdziesiąt klatek cała siatka jest ponownie triangulowana z ujemnym Delaunayem, aby utrzymać czystość geometryi granic — ale tylko wierzchołki bliskie aktywnym krawędziom potrzebują tego, ponieważ większość wnętrza każdej płyty pozostaje niewzmocniona. Nic z tego nie jest geologicznie dokładne, ale pokrywa właściwą łańcuchowość przyczyn: konwekcja prowadzi do ruchu, ruch na granicach buduje lub zniszcza pokrycie, a kształt pokrycia wpływa na to, co widzisz na ekranie.

function updateBoundaryMesh(plateA, plateB, velocity) {
  const relSpeed = velocity.dot(boundary.normal);  // + convergent, - divergent
  if (relSpeed > 0) {
    // convergent: raise terrain along the boundary, then smooth
    boundary.vertices.forEach(v => v.y = Math.min(v.y + relSpeed * dt * MOUNTAIN_SCALE, MAX_HEIGHT));
    smoothTerrain(plateA, 3); smoothTerrain(plateB, 3);
  } else {
    // divergent: insert new vertices as young, hot crust
    insertVertices(plateA, interpolateBoundary(boundary, -relSpeed * dt), HOT_CRUST_COLOR);
  }
}

Często zadawane pytania

Wypróbuj na żywo

Wszystko powyżej działa bezpośrednio w Twojej przeglądarce — otwórz Shifting Continental Plates – A Dynamic Model i zmieniaj parametry podczas działania. Nic nie jest instalowane ani przesyłane na serwer, cały model działa w jednej karcie.