PPM_DB3 · climatol 4.5 · 250 estaciones · 2000–2025
Desierto de Atacama antes y después de climatol
El desierto no es un régimen sino cuatro, y tratarlos juntos era el error. Separados, la homogeneización deja de inventar saltos y — lo que más importa aquí — deja de borrar la lluvia que sí cayó.
Una sola corrida, sin separar regímenes, marcaba 44 saltos. Segmentadas por régimen quedan 4: Visviri, enero de 2005; El Corral, enero de 2010; Los Loros, enero de 2016; Rio Lluta En Alcerreca, enero de 2021. La diferencia no es indulgencia sino aritmética de la referencia: entre estaciones del mismo régimen la correlación mensual mediana es 0,78; entre regímenes distintos, 0,16. Comparar el altiplano con la costa fabrica los saltos que luego se «corrigen».
Que queden tan pocos depende de una decisión que conviene declarar. El diseño de dos pasos dice que la escala anual detecta los saltos y la mensual los corrige. climatol desactiva la detección mensual cuando recibe metadatos, así que si esa entrega se activa solo cuando el paso anual encontró algo, unos regímenes buscan y otros no. Aquí se entrega siempre, vacía si hace falta: los cuatro reciben el mismo trato. Con el criterio anterior aparecían seis saltos mensuales más, que eran artefacto del procedimiento y no del dato.
El de Cancosa vale la pena por separado: climatol lo situó en marzo de 2019 sin que se le dijera nada, a dos meses de la fecha en que la limpieza ya había cortado esa serie por dejar de ser creíble. Al mirarlo de cerca resultó no ser una inhomogeneidad sino un cambio de variable medida: hasta 2018 concentra el 96 % de su lluvia en diciembre-marzo, como cualquier estación de altiplano; desde 2019 sus máximos son julio (446 mm), junio (432) y mayo (421). Y los valores no vienen del portal primario sino de CR2, el redistribuidor. Se conserva marcada, no borrada.
que acaban en cero (eran 110)
(la corrida única daba 44)
es observación
corridos por separado
Las seis pruebas
Una estación por localidad
Totales anuales. La barra gris ancha es el año ya homogeneizado y rellenado; el núcleo de color es lo que midió el pluviómetro. Si el núcleo llega arriba del gris, la homogeneización bajó ese año; si el gris asoma, lo subió. El año observado solo se dibuja cuando tiene los doce meses medidos —un mes sin datos no es un mes seco—, así que su ausencia marca los años que existen únicamente gracias al relleno.
Mes a mes, arriba lo medido y abajo lo homogeneizado
Las bandas grises son los meses sin observación. En la fila ORIGINAL la línea desaparece ahí; en HOMOG. continúa, y esa continuación es estimación, no medida. A la derecha, el total anual medio de cada versión.
La validación que decide. Dos pluviómetros de instituciones distintas a menos de 1 km miden la misma lluvia, y ninguno le fue declarado a climatol como pareja del otro: si la homogeneización estropeara algo, se notaría aquí. Hay 15 parejas con al menos doce meses medidos en ambas.
Sobre los meses que ambas midieron, no cambia nada. El RMSE entre parejas es 11,33 mm antes y 11,33 mm después, y la correlación 0,891 en los dos casos: la homogeneización no toca la celda medida. Lo que sí cambia es qué meses se comparan, porque la serie homogeneizada está completa: suma 2.409 pares de meses en que una de las dos es una estimación, y ahí el RMSE es 9,78 mm. Una versión anterior de esta figura leyó esa diferencia como «el relleno degrada a las vecinas»; era aritmética de conjuntos de meses distintos.
Lo que corrigió la v5: no rellenar lo que no se puede. En la v4, Río Copiapó en Pastillo midió 278,4 mm en marzo de 2014 y a Altar de la Virgen, a 522 m, el relleno le puso 0,0. Primero se probó reasignarle el régimen de su vecina: ese mes pasó a 227,2 mm, pero otros se fueron al lado contrario —307 mm donde Pastillo midió 81—. climatol rellena escalando por la media propia de cada serie, y la de Altar salía de cuatro años, 2022-2026, con un julio de 188 mm: cualquier media así sesga todos los meses que rellena. Desde la v5, una serie con menos de 48 meses medidos no se homogeneiza ni se rellena: sigue en la base observada, marcada, y salen 25 series que antes se rellenaban.
Y lo que no corrigió. En los pares que quedan, el relleno
se comporta casi igual que en la v4: la mejora del conjunto viene de no
rellenar donde no se podía, no de rellenar mejor. Las parejas que más se
separan siguen siendo del altiplano —Putre, Alcérreca, Chilcaya, Tacora—,
donde una tormenta convectiva cae sobre un pluviómetro y no sobre el de al
lado, y ninguna estimación regional puede reproducirla. Para estadística
espacial o comparación entre estaciones, usar solo bandera = 0.
Lo que este trabajo no resuelve. Un 39,0 % de la matriz es relleno. Sirve para detectar saltos, que era el objetivo. Para calcular tendencias, o para cualquier comparación entre estaciones, conviene usar solo las celdas observadas.