পৃথিবী ও মহাকাশ বিজ্ঞান · বই 5 · স্নাতক তৃতীয় বর্ষ

Universitaire aard- en ruimtewetenschappen — jaar 3

বইগুলো এখনো বাংলা ভাষায় পাওয়া যায় না — Nederlands সংস্করণ দেখানো হচ্ছে।

Universitaire aard- en ruimtewetenschappen — jaar 3 · স্নাতক তৃতীয় বর্ষ

4Geomagnetisme en paleomagnetisme

Een basaltstroom die zestig miljoen jaar geleden uitbarstte, koelde af met zijn piepkleine magnetietkorrels vastgelegd evenwijdig aan het veld van die tijd. Vandaag geboord, ter plaatse georiënteerd, in kleine cilinders gezaagd en stap voor stap verhit in een kamer die is afgeschermd van het huidige veld, vertelt hij nog steeds hoe ver hij toen van de pool lag en welke kant hij op wees. Duizenden zulke metingen, van elk continent en uit elke tijd, hebben laten zien hoe de continenten bewogen voordat de oceaanbodems ontstonden die hun beweging registreren, en in welke volgorde het veld is omgekeerd door de lagen van een sedimentaire sectie heen. Dit hoofdstuk beschrijft het veld zoals het er vandaag uitziet, hoe gesteenten het vastleggen, hoe een richting wordt bepaald en gemiddeld, en wat een paleomagnetische pool zegt over de beweging van een plaat.

Wat je al weet

De Aarde als magneet en de magnetische declinatie staan in de schoolboeken van de natuurkundeboeken; magnetische domeinen, hysterese, remanentie en de Curietemperatuur in hun Boek 3; het veld van een magnetische dipool in hun Boek 1. Magnetische omkeringen en de normale en omgekeerde polariteiten kwamen aan bod in leerjaar 11. De grootcirkel en de bolvormige driehoek werden gedefinieerd in Boek 1, in het hoofdstuk over bolastronomie; het polariteitschron, de geomagnetische polariteitsschaal, de marien magnetische anomalie en eindige rotaties om een Eulerpool in het hoofdstuk over plaatkinematica op een bol; het hotspotspoor in het hoofdstuk over hotspots. Sferische harmonischen en de Schmidt-genormaliseerde Legendre-functies PnmP_n^m werden voor de zwaartekracht gebruikt in het hoofdstuk over zwaartekracht en de vorm van de Aarde (zie Hoofdstuk 3).

4.1 Het hoofdveld

Definitie 4.1 (Geomagnetisch hoofdveld; geocentrische axiale dipool)

Het geomagnetisch hoofdveld is het deel van het magnetisch veld van de Aarde dat binnen de Aarde wordt opgewekt, in de vloeibare ijzeren buitenkern; het verandert over jaren tot miljoenen jaren, anders dan de velden van de ionosfeer en de magnetosfeer, die over seconden tot dagen veranderen. De geocentrische axiale dipool is het veld van een magnetische dipool die in het middelpunt van de Aarde is geplaatst en met zijn rotatieas is uitgelijnd, wijzend naar het zuiden, zodat de veldlijnen het zuidelijk halfrond verlaten en het noordelijk halfrond binnengaan.

Op elk punt van het oppervlak is het veld F\mathbf F een vector. Het wordt beschreven in een lokaal assenstelsel (figuur hieronder): XX naar geografisch noord, YY naar het oosten, ZZ naar beneden. Zijn horizontale component is H=X2+Y2H = \sqrt{X^2 + Y^2}, zijn totale intensiteit F=∣F∣F = |\mathbf F|, zijn magnetische declinatie DD de hoek van geografisch noord naar H\mathbf H, rechtsom geteld, en zijn inclinatie II de hoek van de horizontaal naar F\mathbf F.

Definitie 4.2 (Magnetische inclinatie)

De magnetische inclinatie II van het veld op een punt is de hoek tussen de veldvector en het horizontale vlak, positief wanneer het veld naar beneden wijst: tan⁡I=Z/H\tan I = Z/H, met −90∘≤I≤90∘-90{}^{\circ} \le I \le 90{}^{\circ}.

De elementen van het veld op een punt: de noord-X-, oost-Y- en neer-Z-componenten, de horizontale component H, het totale veld F, de magnetische declinatie D (in het horizontale vlak, grijs) en de inclinatie I (in het verticale vlak dat F).
De elementen van het veld op een punt: de noord-XX-, oost-YY- en neer-ZZ-componenten, de horizontale component H\mathbf H, het totale veld F\mathbf F, de magnetische declinatie DD (in het horizontale vlak, grijs) en de inclinatie II (in het verticale vlak dat F\mathbf F).

Stelling 4.3 (De dipoolformule)

Aan het oppervlak van een bol met straal aa geldt voor het veld van een geocentrische axiale dipool met moment mm op breedtegraad λ\lambda:

H=B0cos⁡λ,Z=2B0sin⁡λ,F=B01+3sin⁡2λ,tan⁡I=2tan⁡λ,D=0,H = B_0\cos\lambda,\qquad Z = 2B_0\sin\lambda,\qquad F = B_0\sqrt{1 + 3\sin^2\lambda},\qquad \tan I = 2\tan\lambda,\qquad D = 0,

waarbij B0=μ0m/(4πa3)B_0 = \mu_0 m/(4\pi a^3) het veld aan de evenaar is. Dit is de dipoolformule.

Bewijs. Buiten zijn bronnen wordt het veld afgeleid van een magnetische scalaire potentiaal, B=−∇V\mathbf B = -\nabla V. Voor een dipool met moment mm die naar −z-z wijst, het geografische zuiden, geldt V=−μ0mcos⁡θ/(4πr2)V = -\mu_0 m\cos\theta/(4\pi r^2), met θ\theta de colatitude. In bolcoördinaten is Br=−∂V/∂r=−2μ0mcos⁡θ/(4πr3)B_r = -\partial V/\partial r = -2\mu_0 m\cos\theta/(4\pi r^3) en Bθ=−(1/r) ∂V/∂θ=−μ0msin⁡θ/(4πr3)B_\theta = -(1/r)\,\partial V/\partial\theta = -\mu_0 m\sin\theta/(4\pi r^3), en Bφ=0B_\varphi = 0. Bij r=ar = a is de neerwaartse component Z=−Br=2B0cos⁡θ=2B0sin⁡λZ = -B_r = 2B_0\cos\theta = 2B_0\sin\lambda, en omdat eθ\mathbf e_\theta naar het zuiden wijst, is de naar het noorden gerichte component X=−Bθ=B0sin⁡θ=B0cos⁡λX = -B_\theta = B_0\sin\theta = B_0\cos\lambda, met Y=0Y = 0. Dan geldt H=XH = X, D=0D = 0, tan⁡I=Z/H=2tan⁡λ\tan I = Z/H = 2\tan\lambda en F2=B02(cos⁡2λ+4sin⁡2λ)=B02(1+3sin⁡2λ)F^2 = B_0^2(\cos^2\lambda + 4\sin^2\lambda) = B_0^2(1 + 3\sin^2\lambda). ∎

De inclinatie is nul aan de evenaar, 90∘90{}^{\circ} aan de noordpool, en negatief op het zuidelijk halfrond; de intensiteit is aan de polen twee keer zo groot als aan de evenaar. Het werkelijke veld is geen perfecte dipool, en het wordt, net als de zwaartekracht, beschreven met bolfuncties.

Definitie 4.4 (International Geomagnetic Reference Field)

Het International Geomagnetic Reference Field (IGRF) is het internationaal overeengekomen model van het hoofdveld: de potentiaal, op en boven het oppervlak,

V(r,θ,φ)=a∑n=1N(ar)n+1∑m=0n(gnmcos⁡mφ+hnmsin⁡mφ)Pnm(cos⁡θ),V(r,\theta,\varphi) = a\sum_{n=1}^{N}\Bigl(\frac ar\Bigr)^{n+1}\sum_{m=0}^{n} \bigl(g_n^m\cos m\varphi + h_n^m\sin m\varphi\bigr)P_n^m(\cos\theta),

met a=6371.2 kma = 6371.2\,\mathrm{km}, B=−∇V\mathbf B = -\nabla V, en Gauss-coëfficiënten gnmg_n^m, hnmh_n^m (in nT\mathrm{nT}) die sinds 1900 elke vijf jaar worden gegeven, samen met hun veranderingssnelheid over de komende vijf jaar.

Het model wordt onderhouden door de International Association of Geomagnetism and Aeronomy; de veertiende generatie werd in november 2024 voltooid en is afgekapt bij graad N=13N = 13. De eerste term, g10g_1^0, is de axiale dipool: vergeleken met het bewijs hierboven geldt B0=−g10B_0 = -g_1^0, en g10=−29 350 nTg_1^0 = -29\,350\,\mathrm{nT} in 2025. Met de twee andere graad-1-termen is de dipool gekanteld; zijn as snijdt het oppervlak in de noordelijke geomagnetische pool, op 80.79∘80.79{}^{\circ} N, 72.76∘72.76{}^{\circ} W in 2025, en zijn moment is m=4πa3(g10)2+(g11)2+(h11)2/μ0=7.69×1022 A m2m = 4\pi a^3\sqrt{(g_1^0)^2 + (g_1^1)^2 + (h_1^1)^2}/\mu_0 = 7.69 \times 10^{22}\,\mathrm{A}\,\mathrm{m}^{2}.

