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

Ciencias de la Tierra y del espacio universitarias — tercer año

书籍暂未推出简体中文版——当前显示 Español 版。

Ciencias de la Tierra y del espacio universitarias — tercer año · 本科三年级

13Geocronología avanzada

Un cristal de circón de una décima de milímetro de largo ha mantenido en marcha dos relojes desde que cristalizó su granito: el uranio-238 que se desintegra en plomo-206 y el uranio-235 que se desintegra en plomo-207. Cuando los dos relojes coinciden, la edad es indudable. Cuando discrepan, la propia discrepancia puede datarse: el cristal perdió parte de su plomo en un evento posterior, y los dos relojes, leídos en conjunto, dan la edad del granito y la del evento. Este capítulo lee varios relojes a la vez, en uranio y plomo, en argón, en neodimio, en trazas de fisión y helio, y en los átomos raros producidos por los rayos cósmicos, y indica hasta qué punto cabe fiarse de cada edad.

Lo que ya sabes

La desintegración radiactiva, la ecuación de la edad, la isócrona y su razón inicial, la recta de mínimos cuadrados, la desintegración ramificada, la edad potasio–argón, la temperatura de cierre y la edad de enfriamiento están en el volumen de primer año, en su capítulo sobre datación radiométrica; la fórmula de erosión estacionaria de los nucleidos cosmogénicos está en su capítulo sobre erosión y transporte. El reservorio condrítico uniforme, la notación épsilon y el manto empobrecido están en el volumen de segundo año, en su capítulo sobre geoquímica isotópica; la ionización térmica y la espectrometría de masas multicolector, en su capítulo sobre métodos de campo y de laboratorio. La termocronología y la relación edad–altitud están en el Capítulo 10.

13.1 Uranio, plomo y la concordia

El uranio tiene dos isótopos de vida larga. El uranio-238 se desintegra a través de una cadena de hijas de vida corta hasta el plomo-206, con λ238=1.551 25×10−10 yr−1\lambda_{238} = 1.551\,25 \times 10^{-10}\,\mathrm{yr}^{-1}; el uranio-235 se desintegra en plomo-207 con λ235=9.8485×10−10 yr−1\lambda_{235} = 9.8485 \times 10^{-10}\,\mathrm{yr}^{-1}, seis veces más rápido. Su razón actual es 235U/238U=0.0072567^{235}\mathrm{U}/^{238}\mathrm{U} = 0.0072567, o bien 238U/235U=137.8^{238}\mathrm{U}/^{235}\mathrm{U} = 137.8, la misma en toda roca terrestre. El circón, ZrSiO4\mathrm{ZrSiO_4}, incorpora mucho uranio y casi nada de plomo al crecer, y resiste la meteorización y sobrevive al metamorfismo: casi todo su plomo es radiogénico.

Se definen x=207Pb∗/235Ux = {}^{207}\mathrm{Pb}^*/{}^{235}\mathrm{U} e y=206Pb∗/238Uy = {}^{206}\mathrm{Pb}^*/{}^{238}\mathrm{U} como razones actuales del plomo radiogénico respecto a su progenitor. Cada una da una edad, t206=ln⁡(1+y)/λ238t_{206} = \ln(1 + y)/\lambda_{238} y t207=ln⁡(1+x)/λ235t_{207} = \ln(1 + x)/\lambda_{235}.

Definición 13.1 (Concordia, discordia, intersecciones, diagrama de Tera–Wasserburg)

La concordia es la curva del plano (x,y)(x, y) en la que las dos edades uranio–plomo son iguales. Una discordia es una recta a lo largo de la cual se alinean los análisis discordantes de una misma población de cristales; sus intersecciones con la concordia son la intersección superior (la más antigua) y la intersección inferior (la más joven). El diagrama de Tera–Wasserburg representa X=238U/206Pb∗X = {}^{238}\mathrm{U}/{}^{206}\mathrm{Pb}^* frente a Y=207Pb∗/206Pb∗Y = {}^{207}\mathrm{Pb}^*/{}^{206}\mathrm{Pb}^*.

Teorema 13.2 (La concordia)

Un cristal que ha permanecido cerrado al uranio y al plomo desde que se formó, hace un tiempo tt, se sitúa en

C(t)=(eλ235t−1,  eλ238t−1).C(t) = \bigl(e^{\lambda_{235}t} - 1,\; e^{\lambda_{238}t} - 1\bigr).

A medida que tt crece, la curva sube cada vez más despacio frente a xx: su pendiente, dy/dx=(λ238/λ235) e(λ238−λ235)t\mathrm dy/\mathrm dx = (\lambda_{238}/\lambda_{235})\,e^{(\lambda_{238} - \lambda_{235})t}, decrece desde λ238/λ235=0.1575\lambda_{238}/\lambda_{235} = 0.1575 en el origen, de modo que la concordia es cóncava hacia el eje xx.

Demostración. Para cada progenitor, el número de hijas radiogénicas por átomo progenitor que queda es eλt−1e^{\lambda t} - 1, la ecuación de la edad del volumen de primer año. Derivando ambas coordenadas respecto de tt y dividiendo se obtiene la pendiente; dado que λ235>λ238\lambda_{235} > \lambda_{238}, decrece con tt. ∎

El plomo es un ion grande que encaja mal en la red del circón. En las zonas de un cristal dañadas por su propia radiactividad, puede perderse en un evento posterior de calentamiento o circulación de fluidos.

Teorema 13.3 (La pérdida episódica de plomo da una cuerda)

Un cristal que se formó hace un tiempo t1t_1 y que perdió, hace un tiempo t2t_2, una fracción ff de su plomo radiogénico, y que por lo demás permaneció cerrado, se sitúa hoy en

(1−f) C(t1)+f C(t2),(1 - f)\,C(t_1) + f\,C(t_2),

sobre la cuerda que une los puntos t1t_1 y t2t_2 de la concordia. Los cristales que perdieron fracciones distintas se alinean sobre esa cuerda, una discordia cuya intersección superior es t1t_1 y la inferior es t2t_2.

Demostración. Se considera un sistema de desintegración, con constante λ\lambda, y se cuenta por átomo de progenitor presente hoy. Hace un tiempo t2t_2 había eλt2e^{\lambda t_2} átomos progenitores, y hijas radiogénicas eλt2(eλ(t1−t2)−1)=eλt1−eλt2e^{\lambda t_2}(e^{\lambda(t_1 - t_2)} - 1) = e^{\lambda t_1} - e^{\lambda t_2}. El evento dejó (1−f)(eλt1−eλt2)(1 - f)(e^{\lambda t_1} - e^{\lambda t_2}); desde entonces el progenitor generó eλt2−1e^{\lambda t_2} - 1 más. La suma es (1−f)(eλt1−1)+f(eλt2−1)(1 - f)(e^{\lambda t_1} - 1) + f(e^{\lambda t_2} - 1). Los dos isótopos de plomo se comportan igual químicamente, así que ff es la misma para ambos, y las dos coordenadas se combinan en el punto enunciado. ∎

(a) La concordia, con sus edades en millones de años, y cinco circones discordantes inventados de un mismo granito (en rojo; sus elipses de error son más pequeñas que los símbolos), alineados sobre una discordia (discontinua) que corta a la concordia cerca de 311 y 100\, Ma. (b) Los mismos circones y la misma recta en el diagrama de Tera–Wasserburg, donde una discordia sigue siendo una recta.
(a) La concordia, con sus edades en millones de años, y cinco circones discordantes inventados de un mismo granito (en rojo; sus elipses de error son más pequeñas que los símbolos), alineados sobre una discordia (discontinua) que corta a la concordia cerca de 311 y 100 Ma100\,\mathrm{Ma}. (b) Los mismos circones y la misma recta en el diagrama de Tera–Wasserburg, donde una discordia sigue siendo una recta.

Las intersecciones se hallan donde a+b x(t)=y(t)a + b\,x(t) = y(t), una ecuación en tt que se resuelve por el método de Newton: con g(t)=eλ238t−1−a−b(eλ235t−1)g(t) = e^{\lambda_{238}t} - 1 - a - b(e^{\lambda_{235}t} - 1) y g′(t)=λ238eλ238t−bλ235eλ235tg'(t) = \lambda_{238}e^{\lambda_{238}t} - b\lambda_{235}e^{\lambda_{235}t}, iterando t←t−g/g′t \leftarrow t - g/g' a partir de una estimación cercana a cada intersección. Tres propiedades del resultado importan en la práctica. En primer lugar, la intersección superior solo es precisa si la discordia corta a la concordia con un ángulo grande. Una discordia joven, como en la figura anterior, discurre casi a lo largo de la concordia, y un pequeño cambio en su pendiente desplaza la intersección mucho. En segundo lugar, si la pérdida de plomo es reciente, t2=0t_2 = 0, la discordia pasa por el origen, y todos los cristales tienen la misma razón 207Pb∗/206Pb∗^{207}\mathrm{Pb}^*/^{206}\mathrm{Pb}^*,

