地球与空间科学 · 第 3 册 · 本科一年级

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

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

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

10El balance térmico de la Tierra

Bajo una montaña de los Alpes japoneses, en una antigua cavidad minera, un tanque con mil toneladas de centelleador líquido vigila destellos de luz muy débiles. La mayoría procede de antineutrinos producidos en los reactores nucleares del país de alrededor, pero algunos vienen de abajo: el uranio y el torio del interior de la Tierra emiten antineutrinos al desintegrarse, y estos atraviesan todo el planeta como si estuviera vacío. En 2135 días de recuento el detector registró 841 sucesos candidatos, de los que 106 se atribuyeron a la propia Tierra. A partir de ese puñado de destellos, su equipo estimó en 2011 que la desintegración del uranio y el torio calienta la Tierra con unos 20 teravatios. El planeta en conjunto pierde unos 47. Este capítulo explica la diferencia: de dónde procede el calor de la Tierra, cómo se mide el calor que sale de ella, cómo sube la temperatura con la profundidad cuando el calor se transporta por conducción, por qué el manto debe convectar, y por qué las placas son la frontera superior y fría de esa convección.

Lo que ya sabes

El calor fluye de lo caliente a lo frío; la conductividad térmica kk de una roca mide lo fácilmente que el calor la atraviesa, y la difusividad térmica κ=k/(ρc)\kappa = k/(\rho c), lo deprisa que se extiende un cambio de temperatura (física, segundo año universitario). Los núcleos radiactivos se desintegran con una semivida fija (física). La temperatura sube con la profundidad a lo largo del geotermo, y el manto convecta (libro 2, años 11 y 12). La litosfera es la cáscara exterior resistente de la Tierra (Capítulo 9).

10.1 De dónde viene el calor

Definición 10.1 (Flujo de calor)

El flujo de calor qq en un punto de la superficie de la Tierra, o densidad de flujo de calor superficial, es la potencia de origen interno que atraviesa la superficie allí, por unidad de área, contada positiva hacia fuera. Se mide en W m−2\mathrm{W}\,\mathrm{m}^{-2}, en la práctica en mW m−2\mathrm{mW}\,\mathrm{m}^{-2}. Integrado sobre toda la superficie, da la pérdida de calor de la Tierra Q=∮q dSQ = \oint q\,\mathrm dS.

El flujo de calor se ha medido en decenas de miles de lugares. Una recopilación de 38 34738\,347 medidas, publicada en 2010, da los valores medios siguientes, contando por separado los océanos y los continentes:

superficie (m2\mathrm{m}^{2})pérdida de calor (TW\mathrm{TW})flujo de calor medio (mW m−2\mathrm{mW}\,\mathrm{m}^{-2})
continentes2.073×10142.073 \times 10^{14}14.770.9
océanos3.028×10143.028 \times 10^{14}31.9105.4
Tierra entera5.101×10145.101 \times 10^{14}46.791.6

El total, redondeado a su incertidumbre, es Q=(47±2) TWQ = (47 \pm 2)\,\mathrm{TW}. Dos tercios de él salen por el fondo oceánico, que es más joven y más delgado que los continentes. Es un flujo pequeño: el Sol entrega a la Tierra unos 173 000 TW173\,000\,\mathrm{TW}, y el calor interno representa el 0.03 %0.03\,\% del balance de energía en la superficie. Pero es el flujo que impulsa los volcanes, las placas y la formación de montañas, y tiene dos fuentes.

Definición 10.2 (Producción de calor radiogénico)

La producción de calor radiogénico AA de una roca es el calor liberado por unidad de volumen y de tiempo por la desintegración de los núcleos radiactivos que contiene, en W m−3\mathrm{W}\,\mathrm{m}^{-3}, en la práctica en µW m−3\text{µ}\mathrm{W}\,\mathrm{m}^{-3}. Una roca de densidad ρ\rho que contiene fracciones másicas CiC_i de los elementos productores de calor, cada uno de los cuales libera HiH_i vatios por kilogramo de elemento, tiene

A=ρ∑iCiHi.A = \rho \sum_i C_i H_i .

Cuatro isótopos producen casi todo este calor, porque son de vida larga, suficientemente abundantes y energéticos: X238X22238U\ce{^{238}U}, X235X22235U\ce{^{235}U}, X232X22232Th\ce{^{232}Th} y X40X2240K\ce{^{40}K}. El calor por kilogramo de cada elemento es el calor por kilogramo de su isótopo radiactivo multiplicado por la proporción de ese isótopo en el elemento (Ejercicio 10.3):

isótoposemividacalor por kg de isótopoproporción delcalor por kg de elemento
(Gyr\mathrm{Gyr})(W kg−1\mathrm{W}\,\mathrm{kg}^{-1})elemento (átomos)(W kg−1\mathrm{W}\,\mathrm{kg}^{-1})
X238X22238U\ce{^{238}U}4.4689.513×10−59.513 \times 10^{-5}99.274 %99.274\,\%U: 9.85×10−59.85 \times 10^{-5}
X235X22235U\ce{^{235}U}0.7045.685×10−45.685 \times 10^{-4}0.720 %0.720\,\%
X232X22232Th\ce{^{232}Th}14.02.628×10−52.628 \times 10^{-5}casi todoTh: 2.63×10−52.63 \times 10^{-5}
X40X2240K\ce{^{40}K}1.2482.847×10−52.847 \times 10^{-5}0.0117 %0.0117\,\%K: 3.40×10−93.40 \times 10^{-9}

El potasio libera treinta mil veces menos calor por kilogramo que el uranio, pero es un elemento principal de muchas rocas, contado en tanto por ciento, mientras que el uranio y el torio se cuentan en partes por millón. Las rocas ricas en los tres son los granitos y sus parientes; las más pobres son las peridotitas del manto. Una recopilación de 21882188 análisis da estas medianas, en µW m−3\text{µ}\mathrm{W}\,\mathrm{m}^{-3}:

riolita2.55lutitas1.44gabro0.35
granito2.43areniscas0.82basalto0.21
granodiorita1.88rocas carbonatadas0.42rocas ultramáficas0.11

La dispersión dentro de cada grupo es grande: un cuarto de los granitos produce menos de 1.74 µW m−31.74\,\text{µ}\mathrm{W}\,\mathrm{m}^{-3}, un cuarto más de 3.23 µW m−33.23\,\text{µ}\mathrm{W}\,\mathrm{m}^{-3}. La corteza continental está enriquecida en estos elementos, y la corteza oceánica y el manto empobrecidos; pero el manto es tan grande que aún contiene una gran parte de las fuentes de calor del planeta.

Definición 10.3 (Enfriamiento secular)

El enfriamiento secular es la parte de la pérdida de calor de la Tierra suministrada por la lenta disminución de la temperatura de su interior, manto y núcleo: la liberación del calor almacenado desde la formación del planeta, a un ritmo −∑jMjcj dTj/dt-\sum_j M_j c_j\,\mathrm dT_j/\mathrm dt para depósitos de masa MjM_j, calor específico cjc_j y temperatura media TjT_j.

La conservación de la energía da el balance de toda la Tierra: el calor perdido en la superficie es la producción radiogénica más el enfriamiento secular,

Q=Hrad+Qsc.Q = H_{\text{rad}} + Q_{\text{sc}} .

Los geoneutrinos miden una parte del primer término: el resultado de 2011, que combinó dos detectores, atribuye 20 TW20\,\mathrm{TW} al uranio y al torio, y se sabe que el X40X2240K\ce{^{40}K}, cuyos antineutrinos son demasiado débiles para detectarse así, añade unos 4 TW4\,\mathrm{TW}.

Alrededor de la mitad de la pérdida de calor es, por tanto, radiogénica, y la otra mitad, unos 23 TW23\,\mathrm{TW}, es enfriamiento secular. Parte de ese enfriamiento es del núcleo: el calor que atraviesa la frontera núcleo–manto se estima en 5 a 15 TW5\text{ a }15\,\mathrm{TW}. Estas cifras son inciertas: la medida misma tiene una incertidumbre de unos 9 TW9\,\mathrm{TW}, y los resultados de los geoneutrinos hasta ahora solo acotan la parte radiogénica entre 15 a 41 TW15\text{ a }41\,\mathrm{TW}, mientras que los modelos de Tierra de distintos tipos varían entre 14 a 46 TW14\text{ a }46\,\mathrm{TW}. El total, medido en la superficie, es el término mejor conocido del balance (figura siguiente).

El presupuesto térmico de la Tierra. Arriba: el calor perdido a través del fondo oceánico y de los continentes, (47 ± 2)\, TW en total. Abajo: sus fuentes, con la parte radiogénica del uranio y del torio medida con geoneutrinos y la del potasio estimada; el resto, por diferencia, es enfriamiento secular, del que el núcleo aporta una parte mal conocida. La barra naranja muestra el intervalo de la parte radiogénica permitido por las medidas de geoneutrinos.
El presupuesto térmico de la Tierra. Arriba: el calor perdido a través del fondo oceánico y de los continentes, (47±2) TW(47 \pm 2)\,\mathrm{TW} en total. Abajo: sus fuentes, con la parte radiogénica del uranio y del torio medida con geoneutrinos y la del potasio estimada; el resto, por diferencia, es enfriamiento secular, del que el núcleo aporta una parte mal conocida. La barra naranja muestra el intervalo de la parte radiogénica permitido por las medidas de geoneutrinos.