Het hoofdveld in 2025 volgens het veertiende IGRF, aan het oppervlak van de bol met straal 6371.2\, km. (a) Inclinatie in graden; de dikke lijn, I = 0, is de dip-evenaar. (b) Totale intensiteit in nT; het kruisje markeert het zwakste veld, ongeveer 22\,100\, nT, boven de Zuid-Atlantische Oceaan, tegenover ongeveer 66\,600\, nT op het sterkste punt, ten zuiden van Australië.
Het hoofdveld in 2025 volgens het veertiende IGRF, aan het oppervlak van de bol met straal 6371.2 km6371.2\,\mathrm{km}. (a) Inclinatie in graden; de dikke lijn, I=0I = 0, is de dip-evenaar. (b) Totale intensiteit in nT\mathrm{nT}; het kruisje markeert het zwakste veld, ongeveer 22 100 nT22\,100\,\mathrm{nT}, boven de Zuid-Atlantische Oceaan, tegenover ongeveer 66 600 nT66\,600\,\mathrm{nT} op het sterkste punt, ten zuiden van Australië.

De kaarten hierboven laten zien hoe sterk het veld van een dipool afwijkt. De dip-evenaar wijkt tot ongeveer 12∘12{}^{\circ} af van de geografische evenaar, en de intensiteit, die een dipool langs een breedtecirkel uniform zou maken, varieert met een factor drie, met een brede lage zone boven de Zuid-Atlantische Oceaan.

Definitie 4.5 (Seculaire variatie)

De seculaire variatie is de langzame verandering van het hoofdveld in de tijd, over jaren tot eeuwen, gemeten aan de hand van de veranderingssnelheid van de Gauss-coëfficiënten of van de elementen op één plaats.

In de afgelopen 125 jaar is het dipoolmoment gedaald van 8.32×1022 A m28.32 \times 10^{22}\,\mathrm{A}\,\mathrm{m}^{2} naar 7.69×1022 A m27.69 \times 10^{22}\,\mathrm{A}\,\mathrm{m}^{2}, met 7.6 %7.6\,\%, en is de noordelijke geomagnetische pool verschoven van 78.61∘78.61{}^{\circ} N, 68.79∘68.79{}^{\circ} W naar zijn huidige positie (figuur hieronder); g10g_1^0 neemt nog steeds toe met 12.6 nT12.6\,\mathrm{nT} per jaar, wat wijst op een verzwakking van de dipool.

Het dipoolgedeelte van het hoofdveld, 1900–2025, elke vijf jaar, volgens het IGRF: het moment (links) en de breedte van de noordelijke geomagnetische pool (rechts).
Het dipoolgedeelte van het hoofdveld, 1900–2025, elke vijf jaar, volgens het IGRF: het moment (links) en de breedte van de noordelijke geomagnetische pool (rechts).

Waar komt het veld vandaan? Een permanente magneet valt af: de Curietemperatuur van magnetiet is ongeveer 580 ∘C580\,{}^{\circ}\mathrm{C} en die van hematiet 675 ∘C675\,{}^{\circ}\mathrm{C}, en de diepe Aarde is veel heter. Het veld wordt voortdurend opgewekt door de beweging van de geleidende vloeistof van de buitenkern.

Definitie 4.6 (Geodynamo)

De geodynamo is het zichzelf in stand houdende proces waarbij convectie van de elektrisch geleidende vloeibare ijzerlegering van de buitenkern, afgebogen door de rotatie van de Aarde, elektrische stromen opwekt waarvan het magnetische veld zichzelf handhaaft tegen ohmse vervaging.

Propositie 4.7 (Inductievergelijking)

In een geleider met magnetische diffusiviteit η=1/(μ0σ)\eta = 1/(\mu_0\sigma) die beweegt met de snelheid u\mathbf u, gehoorzaamt het magnetisch veld aan

∂B∂t=∇×(u×B)+η∇2B.\frac{\partial\mathbf B}{\partial t} = \nabla\times(\mathbf u\times\mathbf B) + \eta\nabla^2\mathbf B .

Bewijs. Op dit niveau zonder bewijs aangenomen. ∎

De eerste term rekt de veldlijnen uit die door de stroming worden meegevoerd en draait ze, de tweede laat het veld vervallen, in een lichaam van grootte LL, in een tijd van orde L2/ηL^2/\eta. Een dynamo werkt wanneer de stroming sneller veld aanmaakt dan het wegdiffundeert. Computersimulaties van een convecterende, roterende schil van vloeistof leveren een overheersend axiaal dipool op en af en toe een omkering; maar de vergelijkingen van de geodynamo, gekoppeld aan die van de stroming, vallen buiten dit boek. Naast de volledige omkeringen bevat het archief kortere gebeurtenissen.

Definitie 4.8 (Geomagnetische excursie)

Een geomagnetische excursie is een grote afwijking van de veldrichting ten opzichte van die van het axiale dipool, waarbij de virtuele geomagnetische pool (Definitie 4.16) zich ver van de geografische pool wegbeweegt, en die terugkeert naar de oorspronkelijke polariteit in plaats van een omkering te voltooien.

4.2 Hoe gesteenten het veld vastleggen

Ferromagnetische mineralen in de brede zin, vooral magnetiet in stollingsgesteenten en hematiet in rode sedimenten, dragen een remanentie. Boven hun Curietemperatuur heeft een korrel er geen; daaronder is een kleine korrel een enkel magnetisch domein waarvan het moment thermisch geactiveerd tussen zijn voorkeursrichtingen omklapt, na een relaxatietijd

τ=1Cexp⁡(KVkBT),\tau = \frac1C\exp\Bigl(\frac{K V}{k_B T}\Bigr),

waarbij KK de anisotropie-energie per volume-eenheid is, VV het volume van de korrel en CC een frequentiefactor. Dit is de relaxatiewet van Néel, hier aangenomen. Omdat τ\tau exponentieel groeit naarmate TT daalt, gaat een korrel binnen een smal temperatuurinterval over van vrij omklappen naar het vasthouden van zijn moment gedurende miljarden jaren.

Definitie 4.9 (Blokkeringstemperatuur)

De blokkeringstemperatuur van een magnetische korrel is de temperatuur waaronder zijn relaxatietijd groter is dan de duur van het experiment, zodat zijn moment vastligt: bij afkoeling door deze temperatuur legt de korrel het veld vast; bij verhitting door deze temperatuur in een nulveld verliest de korrel die registratie.

Definitie 4.10 (Natuurlijke, thermische, detritische en chemische remanente magnetisatie)

De natuurlijke remanente magnetisatie (NRM) van een gesteente is de remanentie die het draagt op het moment dat het wordt verzameld, de som van componenten die op verschillende tijdstippen zijn verkregen. Een thermoremanente magnetisatie (TRM) ontstaat door afkoeling door de blokkeringstemperaturen in een veld, zoals in een lavastroom of intrusie. Een detritische remanente magnetisatie (DRM) ontstaat wanneer magnetische korrels die in water bezinken, of in het natte sediment net onder het oppervlak, zich in het veld richten voordat ze worden vastgelegd. Een chemische remanente magnetisatie (CRM) ontstaat wanneer een magnetisch mineraal groeit, door een blokkeringsvolume heen, in een veld, zoals hematiet doet wanneer een rode laag zich vormt.

Een TRM legt het veld vast op het moment van afkoeling; een DRM binnen de tijd die het sediment nodig heeft om te compacteren; een CRM, telkens wanneer het mineraal groeide, wat lang na de afzetting kan zijn. Een latere verhitting, of het huidige veld dat gedurende lange tijd inwerkt op korrels met een korte relaxatietijd, voegt een secundaire component toe bovenop: een viskeuze remanentie langs het recente veld.

4.3 Van monster naar richting

Georiënteerde boorkernen worden op meerdere locaties in een formatie geboord, enkele kernen per locatie, en gesneden tot monsters van enkele kubieke centimeters. Van elk monster wordt de NRM gemeten en vervolgens stapsgewijs verwijderd.

Definitie 4.11 (Stapsgewijze demagnetisatie)

Stapsgewijze demagnetisatie is het verwijderen van de remanentie van een monster in stappen van toenemende intensiteit, waarbij na elke stap wordt de resterende vector gemeten: door verhitting tot steeds hogere temperaturen en afkoeling in een nulveld (thermische demagnetisatie), waarbij eerst de componenten met een lagere blokkeringstemperatuur verdwijnen, of door wisselvelden met toenemende piekwaarde die naar nul terugvallen (wisselveld-demagnetisatie), waarbij eerst die met een lagere coërciviteit verdwijnen.

Definitie 4.12 (Zijderveld-diagram; karakteristieke remanente magnetisatie)

Een Zijderveld-diagram zet de eindpunten van de resterende remanentievector na elke stap van demagnetisatie tegelijk uit op twee vlakken, het horizontale (noord tegen oost) en een verticaal vlak (noord tegen omlaag), met een gemeenschappelijke noordas. Een component die afzonderlijk over een reeks stappen wordt verwijderd, geeft een rechte lijn; een lijnstuk dat in de oorsprong eindigt is de karakteristieke remanente magnetisatie (ChRM), de meest stabiele component, die als de primaire wordt geïnterpreteerd wanneer de toetsen dat toelaten.