207Pb∗206Pb∗=1137.8 eλ235t1−1eλ238t1−1,\frac{^{207}\mathrm{Pb}^*}{^{206}\mathrm{Pb}^*} = \frac{1}{137.8}\, \frac{e^{\lambda_{235}t_1} - 1}{e^{\lambda_{238}t_1} - 1},

lo que da t1t_1 sin necesidad del uranio: la edad plomo–plomo. En tercer lugar, una recta en el diagrama de la concordia, x=(y−a)/bx = (y - a)/b, se convierte en las coordenadas de Tera–Wasserburg Y=(1−aX)/(137.8 b)Y = (1 - aX)/(137.8\,b), de nuevo una recta. El diagrama de Tera–Wasserburg usa razones medidas directamente, sin necesidad de calcular el uranio-235 a partir del uranio-238, y una mezcla con plomo común, no radiogénico y rico en plomo-207, se representa sobre una recta que corta el eje YY en la razón del plomo común.

Las edades uranio–plomo se miden de dos maneras. Cristales sueltos o fragmentos, disueltos tras eliminar las zonas dañadas mediante disolución parcial y separados su plomo y su uranio, se miden por ionización térmica en espectrometría de masas, con las edades más precisas de la geología. Puntos de unas pocas decenas de micrómetros de diámetro se miden in situ, con una microsonda iónica o por ablación láser, menos precisamente, pero zona a zona a través de un cristal que creció en varios episodios.

En el laboratorio — El laboratorio de geocronología

Un espectrómetro de masas en un laboratorio limpio: la torre de muestras a la izquierda, el tubo de vuelo y el electroimán (ilustración generada por IA).

Una datación de circones empieza con una roca triturada y termina en un laboratorio limpio. Los geocronólogos separan los circones por densidad y magnetismo, escogen los cristales más limpios al microscopio binocular y obtienen imágenes de su zonado. Los cristales se disuelven en parte para eliminar sus zonas dañadas, y luego se disuelven por completo con un trazador de isótopos conocidos, y el uranio y el plomo se separan en diminutas columnas de resina, en aire filtrado para que no entre plomo del polvo. Un espectrómetro de masas mide las razones isotópicas. Las muestras de argón siguen un camino distinto: un reactor, después una línea de gas en la que un láser calienta cada paso y el argón se depura antes del espectrómetro.

13.2 Datación argón–argón

El método potasio–argón del volumen de primer año mide el potasio y el argón en dos alícuotas separadas de un mineral. El método argón–argón mide solo argón, en una alícuota, como razones de isótopos. La muestra se irradia con neutrones rápidos en el núcleo de un reactor nuclear, donde la reacción 39K(n,p)39Ar^{39}\mathrm{K}(n,p)^{39}\mathrm{Ar} convierte una fracción fija de su potasio-39 en argón-39, un indicador del potasio. Un mineral de edad conocida, el monitor, se irradia junto a ella.

Definición 13.4 (Datación argón–argón, factor J, espectro de edades, edad meseta)

La datación 40Ar/39Ar^{40}\mathrm{Ar}/^{39}\mathrm{Ar} data un mineral portador de potasio mediante la razón R=40Ar∗/39ArKR = {}^{40}\mathrm{Ar}^*/{}^{39}\mathrm{Ar}_{\mathrm K} de su argón-40 radiogénico al argón-39 formado a partir de su potasio por irradiación neutrónica. El factor J de una irradiación es J=(eλtm−1)/RmJ = (e^{\lambda t_m} - 1)/R_m, hallado a partir de un monitor de edad conocida tmt_m y de la razón medida RmR_m, siendo λ\lambda la constante de desintegración total del potasio-40. Cuando la muestra se desgasifica en etapas de temperatura creciente, su espectro de edades representa la edad de cada etapa frente a la fracción acumulada de argón-39 liberado. Una edad meseta es la media ponderada de las edades de etapas consecutivas que concuerdan dentro de sus incertidumbres y liberan una gran parte del argón-39.

Teorema 13.5 (La edad argón–argón)

Una muestra irradiada junto a un monitor tiene la edad

t=1λln⁡(1+JR).t = \frac1\lambda\ln(1 + JR).

Demostración. El argón radiogénico es 40Ar∗=(λe/λ) 40K (eλt−1)^{40}\mathrm{Ar}^* = (\lambda_e/\lambda)\, {}^{40}\mathrm{K}\,(e^{\lambda t} - 1), por la desintegración ramificada del volumen de primer año. El argón-39 formado es 39ArK=39K Φ^{39}\mathrm{Ar}_{\mathrm K} = {}^{39}\mathrm{K}\,\Phi, donde Φ\Phi recoge la dosis neutrónica y la sección eficaz de la reacción. Así, R=(eλt−1)/JR = (e^{\lambda t} - 1)/J con J=λ Φ 39K/(λe 40K)J = \lambda\,\Phi\,{}^{39}\mathrm{K}/(\lambda_e\,{}^{40}\mathrm{K}). La razón 40K/39K^{40}\mathrm{K}/^{39}\mathrm{K} es la misma en todo el potasio, y Φ\Phi la misma para dos muestras irradiadas lado a lado, así que JJ es el mismo para la muestra y el monitor; la edad y la razón del monitor lo determinan, e invertir RR da tt. ∎

Las constantes son las del método potasio–argón, λ=5.543×10−10 yr−1\lambda = 5.543 \times 10^{-10}\,\mathrm{yr}^{-1}. Un monitor habitual es una sanidina de una toba volcánica datada en 28.201 Ma28.201\,\mathrm{Ma}; una razón del monitor Rm=1.3129R_m = 1.3129 da J=0.01200J = 0.01200. El argón atmosférico, con 40Ar/36Ar=295.5^{40}\mathrm{Ar}/^{36}\mathrm{Ar} = 295.5, se resta con el argón-36 medido en cada etapa, y el argón formado por neutrones a partir del calcio y del propio potasio se corrige con sales irradiadas en el mismo paquete.

El calentamiento por etapas es la gran fortaleza del método. La muestra se calienta en un horno o con un láser en etapas de temperatura creciente, y cada etapa proporciona una edad. Si el mineral permaneció cerrado desde que se enfrió, todas las etapas coinciden: el espectro es plano. Si perdió argón de los bordes de sus granos, las primeras etapas, que proceden de los bordes, son más jóvenes, y las edades suben en forma de escalera hacia una meseta dada por los núcleos. Los criterios recomendados para una meseta son al menos cinco etapas consecutivas que acumulan más de 50 %50\,\% del argón-39, sin tendencia entre ellas; las etapas de baja temperatura de una escalera suelen interpretarse como pérdida de argón radiogénico. El argón en exceso, atrapado a partir de un fluido cuando el mineral creció, produce el efecto contrario: hace que algunas etapas sean demasiado viejas, a menudo la primera y la última, en un espectro con forma de silla de montar. La suma de todas las etapas, la edad del gas total, es lo que habría dado una edad potasio–argón: un promedio que oculta la historia.

El espectro de edades inventado de una biotita del granito de la primera figura: cada etapa se dibuja a su edad, y su anchura es la fracción de argón-39 que liberó. (a) Tres etapas de baja temperatura suben en escalera hasta una meseta. (b) Las siete etapas de meseta, que concentran 82\,\% del argón-39, con sus intervalos de dos desviaciones típicas (sombreados) y la edad de meseta (en rojo).
El espectro de edades inventado de una biotita del granito de la primera figura: cada etapa se dibuja a su edad, y su anchura es la fracción de argón-39 que liberó. (a) Tres etapas de baja temperatura suben en escalera hasta una meseta. (b) Las siete etapas de meseta, que concentran 82 %82\,\% del argón-39, con sus intervalos de dos desviaciones típicas (sombreados) y la edad de meseta (en rojo).

13.3 Edades modelo del neodimio