Las fuentes radiogénicas fueron mayores en el pasado, porque los isótopos eran más abundantes: X235X22235U\ce{^{235}U} y X40X2240K\ce{^{40}K}, los de vida corta, contribuían entonces mucho más (libro 2, año 12). El motor térmico de la Tierra se apaga despacio, por ambas causas.

10.2 Medir el flujo de calor

El flujo de calor no se mide directamente: se calcula a partir de un gradiente de temperatura y de una conductividad.

Proposición 10.4 (Ley de Fourier)

En un medio de conductividad térmica kk, el calor se conduce cuesta abajo del gradiente de temperatura, a una tasa por unidad de área proporcional a él. Con zz la profundidad, contada positiva hacia abajo, y qq contado positivo hacia arriba,

q=k dTdz.q = k\,\frac{\mathrm dT}{\mathrm dz} .

Demostración. Admitido a este nivel. ∎

La ley está establecida en el volumen de segundo año universitario de los libros de física. La mayoría de las rocas conducen el calor a 1.5 a 3 W m−1 K−11.5\text{ a }3\,\mathrm{W}\,\mathrm{m}^{-1}\,\mathrm{K}^{-1} por debajo de 200 ∘C200\,{}^{\circ}\mathrm{C}, y menos cuando están calientes; los modelos del flujo de calor continental dan a la corteza superior granítica unos 2.7 W m−1 K−12.7\,\mathrm{W}\,\mathrm{m}^{-1}\,\mathrm{K}^{-1};

un gradiente de 25 K km−125\,\mathrm{K}\,\mathrm{km}^{-1} en ella significa un flujo de calor de unos 70 mW m−270\,\mathrm{mW}\,\mathrm{m}^{-2}.

Cuando las rocas que atraviesa un sondeo tienen conductividades distintas, el gradiente cambia de una a otra, pero el flujo de calor no, mientras no se produzca ni se almacene calor entre ellas. La proposición siguiente convierte eso en una prueba.

Proposición 10.5 (Diagrama de Bullard)

Sea R(z)=∫0zdζ/k(ζ)R(z) = \int_0^z \mathrm d\zeta/k(\zeta) la resistencia térmica, por unidad de área, de las rocas entre la superficie y la profundidad zz, en m2 K W−1\mathrm{m}^{2}\,\mathrm{K}\,\mathrm{W}^{-1}. Si el flujo de calor qq es el mismo a toda profundidad de un sondeo, la temperatura es una función lineal de la resistencia térmica: T(z)=T0+q R(z)T(z) = T_0 + q\,R(z).

Demostración. Por la ley de Fourier, dT/dz=q/k(z)\mathrm dT/\mathrm dz = q/k(z). Integrando desde 00 hasta zz con qq constante, T(z)−T0=q∫0zdζ/k(ζ)=q R(z)T(z) - T_0 = q\int_0^z \mathrm d\zeta/k(\zeta) = q\,R(z). ∎

Método 10.6 (El flujo de calor a partir de un sondeo)

  1. Registrese la temperatura a lo largo del sondeo mucho después de la perforación, cuando la perturbación del fluido de perforación se ha disipado, y descártense los primeros decenas de metros, que siguen las estaciones.
  2. Mídase en laboratorio la conductividad de las muestras de testigo de cada roca atravesada, y calcúlese la resistencia térmica R(z)R(z) a cada profundidad registrada.
  3. Represéntese TT frente a RR: la pendiente de la recta es el flujo de calor. Una curvatura revela agua que circula por las rocas, un cambio pasado de la temperatura superficial o un error en una conductividad.
  4. Corríjase, donde haga falta, por el relieve, que curva las isotermas bajo colinas y valles, y por el calentamiento lento desde la última glaciación, que todavía perturba el primer kilómetro aproximadamente.

En el mar no hace falta un sondeo. La temperatura en el fondo del océano profundo es muy estable, de modo que unos pocos metros de sedimento blando ya conservan el gradiente de equilibrio. Una lanza de acero con una hilera de sensores de temperatura se deja caer desde un barco y se clava en el sedimento con un peso pesado (figura siguiente); los sensores registran el gradiente, y un pulso de calor enviado a lo largo de la lanza da la conductividad del sedimento in situ. El fondo oceánico tiene una trampa. Sobre la corteza joven, el agua de mar circula por el basalto agrietado y se lleva gran parte del calor por advección, de modo que el flujo de calor conductivo medido queda muy por debajo de la pérdida real. La pérdida a través del fondo oceánico más joven se calcula, por tanto, a partir del modelo de enfriamiento de Sección 10.5, calibrado en corteza antigua donde la circulación se ha detenido: esta parte calculada, 23.1 TW23.1\,\mathrm{TW}, es aproximadamente la mitad del total mundial.

Medición del flujo de calor (esquema, no a escala). (a) En tierra, la temperatura se registra a lo largo de un sondeo una vez que ha vuelto al equilibrio, y la conductividad de cada roca se mide en testigos: el gradiente es más empinado donde la conductividad es menor, y el flujo de calor, su producto, es el mismo. (b) En el mar, una lanza con sensores se clava unos metros en el sedimento del fondo oceánico.
Medición del flujo de calor (esquema, no a escala). (a) En tierra, la temperatura se registra a lo largo de un sondeo una vez que ha vuelto al equilibrio, y la conductividad de cada roca se mide en testigos: el gradiente es más empinado donde la conductividad es menor, y el flujo de calor, su producto, es el mismo. (b) En el mar, una lanza con sensores se clava unos metros en el sedimento del fondo oceánico.

Sobre el terreno — Flujo de calor en el mar

Una sonda de flujo de calor colgada de la grúa de popa de un buque de investigación (ilustración generada por IA).

En un buque de investigación, la sonda de flujo de calor es una lanza pesada suspendida de la grúa de popa. Se baja con el cable del barco hasta justo encima del fondo marino, y luego se suelta: el peso la clava en el sedimento, y la tripulación espera en cubierta a que los sensores se estabilicen a la temperatura del fango, después la saca y se desplaza a la siguiente estación, a poca distancia. Una línea de estas estaciones a través del flanco de una dorsal da la disminución del flujo de calor con la edad del fondo, y la dispersión entre estaciones vecinas muestra dónde circula agua por la corteza subyacente.

10.3 El geotermo conductivo estacionario

Definición 10.7 (Geotermo)

El geotermo bajo un punto de la superficie es la temperatura en función de la profundidad allí, T(z)T(z). Es un geotermo conductivo estacionario cuando no cambia con el tiempo y el calor se transporta a lo largo de él solo por conducción.

Dentro de la litosfera, lejos de los volcanes activos y de la circulación de agua, el calor se transporta por conducción, y a escala de millones de años el geotermo es estacionario. Su forma se deduce entonces de un balance de energía.

Teorema 10.8 (Conducción estacionaria con producción de calor)

A lo largo de un geotermo conductivo estacionario, en un medio de conductividad uniforme kk con producción de calor A(z)A(z) y sin flujo horizontal de calor,

dqdz=−A(z),es decirk d2Tdz2=−A(z).\frac{\mathrm dq}{\mathrm dz} = -A(z), \qquad\text{es decir}\qquad k\,\frac{\mathrm d^2T}{\mathrm dz^2} = -A(z).

Demostración. Consideremos una losa horizontal de área SS entre las profundidades zz y z+dzz + \mathrm dz. El calor entra por su base a razón de q(z+dz) Sq(z + \mathrm dz)\,S, sale por su techo a razón de q(z) Sq(z)\,S, y se produce en su interior a razón de A(z) S dzA(z)\,S\,\mathrm dz. En un régimen estacionario no se almacena calor, así que q(z) S=q(z+dz) S+A(z) S dzq(z)\,S = q(z + \mathrm dz)\,S + A(z)\,S\,\mathrm dz, lo que da dq/dz=−A\mathrm dq/\mathrm dz = -A. Por la ley de Fourier q=k dT/dzq = k\,\mathrm dT/\mathrm dz, y con kk uniforme, k T′′=−Ak\,T'' = -A. ∎

Leído desde la superficie hacia abajo, el flujo de calor disminuye con la profundidad en el calor producido sobre cada nivel: lo que atraviesa la superficie es el calor que viene de abajo más todo el calor producido por el camino. Dos distribuciones de AA dan soluciones cerradas.

Proposición 10.9 (Una capa radiogénica uniforme)

Sea la producción de calor AA desde la superficie hasta la profundidad hh y cero debajo, la superficie a la temperatura T0T_0, y el flujo de calor en la superficie q0q_0. Entonces, para 0≤z≤h0 \le z \le h

q(z)=q0−Az,T(z)=T0+q0zk−Az22k,q(z) = q_0 - A z, \qquad T(z) = T_0 + \frac{q_0 z}{k} - \frac{A z^2}{2k},