Het monster in de onderstaande figuur draagt twee componenten: een viskeuze overprint langs het huidige veld, verwijderd bij 300 ∘C300\,{}^{\circ}\mathrm{C}, en een ChRM die vrijkomt van 300 tot 580 ∘C580\,{}^{\circ}\mathrm{C}, de Curietemperatuur van magnetiet. Boven 300 ∘C300\,{}^{\circ}\mathrm{C} lopen beide projecties recht naar de oorsprong; de richting van dat lijnstuk is de ChRM.

(a) Thermische demagnetisatie van een synthetisch monster, van 20\, C tot 580\, C: gevulde stippen, de horizontale projectie; open stippen, de verticale. (b) De zes ChRM-richtingen van vindplaats A van het weekendprobleem (stippen), hun gemiddelde (ster) en hun kegel van 95 procent betrouwbaarheid (cirkel).
(a) Thermische demagnetisatie van een synthetisch monster, van 20 ∘C20\,{}^{\circ}\mathrm{C} tot 580 ∘C580\,{}^{\circ}\mathrm{C}: gevulde stippen, de horizontale projectie; open stippen, de verticale. (b) De zes ChRM-richtingen van vindplaats A van het weekendprobleem (stippen), hun gemiddelde (ster) en hun kegel van 95 procent betrouwbaarheid (cirkel).

De richtingen van de monsters van een vindplaats spreiden zich rond een gemiddelde, en de statistiek van die spreiding is die van een verdeling op de bol.

Propositie 4.13 (Fisher-verdeling)

De Fisher-verdeling met de precisieparameter κ>0\kappa > 0 rond een gemiddelde richting geeft aan een eenheidsvector de kansdichtheid, per eenheid van ruimtehoek,

f(θ)=κ4πsinh⁡κ eκcos⁡θ,f(\theta) = \frac{\kappa}{4\pi\sinh\kappa}\,e^{\kappa\cos\theta},

waarbij θ\theta de hoek met het gemiddelde is. Gegeven NN eenheidsvectoren ui\mathbf u_i die eruit getrokken zijn, met resultante R=∑ui\mathbf R = \sum\mathbf u_i van lengte RR, is de beste schatting van de gemiddelde richting die van R\mathbf R; de precisie wordt geschat door k=(N−1)/(N−R)k = (N - 1)/(N - R); en het ware gemiddelde ligt met kans 1−p1 - p binnen de hoek α1−p\alpha_{1-p} van het geschatte gemiddelde, waarbij

cos⁡α1−p=1−N−RR[(1p)1/(N−1)−1].\cos\alpha_{1-p} = 1 - \frac{N - R}{R}\Bigl[\Bigl(\frac1p\Bigr)^{1/(N-1)} - 1\Bigr].

Bewijs. Op dit niveau zonder bewijs aangenomen. ∎

Definitie 4.14 (Alpha-95)

De alpha-95 (α95\alpha_{95}) van een gemiddelde richting is de halve hoek van zijn kegel van 95 procent betrouwbaarheid, de hoek α1−p\alpha_{1-p} van Propositie 4.13 voor p=0.05p = 0.05.

De dichtheid integreert tot één: met dΩ=2πsin⁡θ dθ\mathrm d\Omega = 2\pi\sin\theta\,\mathrm d\theta volgt ∫0πf 2πsin⁡θ dθ=κ2sinh⁡κ⋅eκ−e−κκ=1\int_0^\pi f\,2\pi\sin\theta\,\mathrm d\theta = \frac{\kappa}{2\sinh\kappa}\cdot\frac{e^\kappa - e^{-\kappa}}{\kappa} = 1. Richtingen worden als eenheidsvectoren gemiddeld, nooit als hoeken: het gemiddelde van de declinaties 350∘350{}^{\circ} en 10∘10{}^{\circ} is noord, niet zuid. Een grote kk betekent een strakke cluster, en α95\alpha_{95} krimpt ongeveer evenredig met 1/kN1/\sqrt{kN}.

Methode 4.15 (Veldtests van een paleomagnetische richting)

Men kan met een test in het veld aantonen dat een ChRM ouder is dan een gebeurtenis.

  • Plooitest. Vindplaatsen op beide flanken van een plooi: als de richtingen beter clusteren nadat de lagen teruggedraaid zijn naar horizontaal dan ervoor, is de magnetisatie ouder dan de plooiing.
  • Omkeringstest. Bevat een formatie beide polariteiten, dan moet het gemiddelde van de omgekeerde richtingen, omgekeerd, samenvallen met het gemiddelde van de normale binnen hun betrouwbaarheidskegels; een viskeuze overprint die niet verwijderd is, zou beide in dezelfde zin verschuiven en de antipodaliteit verbreken.
  • Gebakken-contacttest. Een gang verhit het gesteente ernaast, dat een nieuwe TRM verwerft: als het gebakken gesteente de richting van de gang draagt en het verder weg gelegen gesteente een andere, dan is de magnetisatie van de gang sinds het afkoelen niet teruggezet.

4.4 Polen en platen

Als het veld, gemiddeld over enkele duizenden jaren om de seculaire variatie uit te middelen, een geocentrische axiale dipool is, dan zet de dipoolformule de gemiddelde inclinatie om in een breedte, en wijst de gemiddelde declinatie naar de pool.

Definitie 4.16 (Paleobreedte; virtuele geomagnetische pool)

De paleobreedte van een vindplaats is de breedte λ\lambda die volgt uit haar gemiddelde inclinatie via de dipoolformule, tan⁡I=2tan⁡λ\tan I = 2\tan\lambda. De virtuele geomagnetische pool (VGP) van een richting (D,I)(D, I) die op een vindplaats is waargenomen, is de noordpool van de geocentrische dipool die die richting daar zou opwekken.

Stelling 4.17 (De pool van een richting)

Laat een vindplaats breedte λs\lambda_s en lengte φs\varphi_s hebben en een gemiddelde richting (D,I)(D, I). Laat pp de paleocolatitude zijn, cot⁡p=12tan⁡I\cot p = \tfrac12\tan I, met 0≤p≤180∘0 \le p \le 180{}^{\circ}. De virtuele geomagnetische pool heeft breedte λp\lambda_p en lengte φp\varphi_p, gegeven door

sin⁡λp=sin⁡λscos⁡p+cos⁡λssin⁡pcos⁡D,sin⁡β=sin⁡psin⁡Dcos⁡λp,\sin\lambda_p = \sin\lambda_s\cos p + \cos\lambda_s\sin p\cos D,\qquad \sin\beta = \frac{\sin p\sin D}{\cos\lambda_p},

met φp=φs+β\varphi_p = \varphi_s + \beta als cos⁡p≥sin⁡λssin⁡λp\cos p \ge \sin\lambda_s\sin\lambda_p, en φp=φs+180∘−β\varphi_p = \varphi_s + 180{}^{\circ} - \beta anders.

Bewijs. Op de bol vormen de vindplaats SS, de geografische noordpool NN en de VGP PP een bolvormige driehoek (figuur hieronder). De zijde NSNS is de colatitude van de vindplaats, 90∘−λs90{}^{\circ} - \lambda_s; de zijde SPSP is de hoekafstand van de vindplaats tot de pool van de dipool, die de dipoolformule geeft als pp, aangezien de magnetische breedte van de vindplaats 90∘−p90{}^{\circ} - p is; en de hoek bij SS, tussen de richtingen naar NN en naar PP, is de declinatie DD. De cosinusregel voor de zijde NP=90∘−λpNP = 90{}^{\circ} - \lambda_p levert de eerste betrekking; de sinusregel, sin⁡β/sin⁡p=sin⁡D/sin⁡(90∘−λp)\sin\beta/\sin p = \sin D/\sin(90{}^{\circ} - \lambda_p), levert de hoek β\beta bij NN, het verschil in lengte. De sinus legt β\beta alleen vast tot op de supplementaire hoek: de cosinusregel voor de zijde SPSP, cos⁡p=sin⁡λssin⁡λp+cos⁡λscos⁡λpcos⁡β\cos p = \sin\lambda_s \sin\lambda_p + \cos\lambda_s\cos\lambda_p\cos\beta, toont aan dat cos⁡β≥0\cos\beta \ge 0 precies dan wanneer cos⁡p≥sin⁡λssin⁡λp\cos p \ge \sin\lambda_s\sin\lambda_p. ∎

De bolvormige driehoek van de virtuele geomagnetische pool: geografische noordpool N, vindplaats S, pool P. De hoek bij de vindplaats is de declinatie D; de hoek bij de noordpool N is het verschil in lengtegraad .
De bolvormige driehoek van de virtuele geomagnetische pool: geografische noordpool NN, vindplaats SS, pool PP. De hoek bij de vindplaats is de declinatie DD; de hoek bij de noordpool NN is het verschil in lengtegraad β\beta.

De kegel α95\alpha_{95} rond de gemiddelde richting wordt een ovaal rond de pool, met halve assen dp=12α95(1+3cos⁡2p)\mathrm dp = \tfrac12\alpha_{95}(1 + 3\cos^2p) langs de grootcirkel vanaf de vindplaats, en dm=α95sin⁡p/cos⁡I\mathrm dm = \alpha_{95}\sin p/\cos I dwars daarop. Een pool die gemiddeld is over meerdere vindplaatsen van een formatie, elk een VGP, is een paleomagnetische pool met een eigen kegel, A95A_{95}. Gesteenten van dezelfde ouderdom op één starre plaat geven dezelfde pool; gesteenten van verschillende ouderdom geven verschillende polen.