Cuando un fundido se extrae del manto para formar corteza continental, el samario, más compatible, queda en el residuo en mayor proporción que el neodimio: la nueva corteza tiene una razón 147Sm/144Nd^{147}\mathrm{Sm}/^{144}\mathrm{Nd} más baja que su fuente, de unos 0.13 en la mayoría de las rocas corticales, frente a 0.1966 del reservorio condrítico. Su 143Nd/144Nd^{143}\mathrm{Nd}/^{144}\mathrm{Nd} crece entonces más lentamente. Una vez formada, una roca cortical tiende a conservar su razón samario–neodimio, a través de la erosión, la sedimentación e incluso la refusión, de modo que su neodimio recuerda cuándo salió del manto.

Definición 13.6 (Edad modelo)

La edad modelo de una roca en un sistema isotópico es el tiempo en el que su razón, proyectada hacia atrás con la razón padre–hija actual, iguala la de un reservorio modelo. Para el neodimio, TCHURT_{\mathrm{CHUR}} emplea el reservorio condrítico uniforme y TDMT_{\mathrm{DM}} el manto empobrecido.

Teorema 13.7 (Edad modelo)

Una muestra con razones actuales Is=143Nd/144NdI_s = {}^{143}\mathrm{Nd}/^{144}\mathrm{Nd} y Rs=147Sm/144NdR_s = {}^{147}\mathrm{Sm}/^{144}\mathrm{Nd}, comparada con un reservorio de razones actuales ImI_m y RmR_m, tiene la edad modelo

T=1λln⁡(1+Is−ImRs−Rm),λ=λ147.T = \frac1\lambda\ln\left(1 + \frac{I_s - I_m}{R_s - R_m}\right), \qquad \lambda = \lambda_{147}.

Demostración. Un sistema cerrado tenía, hace un tiempo TT, la razón I−R(eλT−1)I - R(e^{\lambda T} - 1): la ecuación de la isócrona leída hacia atrás. Igualando esto para la muestra y el reservorio da Is−Im=(Rs−Rm)(eλT−1)I_s - I_m = (R_s - R_m)(e^{\lambda T} - 1). ∎

Con λ147=6.524×10−12 yr−1\lambda_{147} = 6.524 \times 10^{-12}\,\mathrm{yr}^{-1}, los reservorios son el condrítico, 0.512638 y 0.1966 en la actualidad, y el manto empobrecido, 0.51315 en la actualidad, en torno a εNd=+10\varepsilon_{\mathrm{Nd}} = +10, valor que una evolución cerrada a partir de valores condríticos desde hace 4.55 Ga4.55\,\mathrm{Ga} da la razón 147Sm/144Nd=0.2137^{147}\mathrm{Sm}/^{144}\mathrm{Nd} = 0.2137. Un granito inventado con Is=0.511900I_s = 0.511900 y Rs=0.1100R_s = 0.1100 tiene εNd=−14.4\varepsilon_{\mathrm{Nd}} = -14.4 en la actualidad, TCHUR=1.30 GaT_{\mathrm{CHUR}} = 1.30\,\mathrm{Ga} y TDM=1.84 GaT_{\mathrm{DM}} = 1.84\,\mathrm{Ga}. La figura siguiente muestra ambas edades. La edad del manto empobrecido es la mejor estimación, porque la corteza se extrae de un manto ya empobrecido; es la edad media a la que el neodimio de la muestra salió del manto, su edad de residencia cortical. Dos salvedades la limitan: la corteza no debe haber cambiado su razón samario–neodimio desde entonces, y una roca formada por mezcla de corteza antigua con fundidos de manto nuevo tiene una edad modelo intermedia carente de significado.

_ Nd de un granito inventado proyectado hacia atrás en el tiempo con su propia razón samario–neodimio (rojo). Alcanza el reservorio condrítico en T_ CHUR = 1.30\, Ga y el manto empobrecido (azul) en T_ DM = 1.84\, Ga.
εNd\varepsilon_{\mathrm{Nd}} de un granito inventado proyectado hacia atrás en el tiempo con su propia razón samario–neodimio (rojo). Alcanza el reservorio condrítico en TCHUR=1.30 GaT_{\mathrm{CHUR}} = 1.30\,\mathrm{Ga} y el manto empobrecido (azul) en TDM=1.84 GaT_{\mathrm{DM}} = 1.84\,\mathrm{Ga}.

13.4 Trazas de fisión y helio

El uranio-238 también se desintegra, muy raramente, por fisión espontánea: el núcleo se divide en dos fragmentos que salen despedidos y dejan un rastro de daño en el cristal, de hasta 9 nm9\,\mathrm{nm} de ancho y unos 23 µm23\,\text{µ}\mathrm{m} de largo en apatito. Atacados con ácido, los rastros se hacen visibles al microscopio. Se recuperan al calentarse el cristal, en un intervalo de unos 20 a 350 ∘C350\,{}^{\circ}\mathrm{C} según el mineral; en el apatito se acortan en la zona de recocido parcial de unos 60 a 120 ∘C120\,{}^{\circ}\mathrm{C}, y desaparecen por encima de ella.

Definición 13.8 (Datación por trazas de fisión, recocido de trazas)

La datación por trazas de fisión data un mineral por la densidad de las trazas que deja en él la fisión espontánea del uranio-238, comparada con su contenido en uranio. El recocido de trazas es el acortamiento y la desaparición de las trazas por reparación del cristal cuando se calienta.

Teorema 13.9 (La edad por trazas de fisión)

Con λf\lambda_f la constante de desintegración de la fisión espontánea y λ=λ238\lambda = \lambda_{238} la total, un cristal con NsN_s fisiones espontáneas por unidad de volumen y NN átomos de uranio-238 tiene la edad

t=1λln⁡(1+λλf NsN).t = \frac1\lambda\ln\left(1 + \frac{\lambda}{\lambda_f}\,\frac{N_s}{N}\right).

En el laboratorio, el uranio se mide induciendo la fisión del uranio-235 con neutrones térmicos, y la edad pasa a ser

t=1λln⁡(1+λ ζ g ρd ρsρi),t = \frac1\lambda\ln\left(1 + \lambda\,\zeta\,g\,\rho_d\,\frac{\rho_s}{\rho_i}\right),

con ρs\rho_s y ρi\rho_i las densidades de trazas espontáneas e inducidas, ρd\rho_d la densidad de trazas de un vidrio dosímetro, g=0.5g = 0.5 un factor geométrico y ζ\zeta una constante de calibración obtenida con patrones de edad conocida.

Demostración. El uranio-238 se ha desintegrado N(eλt−1)N(e^{\lambda t} - 1) veces, y la fracción λf/λ\lambda_f/\lambda de estas desintegraciones fueron fisiones: Ns=(λf/λ)N(eλt−1)N_s = (\lambda_f/\lambda)N(e^{\lambda t} - 1), expresión que se invierte en la primera forma. Las densidades de trazas son proporcionales a NsN_s y al número de fisiones inducidas, N σ Φ/137.8N\,\sigma\,\Phi/137.8, donde σ\sigma es la sección eficaz de fisión del uranio-235 y Φ\Phi la fluencia neutrónica, a su vez proporcional a ρd\rho_d. La razón ρs/ρi\rho_s/\rho_i es entonces Ns/NN_s/N por un producto de constantes, y λf\lambda_f, σ\sigma, 137.8137.8 y la eficiencia del ataque químico y del recuento se agrupan con gg en ζ\zeta. ∎

La constante de desintegración por fisión se conoce solo con una incertidumbre de unos 10 %10\,\%, entre 7.9×10−177.9 \times 10^{-17} y 8.7×10−17 yr−18.7 \times 10^{-17}\,\mathrm{yr}^{-1}, por lo que la mayoría de laboratorios prefieren la calibración con ζ\zeta. Las longitudes de las trazas conservan un registro de la trayectoria térmica. Una roca que se enfrió deprisa a través de la zona de recocido parcial tiene trazas largas de longitud casi igual; una que permaneció en ella tiene trazas más cortas de longitudes muy dispersas, y la edad por trazas de fisión es entonces una edad mixta, no el momento de un paso por una única temperatura.

Esquema. (a) Trazas de fisión espontánea, reveladas por ataque químico en un grano de apatito pulido, se cuentan por unidad de área. (b) Una partícula alfa recorre unos 20\, µ m: las emitidas dentro de esa distancia de la superficie (fuera de la línea discontinua) pueden salir del grano, de modo que un grano pequeño retiene menos de su helio.
Esquema. (a) Trazas de fisión espontánea, reveladas por ataque químico en un grano de apatito pulido, se cuentan por unidad de área. (b) Una partícula alfa recorre unos 20 µm20\,\text{µ}\mathrm{m}: las emitidas dentro de esa distancia de la superficie (fuera de la línea discontinua) pueden salir del grano, de modo que un grano pequeño retiene menos de su helio.

