Dónde perforar: la vegetación que no se seca marca el sitio

En roca granítica no hay acuífero cartografiado: hay una franja de roca meteorizada y fisurada, de espesor variable, que a veces tiene agua y casi siempre está donde no se ve. Este trabajo la localiza combinando cuatro años de imagen diaria con dos líneas de tomografía eléctrica, y acaba con una coordenada, una profundidad y un intervalo de captación.

La idea

Al final del estiaje, la mayor parte de la vegetación de una ladera costera del Pacífico ha perdido verdor: vive de la lluvia del año y la lluvia se acabó en octubre. Lo que sigue verde en junio está bebiendo de otra cosa.

Esa otra cosa, en un terreno así, es agua subterránea somera en la roca descompuesta. El rastro no aparece en una imagen: aparece al comparar la mediana de 11 fechas del final de las lluvias con la de 23 fechas del final del estiaje, sobre 32 y 15 escenas de 3 m, entre 2022 y 2026.

La escena entera pierde verdor al secarse. Lo que interesa no es el valor absoluto sino quién pierde menos que el resto, y por eso la escala del panel derecho se centra en la pérdida mediana y no en el cero. Es el mismo problema que una medida sin anclar: si el desplazamiento común se lleva toda la escala, no se ve nada.

Las capas y sus pesos

La capa de verdor no decide sola. Entra en un índice de prospectividad con otras dos, y se lleva el peso mayor:

Capa	Peso	Por qué
Vegetación que resiste la seca	0,45	es la única que responde al agua, no a la forma del terreno
Posición en la red de drenaje	0,35	la escorrentía converge donde se recarga
Espesor de roca descompuesta	0,20	sin espesor alterado no hay dónde almacenar

El resultado delimita una zona objetivo y, dentro de ella, el sitio donde se tienden las dos líneas de tomografía. Aquí es donde el satélite deja de decidir: lo único que ha hecho es reducir una ladera entera a unas hectáreas, para que la geofísica de campo se gaste donde tiene sentido.

El índice se acompaña de su incertidumbre. El panel izquierdo da la probabilidad de que cada punto esté en el 10 % mejor del área cuando se perturban los pesos; el derecho, la desviación típica del índice bajo esa misma perturbación. Una zona objetivo que se deshace al mover los pesos un poco no es una zona objetivo.

El modelo del terreno

En roca granítica el agua se organiza en tres capas: la saprolita, roca descompuesta que almacena pero conduce mal; la capa fisurada debajo, que conduce y es donde se capta; y la roca sana, que no da nada. El pozo tiene que cortar la capa fisurada por debajo del nivel del agua, y ese es el intervalo que hay que identificar antes de perforar.

Lo que vio la geofísica

La línea L1, de 188 m con 48 electrodos cada 4 m, se invirtió combinando dos dispositivos —gradiente y dipolo-dipolo— con un error RMS del 10,9 % sobre 1 555 datos. La sección separa con claridad:

• Saprolita arcillosa húmeda, de 20 a 45 Ω·m, en el extremo noroeste.
• Capa laminada en la base de la saprolita, seguida sin interrupción entre los 76 y los 140 m de la línea: es el marcador que da continuidad al modelo.
• Capa fisurada, de 73 a 114 Ω·m, más somera y más alterada en el sector central. Es el objetivo.
• Bloques y gravas secas del cauce, por encima de 1 400 Ω·m, y una zona húmeda somera de 85 a 150 Ω·m bajo el flanco oeste del cauce, entre 10 y 18 m de profundidad.

Hay además un dato de campo que ancla todo: un afloramiento de agua a 34,5 m de la línea, a la altura del electrodo 28.

L2 cruza a L1 y confirma la geometría desde la otra dirección.

La propuesta

	Sondeo principal	Alternativo
Posición	L1, x = 112 m	L1, x = 124 m
Cota del brocal	19,4 m	18,8 m
Profundidad total	30 m	32 m
Intervalo objetivo	16 – 30 m	22 – 32 m
Espesor de capa fisurada cortada	17,6 m	15,6 m
Resistividad del tramo objetivo	205 Ω·m	378 Ω·m

El sondeo principal corta 4,4 m de suelo y saprolita, 7,0 m de saprolita arenosa y 5,0 m de capa laminada antes de entrar en la capa fisurada a los 16,4 m. La resistividad más baja del tramo objetivo en P1 —205 frente a 378 Ω·m— es lo que lo pone por delante del alternativo: menos resistivo, a igualdad de litología, significa más agua.

Límites de este trabajo

• El satélite no ve agua. Ve vegetación que se comporta como si la hubiera. Es una pista, y por eso se contrasta con geofísica antes de proponer nada.
• El índice de prospectividad es un consenso ponderado, no una medida. Los pesos son un juicio experto, y por eso se publica qué pasa cuando se mueven.
• La tomografía eléctrica resuelve contrastes de resistividad, no litología: la interpretación en saprolita, capa laminada y capa fisurada se apoya en el modelo de meteorización de la roca granítica y en el afloramiento medido.
• Por debajo del límite de investigación de cada línea el dato ya no restringe el modelo, y esa zona se marca como no interpretada en las secciones.
• Ningún estudio indirecto garantiza caudal. Lo que se entrega es la mejor posición, la profundidad y el intervalo a captar; el caudal lo dice el aforo.

Qué se entrega

Producto	Formato
Índice de prospectividad y sus capas	GeoTIFF georreferenciado
Probabilidad e incertidumbre del índice	GeoTIFF
Secciones de resistividad invertidas e interpretadas	imagen y datos del modelo
Columna prevista de cada sondeo propuesto	CSV y figura
Coordenada, profundidad e intervalo de captación	CSV
Informe con la metodología y las comprobaciones	Markdown y PDF

Fuentes y licencias

• PlanetScope Super Dove (Planet Labs PBC), 8 bandas a 3 m, licencia comercial.
• Copernicus DEM GLO-30 (ESA). Uso abierto con atribución.
• Tomografía eléctrica de resistividad levantada en campo, con dispositivos gradiente y dipolo-dipolo, 48 electrodos por línea.

Aviso

Este documento presenta una propuesta de sondeo basada en métodos indirectos. No garantiza la presencia de agua ni el caudal de un pozo, y no sustituye el aforo ni el permiso de aprovechamiento correspondiente.