Definitie 4.18 (Schijnbare poolverschuiving; schijnbaar poolverschuivingspad)

Schijnbare poolverschuiving is de verandering, met de ouderdom van de gesteenten, van de paleomagnetische pool van een plaat gezien vanaf die plaat, die vastgehouden wordt. Een schijnbaar poolverschuivingspad is de kromme die de paleomagnetische polen van een plaat in volgorde van ouderdom verbindt.

De pool beweegt zelf niet: de plaat beweegt, en het pad registreert die beweging in omgekeerde richting. Twee platen die uit elkaar zijn bewogen, hebben verschillende paden voor de tijd vóór hun scheiding, en elk pad wordt door de eindige rotatie die de oceaan ertussen heeft gesloten, op het andere gedraaid. Een paleomagnetische pool legt de breedtegraad van een blok en zijn oriëntatie vast, maar niet zijn lengtegraad: een rotatie van het blok om de draaias verandert geen van beide.

Een verplaatst blok, een terrein, wordt getoetst aan het pad van de stabiele plaat waar het nu deel van uitmaakt. Uit de referentiepool van de juiste ouderdom berekent men de verwachte richting (Dx,Ix)(D_x, I_x) op de vindplaats, door Stelling 4.17 achterstevoren te lezen. Een rotatie van het blok met de klok mee om een verticale as wordt gemeten met R=Do−DxR = D_o - D_x, de verplaatsing naar de pool toe met P=λx−λoP = \lambda_x - \lambda_o, het verschil tussen verwachte en waargenomen paleobreedten; de fout van een declinatie is ongeveer ΔD=arcsin⁡(sin⁡α95/cos⁡I)\Delta D = \arcsin(\sin\alpha_{95}/\cos I), en fouten worden gecombineerd als de wortel uit de som van hun kwadraten. De Hawaï–Keizerketen werd in Boek 3 aan deze toets onderworpen, in het hoofdstuk over hotspots: de paleobreedten gemeten op de Keizerzeebergen toonden de hotspot zelf naar het zuiden bewegen, met meer dan 40 mm/yr40\,\mathrm{mm}/\mathrm{yr} van 81 tot 47 Ma47\,\mathrm{Ma}.

Poolprojectie met gelijke oppervlakten van het noordpoolgebied, tot 45 NB (buitenste cirkel), breedtecirkels elke 10 vanaf 50 NB, meridianen elke 30 (aangegeven in graden oost). Vierkanten: het schijnbare poolverschuivingspad van een verzonnen stabiel continent, met de A_95-cirkels van zijn polen; driehoeken: de virtuele geomagnetische polen van de drie formaties van het verzonnen Kesta-terrein uit het weekendprobleem, met hun _95-cirkels (gestippeld). De geografische pool is de stip in het midden.
Poolprojectie met gelijke oppervlakten van het noordpoolgebied, tot 45∘45{}^{\circ} NB (buitenste cirkel), breedtecirkels elke 10∘10{}^{\circ} vanaf 50∘50{}^{\circ} NB, meridianen elke 30∘30{}^{\circ} (aangegeven in graden oost). Vierkanten: het schijnbare poolverschuivingspad van een verzonnen stabiel continent, met de A95A_{95}-cirkels van zijn polen; driehoeken: de virtuele geomagnetische polen van de drie formaties van het verzonnen Kesta-terrein uit het weekendprobleem, met hun α95\alpha_{95}-cirkels (gestippeld). De geografische pool is de stip in het midden.

Een pool die beweegt ten opzichte van de draaias kan betekenen dat de platen bewogen zijn, of dat de hele vaste Aarde onder de draaias is gedraaid.

Definitie 4.19 (Ware poolverschuiving)

Ware poolverschuiving is de rotatie van de hele vaste Aarde, mantel en korst samen, ten opzichte van haar draaias, aangedreven door veranderingen in de verdeling van massa in haar binnenste. Zij verplaatst alle paleomagnetische polen van alle platen samen, met dezelfde rotatie, terwijl plaatbeweging de polen van elke plaat afzonderlijk verplaatst.

Zij wordt alleen gedetecteerd door de paden van alle platen te vergelijken binnen een referentiekader dat onafhankelijk is van de draaias, zoals dat van de hotspots; hoe groot ze is geweest, is nog steeds omstreden, en er wordt hier geen snelheid genoemd.

4.5 Magnetostratigrafie

Definitie 4.20 (Magnetostratigrafie)

Magnetostratigrafie is de onderverdeling van een reeks gesteenten in polariteitszones, intervallen van uniforme normale of omgekeerde polariteit, en de datering van de reeks door het patroon van zones te laten overeenkomen met dat van de geomagnetische polariteitsschaal.

Een omkering wordt, binnen de resolutie van de bemonstering, overal op hetzelfde moment vastgelegd: een omkeringsgrens in een profiel is een isochroon oppervlak. Het patroon van zones, hun aantal en hun relatieve dikten, is wat wordt vergeleken; een enkele zone, op zichzelf, kan elk chron van die polariteit zijn. Eén of twee onafhankelijke ouderdommen, uit een vulkanische as of een fossiel, verankeren de overeenkomst, waarna elke omkeringsgrens het profiel dateert. In de figuur hieronder toont een profiel van 64 m64\,\mathrm{m} toont drie normale zones; gepast op de normale subchrons van chron C2An, van 3.596 tot 2.581 Ma2.581\,\mathrm{Ma}, geven de grenzen een sedimentatiesnelheid van ongeveer 40 m/Myr40\,\mathrm{m}/\mathrm{Myr} onder 20 m20\,\mathrm{m} en 55 m/Myr55\,\mathrm{m}/\mathrm{Myr} boven.

Een magnetostratigrafische correlatie: de polariteitszones van een verzonnen sectie (zwart, normaal; wit, omgekeerd) gekoppeld aan de polariteitstijdschaal tussen 3.8 en 2.3\, Ma. De gestippelde lijnen verbinden elke omkeringsgrens met zijn ouderdom.
Een magnetostratigrafische correlatie: de polariteitszones van een verzonnen sectie (zwart, normaal; wit, omgekeerd) gekoppeld aan de polariteitstijdschaal tussen 3.8 en 2.3 Ma2.3\,\mathrm{Ma}. De gestippelde lijnen verbinden elke omkeringsgrens met zijn ouderdom.

In het laboratorium — Een afgeschermd paleomagnetisch laboratorium

Het laden van monsters op de monsterlade van een magnetometer in een afgeschermde ruimte (door AI gegenereerde illustratie).

Het laboratorium is een ruimte binnen een magnetische afscherming, lagen van metaal met een hoge permeabiliteit die het aardveld terugbrengen tot een klein deel van zijn waarde, zodat de monsters tijdens het hanteren geen nieuwe remanentie verwerven. Een paleomagnetist laadt de monsters, kleine cilinders van gesteente die elk zijn voorzien van een oriëntatiepijl, op een lange lade die ze een voor een meeneemt naar een cryogene magnetometer, waarvan de supergeleidende sensoren, gekoeld door vloeibaar helium, de drie componenten van remanenties meten die veel zwakker zijn dan een kompas kan waarnemen. Tussen de metingen gaan de monsters naar een afgeschermde oven voor de volgende verhittingsstap, of door de spoelen van een wisselveld-demagnetiseerder die op dezelfde lijn is aangesloten, en de software tekent elk Zijderveld-diagram terwijl de stappen zich opstapelen.

4.6 Oefeningen

Oefening 4.1 ★

Bereken in een geocentrisch axiaal dipoolveld de inclinatie op breedtegraad 30∘30{}^{\circ} NB, en de breedtegraad van een locatie waar de inclinatie −20∘-20{}^{\circ} is.

Oplossing

Oplossing van Oefening 4.1.

tan⁡I=2tan⁡30∘=1.1547\tan I = 2\tan30{}^{\circ} = 1.1547, dus I=49.1∘I = 49.1{}^{\circ}. Omgekeerd, tan⁡λ=12tan⁡(−20∘)=−0.1820\tan\lambda = \tfrac12\tan(-20{}^{\circ}) = -0.1820, dus λ=−10.3∘\lambda = -10.3{}^{\circ}, dat is 10.3∘10.3{}^{\circ} Z.

Oefening 4.2 ★

Op een station heeft het veld de componenten X=20 000 nTX = 20\,000\,\mathrm{nT}, Y=−2000 nTY = -2000\,\mathrm{nT} en Z=40 000 nTZ = 40\,000\,\mathrm{nT}. Bereken HH, FF, de magnetische declinatie en de inclinatie.

Oplossing

Oplossing van Oefening 4.2.

H=200002+20002=20 100 nTH = \sqrt{20000^2 + 2000^2} = 20\,100\,\mathrm{nT}; F=200002+20002+400002=44 770 nTF = \sqrt{20000^2 + 2000^2 + 40000^2} = 44\,770\,\mathrm{nT}. De magnetische declinatie is arctan⁡(−2000/20000)=−5.7∘\arctan(-2000/20000) = -5.7{}^{\circ}, 5.7∘5.7{}^{\circ} naar het westen van het noorden; de inclinatie arctan⁡(40000/20100)=63.3∘\arctan(40000/20100) = 63.3{}^{\circ}, naar beneden.

Oefening 4.3 ★