y debajo de hh el flujo de calor es uniforme, qm=q0−Ahq_m = q_0 - A h, y el geotermo una recta de pendiente qm/kq_m/k.

Demostración. Integrando dq/dz=−A\mathrm dq/\mathrm dz = -A desde 00 con q(0)=q0q(0) = q_0: q=q0−Azq = q_0 - Az. Entonces dT/dz=(q0−Az)/k\mathrm dT/\mathrm dz = (q_0 - Az)/k; integrando con T(0)=T0T(0) = T_0 se obtiene la parábola. Debajo de hh, dq/dz=0\mathrm dq/\mathrm dz = 0 y qq conserva su valor en hh. ∎

Teorema 10.10 (Producción de calor exponencial y relación lineal del flujo de calor)

Sea la producción de calor que disminuye con la profundidad como A(z)=A0e−z/DA(z) = A_0 e^{-z/D}, el flujo de calor tiende a qmq_m a gran profundidad, y la superficie está a T0T_0. Entonces

q(z)=qm+A0D e−z/D,T(z)=T0+qmzk+A0D2k(1−e−z/D),q(z) = q_m + A_0 D\,e^{-z/D}, \qquad T(z) = T_0 + \frac{q_m z}{k} + \frac{A_0 D^2}{k}\bigl(1 - e^{-z/D}\bigr),

y el flujo de calor en la superficie es

q0=qm+A0D.q_0 = q_m + A_0 D .

Demostración. Por Teorema 10.8, dq/dz=−A0e−z/D\mathrm dq/\mathrm dz = -A_0e^{-z/D}, así que q(z)=A0De−z/D+Cq(z) = A_0 D e^{-z/D} + C; cuando z→∞z \to \infty, q→Cq \to C, luego C=qmC = q_m. En z=0z = 0, q0=qm+A0Dq_0 = q_m + A_0 D. Entonces dT/dz=q/k=qm/k+(A0D/k) e−z/D\mathrm dT/\mathrm dz = q/k = q_m/k + (A_0D/k)\,e^{-z/D}, e integrando desde 00 con T(0)=T0T(0) = T_0, T=T0+qmz/k+(A0D2/k)(1−e−z/D)T = T_0 + q_m z/k + (A_0 D^2/k)(1 - e^{-z/D}). ∎

La última línea es una predicción que puede contrastarse en la superficie. Los plutones graníticos de una región comparten el mismo flujo profundo qmq_m y, si solo difieren en lo ricos que son en uranio, torio y potasio, la misma escala DD; entonces el flujo de calor medido sobre cada plutón debe ser una función lineal de la producción de calor medida en sus rocas, con pendiente DD, una longitud, y ordenada en el origen qmq_m, el flujo reducido, que procede de abajo y no del radiogénico. Francis Birch observó por primera vez una recta así, publicada en 1968, y el mismo año se definieron tres provincias por sus rectas:

provinciaordenada en el origen qmq_m (mW m−2\mathrm{mW}\,\mathrm{m}^{-2})pendiente DD (km\mathrm{km})
Estados Unidos orientales337.5
provincia de Basin and Range599.4
Sierra Nevada1710

Las ordenadas en el origen, publicadas como 0.8, 1.4 y 0.4×10−6 cal cm−2 s−10.4\times10^{-6}\,\mathrm{cal}\,\mathrm{cm}^{-2}\,\mathrm{s}^{-1}, se convierten con 1 µcal cm−2 s−1=41.84 mW m−21\,\text{µ}\mathrm{cal}\,\mathrm{cm}^{-2}\,\mathrm{s}^{-1} = 41.84\,\mathrm{mW}\,\mathrm{m}^{-2} (Ejercicio 10.4). Las pendientes varían poco, de 7 a 10 km10\,\mathrm{km}; el flujo de calor reducido varía mucho, y sus mapas muestran dónde el manto de debajo está caliente, como bajo la corteza adelgazada de Basin and Range. La capa uniforme de Proposición 10.9 da también una recta, de pendiente hh: ambos modelos ajustan los datos, pero solo el exponencial sigue siendo una recta de la misma pendiente cuando los plutones se erosionan hasta profundidades distintas (Ejercicio 10.10), por lo que se prefiere.

Geotermos conductivos estacionarios de un continente con producción de calor A_0 = 2.5\, µ W\, m-3 que disminuye sobre D = 10\, km, k = 2.5\, W\, m-1\, K-1, para tres flujos de calor en superficie (flujos reducidos de 15, 30 y 45\, mW\, m-2). Cada uno sigue su curva conductiva (discontinua más allá) hasta la adiabática del manto, a la que llega en el punto: la base de la litosfera térmica. El geotermo del fondo oceánico es el de un semiespacio enfriado durante 60\, Myr. La conductividad uniforme es una simplificación.
Geotermos conductivos estacionarios de un continente con producción de calor A0=2.5 µW m−3A_0 = 2.5\,\text{µ}\mathrm{W}\,\mathrm{m}^{-3} que disminuye sobre D=10 kmD = 10\,\mathrm{km}, k=2.5 W m−1 K−1k = 2.5\,\mathrm{W}\,\mathrm{m}^{-1}\,\mathrm{K}^{-1}, para tres flujos de calor en superficie (flujos reducidos de 15, 30 y 45 mW m−245\,\mathrm{mW}\,\mathrm{m}^{-2}). Cada uno sigue su curva conductiva (discontinua más allá) hasta la adiabática del manto, a la que llega en el punto: la base de la litosfera térmica. El geotermo del fondo oceánico es el de un semiespacio enfriado durante 60 Myr60\,\mathrm{Myr}. La conductividad uniforme es una simplificación.

Método 10.11 (La temperatura a profundidad bajo un plutón)

  1. Mídase el flujo de calor en superficie q0q_0 sobre el plutón (método del sondeo) y su producción de calor A0A_0 a partir de sus contenidos en U, Th y K.
  2. A partir de la recta de la provincia, o con un valor de DD de unos 10 km10\,\mathrm{km}, calcúlese qm=q0−A0Dq_m = q_0 - A_0 D.
  3. Aplíquese Teorema 10.10. Por debajo de unas pocas veces DD, e−z/De^{-z/D} es despreciable y el geotermo es la recta T0+A0D2/k+qmz/kT_0 + A_0D^2/k + q_m z/k: la capa radiogénica ha elevado toda temperatura profunda en A0D2/kA_0D^2/k.

En la figura, la temperatura en el Moho va de unos 320∘C320{}^{\circ}\mathrm{C} para el flujo más bajo a más de 700∘C700{}^{\circ}\mathrm{C} para el más alto: el flujo de calor en superficie lleva una lectura directa de la temperatura de la corteza profunda, que gobierna su resistencia (Capítulo 16).

10.4 Conducción frente a convección

La conducción es lenta. Lo lenta que es se deduce de la ecuación del calor, que los libros de física deducen de la ley de Fourier y del balance de Teorema 10.8 con el almacenamiento añadido: en una dimensión, sin producción de calor, ∂T/∂t=κ ∂2T/∂z2\partial T/\partial t = \kappa\,\partial^2 T/\partial z^2.

Proposición 10.12 (La escala de tiempo conductiva)

Una variación de temperatura impuesta en la frontera de un cuerpo conductor se extiende una distancia LL en un tiempo del orden de

τ=L2κ.\tau = \frac{L^2}{\kappa} .

Demostración. Escríbase z=Lz∗z = L z^*, t=τt∗t = \tau t^*, con z∗z^* y t∗t^* adimensionales. La ecuación del calor queda ∂T/∂t∗=(κτ/L2) ∂2T/∂z∗2\partial T/\partial t^* = (\kappa\tau/L^2)\,\partial^2 T/\partial z^{*2}. Sus soluciones varían en la escala z∗∼1z^* \sim 1 cuando t∗∼1t^* \sim 1 solo si el coeficiente es de orden uno, es decir, si τ∼L2/κ\tau \sim L^2/\kappa. ∎

Las rocas tienen κ≈10−6 m2 s−1\kappa \approx 10^{-6}\,\mathrm{m}^{2}\,\mathrm{s}^{-1}. El ciclo térmico diario llega a unas decenas de centímetros bajo el suelo, el anual a unos metros, y el calentamiento desde la última glaciación a un kilómetro más o menos. En una litosfera de 100 km100\,\mathrm{km} de espesor, τ=1016 s\tau = 10^{16}\,\mathrm{s}, unos 300 Myr300\,\mathrm{Myr}. A través del manto, de 2900 km2900\,\mathrm{km} de espesor, τ=8.4×1018 s\tau = 8.4 \times 10^{18}\,\mathrm{s}, o 270 Gyr270\,\mathrm{Gyr}: sesenta veces la edad de la Tierra, 4.54 Gyr4.54\,\mathrm{Gyr}. La conducción sola no podría haber llevado el calor del manto profundo a ningún punto cercano a la superficie. El manto convecta, y en su interior el calor lo transporta el movimiento de la propia roca, a las velocidades de las placas, centímetros por año.

Definición 10.13 (Adiabática del manto)

