Chemistry Labs

Chimie nucléaire et radiochimie

Datation radiométrique

Comment les constantes de décroissance et les rapports isotopes père-fils fournissent des horloges géochronologiques, du radiocarbone aux isochrones et à l'âge de la Terre.

IntuitionIntuition : des sabliers radioactifs scellés dans la roche

Un sablier mesure le temps grâce au flux régulier de grains de sable s'écoulant du compartiment supérieur vers le compartiment inférieur. Dans la nature, les noyaux radioactifs agissent comme des sabliers atomiques : les atomes pères instables se transforment en atomes fils radiogéniques selon des lois quantiques invariables, insensibles à la température, à la pression ou aux liaisons chimiques. Lorsqu'un magma cristallise ou qu'un organisme cesse sa respiration métabolique, le sablier se retourne : les atomes pères sont scellés dans le réseau cristallin, et l'accumulation continue des descendants enregistre le temps profond.

Sélectionnez un géochronomètre pour suivre sa trajectoire de décroissance : la fraction restante de 14C sur 50 ka (jalons à 5,73, 11,5 et 50 ka), l'accumulation radiogénique 40Ar/40K jusqu'à l'âge de 4,54 Ga de la Terre, et le rapport 207Pb/206Pb* traçant des milliards d'années d'évolution planétaire.

ScolaireNiveau scolaire : équation fondamentale de l'âge et datation au radiocarbone

Définition: L'équation fondamentale de l'âge et le système fermé