In 2025 geldt g10=−29 350 nTg_1^0 = -29\,350\,\mathrm{nT}. Neem het veld als een axiale dipool en bereken het veld aan de evenaar en aan de polen, en het dipoolmoment (a=6371.2 kma = 6371.2\,\mathrm{km}, μ0=4π×10−7 H m−1\mu_0 = 4\pi\times10^{-7}\,\mathrm{H}\,\mathrm{m}^{-1}).

Oplossing

Oplossing van Oefening 4.3.

B0=−g10=29 350 nTB_0 = -g_1^0 = 29\,350\,\mathrm{nT} aan de evenaar, en 2B0=58 700 nT2B_0 = 58\,700\,\mathrm{nT} bij de polen. m=4πa3B0/μ0=a3B0/10−7=2.5862×1020×2.935×10−5/10−7=7.59×1022 A m2m = 4\pi a^3B_0/\mu_0 = a^3B_0/10^{-7} = 2.5862 \times 10^{20}\times 2.935 \times 10^{-5}/10^{-7} = 7.59 \times 10^{22}\,\mathrm{A}\,\mathrm{m}^{2}, iets minder dan de 7.69×1022 A m27.69 \times 10^{22}\,\mathrm{A}\,\mathrm{m}^{2} van de volledige dipool, die g11g_1^1 en h11h_1^1 toevoegt.

Oefening 4.4 ★

Benoem het soort remanentie dat wordt verworven en wanneer het het veld vastlegt, voor: (a) een lavastroom die afkoelt; (b) silt die bezinkt in een rustig meer; (c) een zandsteen waarvan de ijzeroxidecoatings tijdens de vroege begraving in hematiet veranderen; (d) een basalt die een miljoen jaar in het huidige veld ligt.

Oplossing

Oplossing van Oefening 4.4.

(a) Een thermoremanente magnetisatie, die het veld vastlegt toen de lava afkoelde door de blokkeringstemperaturen van de korrels heen. (b) Een detritische remanente magnetisatie, die het veld vastlegt binnen de tijd die het silt nodig had om te bezinken en door compactie vergrendeld te raken. (c) Een chemische remanente magnetisatie, die het veld vastlegt tijdens de groei van het hematiet, tijdens de vroege begraving en dus enigszins na de afzetting. (d) Een viskeuze remanente magnetisatie langs het huidige veld, toegevoegd aan de oorspronkelijke TRM door de korrels met de kortste relaxatietijd; zij wordt verwijderd door de eerste stappen van de demagnetisatie.

Oefening 4.5 ★★

Vier monsters van een locatie geven de richtingen (D,I)(D, I): (12∘,48∘)(12{}^{\circ}, 48{}^{\circ}), (22∘,44∘)(22{}^{\circ}, 44{}^{\circ}), (16∘,57∘)(16{}^{\circ}, 57{}^{\circ}) en (6∘,53∘)(6{}^{\circ}, 53{}^{\circ}). Bereken hun eenheidsvectoren, de lengte van de resultante RR, de gemiddelde richting, kk en α95\alpha_{95}.

Oplossing

Oplossing van Oefening 4.5.

De eenheidsvectoren (cos⁡Icos⁡D, cos⁡Isin⁡D, sin⁡I)(\cos I\cos D,\ \cos I\sin D,\ \sin I) zijn

(0.6545,0.1391,0.7431),(0.6670,0.2695,0.6947),(0.5235,0.1501,0.8387),(0.5985,0.0629,0.7986).\begin{gather*} (0.6545, 0.1391, 0.7431),\quad (0.6670, 0.2695, 0.6947),\\ (0.5235, 0.1501, 0.8387),\quad (0.5985, 0.0629, 0.7986). \end{gather*}

De som is (2.4435,0.6216,3.0751)(2.4435, 0.6216, 3.0751), met lengte R=3.9766R = 3.9766. Gemiddelde: D=arctan⁡(0.6216/2.4435)=14.3∘D = \arctan(0.6216/2.4435) = 14.3{}^{\circ}, I=arcsin⁡(3.0751/3.9766)=50.7∘I = \arcsin(3.0751/3.9766) = 50.7{}^{\circ}. k=3/(4−3.9766)=128k = 3/(4 - 3.9766) = 128. Met 201/3−1=1.714420^{1/3} - 1 = 1.7144, cos⁡α95=1−(0.0234/3.9766)×1.7144=0.98991\cos\alpha_{95} = 1 - (0.0234/3.9766)\times1.7144 = 0.98991, dus α95=8.1∘\alpha_{95} = 8.1{}^{\circ}: met vier monsters blijft de kegel breed, hoe strak de groep ook is.

Oefening 4.6 ★★

Een locatie op 40.0∘40.0{}^{\circ} NB, 10.0∘10.0{}^{\circ} OL geeft de gemiddelde richting D=20∘D = 20{}^{\circ}, I=50∘I = 50{}^{\circ}. Bereken de paleobreedte en de virtuele geomagnetische pool.

Oplossing

Oplossing van Oefening 4.6.

λ=arctan⁡(12tan⁡50∘)=30.79∘\lambda = \arctan(\tfrac12\tan50{}^{\circ}) = 30.79{}^{\circ}, dus p=59.21∘p = 59.21{}^{\circ}. sin⁡λp=sin⁡40∘cos⁡59.21∘+cos⁡40∘sin⁡59.21∘cos⁡20∘=0.9474\sin\lambda_p = \sin40{}^{\circ}\cos59.21{}^{\circ} + \cos40{}^{\circ}\sin59.21{}^{\circ} \cos20{}^{\circ} = 0.9474, λp=71.34∘\lambda_p = 71.34{}^{\circ} NB. sin⁡β=sin⁡59.21∘sin⁡20∘/cos⁡71.34∘=0.9184\sin\beta = \sin59.21{}^{\circ}\sin20{}^{\circ}/\cos71.34{}^{\circ} = 0.9184, β=66.66∘\beta = 66.66{}^{\circ}. Omdat cos⁡p=0.5119<sin⁡40∘×0.9474=0.6090\cos p = 0.5119 < \sin40{}^{\circ}\times0.9474 = 0.6090, is het tweede geval van toepassing: φp=10+180−66.66=123.34∘\varphi_p = 10 + 180 - 66.66 = 123.34{}^{\circ} OL. De pool ligt op 71.3∘71.3{}^{\circ} NB, 123.3∘123.3{}^{\circ} OL.

Oefening 4.7 ★★

Het IGRF geeft een dipoolmoment van 8.32×1022 A m28.32 \times 10^{22}\,\mathrm{A}\,\mathrm{m}^{2} in 1900 en 7.69×1022 A m27.69 \times 10^{22}\,\mathrm{A}\,\mathrm{m}^{2} in 2025, en een veranderingssnelheid van g10g_1^0 van +12.6 nT/yr+12.6\,\mathrm{nT}/\mathrm{yr} in 2025. (a) Bereken de relatieve verandering van het moment per eeuw. (b) Als g10g_1^0 die snelheid zou behouden, na hoeveel jaar zou het nul bereiken? Waarom is dit geen voorspelling?

Oplossing

Oplossing van Oefening 4.7.

(a) (7.69−8.32)/8.32=−7.6 %(7.69 - 8.32)/8.32 = -7.6\,\% in 125 jaar, of −6.1 %-6.1\,\% per eeuw. (b) 29350/12.6=232929350/12.6 = 2329, ongeveer 2300 jaar. Het is geen voorspelling: de snelheid zelf verandert van het ene vijfjaarmodel naar het volgende; het moment is in het verleden zowel hoger als lager geweest zonder om te keren, en het veld keert op onregelmatige tijdstippen om die geen enkele snelheid voorspelt.

Oefening 4.8 ★★

De resterende remanentie van een monster, in mA/m\mathrm{mA}/\mathrm{m} (noord, oost, neer), is (32.2,18.1,38.2)(32.2, 18.1, 38.2) na verhitting tot 20 ∘C20\,{}^{\circ}\mathrm{C}, (25.8,18.1,24.6)(25.8, 18.1, 24.6) na 300 ∘C300\,{}^{\circ}\mathrm{C}, en (0.0,0.0,0.0)(0.0, 0.0, 0.0) na 580 ∘C580\,{}^{\circ}\mathrm{C}; het Zijderveld-diagram is recht tussen elk paar. Bereken de richting en intensiteit van de component die tussen 20 en 300 ∘C300\,{}^{\circ}\mathrm{C} wordt verwijderd, en van die welke boven 300 ∘C300\,{}^{\circ}\mathrm{C} wordt verwijderd. Welke is de ChRM, en wat is de andere?

Oplossing

Oplossing van Oefening 4.8.

Tussen 20 en 300 ∘C300\,{}^{\circ}\mathrm{C} verloor de vector (6.4,0.0,13.6)(6.4, 0.0, 13.6): een component van 15.0 mA/m15.0\,\mathrm{mA}/\mathrm{m}, met D=0.0∘D = 0.0{}^{\circ} en I=arctan⁡(13.6/6.4)=64.8∘I = \arctan(13.6/6.4) = 64.8{}^{\circ}. Boven 300 ∘C300\,{}^{\circ}\mathrm{C} verloor hij (25.8,18.1,24.6)(25.8, 18.1, 24.6): 40.0 mA/m40.0\,\mathrm{mA}/\mathrm{m}, D=arctan⁡(18.1/25.8)=35.1∘D = \arctan(18.1/25.8) = 35.1{}^{\circ}, I=arctan⁡(24.6/31.5)=38.0∘I = \arctan(24.6/31.5) = 38.0{}^{\circ}. De tweede, die naar de oorsprong gaat, is de ChRM; de eerste, langs een steile, noordwaarts gerichte richting, bij lage temperatuur verwijderd, is een viskeuze overprint langs het recente veld.