La adiabática del manto es el geotermo del interior bien mezclado del manto en convección: la temperatura que toma la roca del manto al subir o bajar con la rapidez suficiente para no intercambiar calor con su entorno, enfriándose por expansión al subir y calentándose por compresión al bajar.

La termodinámica da su gradiente, dT/dz=αgT/cp\mathrm dT/\mathrm dz = \alpha g T/c_p, con α\alpha la expansividad térmica y cpc_p el calor específico (física, admitido aquí). Es pequeño. Una evaluación reciente de la física mineral da 0.54 K km−10.54\,\mathrm{K}\,\mathrm{km}^{-1} a 50 km50\,\mathrm{km} de profundidad, que baja hasta 0.36 K km−10.36\,\mathrm{K}\,\mathrm{km}^{-1} a 410 km410\,\mathrm{km}, y luego de 0.41 a 0.23 K km−10.23\,\mathrm{K}\,\mathrm{km}^{-1} en el manto inferior, con 1646 K1646\,\mathrm{K} a 50 km50\,\mathrm{km} y 2587 K2587\,\mathrm{K} a 2800 km2800\,\mathrm{km}.

Extrapolada a la superficie, la adiabática da la temperatura potencial del manto, 1619 K1619\,\mathrm{K} o unos 1350∘C1350{}^{\circ}\mathrm{C}: la temperatura que tendría el manto si se subiera a la superficie sin intercambiar calor. La noción de temperatura potencial se define, para la atmósfera y el océano, en Capítulo 26. Una temperatura que aumenta solo 0.5 K0.5\,\mathrm{K} por kilómetro, mientras el geotermo conductivo cerca de la superficie sube 25 K25\,\mathrm{K}, es la firma de la convección.

Historia — John Perry y la Tierra que se enfría

En 1862–1863 William Thomson, después Lord Kelvin, calculó la edad de la Tierra a partir de su enfriamiento. La trató como un sólido que se había enfriado por conducción desde una temperatura uniforme de unos 3900∘C3900{}^{\circ}\mathrm{C}; con el gradiente medio de su época, unos 36 K km−136\,\mathrm{K}\,\mathrm{km}^{-1}, que eligió a partir de las medidas de entonces, y una difusividad de 1.2×10−6 m2 s−11.2 \times 10^{-6}\,\mathrm{m}^{2}\,\mathrm{s}^{-1}, el tiempo necesario para alcanzar el gradiente actual fue de unos 96 Myr96\,\mathrm{Myr}, entre cotas de 24 y 400 Myr400\,\mathrm{Myr} (Ejercicio 10.9). Los geólogos lo juzgaron demasiado corto. En 1895, un año antes del descubrimiento de la radiactividad, su antiguo ayudante John Perry señaló en cambio una suposición: si el interior era en parte fluido y convectaba, transportaría el calor como si tuviera una conductividad enorme, y todo el interior, no solo una cáscara exterior delgada, alimentaría el flujo de calor en superficie. Su argumento fue ignorado y olvidado; la convección del manto, establecida mucho después, le dio la razón: el manto profundo se enfría por convección, y eso es el enfriamiento secular del balance.

Proposición 10.14 (La Tierra bien agitada de Perry)

Sea una esfera de radio RR con un interior bien agitado a la temperatura uniforme Ti(t)T_i(t), con capacidad calorífica ρc\rho c por unidad de volumen, bajo una tapa conductora de espesor L≪RL \ll R y conductividad kk cuya superficie está a 0∘C0{}^{\circ}\mathrm{C}, y sin producción de calor. Entonces

Ti(t)=Ti(0) e−t/τP,τP=RL3κ,T_i(t) = T_i(0)\,e^{-t/\tau_P}, \qquad \tau_P = \frac{RL}{3\kappa},

y el flujo de calor en superficie kTi/LkT_i/L decae con la misma constante de tiempo.

Demostración. Como la tapa es delgada, está en estado estacionario y transporta el flujo de calor q=kTi/Lq = kT_i/L. El interior pierde su calor a través del área 4πR24\pi R^2: 43πR3ρc dTi/dt=−4πR2kTi/L\tfrac43\pi R^3\rho c\,\mathrm dT_i/\mathrm dt = -4\pi R^2 kT_i/L, que es dTi/dt=−3κTi/(RL)\mathrm dT_i/\mathrm dt = -3\kappa T_i/(RL) con κ=k/(ρc)\kappa = k/(\rho c), cuya solución es la exponencial de τP=RL/(3κ)\tau_P = RL/(3\kappa). ∎

Para R=6400 kmR = 6400\,\mathrm{km}, L=100 kmL = 100\,\mathrm{km} y κ=10−6 m2 s−1\kappa = 10^{-6}\,\mathrm{m}^{2}\,\mathrm{s}^{-1}, τP=6.8 Gyr\tau_P = 6.8\,\mathrm{Gyr}: una Tierra que convecta mantiene un flujo de calor alto durante miles de millones de años, donde la Tierra conductora de Kelvin lo perdió en decenas de millones (Ejercicio 10.11).

10.5 La litosfera como capa límite térmica

Definición 10.15 (Capa límite térmica, litosfera térmica)

En un fluido que convecta, una capa límite térmica es una capa a lo largo de un contorno en la que el calor pasa por conducción y la temperatura cambia del valor del contorno al del interior bien mezclado. La litosfera térmica es la capa límite térmica superior del manto: la cáscara exterior de la Tierra, desde la superficie hasta la profundidad a la que el geotermo conductivo se encuentra con la adiabática del manto.

Las placas son la parte fría superior de la convección del manto, su capa límite superior. Su base, donde el geotermo conductivo empinado se curva hacia la adiabática suave, está a 70 a 230 km70\text{ a }230\,\mathrm{km} en la figura de los geotermos para flujos de calor continentales de 70 a 40 mW m−240\,\mathrm{mW}\,\mathrm{m}^{-2}, en el rango medido por sismología (Capítulo 9): la litosfera definida por su resistencia y la definida por su temperatura son aproximadamente la misma cáscara, porque la resistencia depende de la temperatura. En el fondo del manto, la capa D′′'' es la capa límite inferior, calentada desde abajo por el núcleo.

Bajo los océanos, la capa límite nace en la dorsal, donde el manto caliente sube a la superficie, y se engrosa a medida que la placa se aleja y se enfría. El enfriamiento de un semiespacio desde una temperatura uniforme, resuelto por Fourier y usado por Kelvin para la Tierra entera, lo describe bien (Capítulo 13 y el libro de tercer año universitario lo deducen): las isotermas se hunden como la raíz cuadrada de la edad del fondo (figura de abajo), y el flujo de calor disminuye como

q(t)=Ct,C≈490 mW m−2 Myr1/2,q(t) = \frac{C}{\sqrt t}, \qquad C \approx 490\,\mathrm{mW}\,\mathrm{m}^{-2}\,\mathrm{Myr}^{1/2},

el valor empleado en el balance global. A 100 Myr100\,\mathrm{Myr} el semiespacio se ha enfriado hasta unos 130 km130\,\mathrm{km} (Ejercicio 10.12). La placa real deja de engrosarse antes: las ondas sísmicas que cruzan el Pacífico se ajustan mejor con una placa enfriada que tiende a un espesor de 85 a 95 km85\text{ a }95\,\mathrm{km}, el aplanamiento del fondo oceánico antiguo que se vuelve a encontrar en Capítulo 13.

El balance térmico, el flujo de calor y la estructura térmica de las placas se reúnen en los capítulos siguientes: la placa enfriada se hunde (Capítulo 13), y la litosfera fría que se hunde en las zonas de subducción es la corriente descendente de la convección a la que pertenece.

Isotermas en una placa oceánica que se enfría: un semiespacio inicialmente a 1350 C, enfriado desde su superficie, con = 10-6\, m2\, s-1. Cada isoterma se hunde como la raíz cuadrada de la edad; la capa sombreada sobre la isoterma de 1200 C es la litosfera térmica. Las flechas son esquemáticas.
Isotermas en una placa oceánica que se enfría: un semiespacio inicialmente a 1350∘C1350{}^{\circ}\mathrm{C}, enfriado desde su superficie, con κ=10−6 m2 s−1\kappa = 10^{-6}\,\mathrm{m}^{2}\,\mathrm{s}^{-1}. Cada isoterma se hunde como la raíz cuadrada de la edad; la capa sombreada sobre la isoterma de 1200∘C1200{}^{\circ}\mathrm{C} es la litosfera térmica. Las flechas son esquemáticas.

10.6 Ejercicios

Ejercicio 10.1 ★

Un sondeo en granito de conductividad 2.5 W m−1 K−12.5\,\mathrm{W}\,\mathrm{m}^{-1}\,\mathrm{K}^{-1} muestra un gradiente estacionario de 22 K km−122\,\mathrm{K}\,\mathrm{km}^{-1}. Calcúlese el flujo de calor y la potencia que atraviesa un kilómetro cuadrado de superficie en ese punto.

Solución

Solución de Ejercicio 10.1.