Pour un isotope père radioactif PP se désintégrant en descendant radiogénique D∗D^* avec la constante λ\lambda, l'accumulation de descendants s'écrit D∗=P(eλt−1)D^* = P(e^{\lambda t} - 1). En isolant l'âge de cristallisation tt, on obtient t=1λln⁡(1+D∗P)t = \frac{1}{\lambda}\ln\left(1 + \frac{D^*}{P}\right). Cette relation ne reflète l'âge géologique réel que si le minéral est resté un système fermé strict (sans perte ni gain d'atomes pères ou fils) depuis son refroidissement sous sa température de fermeture (TcT_c).

La datation au radiocarbone repose sur la spallation cosmique dans la haute atmosphère, où des neutrons thermiques réagissent avec l'azote : X714X27214N+X01X2021n→X614X26214C+X11X2121p\ce{^{14}_7N + ^1_0n -> ^{14}_6C + ^1_1p}. Le carbone-14 formé est oxydé en X14X2214COX2\ce{^{14}CO2} et assimilé dans la biosphère par photosynthèse, maintenant un rapport stationnaire X14X2214C/X12X2212C≈1,2×10−12\ce{^{14}C/^{12}C} \approx 1{,}2 \times 10^{-12} chez les êtres vivants. À la mort, les échanges cessent et X14X2214C\ce{^{14}C} se désintègre en azote-14 par émission β−\beta^- avec une demi-vie de t1/2=5730±40 anst_{1/2} = 5730\pm 40\text{ ans} (valeur de Cambridge). Par convention internationale, les âges bruts sont exprimés en 'années radiocarbone BP' (référence 1950) avec la demi-vie de Libby (5568 ans5568\text{ ans}). Les variations du champ géomagnétique et de l'activité solaire imposent des courbes d'étalonnage comme IntCal20, basées sur la dendrochronologie. La spectrométrie de masse par accélérateur (SMA) compte directement les atomes sur des échantillons milligrammiques jusqu'à ∼50000 ans\sim 50000\text{ ans}.

Principaux systèmes de datation géochronologique et minéraux cibles
SystèmeMode de désintégrationDemi-viePlage de datationMatériaux et minéraux types
X14X2214C→X14X2214N\ce{^{14}C -> ^{14}N}Désintégration bêta-5730 ans300 à 50 000 ansCharbon de bois, os, coquilles, tourbe
X40X2240K→X40X2240Ar\ce{^{40}K -> ^{40}Ar}Capture électronique (10,7%)1,248 Ga100 ka à 4,5 GaSanidine, biotite, muscovite, hornblende
X87X2287Rb→X87X2287Sr\ce{^{87}Rb->^{87}Sr}Désintégration bêta-48,8 Ga10 Ma à 4,6 GaMicas, feldspaths potassiques, roches totales
X238X22238U→X206X22206Pb\ce{^{238}U->^{206}Pb}Chaîne alpha et bêta4,468 Ga1 Ma à 4,6 GaZircon, baddeleyite, monazite, titanite
X235X22235U→X207X22207Pb\ce{^{235}U->^{207}Pb}Chaîne alpha et bêta0,7038 Ga1 Ma à 4,6 GaCouplé à 238U sur diagramme concordia
X147X22147Sm→X143X22143Nd\ce{^{147}Sm->^{143}Nd}Désintégration alpha106 Ga100 Ma à 4,6 GaBasaltes, anorthosites lunaires, météorites

Exemple: Calcul de l'âge radiocarbone non étalonné

Une poutre de cèdre antique découverte lors d'une fouille archéologique présente une activité spécifique en X14X2214C\ce{^{14}C} égale à exactement 12,5%12{,}5\% (1/81/8) de celle d'un étalon de bois moderne préindustriel. En utilisant la demi-vie physique de Cambridge t1/2=5730 anst_{1/2} = 5730\text{ ans}, calculez l'âge radiométrique brut de cet artefact avant étalonnage dendrochronologique.

Solution

La fraction d'activité restante est N/N0=0,125=(1/2)3N/N_0 = 0{,}125 = (1/2)^3. Cela correspond à exactement trois demi-vies écoulées : t=3×t1/2=3×5730=17190 anst = 3 \times t_{1/2} = 3 \times 5730 = 17190\text{ ans}. Selon la loi exponentielle : t=−t1/2ln⁡2ln⁡(NN0)=−57300,69315×ln⁡(0,125)≈17190 anst = -\frac{t_{1/2}}{\ln 2}\ln\left(\frac{N}{N_0}\right) = -\frac{5730}{0{,}69315}\times \ln(0{,}125) \approx 17190\text{ ans}. Les archéologues calibrent ensuite cette date brute sur la courbe IntCal20 pour obtenir l'âge calendaire réel en années cal BP.

UniversitaireNiveau universitaire : décroissance branchée K-Ar, 40Ar/39Ar et isochrones Rb-Sr

t=1λln⁡(1+λλe40Ar∗40K)withλ=5.543×10−10 y−1, λe=0.580×10−10 y−1t = \frac{1}{\lambda}\ln\left(1 + \frac{\lambda}{\lambda_e}\frac{^{40}\text{Ar}^*}{^{40}\text{K}}\right) \quad \text{with} \quad \lambda = 5.543 \times 10^{-10}\ \text{y}^{-1},\ \lambda_e = 0.580 \times 10^{-10}\ \text{y}^{-1}

Le potassium-40 présente une décroissance branchée : capture électronique vers le gaz rare X40X2240Ar\ce{^{40}Ar} (10,72%10{,}72\%, λe=0,580×10−10 an−1\lambda_e = 0{,}580 \times 10^{-10}\ \text{an}^{-1}) et désintégration β−\beta^- vers X40X2240Ca\ce{^{40}Ca} (89,28%89{,}28\%, λβ=4,963×10−10 an−1\lambda_\beta = 4{,}963 \times 10^{-10}\ \text{an}^{-1}), avec une constante totale λ=5,543×10−10 an−1\lambda = 5{,}543 \times 10^{-10}\ \text{an}^{-1} (t1/2=1,248 Gat_{1/2} = 1{,}248\text{ Ga}). Le magma liquide dégazant totalement l'argon, le chronomètre est remis à zéro (X40X2240Ar0≈0\ce{^{40}Ar}_0 \approx 0). Après cristallisation, X40X2240Ar∗\ce{^{40}Ar}^* s'accumule dans le réseau des sanidines ou biotites. La méthode moderne X40X2240Ar/X39X2239Ar\ce{^{40}Ar/^{39}Ar} irradie l'échantillon par neutrons en réacteur pour convertir X39X2239K(n,p)X39X2239Ar\ce{^{39}K}(n,p)\ce{^{39}Ar} ; un chauffage par paliers laser en spectrométrie de masse fournit un plateau d'âge distinguant les âges primaires des fuites thermiques d'argon.

87Sr86Sr=(87Sr86Sr)0+87Rb86Sr(eλt−1)\frac{^{87}\text{Sr}}{^{86}\text{Sr}} = \left(\frac{^{87}\text{Sr}}{^{86}\text{Sr}}\right)_0 + \frac{^{87}\text{Rb}}{^{86}\text{Sr}}\left(e^{\lambda t} - 1\right)

Lorsque des descendants initiaux D0D_0 sont présents dès la cristallisation, la datation sur échantillon unique est indéterminée. La méthode isochrone résout ce défi en normalisant chaque abondance par un isotope stable non radiogénique du fils, tel que X86X2286Sr\ce{^{86}Sr} pour le couple rubidium-strontium (X87X2287Rb→X87X2287Sr\ce{^{87}Rb -> ^{87}Sr}, t1/2=48,8 Gat_{1/2} = 48{,}8\text{ Ga}, λ=1,42×10−11 an−1\lambda = 1{,}42 \times 10^{-11}\ \text{an}^{-1}). Les différents minéraux cristallisant à partir du même magma présentent des rapports élémentaires Rb/Sr\ce{Rb/Sr} contrastés mais un rapport isotopique initial (X87X2287Sr/X86X2286Sr)0(\ce{^{87}Sr/^{86}Sr})_0 strictement identique. Le diagramme reportant X87X2287Sr/X86X2286Sr\ce{^{87}Sr/^{86}Sr} en fonction de X87X2287Rb/X86X2286Sr\ce{^{87}Rb/^{86}Sr} donne une droite isochrone dont la pente m=eλt−1≈λtm = e^{\lambda t} - 1 \approx \lambda t fournit directement l'âge de cristallisation tt.

Exemple: Datation potassium-argon d'un cristal de sanidine volcanique

L'analyse par spectrométrie de masse d'un monocristal de sanidine extrait d'un tuf volcanique donne un rapport atomique d'argon radiogénique au potassium de X40X2240Ar∗X40X2240K=0,0950\frac{\ce{^{40}Ar}^*}{\ce{^{40}K}} = 0{,}0950. Avec les constantes λ=5,543×10−10 an−1\lambda = 5{,}543 \times 10^{-10}\ \text{an}^{-1} et λe=0,580×10−10 an−1\lambda_e = 0{,}580 \times 10^{-10}\ \text{an}^{-1}, calculez l'âge de cristallisation de cette formation volcanique en millions d'années (Ma).

Solution

En appliquant la formule branchée K-Ar : t=1λln⁡(1+λλeX40X2240Ar∗X40X2240K)t = \frac{1}{\lambda}\ln\left(1 + \frac{\lambda}{\lambda_e}\frac{\ce{^{40}Ar}^*}{\ce{^{40}K}}\right). Le facteur de branchement vaut λλe=5,543×10−100,580×10−10≈9,5569\frac{\lambda}{\lambda_e} = \frac{5{,}543 \times 10^{-10}}{0{,}580 \times 10^{-10}} \approx 9{,}5569. L'argument du logarithme vaut 1+9,5569×0,0950=1,90791 + 9{,}5569 \times 0{,}0950 = 1{,}9079. Son logarithme naturel est ln⁡(1,9079)≈0,6460\ln(1{,}9079) \approx 0{,}6460. En divisant par la constante totale λ\lambda : t=0,64605,543×10−10 an−1≈1,1654×109 ans=1165,4 Mat = \frac{0{,}6460}{5{,}543 \times 10^{-10}\ \text{an}^{-1}} \approx 1{,}1654 \times 10^9\text{ ans} = 1165{,}4\text{ Ma} (1,165 Ga1{,}165\text{ Ga}).

AvancéNiveau avancé : diagramme concordia U-Pb, zircons de Jack Hills et âge de la Terre

Le zircon (ZrSiOX4\ce{ZrSiO4}) est la référence absolue en géochronologie. Lors de la cristallisation, le site du zirconium (r=84 pmr = 84\text{ pm}) accueille UX4+\ce{U^4+} (r=100 pmr = 100\text{ pm}) mais exclut totalement l'ion PbX2+\ce{Pb^2+} (r=119 pmr = 119\text{ pm}) avec un facteur d'incompatibilité supérieur à 10510^5, garantissant une teneur initiale en plomb radiogénique quasi nulle (Pb0≈0\ce{Pb}_0 \approx 0). Comme le zircon abrite deux horloges isotopiques couplées de demi-vies contrastées—X238X22238U→X206X22206Pb\ce{^{238}U -> ^{206}Pb} (t1/2=4,468 Gat_{1/2} = 4{,}468\text{ Ga}, λ238=1,55125×10−10 an−1\lambda_{238} = 1{,}55125 \times 10^{-10}\ \text{an}^{-1}) et X235X22235U→X207X22207Pb\ce{^{235}U -> ^{207}Pb} (t1/2=0,7038 Gat_{1/2} = 0{,}7038\text{ Ga}, λ235=9,8485×10−10 an−1\lambda_{235} = 9{,}8485 \times 10^{-10}\ \text{an}^{-1})—le diagramme reportant X206X22206Pb∗/X238X22238U\ce{^{206}Pb}^*/\ce{^{238}U} contre X207X22207Pb∗/X235X22235U\ce{^{207}Pb}^*/\ce{^{235}U} engendre la courbe concordia de Wetherill. Les minéraux fermés se placent sur la courbe ; ceux ayant subi une fuite de plomb définissent une droite discordia dont l'intersection supérieure date la cristallisation et l'inférieure l'événement thermique. Des zircons détritiques de Jack Hills en Australie occidentale atteignent 4,404 Ga4{,}404\text{ Ga}, témoins de la croûte à l'Hadéen.

Prendre le rapport des deux équations radiogéniques élimine la mesure de la concentration absolue en uranium : X207X22207Pb∗X206X22206Pb∗=X235X22235UX238X22238U(eλ235t−1eλ238t−1)=1137,88(eλ235t−1eλ238t−1)\frac{\ce{^{207}Pb}^*}{\ce{^{206}Pb}^*} = \frac{\ce{^{235}U}}{\ce{^{238}U}}\left(\frac{e^{\lambda_{235} t} - 1}{e^{\lambda_{238} t} - 1}\right) = \frac{1}{137{,}88}\left(\frac{e^{\lambda_{235} t} - 1}{e^{\lambda_{238} t} - 1}\right), où 137,88 est le rapport X238X22238U/X235X22235U\ce{^{238}U/^{235}U} actuel. En 1956, Clair Patterson a mesuré la composition isotopique du plomb primordial dans la troïlite (FeS\ce{FeS}) sans uranium de la météorite ferreuse de Canyon Diablo. En reliant ces mesures au plomb des sédiments océaniques terrestres sur une isochrone Pb-Pb commune, Patterson a établi l'âge de la Terre à 4,55±0,07 Ga4{,}55 \pm 0{,}07\text{ Ga} (4,54 Ga4{,}54\text{ Ga}). L'analyse U-Pb de pointe des inclusions riches en calcium et aluminium (CAI) des chondrites date la condensation de la nébuleuse solaire à 4,567 Ga4{,}567\text{ Ga}.

Références

  • Atmospheric Helium Three and Radiocarbon from Cosmic Radiation · W. F. Libby, 1946
  • Precision Measurement of Half-Lives and Specific Activities of 235U and 238U · A. H. Jaffey, K. F. Flynn, L. E. Glendenin, W. C. Bentley, A. M. Essling, 1971
  • The IntCal20 Northern Hemisphere Radiocarbon Age Calibration Curve (0-55 cal kBP) · P. J. Reimer et al., 2020