Oefening 4.9 ★★

Een formatie geeft I=−35∘I = -35{}^{\circ} met α95=5∘\alpha_{95} = 5{}^{\circ}. Bereken de paleobreedte, het bereik van paleobreedtes dat is toegestaan door I±α95I \pm \alpha_{95}, en de halve assen dp\mathrm dp en dm\mathrm dm van de betrouwbaarheidsovaal van de pool.

Oplossing

Oplossing van Oefening 4.9.

λ=arctan⁡(12tan⁡(−35∘))=−19.3∘\lambda = \arctan(\tfrac12\tan(-35{}^{\circ})) = -19.3{}^{\circ}. Voor I=−40∘I = -40{}^{\circ} en −30∘-30{}^{\circ}: −22.8∘-22.8{}^{\circ} en −16.1∘-16.1{}^{\circ}, een interval van 6.7∘6.7{}^{\circ}. p=90∘−λ=109.3∘p = 90{}^{\circ} - \lambda = 109.3{}^{\circ}, cos⁡2p=0.1096\cos^2p = 0.1096: dp=2.5×(1+0.329)=3.3∘\mathrm dp = 2.5\times(1 + 0.329) = 3.3{}^{\circ} en dm=5×sin⁡109.3∘/cos⁡35∘=5.8∘\mathrm dm = 5\times\sin109.3{}^{\circ}/\cos35{}^{\circ} = 5.8{}^{\circ}.

Oefening 4.10 ★★★

(a) Bepaal uit Stelling 4.3 de breedtegraad waarop de inclinatie 45∘45{}^{\circ} is, en de verhouding van het veld aan de polen tot dat aan de evenaar. (b) Toon aan dat dI/dλ=2/(1+3sin⁡2λ)\mathrm dI/\mathrm d\lambda = 2/(1 + 3\sin^2\lambda), en leid af waar een fout in de inclinatie de grootste fout in de paleobreedte oplevert. (c) Waarom zijn paleobreedtes die uit lage inclinaties volgen het nauwkeurigst uitgedrukt in graden?

Oplossing

Oplossing van Oefening 4.10.

(a) tan⁡λ=12\tan\lambda = \tfrac12, λ=26.6∘\lambda = 26.6{}^{\circ}. FF is B0B_0 aan de evenaar en 2B02B_0 bij de polen: een verhouding van 2. (b) tan⁡I=2tan⁡λ\tan I = 2\tan\lambda levert (1+tan⁡2I) dI=2(1+tan⁡2λ) dλ(1 + \tan^2I)\,\mathrm dI = 2(1 + \tan^2\lambda)\,\mathrm d\lambda, dus

dIdλ=2/cos⁡2λ1+4tan⁡2λ=2cos⁡2λ+4sin⁡2λ=21+3sin⁡2λ.\frac{\mathrm dI}{\mathrm d\lambda} = \frac{2/\cos^2\lambda}{1 + 4\tan^2\lambda} = \frac{2}{\cos^2\lambda + 4\sin^2\lambda} = \frac{2}{1 + 3\sin^2\lambda}.

Een fout δI\delta I levert δλ=12(1+3sin⁡2λ) δI\delta\lambda = \tfrac12(1 + 3\sin^2\lambda)\,\delta I, van 12δI\tfrac12\delta I aan de evenaar tot 2 δI2\,\delta I bij de polen: de grootste fout ligt nabij de polen. (c) Nabij de evenaar verandert de inclinatie twee keer zo snel als de breedtegraad, dus een gegeven fout in de inclinatie kost de helft zoveel breedtegraad; nabij de polen verandert de inclinatie nauwelijks met de breedtegraad.

Oefening 4.11 ★★★

Een locatie ligt op 35.0∘35.0{}^{\circ} Z, 150.0∘150.0{}^{\circ} O op een stabiele plaat waarvan de referentiepool, voor de ouderdom van een formatie die daar is bemonsterd, ligt op 70.0∘70.0{}^{\circ} NB, 300.0∘300.0{}^{\circ} OL. (a) Bereken de hoekafstand pp van de locatie tot de pool, de verwachte inclinatie en declinatie, en de verwachte paleobreedte. (b) De formatie geeft Do=30.0∘D_o = 30.0{}^{\circ}, Io=−60.0∘I_o = -60.0{}^{\circ}. Bereken de rotatie RR en de poolwaartse verschuiving PP, en zeg wat ze betekenen.

Oplossing

Oplossing van Oefening 4.11.

(a) cos⁡p=sin⁡(−35∘)sin⁡70∘+cos⁡35∘cos⁡70∘cos⁡150∘=−0.5390−0.2426=−0.7816\cos p = \sin(-35{}^{\circ})\sin70{}^{\circ} + \cos35{}^{\circ}\cos70{}^{\circ}\cos150{}^{\circ} = -0.5390 - 0.2426 = -0.7816, p=141.41∘p = 141.41{}^{\circ}. Dan tan⁡Ix=2cot⁡p=−2.506\tan I_x = 2\cot p = -2.506, Ix=−68.2∘I_x = -68.2{}^{\circ}, en λx=90∘−p=−51.4∘\lambda_x = 90{}^{\circ} - p = -51.4{}^{\circ}. cos⁡Dx=(sin⁡70∘−sin⁡(−35∘)cos⁡p)/(cos⁡35∘sin⁡p)=0.9617\cos D_x = (\sin70{}^{\circ} - \sin(-35{}^{\circ})\cos p)/(\cos35{}^{\circ}\sin p) = 0.9617, en omdat de pool oostelijk ligt van de meridiaan van de locatie (sin⁡150∘>0\sin150{}^{\circ} > 0), Dx=15.9∘D_x = 15.9{}^{\circ}. (b) λo=arctan⁡(12tan⁡(−60∘))=−40.9∘\lambda_o = \arctan(\tfrac12\tan(-60{}^{\circ})) = -40.9{}^{\circ}. R=30.0−15.9=14.1∘R = 30.0 - 15.9 = 14.1{}^{\circ}: een rotatie van het blok met de klok mee rond een verticale as. P=λx−λo=−51.4+40.9=−10.5∘P = \lambda_x - \lambda_o = -51.4 + 40.9 = -10.5{}^{\circ}: de formatie lag 10.5∘10.5{}^{\circ} verder naar het noorden dan de stabiele plaat aangeeft dat ze zou moeten, dus het blok is sindsdien ongeveer 10.5∘10.5{}^{\circ} verschoven, zo’n 1170 km1170\,\mathrm{km} naar het zuiden ten opzichte van de plaat.

Oefening 4.12 ★★★

(a) Toon uit Propositie 4.13 aan dat de kans dat een richting binnen θ\theta van het gemiddelde ligt, is gelijk aan (eκ−eκcos⁡θ)/(eκ−e−κ)(e^\kappa - e^{\kappa\cos\theta})/(e^\kappa - e^{-\kappa}). (b) Toon voor grote κ\kappa en kleine θ\theta aan dat dit ongeveer 1−e−κθ2/21 - e^{-\kappa\theta^2/2} is, en leid af dat 63 %63\,\% van de richtingen binnen θ63≈81∘/κ\theta_{63} \approx 81{}^{\circ}/\sqrt\kappa ligt. (c) Bereken θ63\theta_{63} voor κ=50\kappa = 50, en controleer de 63 %63\,\% met de exacte formule.

Oplossing

Oplossing van Oefening 4.12.

(a) Met dΩ=2πsin⁡θ′ dθ′\mathrm d\Omega = 2\pi\sin\theta'\,\mathrm d\theta',