q=k dT/dz=2.5×0.022 K m−1=0.055 W m−2=55 mW m−2q = k\,\mathrm dT/\mathrm dz = 2.5\times0.022\,\mathrm{K}\,\mathrm{m}^{-1} = 0.055\,\mathrm{W}\,\mathrm{m}^{-2} = 55\,\mathrm{mW}\,\mathrm{m}^{-2}. A través de 1 km21\,\mathrm{km}^{2} =106 m2= 10^{6}\,\mathrm{m}^{2}: 5.5×104 W5.5 \times 10^{4}\,\mathrm{W}, es decir 55 kW55\,\mathrm{kW}.

Ejercicio 10.2 ★

A partir de las áreas y los flujos de calor medios de los continentes (2.073×1014 m22.073 \times 10^{14}\,\mathrm{m}^{2}, 70.9 mW m−270.9\,\mathrm{mW}\,\mathrm{m}^{-2}) y de los océanos (3.028×1014 m23.028 \times 10^{14}\,\mathrm{m}^{2}, 105.4 mW m−2105.4\,\mathrm{mW}\,\mathrm{m}^{-2}), calcúlese la pérdida de calor de cada uno, la total y el flujo de calor medio sobre la superficie de la Tierra (5.101×1014 m25.101 \times 10^{14}\,\mathrm{m}^{2}). ¿Qué fracción de los 173 000 TW173\,000\,\mathrm{TW} que recibe del Sol representa?

Solución

Solución de Ejercicio 10.2.

Continentes 2.073×1014×0.0709=1.470×1013 W=14.70 TW2.073 \times 10^{14}\times0.0709 = 1.470 \times 10^{13}\,\mathrm{W} = 14.70\,\mathrm{TW}; océanos 3.028×1014×0.1054=31.92 TW3.028 \times 10^{14}\times0.1054 = 31.92\,\mathrm{TW}; total 46.61 TW46.61\,\mathrm{TW}, y el flujo de calor medio 4.661×1013/5.101×1014=91.4 mW m−24.661 \times 10^{13}/5.101 \times 10^{14} = 91.4\,\mathrm{mW}\,\mathrm{m}^{-2} (los valores publicados, 46.7 y 91.6, difieren solo por el redondeo de las dos medias). La fracción de la entrada solar es 46.6/173000=2.7×10−446.6/173000 = 2.7 \times 10^{-4}, o 0.027 %0.027\,\%.

Ejercicio 10.3 ★

A partir de la tabla de los cuatro isótopos, calcúlese el calor liberado por kilogramo de uranio natural (99.274 %99.274\,\% X238X22238U\ce{^{238}U} y 0.720 %0.720\,\% X235X22235U\ce{^{235}U} de los átomos; sus masas difieren en unos 1 %1\,\%, así que tómense estas como fracciones másicas) y por kilogramo de potasio natural, cuyo X40X2240K\ce{^{40}K} constituye el 0.0117 %0.0117\,\% de los átomos (tómese la masa de un átomo de X40X2240K\ce{^{40}K} como 39.96/39.1039.96/39.10 veces la masa media de un átomo de potasio).

Solución

Solución de Ejercicio 10.3.

Uranio: 0.99274×9.513×10−5+0.00720×5.685×10−4=9.444×10−5+4.09×10−6=9.85×10−5 W kg−10.99274\times9.513 \times 10^{-5} + 0.00720\times5.685 \times 10^{-4} = 9.444 \times 10^{-5} + 4.09 \times 10^{-6} = 9.85 \times 10^{-5}\,\mathrm{W}\,\mathrm{kg}^{-1}; el raro X235X22235U\ce{^{235}U} aporta el 4 %4\,\% de ella. Potasio: la fracción másica de X40X2240K\ce{^{40}K} es 1.17×10−4×39.96/39.10=1.196×10−41.17 \times 10^{-4}\times39.96/39.10 = 1.196 \times 10^{-4}, así que 1.196×10−4×2.847×10−5=3.40×10−9 W kg−11.196 \times 10^{-4}\times2.847 \times 10^{-5} = 3.40 \times 10^{-9}\,\mathrm{W}\,\mathrm{kg}^{-1}.

Ejercicio 10.4 ★

Conviértanse las ordenadas en el origen de las tres provincias, 0.8, 1.4 y 0.4×10−6 cal cm−2 s−10.4\times10^{-6}\,\mathrm{cal}\,\mathrm{cm}^{-2}\,\mathrm{s}^{-1}, a mW m−2\mathrm{mW}\,\mathrm{m}^{-2} (1 cal=4.184 J1\,\mathrm{cal} = 4.184\,\mathrm{J}). Prevéase el flujo de calor en superficie sobre un granito de producción de calor 3.0 µW m−33.0\,\text{µ}\mathrm{W}\,\mathrm{m}^{-3} en cada provincia.

Solución

Solución de Ejercicio 10.4.

Conversión: 1 µcal cm−2 s−1=4.184×10−6 J s−1/10−4 m2=41.84 mW m−21\,\text{µ}\mathrm{cal}\,\mathrm{cm}^{-2}\,\mathrm{s}^{-1} = 4.184 \times 10^{-6}\,\mathrm{J}\,\mathrm{s}^{-1}/10^{-4}\,\mathrm{m}^{2} = 41.84\,\mathrm{mW}\,\mathrm{m}^{-2}: Estados Unidos orientales 33.5 mW m−233.5\,\mathrm{mW}\,\mathrm{m}^{-2}, Basin and Range 58.6 mW m−258.6\,\mathrm{mW}\,\mathrm{m}^{-2}, Sierra Nevada 16.7 mW m−216.7\,\mathrm{mW}\,\mathrm{m}^{-2}. Con A0DA_0D en mW m−2\mathrm{mW}\,\mathrm{m}^{-2} cuando A0A_0 está en µW m−3\text{µ}\mathrm{W}\,\mathrm{m}^{-3} y DD en km\mathrm{km}: 33.5+3.0×7.5=56.0 mW m−233.5 + 3.0\times7.5 = 56.0\,\mathrm{mW}\,\mathrm{m}^{-2}, 58.6+3.0×9.4=86.8 mW m−258.6 + 3.0\times9.4 = 86.8\,\mathrm{mW}\,\mathrm{m}^{-2}, 16.7+3.0×10=46.7 mW m−216.7 + 3.0\times10 = 46.7\,\mathrm{mW}\,\mathrm{m}^{-2}. El mismo granito da flujos de calor casi el doble de altos en Basin and Range que en la Sierra Nevada: la diferencia está en el manto de abajo.

Ejercicio 10.5 ★★

Una corteza tiene una producción de calor uniforme de 2.0 µW m−32.0\,\text{µ}\mathrm{W}\,\mathrm{m}^{-3} en sus 10 km10\,\mathrm{km} superiores y ninguna por debajo, una conductividad de 2.5 W m−1 K−12.5\,\mathrm{W}\,\mathrm{m}^{-1}\,\mathrm{K}^{-1}, un flujo de calor en superficie de 60 mW m−260\,\mathrm{mW}\,\mathrm{m}^{-2} y una temperatura en superficie de 10∘C10{}^{\circ}\mathrm{C}. Calcúlense la temperatura a 10 km10\,\mathrm{km}, el flujo de calor bajo la capa radiogénica y la temperatura en el Moho, a 35 km35\,\mathrm{km} de profundidad.

Solución

Solución de Ejercicio 10.5.

T(10 km)=10+0.060×1042.5−2×10−6×1082×2.5=10+240−40=210∘CT(10\,\text{km}) = 10 + \frac{0.060\times10^4}{2.5} - \frac{2 \times 10^{-6}\times10^8}{2\times2.5} = 10 + 240 - 40 = 210{}^{\circ}\mathrm{C}. Debajo, qm=60−2.0×10=40 mW m−2q_m = 60 - 2.0\times10 = 40\,\mathrm{mW}\,\mathrm{m}^{-2}, un gradiente de 16 K km−116\,\mathrm{K}\,\mathrm{km}^{-1}, y T(35 km)=210+16×25=610∘CT(35\,\text{km}) = 210 + 16\times25 = 610{}^{\circ}\mathrm{C}.

Ejercicio 10.6 ★★

Un sondeo atraviesa 300 m300\,\mathrm{m} de lutita (k=1.5 W m−1 K−1k = 1.5\,\mathrm{W}\,\mathrm{m}^{-1}\,\mathrm{K}^{-1}), y luego granito (k=3.0 W m−1 K−1k = 3.0\,\mathrm{W}\,\mathrm{m}^{-1}\,\mathrm{K}^{-1}). Las temperaturas de equilibrio son 12.0∘C12.0{}^{\circ}\mathrm{C} en superficie, 24.0∘C24.0{}^{\circ}\mathrm{C} a 300 m300\,\mathrm{m} y 34.0∘C34.0{}^{\circ}\mathrm{C} a 800 m800\,\mathrm{m}. Calcúlense el gradiente y el flujo de calor en cada roca, la resistencia térmica a 300 m300\,\mathrm{m} y a 800 m800\,\mathrm{m}, y la pendiente del diagrama de Bullard. ¿Es constante el flujo de calor?

Solución

Solución de Ejercicio 10.6.