Cada desintegración alfa emite un núcleo de helio-4. El uranio, el torio y sus cadenas acumulan helio en un cristal, que lo retiene solo cuando está frío.

Definición 13.10 (Datación (U–Th)/He, zona de retención parcial)

La datación (U–Th)/He data un mineral por el helio-4 acumulado a partir de las desintegraciones alfa del uranio, el torio y sus descendientes, comparado con su contenido en uranio y torio. La zona de retención parcial es el intervalo de temperaturas en el que el helio producido se retiene en parte y en parte se pierde por difusión.

Teorema 13.11 (Crecimiento del helio)

Un cristal cerrado al helio durante un tiempo tt contiene

He=8 238U (eλ238t−1)+7 235U (eλ235t−1)+6 232Th (eλ232t−1),\mathrm{He} = 8\,{}^{238}\mathrm U\,(e^{\lambda_{238}t} - 1) + 7\,{}^{235}\mathrm U\,(e^{\lambda_{235}t} - 1) + 6\,{}^{232}\mathrm{Th}\,(e^{\lambda_{232}t} - 1),

en átomos, con los progenitores medidos en la actualidad; el samario-147 añade un término pequeño.

Demostración. Cada cadena pierde su diferencia de masa en partículas alfa, cuatro unidades de masa cada una: (238−206)/4=8(238 - 206)/4 = 8, (235−207)/4=7(235 - 207)/4 = 7, (232−208)/4=6(232 - 208)/4 = 6, mientras que las desintegraciones beta no cambian el número másico. Cada cadena está en equilibrio secular, de modo que a cada desintegración del progenitor le sigue toda su cadena, y el número de desintegraciones de cada progenitor es el de la ecuación de la edad. ∎

Para un cristal joven, eλt−1≈λte^{\lambda t} - 1 \approx \lambda t y la edad se obtiene con una sola división. Una partícula alfa recorre unos 20 µm20\,\text{µ}\mathrm{m} antes de detenerse, así que un cristal pierde el helio producido cerca de su superficie; la edad medida se divide entre la fracción retenida FTF_T, que es menor en granos pequeños (FT=0.70F_T = 0.70 significa que se ha perdido un 30 %30\,\%). El apatito cierra al helio hacia 50 a 80 ∘C80\,{}^{\circ}\mathrm{C}, el circón a de 160 a 200 ∘C200\,{}^{\circ}\mathrm{C}: con las trazas de fisión del apatito, hacia 90 a 120 ∘C120\,{}^{\circ}\mathrm{C}, los cronómetros de baja temperatura datan los últimos kilómetros del ascenso de una roca hasta la superficie, tal como los empleó el Capítulo 10.

13.5 Nucleidos cosmogénicos

Los rayos cósmicos que chocan contra la atmósfera producen lluvias de partículas secundarias. Los neutrones rápidos que llegan al suelo rompen los núcleos de los minerales de los primeros metros de roca, por espalación: el oxígeno del cuarzo se convierte en berilio-10, y el silicio en aluminio-26. Los muones producen algunos nucleidos a mayor profundidad, y dominan por debajo de unos 3 m3\,\mathrm{m}.

Definición 13.12 (Nucleido cosmogénico, edad de exposición, tasa de producción)

Un nucleido cosmogénico es un nucleido producido en un mineral en la superficie o cerca de ella por las partículas secundarias de los rayos cósmicos. Su tasa de producción PP es el número de átomos producidos por unidad de masa de mineral y por año. La edad de exposición de una superficie es el tiempo transcurrido desde que quedó expuesta por primera vez a los rayos cósmicos, y se obtiene a partir de la concentración de un nucleido cosmogénico.

La tasa de producción disminuye con la profundidad como P0e−ρz/ΛP_0e^{-\rho z/\Lambda}, donde la longitud de atenuación se toma como Λ=160 g cm−2\Lambda = 160\,\mathrm{g}\,\mathrm{cm}^{-2}, unos 60 cm60\,\mathrm{cm} de granito de densidad 2.7 g cm−32.7\,\mathrm{g}\,\mathrm{cm}^{-3}. En la superficie aumenta con la altitud, porque hay menos aire que proteja la roca, y con la latitud geomagnética, porque el campo desvía menos partículas cerca de los polos. Las tasas se dan para una referencia a nivel del mar y latitud alta, y se escalan a cada lugar: una calibración da P0=4.96 aˊtomos g−1 yr−1P_0 = 4.96\,\text{átomos}\,\mathrm{g}^{-1}\,\mathrm{yr}^{-1} de berilio-10 en cuarzo para esa referencia. El berilio-10 tiene una vida media de 1.388 Myr1.388\,\mathrm{Myr}, λ=4.994×10−7 yr−1\lambda = 4.994 \times 10^{-7}\,\mathrm{yr}^{-1}.

Teorema 13.13 (Exposición con desintegración y erosión)

Una superficie que desciende a la tasa constante ε\varepsilon, expuesta durante un tiempo tt sin nucleidos heredados, conserva en su superficie la concentración

N(t)=P0μ(1−e−μt),μ=λ+ρεΛ.N(t) = \frac{P_0}{\mu}\bigl(1 - e^{-\mu t}\bigr), \qquad \mu = \lambda + \frac{\rho\varepsilon}{\Lambda}.

Sin erosión, N=(P0/λ)(1−e−λt)N = (P_0/\lambda)(1 - e^{-\lambda t}); tras un tiempo largo, N=P0/μN = P_0/\mu, y sin desintegración esta es la fórmula P0Λ/(ρε)P_0\Lambda/(\rho\varepsilon) del volumen de primer año.

Demostración. Sígase el grano que hoy está en la superficie. Hace un tiempo ss estaba a la profundidad εs\varepsilon s, donde ganó átomos a la tasa P0e−ρεs/ΛP_0e^{-\rho\varepsilon s/\Lambda}; de esos átomos queda hoy la fracción e−λse^{-\lambda s}. Sumando a lo largo de la exposición,

N=∫0tP0e−ρεs/Λe−λs ds=P0μ(1−e−μt).N = \int_0^t P_0e^{-\rho\varepsilon s/\Lambda}e^{-\lambda s}\,\mathrm ds = \frac{P_0}{\mu}\bigl(1 - e^{-\mu t}\bigr).

Los límites se obtienen haciendo ε=0\varepsilon = 0, dejando que t→∞t \to \infty, y después λ=0\lambda = 0. ∎

Sin erosión la concentración se satura en P0/λP_0/\lambda tras unas pocas semividas: una superficie alcanza el 90 %90\,\% de ese valor tras ln⁡10/λ=4.6 Myr\ln 10/\lambda = 4.6\,\mathrm{Myr}, y a partir de ahí ya no puede datarse. Con erosión, el estado estacionario se alcanza antes, tras unas pocas veces 1/μ1/\mu, y su concentración da ε\varepsilon; la tasa se promedia sobre el tiempo Λ/(ρε)\Lambda/(\rho\varepsilon) que el grano tardó en ascender una longitud de atenuación, unas pocas decenas de miles de años a las tasas de unas pocas decenas de metros por millón de años. Los afloramientos de lecho rocoso se erosionan a unos pocos metros a unas pocas decenas de metros por millón de años; las cuencas fluviales, medidas a partir de la arena de sus ríos, más deprisa.

Dos nucleidos con vidas medias distintas datan el enterramiento. El berilio-10 y el aluminio-26 se producen juntos en el cuarzo en una razón fija; cuando un grano queda enterrado a gran profundidad en una cueva o bajo sedimento grueso, la producción se detiene, y el aluminio-26, que se desintegra más deprisa, se rezaga: la razón entre ambos varía con el tiempo y fecha el enterramiento.

13.6 La estadística de una edad

Una edad lleva asociada su incertidumbre, y un conjunto de edades plantea la cuestión de si concuerdan. Sean nn medidas xix_i de una misma magnitud con desviaciones típicas σi\sigma_i. Su media ponderada es xˉ=∑wixi/W\bar x = \sum w_ix_i/W, con wi=1/σi2w_i = 1/\sigma_i^2 y W=∑wiW = \sum w_i, y su desviación típica es 1/W1/\sqrt W.

Definición 13.14 (Media de los cuadrados de las desviaciones ponderadas)