∫0θκeκcos⁡θ′4πsinh⁡κ 2πsin⁡θ′ dθ′=eκ−eκcos⁡θ2sinh⁡κ=eκ−eκcos⁡θeκ−e−κ.\int_0^\theta\frac{\kappa e^{\kappa\cos\theta'}}{4\pi\sinh\kappa}\,2\pi\sin\theta'\, \mathrm d\theta' = \frac{e^\kappa - e^{\kappa\cos\theta}}{2\sinh\kappa} = \frac{e^\kappa - e^{\kappa\cos\theta}}{e^\kappa - e^{-\kappa}} .

(b) Voor grote κ\kappa is e−κe^{-\kappa} verwaarloosbaar; delen door eκe^\kappa levert 1−e−κ(1−cos⁡θ)≈1−e−κθ2/21 - e^{-\kappa(1 - \cos\theta)} \approx 1 - e^{-\kappa\theta^2/2}. Dit is gelijk aan 1−e−1=63 %1 - e^{-1} = 63\,\% wanneer κθ2/2=1\kappa\theta^2/2 = 1, θ=2/κ\theta = \sqrt{2/\kappa} radialen, oftewel 81∘/κ81{}^{\circ}/\sqrt\kappa. (c) 81/50=11.5∘81/\sqrt{50} = 11.5{}^{\circ}; de exacte formule geeft 63.1 %63.1\,\%.

4.7 Probleem: De roodbedden van Kesta

Probleem 4.1

Weekendprobleem — drie formaties van een verplaatst terrein, hun richtingen, polen en rotaties, en de ouderdom van de jongste

Kesta is een verzonnen terrein, nu op 47.0∘47.0{}^{\circ} NB, 15.0∘15.0{}^{\circ} OL, aan de rand van een verzonnen stabiel continent. Het ontsluit drie roodbedformaties: A, 60 Ma60\,\mathrm{Ma}; B, 30 Ma30\,\mathrm{Ma}; C, van ongeveer 3.7 tot 2.4 Ma2.4\,\mathrm{Ma}. Elk werd op meerdere locaties bemonsterd, en elk specimen werd thermisch gedemagnetiseerd; de ChRM wordt gedragen door hematiet, en de formaties liggen horizontaal.

Deel I — Gemiddelde richtingen. De zes locaties van formatie A geven (D,I)(D, I): (50∘,37∘)(50{}^{\circ}, 37{}^{\circ}), (28∘,37∘)(28{}^{\circ}, 37{}^{\circ}), (21∘,42∘)(21{}^{\circ}, 42{}^{\circ}), (24∘,27∘)(24{}^{\circ}, 27{}^{\circ}), (37∘,39∘)(37{}^{\circ}, 39{}^{\circ}), (46∘,38∘)(46{}^{\circ}, 38{}^{\circ}).

  1. Waarom kan een ChRM die door hematiet wordt gedragen in een roodbed jonger zijn dan de laag? Wat zou aantonen dat hij niet veel jonger is?
  2. Bereken de eenheidsvector van de eerste richting.
  3. De som van de zes eenheidsvectoren is (3.8943,2.6500,3.5717)(3.8943, 2.6500, 3.5717). Bereken RR en de gemiddelde richting van formatie A.
  4. Bereken kk en α95\alpha_{95} voor formatie A.
  5. Formatie B geeft, uit 7 locaties, D=23.9∘D = 23.9{}^{\circ}, I=51.5∘I = 51.5{}^{\circ} en R=6.8979R = 6.8979; formatie C, uit 8 locaties, D=3.4∘D = 3.4{}^{\circ}, I=63.7∘I = 63.7{}^{\circ} en R=7.8707R = 7.8707. Bereken kk en α95\alpha_{95} van elk.
  6. Welke formatie is het best gegroepeerd? Welk gemiddelde is het best bepaald?

Deel II — Paleobreedtes en polen.

  1. Bereken de paleobreedtes van A, B en C.
  2. Bereken voor formatie A de paleocolatitude pp en de breedtegraad van de virtuele geomagnetische pool.
  3. Bereken de lengtegraad van die pool, en controleer welk geval uit Stelling 4.17 van toepassing is.
  4. Bereken dp\mathrm dp en dm\mathrm dm voor de pool van A.
  5. De polen van B en C liggen op 66.5∘66.5{}^{\circ} NB, 135.6∘135.6{}^{\circ} OL en 87.1∘87.1{}^{\circ} NB, 139.1∘139.1{}^{\circ} OL. Beschrijf het pad van de drie polen als functie van de ouderdom.

Deel III — De ouderdom van formatie C. In formatie C, vanaf de basis zijn de specimens omgekeerd tot 4.2 m4.2\,\mathrm{m}, normaal tot 14.8 m14.8\,\mathrm{m}, omgekeerd tot 19.7 m19.7\,\mathrm{m}, normaal tot 24.6 m24.6\,\mathrm{m}, omgekeerd tot 29.2 m29.2\,\mathrm{m} en normaal tot 54.0 m54.0\,\mathrm{m}, daarna omgekeerd tot de top bij 64.0 m64.0\,\mathrm{m}. Een asslaag op 22 m22\,\mathrm{m} is gedateerd op ongeveer 3.2 Ma3.2\,\mathrm{Ma}. De normale subchronen van polariteitschron C2An omvatten 3.596–3.330, 3.207–3.116 en 3.032–2.581 Ma2.581\,\mathrm{Ma}.

  1. Koppel de polariteitszones aan de tijdschaal, met behulp van de asslaag.
  2. Bereken de accumulatiesnelheid tussen de eerste en de derde omkeringsgrens, en tussen de vierde en de zesde.
  3. Schat de ouderdom van de basis en de top van de sectie.
  4. Waarom is de snelheid een gemiddelde over honderdduizenden jaren, en wat zou haar kunnen doen veranderen?

Deel IV — Tegen het stabiele continent, en de drift. De referentiepolen van het stabiele continent liggen op 78.0∘78.0{}^{\circ} NB, 190.0∘190.0{}^{\circ} OL (A95=3.0∘A_{95} = 3.0{}^{\circ}) bij 60 Ma60\,\mathrm{Ma}; 83.5∘83.5{}^{\circ} NB, 175.0∘175.0{}^{\circ} OL (2.5∘2.5{}^{\circ}) bij 30 Ma30\,\mathrm{Ma}; en 87.5∘87.5{}^{\circ} NB, 150.0∘150.0{}^{\circ} OL (2.0∘2.0{}^{\circ}) bij 4 Ma4\,\mathrm{Ma}.

  1. Bereken de hoekafstand van Kesta tot de referentiepool van 60 Ma60\,\mathrm{Ma}, en de verwachte inclinatie en paleobreedte.
  2. Bereken de verwachte declinatie bij 60 Ma60\,\mathrm{Ma}.
  3. De verwachte richtingen bij 30 en 4 Ma4\,\mathrm{Ma} zijn (2.9∘,60.0∘)(2.9{}^{\circ}, 60.0{}^{\circ}) en (2.5∘,63.6∘)(2.5{}^{\circ}, 63.6{}^{\circ}). Bereken de poolwaartse verschuiving PP van A, B en C.
  4. Bereken de rotatie RR van A, B en C.
  5. Bereken de fouten ΔDo\Delta D_o en ΔDx\Delta D_x voor A, en de fout van de rotatie.
  6. De fouten van RR voor B en C zijn 13.8∘13.8{}^{\circ} en 17.9∘17.9{}^{\circ}. Welke rotaties zijn significant op 95 procent?
  7. Beschrijf de geschiedenis van Kesta ten opzichte van het stabiele continent.
  8. Bereken de noordwaartse drift van Kesta sinds formatie A, ten opzichte van de rotatie-as, in graden en in kilometer (111.19 km111.19\,\mathrm{km} per graad).
  9. Bereken de paleobreedtes die zijn toegestaan door I±α95I \pm \alpha_{95} voor formatie A, en neem als één standaardafwijking de halve breedte van dat interval gedeeld door 1.96. Geef de noordwaartse drift van Kesta sinds formatie A, in kilometer, met de bijbehorende standaardafwijking.
Oplossing

Oplossing van Probleem 4.1.

1. Hematiet in een rode laag kan groeien uit ijzeroxiden en hydroxiden lang nadat het sediment is afgezet, waardoor een chemische remanente magnetisatie van die latere ouderdom ontstaat. Een positieve plooitest, als de lagen kort na de afzetting zijn geplooid, of een positieve omkeringstest over dicht bij elkaar liggende lagen van beide polariteiten zou aantonen dat de remanente magnetisatie vroeg is.

2. (cos⁡37∘cos⁡50∘,cos⁡37∘sin⁡50∘,sin⁡37∘)=(0.5134,0.6118,0.6018)(\cos37{}^{\circ}\cos50{}^{\circ}, \cos37{}^{\circ}\sin50{}^{\circ}, \sin37{}^{\circ}) = (0.5134, 0.6118, 0.6018).

3. R=3.89432+2.65002+3.57172=5.9114R = \sqrt{3.8943^2 + 2.6500^2 + 3.5717^2} = 5.9114. D=arctan⁡(2.6500/3.8943)=34.2∘D = \arctan(2.6500/3.8943) = 34.2{}^{\circ}, I=arcsin⁡(3.5717/5.9114)=37.2∘I = \arcsin(3.5717/5.9114) = 37.2{}^{\circ}.

4. k=5/(6−5.9114)=56.5k = 5/(6 - 5.9114) = 56.5. Met 201/5−1=0.820620^{1/5} - 1 = 0.8206, cos⁡α95=1−(0.0886/5.9114)×0.8206=0.98771\cos\alpha_{95} = 1 - (0.0886/5.9114)\times0.8206 = 0.98771, α95=9.0∘\alpha_{95} = 9.0{}^{\circ}.

5. B: k=6/(7−6.8979)=58.8k = 6/(7 - 6.8979) = 58.8; 201/6−1=0.647520^{1/6} - 1 = 0.6475, cos⁡α95=1−(0.1021/6.8979)×0.6475=0.99042\cos\alpha_{95} = 1 - (0.1021/6.8979)\times0.6475 = 0.99042, α95=7.9∘\alpha_{95} = 7.9{}^{\circ}. C: k=7/(8−7.8707)=54.1k = 7/(8 - 7.8707) = 54.1; 201/7−1=0.534020^{1/7} - 1 = 0.5340, cos⁡α95=1−(0.1293/7.8707)×0.5340=0.99123\cos\alpha_{95} = 1 - (0.1293/7.8707)\times0.5340 = 0.99123, α95=7.6∘\alpha_{95} = 7.6{}^{\circ}.

6. B heeft de grootste kk, de strakste groepering; C heeft de kleinste α95\alpha_{95}, het best bepaalde gemiddelde, omdat het de meeste vindplaatsen telt. De drie liggen dicht bij elkaar: alle drie formaties zijn goed gegroepeerd.

7. λA=arctan⁡(12tan⁡37.17∘)=20.76∘\lambda_A = \arctan(\tfrac12\tan37.17{}^{\circ}) = 20.76{}^{\circ} N, λB=32.15∘\lambda_B = 32.15{}^{\circ} N, λC=45.33∘\lambda_C = 45.33{}^{\circ} N.

8. p=90−20.76=69.24∘p = 90 - 20.76 = 69.24{}^{\circ}. sin⁡λp=sin⁡47∘cos⁡69.24∘+cos⁡47∘sin⁡69.24∘cos⁡34.23∘=0.2593+0.5272=0.7865\sin\lambda_p = \sin47{}^{\circ}\cos69.24{}^{\circ} + \cos47{}^{\circ}\sin69.24{}^{\circ}\cos34.23{}^{\circ} = 0.2593 + 0.5272 = 0.7865: λp=51.86∘\lambda_p = 51.86{}^{\circ} N.

9. sin⁡β=sin⁡69.24∘sin⁡34.23∘/cos⁡51.86∘=0.8516\sin\beta = \sin69.24{}^{\circ}\sin34.23{}^{\circ}/\cos51.86{}^{\circ} = 0.8516, β=58.40∘\beta = 58.40{}^{\circ}. Omdat cos⁡p=0.3545<sin⁡47∘×0.7865=0.5752\cos p = 0.3545 < \sin47{}^{\circ}\times 0.7865 = 0.5752, geldt het tweede geval, cos⁡β<0\cos\beta < 0: φp=15.00+180−58.40=136.60∘\varphi_p = 15.00 + 180 - 58.40 = 136.60{}^{\circ} O.

10. cos⁡2p=0.1257\cos^2p = 0.1257: dp=4.50×1.377=6.2∘\mathrm dp = 4.50\times1.377 = 6.2{}^{\circ}; dm=8.99×sin⁡69.24∘/cos⁡37.17∘=10.5∘\mathrm dm = 8.99\times\sin69.24{}^{\circ}/\cos37.17{}^{\circ} = 10.5{}^{\circ}.

11. Van 51.9∘51.9{}^{\circ} N bij 60 Ma60\,\mathrm{Ma} naar 66.5∘66.5{}^{\circ} N bij 30 Ma30\,\mathrm{Ma} en 87.1∘87.1{}^{\circ} N bij 4 Ma4\,\mathrm{Ma}, op bijna dezelfde meridiaan (136∘136{}^{\circ}–139∘139{}^{\circ} O): een pad richting de geografische pool, binnen zijn kegel bereikt door de jongste formatie.

12. De aslaag, op 22 m22\,\mathrm{m}, ligt in de tweede normale zone, 19.7–24.6 m24.6\,\mathrm{m}; met een ouderdom nabij 3.2 Ma3.2\,\mathrm{Ma} is die zone C2An.2n (3.207–3.116). Dan is de zone 4.2–14.8 m14.8\,\mathrm{m} C2An.3n, de zone 29.2–54.0 m54.0\,\mathrm{m} C2An.1n, de omgekeerde basis hoort bij het chron onder C2An, en de omgekeerde top hoort bij het chron erboven.

13. (19.7−4.2)/(3.596−3.207)=15.5/0.389=39.8 m/Myr(19.7 - 4.2)/(3.596 - 3.207) = 15.5/0.389 = 39.8\,\mathrm{m}/\mathrm{Myr}; (54.0−24.6)/(3.116−2.581)=29.4/0.535=55.0 m/Myr(54.0 - 24.6)/(3.116 - 2.581) = 29.4/0.535 = 55.0\,\mathrm{m}/\mathrm{Myr}.

14. Basis: 3.596+4.2/39.8=3.70 Ma3.596 + 4.2/39.8 = 3.70\,\mathrm{Ma}. Top: 2.581−10.0/55.0=2.40 Ma2.581 - 10.0/55.0 = 2.40\,\mathrm{Ma}.

15. Alleen omkeringsgrenzen worden gedateerd, een paar honderdduizend jaar uit elkaar; ertussen wordt de snelheid constant verondersteld. Een verandering van de subsidentie, van de sedimenttoevoer of van het klimaat, of een hiaat in de afzetting, zou haar veranderen, en een hiaat zou zelfs een korte zone kunnen wegnemen en de match misleiden.

16. cos⁡p=sin⁡47∘sin⁡78∘+cos⁡47∘cos⁡78∘cos⁡175∘=0.7154−0.1413=0.5741\cos p = \sin47{}^{\circ}\sin78{}^{\circ} + \cos47{}^{\circ}\cos78{}^{\circ} \cos175{}^{\circ} = 0.7154 - 0.1413 = 0.5741, p=54.96∘p = 54.96{}^{\circ}. tan⁡Ix=2cot⁡p=1.4023\tan I_x = 2\cot p = 1.4023, Ix=54.5∘I_x = 54.5{}^{\circ}; λx=35.04∘\lambda_x = 35.04{}^{\circ} N.

17. cos⁡Dx=(sin⁡78∘−sin⁡47∘×0.5741)/(cos⁡47∘sin⁡54.96∘)=0.99976\cos D_x = (\sin78{}^{\circ} - \sin47{}^{\circ}\times0.5741)/ (\cos47{}^{\circ}\sin54.96{}^{\circ}) = 0.99976; de pool ligt oost van de meridiaan van de vindplaats (sin⁡175∘>0\sin175{}^{\circ} > 0), dus Dx=1.3∘D_x = 1.3{}^{\circ}.

18. Verwachte paleobreedtes uit IxI_x = 54.5∘54.5{}^{\circ}, 60.0∘60.0{}^{\circ} en 63.6∘63.6{}^{\circ}: 35.03∘35.03{}^{\circ}, 40.89∘40.89{}^{\circ} en 45.21∘45.21{}^{\circ}. PA=35.03−20.76=14.3∘P_A = 35.03 - 20.76 = 14.3{}^{\circ}, PB=40.89−32.15=8.7∘P_B = 40.89 - 32.15 = 8.7{}^{\circ}, PC=45.21−45.33=−0.1∘P_C = 45.21 - 45.33 = -0.1{}^{\circ}.

19. RA=34.2−1.3=32.9∘R_A = 34.2 - 1.3 = 32.9{}^{\circ}, RB=23.9−2.9=21.0∘R_B = 23.9 - 2.9 = 21.0{}^{\circ}, RC=3.4−2.5=0.9∘R_C = 3.4 - 2.5 = 0.9{}^{\circ}, alle met de klok mee.

20. ΔDo=arcsin⁡(sin⁡8.99∘/cos⁡37.17∘)=11.3∘\Delta D_o = \arcsin(\sin8.99{}^{\circ}/\cos37.17{}^{\circ}) = 11.3{}^{\circ}; ΔDx=arcsin⁡(sin⁡3.0∘/cos⁡54.5∘)=5.2∘\Delta D_x = \arcsin(\sin3.0{}^{\circ}/\cos54.5{}^{\circ}) = 5.2{}^{\circ}; ΔRA=11.32+5.22=12.4∘\Delta R_A = \sqrt{11.3^2 + 5.2^2} = 12.4{}^{\circ}.

21. RA=32.9∘±12.4∘R_A = 32.9{}^{\circ} \pm 12.4{}^{\circ} en RB=21.0∘±13.8∘R_B = 21.0{}^{\circ} \pm 13.8{}^{\circ} zijn significant; RC=0.9∘±17.9∘R_C = 0.9{}^{\circ} \pm 17.9{}^{\circ} is niet.

22. Zestig miljoen jaar geleden lag Kesta ongeveer 14∘14{}^{\circ} breedtegraad ten zuiden van waar het stabiele continent het plaatst, en was het ongeveer 33∘33{}^{\circ} tegen de klok in gedraaid ten opzichte van zijn huidige oriëntatie. Tegen 30 Ma30\,\mathrm{Ma} was het 6∘6{}^{\circ} van die 14∘14{}^{\circ} naar het noorden opgeschoven en 12∘12{}^{\circ} van die 33∘33{}^{\circ} met de klok mee gedraaid; tegen 4 Ma4\,\mathrm{Ma} had het zijn huidige plaats op het continent bereikt, zonder rotatie over. Het dreef langs de marge van het continent naar het noorden en draaide daarbij met de klok mee, en was tegen 4 Ma4\,\mathrm{Ma} op zijn plaats.

23. 47.00−20.76=26.24∘47.00 - 20.76 = 26.24{}^{\circ}, oftewel 26.24×111.19=2918 km26.24\times111.19 = 2918\,\mathrm{km}. Dit omvat de beweging van het stabiele continent zelf; ten opzichte daarvan heeft Kesta zich PA=14.3∘P_A = 14.3{}^{\circ} verplaatst.

24. I=28.18∘I = 28.18{}^{\circ} en 46.16∘46.16{}^{\circ} geven λ=15.00∘\lambda = 15.00{}^{\circ} en 27.50∘27.50{}^{\circ}: een halve breedte van 6.25∘6.25{}^{\circ}, en één standaarddeviatie van 6.25/1.96=3.19∘6.25/1.96 = 3.19{}^{\circ}, oftewel 355 km355\,\mathrm{km}. Sinds formatie A is Kesta naar het noorden gedreven met (2918±355) km(2918 \pm 355)\,\mathrm{km}.

এই অধ্যায়ে সংজ্ঞায়িত পরিভাষা-শব্দ

পরিভাষায় সব 957 টি শব্দ দেখো