Lutita: gradiente 12.0/0.300=40 K km−112.0/0.300 = 40\,\mathrm{K}\,\mathrm{km}^{-1}, q=1.5×0.040=60 mW m−2q = 1.5\times0.040 = 60\,\mathrm{mW}\,\mathrm{m}^{-2}. Granito: 10.0/0.500=20 K km−110.0/0.500 = 20\,\mathrm{K}\,\mathrm{km}^{-1}, q=3.0×0.020=60 mW m−2q = 3.0\times0.020 = 60\,\mathrm{mW}\,\mathrm{m}^{-2}. R(300)=300/1.5=200 m2 K W−1R(300) = 300/1.5 = 200\,\mathrm{m}^{2}\,\mathrm{K}\,\mathrm{W}^{-1}, R(800)=200+500/3.0=366.7 m2 K W−1R(800) = 200 + 500/3.0 = 366.7\,\mathrm{m}^{2}\,\mathrm{K}\,\mathrm{W}^{-1}. Los tres puntos (0,12.0)(0, 12.0), (200,24.0)(200, 24.0), (366.7,34.0)(366.7, 34.0) están alineados, con pendiente 22.0/366.7=0.0600 W m−222.0/366.7 = 0.0600\,\mathrm{W}\,\mathrm{m}^{-2}: el flujo de calor es el mismo, 60 mW m−260\,\mathrm{mW}\,\mathrm{m}^{-2}, en ambas rocas; solo cambia el gradiente, en proporción inversa a la conductividad.

Ejercicio 10.7 ★★

Con κ=10−6 m2 s−1\kappa = 10^{-6}\,\mathrm{m}^{2}\,\mathrm{s}^{-1}, calcúlese la escala de tiempo conductiva para L=1 mL = 1\,\mathrm{m}, 10 m10\,\mathrm{m}, 100 km100\,\mathrm{km} y 2900 km2900\,\mathrm{km}, cada una en una unidad adecuada. ¿Qué ciclos de temperatura en superficie se siguen sintiendo 10 m10\,\mathrm{m} bajo el suelo? Compárese el último tiempo con la edad de la Tierra.

Solución

Solución de Ejercicio 10.7.

τ=L2/κ\tau = L^2/\kappa: 1 m1\,\mathrm{m}, 106 s10^{6}\,\mathrm{s} o 1212 días; 10 m10\,\mathrm{m}, 108 s10^{8}\,\mathrm{s} o 3.23.2 años; 100 km100\,\mathrm{km}, 1016 s10^{16}\,\mathrm{s} o 317 Myr317\,\mathrm{Myr}; 2900 km2900\,\mathrm{km}, 8.4×1018 s8.4 \times 10^{18}\,\mathrm{s} o 267 Gyr267\,\mathrm{Gyr}. A 10 m10\,\mathrm{m} el ciclo diario no se siente en absoluto, y el anual, más corto que los tres años necesarios para atravesar 10 m10\,\mathrm{m}, solo muy débilmente: solo llegan los cambios que duran años o más. El último tiempo es 267/4.54=59267/4.54 = 59 veces la edad de la Tierra.

Ejercicio 10.8 ★★

La adiabática del manto está a 1646 K1646\,\mathrm{K} a 50 km50\,\mathrm{km} de profundidad, con un gradiente de 0.54 K km−10.54\,\mathrm{K}\,\mathrm{km}^{-1} allí, que disminuye hasta 0.36 K km−10.36\,\mathrm{K}\,\mathrm{km}^{-1} a 410 km410\,\mathrm{km}. (a) Extrapólese hasta la superficie con el gradiente a 50 km50\,\mathrm{km}, y dése la temperatura potencial en kelvin y en grados Celsius. (b) Estímese la temperatura a 410 km410\,\mathrm{km} con la media de los dos gradientes, y compárese con los 1799 K1799\,\mathrm{K} publicados justo encima de la discontinuidad.

Solución

Solución de Ejercicio 10.8.

(a) 1646−0.54×50=1619 K1646 - 0.54\times50 = 1619\,\mathrm{K}, es decir 1346∘C1346{}^{\circ}\mathrm{C}, el 1350∘C1350{}^{\circ}\mathrm{C} publicado, redondeado. (b) Gradiente medio (0.54+0.36)/2=0.45 K km−1(0.54 + 0.36)/2 = 0.45\,\mathrm{K}\,\mathrm{km}^{-1}, así que 1646+0.45×360=1808 K1646 + 0.45\times360 = 1808\,\mathrm{K}, 9 K9\,\mathrm{K} (0.5 %0.5\,\%) por encima de 1799 K1799\,\mathrm{K}: el gradiente no disminuye linealmente, pero la estimación es cercana.

Ejercicio 10.9 ★★

Un semiespacio enfriado desde la temperatura uniforme T0T_0, con su superficie mantenida a 0∘C0{}^{\circ}\mathrm{C}, tiene tras un tiempo tt el gradiente superficial G=T0/πκtG = T_0/\sqrt{\pi\kappa t} (admitido). (a) Exprésese tt en función de GG. (b) Con los valores de Kelvin, T0=3900∘CT_0 = 3900{}^{\circ}\mathrm{C}, κ=1.2×10−6 m2 s−1\kappa = 1.2 \times 10^{-6}\,\mathrm{m}^{2}\,\mathrm{s}^{-1} y G=1/50G = 1/50 de grado Fahrenheit por pie; calcúlese GG en K km−1\mathrm{K}\,\mathrm{km}^{-1} y la edad. (c) Rehágase con T0=1300∘CT_0 = 1300{}^{\circ}\mathrm{C}, κ=10−6 m2 s−1\kappa = 10^{-6}\,\mathrm{m}^{2}\,\mathrm{s}^{-1} y G=36 K km−1G = 36\,\mathrm{K}\,\mathrm{km}^{-1}.

Solución

Solución de Ejercicio 10.9.

(a) G2=T02/(πκt)G^2 = T_0^2/(\pi\kappa t) da t=T02/(πκG2)t = T_0^2/(\pi\kappa G^2). (b) G=0.02×59 K/0.3048 m=0.036 45 K m−1=36.45 K km−1G = 0.02\times\tfrac59\,\mathrm{K}/0.3048\,\mathrm{m} = 0.036\,45\,\mathrm{K}\,\mathrm{m}^{-1} = 36.45\,\mathrm{K}\,\mathrm{km}^{-1}; t=39002/(π×1.2×10−6×0.036452)=1.521×107/5.009×10−9=3.04×1015 st = 3900^2/(\pi\times1.2 \times 10^{-6}\times0.03645^2) = 1.521 \times 10^{7}/5.009 \times 10^{-9} = 3.04 \times 10^{15}\,\mathrm{s}, es decir 96 Myr96\,\mathrm{Myr}, la cifra de Kelvin. (c) t=13002/(π×10−6×0.0362)=4.15×1014 s=13 Myrt = 1300^2/(\pi\times10^{-6}\times0.036^2) = 4.15 \times 10^{14}\,\mathrm{s} = 13\,\mathrm{Myr}: una temperatura inicial más realista hace que la Tierra conductora sea aún más joven.

Ejercicio 10.10 ★★★

Una provincia se erosiona: se retira un espesor ee de la parte superior de cada plutón, siendo inalterado el flujo profundo qmq_m. (a) En el modelo exponencial, muéstrese que la nueva producción superficial de calor y el nuevo flujo de calor, A0′A_0' y q0′q_0', siguen cumpliendo q0′=qm+A0′Dq_0' = q_m + A_0'D. (b) En el modelo de capa uniforme de espesor hh, muéstrese que un plutón erosionado en ee se sitúa en una recta de pendiente h−eh - e. (c) Plutones de una misma provincia se erosionaron en cantidades distintas. ¿Qué modelo predice que siguen sobre una sola recta?

Solución

Solución de Ejercicio 10.10.

(a) Tras la erosión, la nueva superficie está a la antigua profundidad ee, donde la producción de calor era A0′=A0e−e/DA_0' = A_0e^{-e/D} y el flujo de calor q(e)=qm+A0De−e/D=qm+A0′Dq(e) = q_m + A_0De^{-e/D} = q_m + A_0'D; medida desde la nueva superficie, la producción sigue disminuyendo como A0′e−z′/DA_0'e^{-z'/D}. El plutón se desplaza hacia abajo a lo largo de la misma recta. (b) La capa tiene ahora un espesor h−eh - e con la misma AA, así que q0′=qm+A(h−e)q_0' = q_m + A(h - e): sobre una recta de pendiente h−eh - e que pasa por qmq_m. (c) Solo el modelo exponencial: en el modelo de capa uniforme, los plutones erosionados a profundidades distintas se situarían sobre rectas de pendientes diferentes, y los puntos de una provincia se dispersarían. Que se observe una sola recta favorece una producción que disminuye aproximadamente de forma exponencial con la profundidad.

Ejercicio 10.11 ★★★