La media de los cuadrados de las desviaciones ponderadas (MSWD) de nn medidas respecto a un modelo ajustado con pp parámetros es

MSWD=1n−p∑i(xi−x^i)2σi2,\mathrm{MSWD} = \frac{1}{n - p}\sum_i \frac{(x_i - \hat x_i)^2}{\sigma_i^2},

donde x^i\hat x_i es el valor del modelo: p=1p = 1 para una media ponderada, p=2p = 2 para una recta como una isócrona o una discordia.

Proposición 13.15 (El MSWD tiene esperanza igual a uno)

Si los xix_i son independientes y normales, con media μ\mu y las desviaciones típicas σi\sigma_i, el MSWD respecto a su media ponderada tiene esperanza 1.

Demostración. Desarrollando xi−μ=(xi−xˉ)+(xˉ−μ)x_i - \mu = (x_i - \bar x) + (\bar x - \mu) y usando ∑wi(xi−xˉ)=0\sum w_i(x_i - \bar x) = 0,

∑wi(xi−μ)2=∑wi(xi−xˉ)2+W(xˉ−μ)2.\sum w_i(x_i - \mu)^2 = \sum w_i(x_i - \bar x)^2 + W(\bar x - \mu)^2.

A la izquierda, cada término tiene esperanza 1, por lo que la suma tiene esperanza nn. En la derecha, xˉ\bar x tiene varianza 1/W1/W, así que el último término tiene esperanza 1. Por tanto, ∑wi(xi−xˉ)2\sum w_i(x_i - \bar x)^2 tiene esperanza n−1n - 1. ∎

El MSWD se dispersa alrededor de 1 de forma más amplia cuando hay pocos grados de libertad, f=n−pf = n - p: su límite superior al 95 %95\,\% es aproximadamente 1+22/f1 + 2\sqrt{2/f}, 2.63 para f=3f = 3 y 2.15 para f=6f = 6. Por encima del límite, las medidas se dispersan más de lo que permiten sus incertidumbres: las muestras no tienen todas la misma edad, o el sistema fue perturbado. Muy por debajo de 1, las incertidumbres estaban sobreestimadas.

Método 13.16 (Ajuste de una isócrona con errores en ambos ejes)

Cuando ambas coordenadas llevan errores σx,i\sigma_{x,i}, σy,i\sigma_{y,i}, correlacionados con el coeficiente rir_i, como en los datos de uranio–plomo e isócronas:

  1. partir de la pendiente bb de la recta de mínimos cuadrados ordinarios;
  2. asignar a cada punto el peso Wi=1/(σy,i2+b2σx,i2−2briσx,iσy,i)W_i = 1/(\sigma_{y,i}^2 + b^2\sigma_{x,i}^2 - 2br_i\sigma_{x,i}\sigma_{y,i}), la inversa de la varianza de su distancia a la recta medida a lo largo de yy;
  3. ajustar la recta ponderada con estos pesos, lo que da una nueva bb, y volver al paso 2 hasta que bb deje de cambiar;
  4. calcular el MSWD con los pesos WiW_i y n−2n - 2 grados de libertad, y la incertidumbre de la pendiente y de la ordenada en el origen a partir de los pesos finales.

Con σx,i=0\sigma_{x,i} = 0 los pesos no dependen de bb, y el método se reduce a los mínimos cuadrados ponderados de yy sobre xx empleados en el problema.

La incertidumbre de una constante de desintegración es la misma para todas las edades de un sistema: no interviene al comparar edades de un mismo sistema, pero debe sumarse al comparar edades de dos sistemas. La constante de desintegración del samario-147, por ejemplo, se conoce con una precisión de unos 0.4 %0.4\,\%; las del potasio-40, de desintegración ramificada y lenta, limitan la precisión de las edades de argón.

13.7 Ejercicios

Ejercicio 13.1 ★

Un circón tiene 206Pb∗/238U=0.04500^{206}\mathrm{Pb}^*/^{238}\mathrm{U} = 0.04500 y 207Pb∗/235U=0.3300^{207}\mathrm{Pb}^*/^{235}\mathrm{U} = 0.3300. Calcúlense sus dos edades uranio–plomo. ¿Es concordante? ¿A qué lado de la concordia se representa?

Solución

Solución de Ejercicio 13.1.

t206=ln⁡(1.04500)/λ238=283.8 Mat_{206} = \ln(1.04500)/\lambda_{238} = 283.8\,\mathrm{Ma}; t207=ln⁡(1.3300)/λ235=289.6 Mat_{207} = \ln(1.3300)/\lambda_{235} = 289.6\,\mathrm{Ma}. Las edades difieren en 6 Ma6\,\mathrm{Ma}: el circón es discordante. Su t207t_{207} es la mayor, por lo que se representa por debajo de la concordia, del lado de la pérdida de plomo.

Ejercicio 13.2 ★

Un monitor de 28.201 Ma28.201\,\mathrm{Ma} da Rm=1.3129R_m = 1.3129. Calcúlese JJ con λ=5.543×10−10 yr−1\lambda = 5.543 \times 10^{-10}\,\mathrm{yr}^{-1}, y después la edad de una muestra de la misma irradiación con R=10.000R = 10.000 (úsese J=0.01200J = 0.01200).

Solución

Solución de Ejercicio 13.2.

J=(e5.543×10−10×2.8201×107−1)/1.3129=0.015754/1.3129=0.01200J = (e^{5.543\times10^{-10}\times2.8201\times10^7} - 1)/1.3129 = 0.015754/1.3129 = 0.01200. Entonces t=ln⁡(1+0.01200×10.000)/λ=ln⁡1.12/(5.543×10−10)=204.5 Mat = \ln(1 + 0.01200\times10.000)/\lambda = \ln 1.12/(5.543\times10^{-10}) = 204.5\,\mathrm{Ma}.

Ejercicio 13.3 ★

Una lutita tiene 143Nd/144Nd=0.51180^{143}\mathrm{Nd}/^{144}\mathrm{Nd} = 0.51180 y 147Sm/144Nd=0.1200^{147}\mathrm{Sm}/^{144}\mathrm{Nd} = 0.1200. Calcúlense sus TCHURT_{\mathrm{CHUR}} y TDMT_{\mathrm{DM}} con los valores del capítulo.

Solución

Solución de Ejercicio 13.3.

TCHUR=ln⁡(1+(0.51180−0.512638)/(0.1200−0.1966))/λ=ln⁡(1.010940)/λ=1.67 GaT_{\mathrm{CHUR}} = \ln\bigl(1 + (0.51180 - 0.512638)/(0.1200 - 0.1966)\bigr)/\lambda = \ln(1.010940)/\lambda = 1.67\,\mathrm{Ga}; TDM=ln⁡(1+(0.51180−0.51315)/(0.1200−0.2137))/λ=ln⁡(1.014408)/λ=2.19 GaT_{\mathrm{DM}} = \ln\bigl(1 + (0.51180 - 0.51315)/(0.1200 - 0.2137)\bigr)/\lambda = \ln(1.014408)/\lambda = 2.19\,\mathrm{Ga}.

Ejercicio 13.4 ★

Un bloque de cuarzo sobre una morrena en el sitio de referencia contiene 2.0×1052.0 \times 10^{5} átomos de berilio-10 por gramo. Suponiendo que no hay erosión ni herencia, calcúlese su edad de exposición sin y después con desintegración, con P0=4.96 aˊtomos g−1 yr−1P_0 = 4.96\,\text{átomos}\,\mathrm{g}^{-1}\,\mathrm{yr}^{-1}.

Solución

Solución de Ejercicio 13.4.

Sin desintegración, t=N/P0=2.0×105/4.96=40.3 kyrt = N/P_0 = 2.0 \times 10^{5}/4.96 = 40.3\,\mathrm{kyr}. Con desintegración, t=−ln⁡(1−Nλ/P0)/λ=−ln⁡(1−0.02014)/λ=40.7 kyrt = -\ln(1 - N\lambda/P_0)/\lambda = -\ln(1 - 0.02014)/\lambda = 40.7\,\mathrm{kyr}, un 1 %1\,\% más antiguo.

Ejercicio 13.5 ★★

Demuéstrese que cerca del origen la concordia es una recta de pendiente λ238/λ235\lambda_{238}/\lambda_{235}, y calcúlese. Explíquese por qué la discordia de cristales que perdieron plomo solo muy recientemente pasa por el origen, y qué mide la edad plomo–plomo de tales cristales.

Solución

Solución de Ejercicio 13.5.