En el modelo de Perry (Proposición 10.14) tómese R=6400 kmR = 6400\,\mathrm{km}, κ=10−6 m2 s−1\kappa = 10^{-6}\,\mathrm{m}^{2}\,\mathrm{s}^{-1}, k=3 W m−1 K−1k = 3\,\mathrm{W}\,\mathrm{m}^{-1}\,\mathrm{K}^{-1} y una temperatura interior inicial de 1300∘C1300{}^{\circ}\mathrm{C}. (a) Calcúlese τP\tau_P para tapas de 50 y 100 km100\,\mathrm{km}. (b) Para la tapa de 50 km50\,\mathrm{km}, calcúlese el flujo de calor en superficie al principio y tras 2 Gyr2\,\mathrm{Gyr}. (c) El semiespacio conductor de Kelvin, con el mismo kk y T0T_0, tiene un flujo de calor en superficie kT0/πκtkT_0/\sqrt{\pi\kappa t}: calcúlese tras 2 Gyr2\,\mathrm{Gyr} y compárese.

Solución

Solución de Ejercicio 10.11.

(a) τP=RL/(3κ)\tau_P = RL/(3\kappa): para 50 km50\,\mathrm{km}, 6.4×106×5×104/3×10−6=1.07×1017 s=3.4 Gyr6.4 \times 10^{6}\times5 \times 10^{4}/3 \times 10^{-6} = 1.07 \times 10^{17}\,\mathrm{s} = 3.4\,\mathrm{Gyr}; para 100 km100\,\mathrm{km}, 6.8 Gyr6.8\,\mathrm{Gyr}. (b) q=kTi/L=3×1300/5×104=78 mW m−2q = kT_i/L = 3\times1300/5 \times 10^{4} = 78\,\mathrm{mW}\,\mathrm{m}^{-2} al principio, y 78 e−2/3.38=43 mW m−278\,e^{-2/3.38} = 43\,\mathrm{mW}\,\mathrm{m}^{-2} tras 2 Gyr2\,\mathrm{Gyr}. (c) 3×1300/π×10−6×6.31×1016=3900/4.45×105=8.8 mW m−23\times1300/\sqrt{\pi\times10^{-6}\times6.31 \times 10^{16}} = 3900/4.45 \times 10^{5} = 8.8\,\mathrm{mW}\,\mathrm{m}^{-2}: la Tierra bien agitada sigue perdiendo cinco veces más calor tras dos mil millones de años. Todo su interior alimenta la superficie, mientras que en una Tierra conductora solo lo hace una cáscara exterior.

Ejercicio 10.12 ★★★

El flujo de calor en superficie de un semiespacio que se enfría es q=k ΔT/πκtq = k\,\Delta T/\sqrt{\pi\kappa t} (Ejercicio 10.9). (a) Con k=3 W m−1 K−1k = 3\,\mathrm{W}\,\mathrm{m}^{-1}\,\mathrm{K}^{-1}, κ=10−6 m2 s−1\kappa = 10^{-6}\,\mathrm{m}^{2}\,\mathrm{s}^{-1} y ΔT=1336 K\Delta T = 1336\,\mathrm{K} (desde la temperatura potencial hasta un fondo marino a unos 10∘C10{}^{\circ}\mathrm{C}), calcúlese la constante CC de q=C/tq = C/\sqrt t, con tt en millones de años, y compárese con los 490 mW m−2 Myr1/2490\,\mathrm{mW}\,\mathrm{m}^{-2}\,\mathrm{Myr}^{1/2} ajustados a los datos. (b) Con C=490C = 490, dense el flujo de calor a 1, 10 y 100 Myr100\,\mathrm{Myr}. (c) La isoterma de 1200∘C1200{}^{\circ}\mathrm{C} de la figura de isotermas en enfriamiento se encuentra donde z=2×1.163κtz = 2\times1.163\sqrt{\kappa t}; calcúlese su profundidad a 100 Myr100\,\mathrm{Myr}.

Solución

Solución de Ejercicio 10.12.

(a) πκ×1 Myr=π×10−6×3.156×1013=9957 m\sqrt{\pi\kappa\times1\,\mathrm{Myr}} = \sqrt{\pi\times10^{-6}\times3.156 \times 10^{13}} = 9957\,\mathrm{m}, así que C=3×1336/9957=0.403 W m−2C = 3\times1336/9957 = 0.403\,\mathrm{W}\,\mathrm{m}^{-2} para tt en Myr, 403 mW m−2 Myr1/2403\,\mathrm{mW}\,\mathrm{m}^{-2}\,\mathrm{Myr}^{1/2}: el 82 %82\,\% de los 490 ajustados. El modelo necesita un k ΔTk\,\Delta T mayor (unos 4900 W m−14900\,\mathrm{W}\,\mathrm{m}^{-1}): una conductividad mayor del manto, o un interior más caliente, como a lo largo de la adiabática. (b) 490/1=490490/\sqrt1 = 490, 490/10=155490/\sqrt{10} = 155 y 490/100=49 mW m−2490/\sqrt{100} = 49\,\mathrm{mW}\,\mathrm{m}^{-2}. (c) z=2.32610−6×3.156×1015=2.326×56 180 m=131 kmz = 2.326\sqrt{10^{-6}\times3.156 \times 10^{15}} = 2.326\times56\,180\,\mathrm{m} = 131\,\mathrm{km}.

10.7 Problema: rocas calientes bajo los plutones

Problema 10.1

Problema de fin de semana — flujo de calor y producción de calor de cinco plutones graníticos, la relación lineal y la temperatura de la corteza profunda

Se han perforado cinco plutones graníticos de una misma provincia. En cada sondeo se registró la temperatura de equilibrio a 100 m100\,\mathrm{m} y a 500 m500\,\mathrm{m}, y la conductividad del granito se midió en testigos; sus contenidos de uranio, torio y potasio se analizaron en muestras de superficie. Los cinco granitos tienen una densidad de 2650 kg m−32650\,\mathrm{kg}\,\mathrm{m}^{-3}.

plutónTT a 100 m100\,\mathrm{m} (∘C{}^{\circ}\mathrm{C})TT a 500 m500\,\mathrm{m} (∘C{}^{\circ}\mathrm{C})kk (W m−1 K−1\mathrm{W}\,\mathrm{m}^{-1}\,\mathrm{K}^{-1})U, Th (ppm)K (%\%)
A13.119.22.92.0, 83.0
B12.420.22.63.5, 143.5
C13.621.43.15.0, 184.0
D12.923.12.77.0, 254.2
E14.225.53.010.0, 354.5

Úsense los calores elementales 9.85×10−5 W kg−19.85 \times 10^{-5}\,\mathrm{W}\,\mathrm{kg}^{-1} (U), 2.63×10−5 W kg−12.63 \times 10^{-5}\,\mathrm{W}\,\mathrm{kg}^{-1} (Th) y 3.40×10−9 W kg−13.40 \times 10^{-9}\,\mathrm{W}\,\mathrm{kg}^{-1} (K). Bajo los plutones, tómese k=2.5 W m−1 K−1k = 2.5\,\mathrm{W}\,\mathrm{m}^{-1}\,\mathrm{K}^{-1} para toda la corteza y una temperatura en superficie de 10∘C10{}^{\circ}\mathrm{C}.

Parte I — Flujo de calor.

  1. Calcúlese el gradiente de temperatura en cada sondeo.
  2. ¿Por qué deben registrarse las temperaturas semanas o meses después de la perforación, y por qué se deja fuera la parte superior de 100 m100\,\mathrm{m}?
  3. Calcúlese el flujo de calor en cada plutón.
  4. Compárense con el flujo de calor continental medio, 70.9 mW m−270.9\,\mathrm{mW}\,\mathrm{m}^{-2}.
  5. Los gradientes de B y de C son iguales. ¿Por qué no lo son sus flujos de calor?
  6. Los plutones están bajo colinas. ¿Elevaría o bajaría el relieve el gradiente medido bajo una cima? (Piénsese en la forma de las isotermas.)

Parte II — Producción de calor y la recta.

  1. Calcúlese la producción de calor debida al uranio en cada plutón.
  2. Lo mismo para el torio.
  3. Lo mismo para el potasio.
  4. Súmense para obtener A0A_0 en cada plutón. ¿Qué elemento contribuye menos?
  5. Represéntese q0q_0 frente a A0A_0. ¿Están los puntos cerca de una recta?
  6. Calcúlense las medias de A0A_0 y de q0q_0, y las sumas ∑(A0−Aˉ0)2\sum (A_0 - \bar A_0)^2 y ∑(A0−Aˉ0)(q0−qˉ0)\sum (A_0 - \bar A_0)(q_0 - \bar q_0).
  7. Dedúzcase la recta de mínimos cuadrados q0=qm+DA0q_0 = q_m + D A_0: su pendiente DD (en km\mathrm{km}) y su ordenada en el origen qmq_m.
  8. Compárense DD y qmq_m con las tres provincias del capítulo. ¿A cuál se parece esta provincia?

Parte III — El geotermo bajo el plutón E. Redondéese DD a 10 km10\,\mathrm{km} y qmq_m a 30 mW m−230\,\mathrm{mW}\,\mathrm{m}^{-2}, y tómese para el plutón E el A0=5.46 µW m−3A_0 = 5.46\,\text{µ}\mathrm{W}\,\mathrm{m}^{-3} redondeado.

  1. Escríbase el geotermo T(z)T(z) de Teorema 10.10 con estos valores.
  2. Calcúlese el flujo de calor a 10, 20 y 35 km35\,\mathrm{km}.
  3. ¿Qué fracción del flujo de calor en superficie sobre E procede de bajo la capa radiogénica?
  4. Calcúlese la temperatura a 10 km10\,\mathrm{km}.
  5. Calcúlese la temperatura a 20 km20\,\mathrm{km}.

Parte IV — La corteza profunda.

  1. Calcúlese la temperatura en el Moho, a 35 km35\,\mathrm{km} de profundidad, bajo el plutón A (A0=1.35 µW m−3A_0 = 1.35\,\text{µ}\mathrm{W}\,\mathrm{m}^{-3}).
  2. ¿En cuánto elevan las capas radiogénicas de A y de E las temperaturas profundas? Explíquese la diferencia entre los dos plutones.
  3. Por debajo de 50 km50\,\mathrm{km}, e−z/De^{-z/D} es despreciable. Hállese la profundidad de la isoterma de 1300∘C1300{}^{\circ}\mathrm{C} bajo el plutón E y bajo el plutón A.
  4. Compárense estas profundidades con el espesor de la litosfera continental que se encuentra en Capítulo 9.
  5. ¿Cuál es la temperatura en el Moho bajo el plutón E?
Solución

Solución de Problema 10.1.

1. Gradiente =(T500−T100)/0.4 km= (T_{500} - T_{100})/0.4\,\mathrm{km}: A 15.25 K km−115.25\,\mathrm{K}\,\mathrm{km}^{-1}, B 19.5019.50, C 19.5019.50, D 25.5025.50, E 28.25 K km−128.25\,\mathrm{K}\,\mathrm{km}^{-1}.

2. Perforar altera la temperatura de la roca (el fluido de perforación enfría o calienta las paredes), y la perturbación tarda semanas o meses en desaparecer. Las decenas de metros superiores siguen el ciclo estacional y los cambios recientes de la temperatura en superficie.

3. q0=k×q_0 = k\timesgradiente: A 2.9×15.25=44.225 mW m−22.9\times15.25 = 44.225\,\mathrm{mW}\,\mathrm{m}^{-2}, B 50.70 mW m−250.70\,\mathrm{mW}\,\mathrm{m}^{-2}, C 60.45 mW m−260.45\,\mathrm{mW}\,\mathrm{m}^{-2}, D 68.85 mW m−268.85\,\mathrm{mW}\,\mathrm{m}^{-2}, E 84.75 mW m−284.75\,\mathrm{mW}\,\mathrm{m}^{-2}.

4. De A a D están por debajo de la media continental, y E por encima en 13.9 mW m−213.9\,\mathrm{mW}\,\mathrm{m}^{-2}, o 20 %20\,\%.

5. El flujo de calor es la conductividad por el gradiente: C conduce mejor (3.1 frente a 2.6), así que el mismo gradiente transporta un 19 %19\,\% más de calor.

6. Las isotermas siguen el relieve solo en parte: bajo una cima se elevan menos que la superficie, así que están más separadas por debajo, y el gradiente medido disminuye (aumenta bajo un valle). Hace falta una corrección por relieve antes de comparar los plutones.

7. ρCUHU\rho C_U H_U: A 2650×2.0×10−6×9.85×10−5=0.522 µW m−32650\times2.0 \times 10^{-6}\times9.85 \times 10^{-5} = 0.522\,\text{µ}\mathrm{W}\,\mathrm{m}^{-3}, B 0.914, C 1.305, D 1.827, E 2.610.

8. Torio: A 0.558, B 0.976, C 1.255, D 1.742, E 2.439.

9. Potasio: A 2650×0.030×3.40×10−9=0.2702650\times0.030\times3.40 \times 10^{-9} = 0.270, B 0.315, C 0.360, D 0.378, E 0.405.

10. A0A_0 = 1.35, 2.20, 2.92, 3.95 y 5.46 µW m−35.46\,\text{µ}\mathrm{W}\,\mathrm{m}^{-3}. El potasio aporta lo menos, desde el 20 %20\,\% del total en A hasta 7 %7\,\% en E; el uranio y el torio se reparten el resto casi por igual.

11. Sí: los cinco puntos suben de forma regular, a unos 1.5 mW m−21.5\,\mathrm{mW}\,\mathrm{m}^{-2} de una recta.

12. Aˉ0=15.88/5=3.176\bar A_0 = 15.88/5 = 3.176, qˉ0=308.975/5=61.795\bar q_0 = 308.975/5 = 61.795; ∑(A0−Aˉ0)2=10.1681\sum(A_0 - \bar A_0)^2 = 10.1681; ∑(A0−Aˉ0)(q0−qˉ0)=101.146\sum(A_0 - \bar A_0)(q_0 - \bar q_0) = 101.146.

13. D=101.146/10.1681=9.95D = 101.146/10.1681 = 9.95 (en mW m−2\mathrm{mW}\,\mathrm{m}^{-2} por µW m−3\text{µ}\mathrm{W}\,\mathrm{m}^{-3}, es decir en km): D=9.95 kmD = 9.95\,\mathrm{km}, y qm=61.795−9.95×3.176=30.20 mW m−2q_m = 61.795 - 9.95\times3.176 = 30.20\,\mathrm{mW}\,\mathrm{m}^{-2}.

14. La pendiente, de unos 10 km10\,\mathrm{km}, es la de las tres provincias; el flujo de calor reducido, 30 mW m−230\,\mathrm{mW}\,\mathrm{m}^{-2}, está cerca de los 33 mW m−233\,\mathrm{mW}\,\mathrm{m}^{-2} de los Estados Unidos orientales: una provincia tranquila con un manto frío debajo, a diferencia de Basin and Range.

15. qm/k=12 K km−1q_m/k = 12\,\mathrm{K}\,\mathrm{km}^{-1} y A0D2/k=5.46×10−6×108/2.5=218.4 KA_0D^2/k = 5.46 \times 10^{-6}\times10^{8}/2.5 = 218.4\,\mathrm{K}: T(z)=10+12z+218.4 (1−e−z/10)T(z) = 10 + 12z + 218.4\,(1 - e^{-z/10}), zz en km, TT en ∘C{}^{\circ}\mathrm{C}.

16. q(z)=30+54.6 e−z/10q(z) = 30 + 54.6\,e^{-z/10}: 50.1 mW m−250.1\,\mathrm{mW}\,\mathrm{m}^{-2} a 10 km10\,\mathrm{km}, 37.437.4 a 20 km20\,\mathrm{km} y 31.7 mW m−231.7\,\mathrm{mW}\,\mathrm{m}^{-2} a 35 km35\,\mathrm{km}.

17. 30/84.6=35 %30/84.6 = 35\,\%: casi dos tercios del flujo de calor sobre E se producen en su capa radiogénica.

18. T(10)=10+120+218.4×0.632=268∘CT(10) = 10 + 120 + 218.4\times0.632 = 268{}^{\circ}\mathrm{C}.

19. T(20)=10+240+218.4×0.865=439∘CT(20) = 10 + 240 + 218.4\times0.865 = 439{}^{\circ}\mathrm{C}.

20. Bajo A, A0D2/k=1.35×100/2.5=54.0 KA_0D^2/k = 1.35\times100/2.5 = 54.0\,\mathrm{K}: T(35)=10+420+54.0×0.970=482∘CT(35) = 10 + 420 + 54.0\times0.970 = 482{}^{\circ}\mathrm{C}.

21. Por A0D2/kA_0D^2/k: 54 K54\,\mathrm{K} bajo A y 218 K218\,\mathrm{K} bajo E. El flujo de calor profundo es el mismo; la corteza de E es cuatro veces más rica en fuentes de calor, y el calor que producen debe conducirse hacia arriba por la misma roca, lo que exige un gradiente más empinado en todo el recorrido.

22. 1300=10+218.4+12z1300 = 10 + 218.4 + 12z da z=89.3 kmz = 89.3\,\mathrm{km} bajo E; 1300=10+54.0+12z1300 = 10 + 54.0 + 12z da 103.0 km103.0\,\mathrm{km} bajo A.

23. Ambas se encuentran dentro del rango de la litosfera continental, desde cien kilómetros o menos bajo regiones jóvenes hasta unos 250 km250\,\mathrm{km} bajo cratones, hacia su lado delgado: con un flujo de calor reducido de 30 mW m−230\,\mathrm{mW}\,\mathrm{m}^{-2} y esta conductividad, el geotermo conductivo alcanza temperaturas del manto a unos 100 km100\,\mathrm{km}; bajo el plutón más radiactivo, 14 km14\,\mathrm{km} antes.

24. T(35)=10+420+218.4×(1−e−3.5)=430+211.8T(35) = 10 + 420 + 218.4\times(1 - e^{-3.5}) = 430 + 211.8: el Moho bajo el plutón E está a 642∘C642{}^{\circ}\mathrm{C}.

本章定义的术语

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