Para tt pequeño, x≈λ235tx \approx \lambda_{235}t e y≈λ238ty \approx \lambda_{238}t, por lo que y≈(λ238/λ235)xy \approx (\lambda_{238}/\lambda_{235})x, pendiente 0.15750.1575. Si el plomo se perdió muy recientemente, t2≈0t_2 \approx 0 y la intersección inferior es C(0)C(0), el origen. A lo largo de una recta que pasa por el origen, y/xy/x es constante, y lo es también 207Pb∗/206Pb∗=x/(137.8 y)^{207}\mathrm{Pb}^*/^{206}\mathrm{Pb}^* = x/(137.8\,y): cada cristal proporciona la edad plomo–plomo de la intersección superior, la edad de formación.

Ejercicio 13.6 ★★

Cuatro edades de escalón de una meseta son 296.4±0.8296.4 \pm 0.8, 295.1±1.0295.1 \pm 1.0, 297.0±0.9297.0 \pm 0.9 y (295.8±0.7) Ma(295.8 \pm 0.7)\,\mathrm{Ma} (1σ1\sigma). Calcúlense su media ponderada, su desviación típica y el MSWD. ¿Concuerdan?

Solución

Solución de Ejercicio 13.6.

Pesos 1.56251.5625, 11, 1.23461.2346, 2.04082.0408, W=5.838W = 5.838; xˉ=296.09 Ma\bar x = 296.09\,\mathrm{Ma}, σ=1/W=0.41 Ma\sigma = 1/\sqrt W = 0.41\,\mathrm{Ma}. Las desviaciones cuadráticas ponderadas suman 2.322.32; con f=3f = 3, MSWD=0.77\mathrm{MSWD} = 0.77, por debajo de 1+22/3=2.631 + 2\sqrt{2/3} = 2.63: las edades concuerdan.

Ejercicio 13.7 ★★

Un apatito tiene ρs=1.20×106 cm−2\rho_s = 1.20 \times 10^{6}\,\mathrm{cm}^{-2} y ρi=2.40×106 cm−2\rho_i = 2.40 \times 10^{6}\,\mathrm{cm}^{-2}, con ρd=1.10×106 cm−2\rho_d = 1.10 \times 10^{6}\,\mathrm{cm}^{-2}, g=0.5g = 0.5 y una calibración de laboratorio ζ=350 yr cm2\zeta = 350\,\mathrm{yr}\,\mathrm{cm}^{2}. Calcúlese su edad de trazas de fisión ¿Qué suceso data?

Solución

Solución de Ejercicio 13.7.

λζgρdρs/ρi=1.55125×10−10×350×0.5×1.10×106×0.5=0.014931\lambda\zeta g\rho_d\rho_s/\rho_i = 1.55125\times10^{-10}\times350\times0.5 \times1.10 \times 10^{6}\times0.5 = 0.014931; t=ln⁡(1.014931)/λ=95.5 Mat = \ln(1.014931)/\lambda = 95.5\,\mathrm{Ma}. Si la roca se enfrió rápidamente, data su paso por la zona de recocido parcial del apatito, entre unos 60 y 120 ∘C120\,{}^{\circ}\mathrm{C}; si se enfrió lentamente, es una edad mixta, y las longitudes de las trazas deben aclararlo.

Ejercicio 13.8 ★★

Un apatito contiene 20 ppm20\,\mathrm{ppm} de uranio (99.2742 %99.2742\,\% uranio-238, 0.7204 %0.7204\,\% uranio-235 en átomos), 40 ppm40\,\mathrm{ppm} torio-232 (semivida 1.40×1010 yr1.40 \times 10^{10}\,\mathrm{yr}) y 1.0 nmol g−11.0\,\mathrm{nmol}\,\mathrm{g}^{-1} de helio. Calcúlese su edad (U–Th)/He con la aproximación lineal, y después la edad corregida por eyección alfa con FT=0.75F_T = 0.75.

Solución

Solución de Ejercicio 13.8.

Por gramo: uranio 20×10−6/238.05×NA=5.06×101620\times10^{-6}/238.05\times N_A = 5.06 \times 10^{16} átomos, de los cuales 5.02×10165.02 \times 10^{16} de uranio-238 y 3.64×10143.64 \times 10^{14} de uranio-235; torio 1.038×10171.038 \times 10^{17}, λ232=ln⁡2/1.40×1010=4.95×10−11 yr−1\lambda_{232} = \ln 2/1.40 \times 10^{10} = 4.95 \times 10^{-11}\,\mathrm{yr}^{-1}; helio 10−9NA=6.02×101410^{-9}N_A = 6.02 \times 10^{14}. La tasa de producción es 8λ238238U+7λ235235U+6λ232Th=(6.23+0.25+3.08)×107=9.57×107 aˊtomos g−1 yr−18\lambda_{238}{}^{238}\mathrm U + 7\lambda_{235}{}^{235}\mathrm U + 6\lambda_{232}\mathrm{Th} = (6.23 + 0.25 + 3.08)\times10^7 = 9.57 \times 10^{7}\,\text{átomos}\,\mathrm{g}^{-1}\,\mathrm{yr}^{-1}, y t=6.02×1014/9.57×107=6.29 Myrt = 6.02 \times 10^{14}/9.57 \times 10^{7} = 6.29\,\mathrm{Myr}. Corregida: 6.29/0.75=8.39 Ma6.29/0.75 = 8.39\,\mathrm{Ma}.

Ejercicio 13.9 ★★

Para el berilio-10 en el emplazamiento de referencia, calcúlese la concentración de saturación sin erosión y el tiempo necesario para alcanzar el 90 %90\,\% de ella. Calcúlese la concentración estacionaria de una superficie que se erosiona a 10 m Myr−110\,\mathrm{m}\,\mathrm{Myr}^{-1} (granito, 2.7 g cm−32.7\,\mathrm{g}\,\mathrm{cm}^{-3}, Λ=160 g cm−2\Lambda = 160\,\mathrm{g}\,\mathrm{cm}^{-2}).

Solución

Solución de Ejercicio 13.9.

P0/λ=4.96/4.994×10−7=9.93×106 aˊtomos g−1P_0/\lambda = 4.96/4.994 \times 10^{-7} = 9.93 \times 10^{6}\,\text{átomos}\,\mathrm{g}^{-1}, alcanzada al 90 %90\,\% tras ln⁡10/λ=4.6 Myr\ln 10/\lambda = 4.6\,\mathrm{Myr}. Con 10 m Myr−110\,\mathrm{m}\,\mathrm{Myr}^{-1} =1×10−3 cm yr−1= 1 \times 10^{-3}\,\mathrm{cm}\,\mathrm{yr}^{-1}: ρε/Λ=2.7×10−3/160=1.69×10−5 yr−1\rho\varepsilon/\Lambda = 2.7\times10^{-3}/160 = 1.69 \times 10^{-5}\,\mathrm{yr}^{-1}, μ=1.74×10−5 yr−1\mu = 1.74 \times 10^{-5}\,\mathrm{yr}^{-1} y N=4.96/1.74×10−5=2.85×105 aˊtomos g−1N = 4.96/1.74 \times 10^{-5} = 2.85 \times 10^{5}\,\text{átomos}\,\mathrm{g}^{-1}, 35 veces menos.

Ejercicio 13.10 ★★★

Se combinan dos medidas x1±σ1x_1 \pm \sigma_1 y x2±σ2x_2 \pm \sigma_2 de una misma edad. Demuéstrese que su MSWD respecto a la media ponderada es (x1−x2)2/(σ12+σ22)(x_1 - x_2)^2/(\sigma_1^2 + \sigma_2^2), y demuéstrese directamente que su esperanza es 1 cuando ambas miden la misma edad.

Solución

Solución de Ejercicio 13.10.

Con W=w1+w2W = w_1 + w_2, x1−xˉ=w2(x1−x2)/Wx_1 - \bar x = w_2(x_1 - x_2)/W y x2−xˉ=−w1(x1−x2)/Wx_2 - \bar x = -w_1(x_1 - x_2)/W, por lo que ∑wi(xi−xˉ)2=(x1−x2)2(w1w22+w2w12)/W2=(x1−x2)2w1w2/W=(x1−x2)2/(σ12+σ22)\sum w_i(x_i - \bar x)^2 = (x_1 - x_2)^2(w_1w_2^2 + w_2w_1^2)/W^2 = (x_1 - x_2)^2w_1w_2/W = (x_1 - x_2)^2/(\sigma_1^2 + \sigma_2^2); con f=1f = 1 esto es el MSWD. Si ambas miden la misma edad, x1−x2x_1 - x_2 tiene media 0 y varianza σ12+σ22\sigma_1^2 + \sigma_2^2, por lo que el valor esperado de su cuadrado dividido entre esa varianza es 1.

Ejercicio 13.11 ★★★

Un circón de 2700 Ma2700\,\mathrm{Ma} perdió plomo a 500 Ma500\,\mathrm{Ma} y hoy presenta 206Pb∗/238U=0.300^{206}\mathrm{Pb}^*/^{238}\mathrm{U} = 0.300. ¿Qué fracción de su plomo perdió? Calcúlese su 207Pb∗/235U^{207}\mathrm{Pb}^*/^{235}\mathrm{U}, sus dos edades uranio–plomo y, por iteración, su edad plomo–plomo. ¿Cuál de las tres edades tiene significado?

Solución

Solución de Ejercicio 13.11.

C(2700)=(13.283,0.52019)C(2700) = (13.283, 0.52019), C(500)=(0.63628,0.08065)C(500) = (0.63628, 0.08065). f=(0.52019−0.300)/(0.52019−0.08065)=0.501f = (0.52019 - 0.300)/(0.52019 - 0.08065) = 0.501, la mitad de su plomo. x=0.499×13.283+0.501×0.63628=6.948x = 0.499\times13.283 + 0.501\times0.63628 = 6.948. t206=1691 Mat_{206} = 1691\,\mathrm{Ma}, t207=2105 Mat_{207} = 2105\,\mathrm{Ma}; la razón x/(137.8 y)=0.1681x/(137.8\,y) = 0.1681 da, por iteración, una edad plomo–plomo de 2537 Ma2537\,\mathrm{Ma}. Ninguna de las tres tiene significado: son mezclas de los dos eventos, y solo las intersecciones las datan. La edad plomo–plomo es un mínimo para la cristalización.

Ejercicio 13.12 ★★★

Una superficie de lecho rocoso en el emplazamiento de referencia contiene 1.0×1061.0 \times 10^{6} átomos de berilio-10 por gramo. Calcúlese su edad de exposición suponiendo que no hay erosión, y después su tasa de erosión suponiendo un estado estacionario. Demuéstrese que cada par (t,ε)(t, \varepsilon) comprendido entre estos dos límites se ajusta a la medición, e indíquese qué segunda medida los separaría.

Solución

Solución de Ejercicio 13.12.

Sin erosión: t=−ln⁡(1−1.0×106×4.994×10−7/4.96)/λ=212 kyrt = -\ln(1 - 1.0 \times 10^{6}\times4.994 \times 10^{-7}/4.96)/\lambda = 212\,\mathrm{kyr}. Estado estacionario: ρε/Λ=P0/N−λ=4.46×10−6 yr−1\rho\varepsilon/\Lambda = P_0/N - \lambda = 4.46 \times 10^{-6}\,\mathrm{yr}^{-1}, ε=2.64×10−4 cm yr−1=2.6 m Myr−1\varepsilon = 2.64 \times 10^{-4}\,\mathrm{cm}\,\mathrm{yr}^{-1} = 2.6\,\mathrm{m}\,\mathrm{Myr}^{-1}. Para cualquier ε\varepsilon inferior a este valor, N=P0(1−e−μt)/μN = P_0(1 - e^{-\mu t})/\mu fija un tt, desde 212 kyr212\,\mathrm{kyr} hasta el infinito: una medida tiene dos incógnitas. Un segundo nucleido con otra semivida, como el aluminio-26 en el mismo cuarzo, o la concentración a profundidad en un testigo de sondeo, los separa.

13.8 Problema: circones del granito de Varden

Problema 13.1

Problema de fin de semana — la cristalización, el enfriamiento y la exhumación de un granito inventado, interpretado con cuatro relojes

El granito de las colinas de Varden, región inventada de la serie, fue exhumado ya en el Triásico, sepultado bajo la cobertera mesozoica y devuelto a la superficie. Cinco de sus circones dan, con errores de una desviación típica en yy (tómese xx como exacto):

circónZ1Z2Z3Z4Z5
x=207Pb∗/235Ux = {}^{207}\mathrm{Pb}^*/^{235}\mathrm U0.34670.32600.30270.27930.2535
y=206Pb∗/238Uy = {}^{206}\mathrm{Pb}^*/^{238}\mathrm U0.047940.045050.042070.038920.03553
σy\sigma_y0.000100.000090.000080.000080.00007

Utilícense las constantes del capítulo. Todos los análisis son datos de ejercicio.

Parte I — Los circones.

  1. ¿Son concordantes? Calcúlense las dos edades uranio–plomo de cada circón.
  2. ¿Por qué disminuyen las edades de Z1 a Z5, con t207>t206t_{207} > t_{206} para cada uno?
  3. Calcúlense las coordenadas Tera–Wasserburg de Z1 y Z5.
  4. Con los pesos wi=1/σy,i2w_i = 1/\sigma_{y,i}^2, calcúlense las medias ponderadas xˉ\bar x y yˉ\bar y, así como Sxx=∑wi(xi−xˉ)2S_{xx} = \sum w_i(x_i - \bar x)^2.
  5. Calcúlense la pendiente b=∑wi(xi−xˉ)(yi−yˉ)/Sxxb = \sum w_i(x_i - \bar x)(y_i - \bar y)/S_{xx} y la ordenada en el origen aa de la discordia.
  6. Compárese bb con la pendiente de la concordia cerca de 310 Ma310\,\mathrm{Ma}. ¿Qué anticipa esto?

Parte II — Las intersecciones.

  1. Búsquese la intersección superior por el método de Newton a partir de 310 Ma310\,\mathrm{Ma}.
  2. Búsquese la intersección inferior a partir de 100 Ma100\,\mathrm{Ma}.
  3. Calcúlense los residuos ponderados y el MSWD. ¿Es aceptable el ajuste?
  4. La desviación típica de la pendiente es 1/Sxx1/\sqrt{S_{xx}}. Calcúlese.
  5. Gírese la recta alrededor de (xˉ,yˉ)(\bar x, \bar y) hasta las pendientes b±2sbb \pm 2s_b y recalcúlese la intersección superior para cada una. ¿Cuál es la incertidumbre de la edad?
  6. ¿Qué fracción de su plomo perdieron Z1 y Z5?

Parte III — La biotita. Una biotita del granito, irradiada con un monitor de 28.201 Ma28.201\,\mathrm{Ma} que da Rm=1.3129R_m = 1.3129, se desgasifica en diez etapas:

etapa12345678910
39Ar^{39}\mathrm{Ar} (%\%)468121415141296
RR7.22511.22613.93714.89114.82014.86914.91814.84714.90714.804
edad (Ma\mathrm{Ma})228.0279.0295.3296.2297.1295.8296.9
σ\sigma (Ma\mathrm{Ma})0.60.91.01.11.11.11.11.11.11.1
  1. Calcúlese JJ; redondéese a 0.01200.
  2. Calcúlense las edades de las etapas 1, 4 y 10.
  3. ¿Qué etapas forman una meseta según los criterios del capítulo?
  4. Calcúlense la edad meseta, su desviación típica y su MSWD.
  5. Calcúlense la razón y la edad del gas total. ¿Por qué es más joven que la meseta?
  6. Suponiendo que el granito cristalizó a 700 ∘C700\,{}^{\circ}\mathrm{C} y que la biotita cierra a 375 ∘C375\,{}^{\circ}\mathrm{C}, calcúlese la velocidad media de enfriamiento.
  7. ¿Qué pudo ocurrir en la intersección inferior? ¿Es el espectro de argón consistente con ella?

Parte IV — La exhumación. Cuatro apatitos, muestreados en una ladera de las colinas, dan edades (U–Th)/He de 43.1, 46.0, 49.2 y 52.0 Ma52.0\,\mathrm{Ma} a 310, 370, 430 y 490 m490\,\mathrm{m}. La cima del granito, en una cresta suave de colina en un emplazamiento donde P0=5.6 aˊtomos g−1 yr−1P_0 = 5.6\,\text{átomos}\,\mathrm{g}^{-1}\,\mathrm{yr}^{-1}, contiene 1.64×1051.64 \times 10^{5} átomos de berilio-10 por gramo de cuarzo.

  1. Calcúlese la pendiente de la relación edad–altitud, en metros por millón de años.
  2. Con una temperatura de cierre de 65 ∘C65\,{}^{\circ}\mathrm{C}, una superficie a 10 ∘C10\,{}^{\circ}\mathrm{C} y una geoterma de 30 ∘C km−130\,{}^{\circ}\mathrm{C}\,\mathrm{km}^{-1}, ¿a qué profundidad cerró el apatito más bajo? ¿Qué velocidad media de exhumación se deduce desde entonces?
  3. Calcúlese la tasa de erosión actual a partir del berilio-10, suponiendo un estado estacionario.
  4. ¿Qué tasa da la fórmula del volumen de primer año, sin desintegración?
  5. ¿Sobre qué intervalo de tiempo se promedia la tasa del berilio-10? Compárense las tres tasas.
  6. Establézcase la edad de cristalización del granito de Varden con su incertidumbre.
Solución

Solución de Problema 13.1.

1. t206t_{206} / t207t_{207}, en Ma\mathrm{Ma}: Z1 301.9 / 302.2, Z2 284.1 / 286.5, Z3 265.7 / 268.5, Z4 246.1 / 250.1, Z5 225.1 / 229.4. Un error de 0.2 %0.2\,\% en yy equivale a unos 0.6 Ma0.6\,\mathrm{Ma}: Z1 es concordante dentro del error; las demás son discordantes, cada vez más discordantes desde Z2 hasta Z5.

2. Los circones perdieron fracciones crecientes de su plomo en un evento posterior: se sitúan cada vez más abajo sobre la discordia. La discordia queda por debajo de la concordia cóncava, donde xx es mayor para un yy dado: t207>t206t_{207} > t_{206}.

3. Z1: X=1/0.04794=20.859X = 1/0.04794 = 20.859, Y=0.3467/(137.8×0.04794)=0.05248Y = 0.3467/(137.8\times0.04794) = 0.05248. Z5: X=28.145X = 28.145, Y=0.05178Y = 0.05178.

4. Pesos 1.0001.000, 1.2351.235, 1.5631.563, 1.5631.563, 2.0412.041 (en 10810^8); xˉ=0.29402\bar x = 0.29402, yˉ=0.040892\bar y = 0.040892, Sxx=7.845×105S_{xx} = 7.845 \times 10^{5}.

5. b=0.13260b = 0.13260, a=yˉ−bxˉ=0.001905a = \bar y - b\bar x = 0.001905.

6. A 310 Ma310\,\mathrm{Ma}, dy/dx=0.1575 e(λ238−λ235)t=0.1218\mathrm dy/\mathrm dx = 0.1575\, e^{(\lambda_{238} - \lambda_{235})t} = 0.1218. La discordia solo es 9 %9\,\% más empinada que la concordia: la corta con un ángulo pequeño, y la intersección superior será imprecisa.

7. Newton desde 310 Ma310\,\mathrm{Ma}: t=311.2 Mat = 311.2\,\mathrm{Ma}.

8. Desde 100 Ma100\,\mathrm{Ma}: t=99.7 Mat = 99.7\,\mathrm{Ma}.

9. (yi−a−bxi)/σi=0.64(y_i - a - bx_i)/\sigma_i = 0.64, −0.91-0.91, 0.350.35, −0.24-0.24, 0.170.17; MSWD=1.435/3=0.48\mathrm{MSWD} = 1.435/3 = 0.48, por debajo del límite 2.632.63 para f=3f = 3: los circones se ajustan a una sola recta.

10. sb=1/7.845×105=0.00113s_b = 1/\sqrt{7.845 \times 10^{5}} = 0.00113.

11. Con b−2sb=0.13034b - 2s_b = 0.13034, la intersección superior es 322.9 Ma322.9\,\mathrm{Ma}; con b+2sb=0.13485b + 2s_b = 0.13485, 302.6 Ma302.6\,\mathrm{Ma}. La edad es 311±10 Ma311 \pm 10\,\mathrm{Ma} a dos desviaciones típicas.

12. f=(x(311.2)−xi)/(x(311.2)−x(99.7))f = (x(311.2) - x_i)/(x(311.2) - x(99.7)): Z1 4.7 %4.7\,\%, Z5 41 %41\,\%.

13. J=0.015754/1.3129=0.0119999J = 0.015754/1.3129 = 0.0119999, 0.012000.01200.

14. Paso 1: 150.0 Ma150.0\,\mathrm{Ma}; paso 4: 296.6 Ma296.6\,\mathrm{Ma}; paso 10: 295.0 Ma295.0\,\mathrm{Ma}.

15. Pasos 4 a 10: siete pasos consecutivos, 82 %82\,\% del argón-39, sin tendencia. Los pasos 1 a 3 suben en escalera.

16. Errores iguales: la media es 2072.9/7=296.1 Ma2072.9/7 = 296.1\,\mathrm{Ma}, σ=1.1/7=0.42 Ma\sigma = 1.1/\sqrt7 = 0.42\,\mathrm{Ma}; la suma de los cuadrados de las desviaciones es 3.833.83, MSWD=3.83/(1.21×6)=0.53\mathrm{MSWD} = 3.83/(1.21\times6) = 0.53, por debajo de 2.152.15.

17. Rˉ=∑fiRi=14.270\bar R = \sum f_iR_i = 14.270, t=285.2 Mat = 285.2\,\mathrm{Ma}: el argón perdido de los bordes, liberado en los primeros pasos, reduce la media, como lo haría con una edad potasio–argón.

18. (700−375)/(311.2−296.1)=21.5 ∘C Myr−1(700 - 375)/(311.2 - 296.1) = 21.5\,{}^{\circ}\mathrm{C}\,\mathrm{Myr}^{-1} (con un intervalo amplio, ya que la edad del circón se conoce con ±10 Ma\pm10\,\mathrm{Ma}).

19. Hacia 100 Ma100\,\mathrm{Ma} el granito estaba enterrado bajo la cobertera mesozoica, y un calentamiento o una circulación de fluidos hizo que los circones perdieran plomo. La biotita perdió argón de sus bordes, lo que registra la escalera; su paso más bajo, 150 Ma150\,\mathrm{Ma}, es una mezcla de argón antiguo y reajustado. La meseta se conservó: la biotita no se calentó lo suficiente para perder todo su argón.

20. La pendiente de mínimos cuadrados de la altitud frente a la edad es 20.1 m Myr−120.1\,\mathrm{m}\,\mathrm{Myr}^{-1}.

21. (65−10)/30=1.83 km(65 - 10)/30 = 1.83\,\mathrm{km}; desde 43.1 Ma43.1\,\mathrm{Ma}, una media de 1830/43.1=42 m Myr−11830/43.1 = 42\,\mathrm{m}\,\mathrm{Myr}^{-1}.

22. ρε/Λ=P0/N−λ=3.415×10−5−0.050×10−5=3.365×10−5 yr−1\rho\varepsilon/\Lambda = P_0/N - \lambda = 3.415 \times 10^{-5} - 0.050 \times 10^{-5} = 3.365 \times 10^{-5}\,\mathrm{yr}^{-1}; ε=3.365×10−5×160/2.7=1.99×10−3 cm yr−1=19.9 m Myr−1\varepsilon = 3.365 \times 10^{-5}\times 160/2.7 = 1.99 \times 10^{-3}\,\mathrm{cm}\,\mathrm{yr}^{-1} = 19.9\,\mathrm{m}\,\mathrm{Myr}^{-1}.

23. P0Λ/(ρN)=20.2 m Myr−1P_0\Lambda/(\rho N) = 20.2\,\mathrm{m}\,\mathrm{Myr}^{-1}, 1.5 %1.5\,\% más: la desintegración apenas influye a esta tasa.

24. Λ/(ρε)=160/(2.7×1.99×10−3)=30 kyr\Lambda/(\rho\varepsilon) = 160/(2.7\times1.99 \times 10^{-3}) = 30\,\mathrm{kyr}. El granito se elevó a unos 20 m Myr−120\,\mathrm{m}\,\mathrm{Myr}^{-1} entre 52 y 43 Ma43\,\mathrm{Ma}, a 42 m Myr−142\,\mathrm{m}\,\mathrm{Myr}^{-1} de media desde entonces, y erosiona a 20 m Myr−120\,\mathrm{m}\,\mathrm{Myr}^{-1} en la actualidad: la exhumación fue más rápida en algún momento del Cenozoico.

25. El granito de Varden cristalizó hace (311±10) Ma(311 \pm 10)\,\mathrm{Ma} (intersección superior, a dos desviaciones típicas).

本章定义的术语

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