Precipitación mensual del Atacama

Cómo se hizo · 9 pasos · 2000–2025

PPM Atacama: de la medición a la publicación

El camino que siguió cada valor de precipitación, de la red que lo midió a esta página: qué se le hizo en cada paso, qué controles tuvo que pasar y qué series se quedaron fuera.

Ver como diagrama interactivo mejor en pantalla grande

La medición

Doce redes, un control de calidad y una base limpia: se marca en vez de borrar y el cero medido es un dato.

  1. 1

    Fuentes

    12 redes

  2. 2

    Control de calidad

    Estévez et al. (2011)

    marca, no borra
  3. 3

    PPM_DB3

    base consolidada

    0 mm es dato

Las redes y sus estaciones

RedPaísCon datosEn climatol
DGAChile156144
DMCChile6231
SENAMHIPerú2424
SENAMHIBolivia1713
INIAChile1615
CEAZAChile1411
APECO-FDFChile83
GHCNvarios76
INDAP-CONADIChile20
SQMChile22
MMA-DMCChile10
SNIHArgentina11
Total310250

«En climatol»: las estaciones con al menos 48 meses medidos, que se homogeneizan.

Control de calidad

La revisión de la estructura del registro y las pruebas de Estévez et al. (2011) para redes meteorológicas, con el rango fijo y el dinámico por separado. Cuatro sirven para un total mensual de lluvia, dos no y una en parte.

Antes de las pruebas: el mes tiene que estar completo

≥ 80 %
de los datos · CEAZA, INIA, SNIH y archivo local
≥ 25 días
SENAMHI Perú
> 20 días
DGA, según su propia marca
sin cuenta
DMC, CR2 y GHCN publican el total del mes sin decir cuántos días lo forman
PruebaQué comprueba¿Aplica?Cómo se aplica
Estructura del registrocódigos, fechas y archivos completos; lo que falla es un erroraplicacódigos de ausencia y negativos → dato ausente · duplicados fuera · fragmentos de un mismo sitio se unen · coordenadas corregidas donde dos registros de la DGA coinciden
Rango fijolímites físicos; lo que falla es un error y no sigueaplicase quita lo imposible: estaciones que no miden lluvia y quiebres de régimen evidentes
Rango dinámicofrente a lo esperable en ese lugar; lo que falla es sospechosoaplicamás de 6 veces el máximo de sus vecinas (≤ 60 km) ese mes: se marca, no se borra
Pasoel cambio entre dos mediciones seguidasno aplicala lluvia es discontinua: Estévez no la aplica
Consistencia internavariables de la misma estación entre síno aplicahay una sola variable
Persistenciaun sensor clavado en un valoren parteun año entero en 0,0 mm en una estación que suele medir lluvia se marca para revisar
Consistencia espacialfrente a las estaciones vecinasaplicaun cero pasa a dato ausente si sus vecinas de altitud parecida midieron más de 20 mm · después, climatol

Además

DGA > CR2
el portal de la DGA manda mes a mes; CR2 solo completa lo que falta
s/p = 0 mm
en la DMC; el «.» es dato ausente
DEM
la altitud sale del DEM donde falta, vale 0 o no calza con el terreno
± 500 m
cada estación se compara con sus vecinas de altitud parecida: OK, ceros sospechosos, sesgo o no evaluable

Estévez, J., Gavilán, P. y Giráldez, J. V. (2011). Guidelines on validation procedures for meteorological data from automatic weather stations. Journal of Hydrology, 402, 144–154.

Régimen y herencia

Menos de 4 años no se rellena. Las fiables estiman su régimen; las que no tienen datos en cada mes lo heredan.

  1. 4

    Régimen

    solo series de 4 años

    fiable ≥3 obs/mes
    saleExcluida < 4 años medidos
    alimenta · no fiableHereda régimen de sus vecinas fiablesy sigue a climatol

Cómo se asigna el régimen

A
altiplano: 60 % o más de la lluvia cae entre diciembre y marzo
B
costa: 25 % o menos en verano; llueve en invierno
C
transición: entre los dos
D
hiperárido: menos de 5 mm al año

Una serie es fiable con 48 meses medidos y 3 datos en cada mes del año. Las que tienen los meses pero no los datos de cada mes heredan el régimen de sus vecinas fiables (≤ 60 km, ≤ 800 m de desnivel).

climatol en dos pasos

El exploratorio fija inht; el anual detecta los saltos y el mensual corrige.

    .dat / .est
  1. 5

    Exploratorio

    fija inht

    entra · referencias *ERA5-Land
  2. inht
  3. 6

    climatol

    detecta y corrige

    año = 12 meses
  4. 7

    Base v5

    valida

    control

Cómo se corre

4 corridas
una por régimen: cada serie se compara solo con las de su mismo clima
ERA5-Land
sus píxeles, en el lugar de las estaciones, son las series de referencia
inht
el umbral de detección de saltos sale del exploratorio de cada régimen
anual → mensual
el total anual (solo años con 12 meses medidos) detecta los saltos; la serie mensual se corrige y se rellena

Publicación

Una base SQLite / D1, con CHIRPS en cada estación, que rehace las figuras y responde por API; los rásters se descargan en /rasters.

  1. 8

    SQLite · D1

    13 tablas

    entraCHIRPS en cada estación
  2. 9

    Sitio + API

    atacama.gfuentes.cl

    API · /rasters
    entra · Zarr en R2Reducción a 7,5″ altitud + NDVI
    desde: GWR + CHIRPS

Cómo se publica

SQLite · D1
la misma base para descargar y para responder la API
CSV
las series y las estaciones, para abrir en una planilla
API
estaciones y series por consulta
versiones
cada versión queda guardada; la anterior se marca como superada

Rásters

GWR por mes con tres protecciones, comparado con ANUSPLIN, CR2MET y CHIRPS; reducción a 7,5″. En R2 solo el Zarr.

  1. Comparación ANUSPLIN · CR2MET · CHIRPS
    desde: GWR + CHIRPS · LOOCV + exclusión
  2. Reducción a 7,5″ altitud + NDVI
    desde: GWR + CHIRPS

Cómo se hacen

GWR por mes
ln(PPM + 1) ~ ln(CHIRPS + 1), con el ancho elegido por validación cruzada
3 protecciones
pendiente mínima, CHIRPS dentro del rango de las vecinas y tope de 1,5 veces su máximo
5 km → 7,5″
cada celda se reparte según la altitud y el NDVI
25–100 km
radios de exclusión con que se valida, frente a ANUSPLIN, CR2MET y CHIRPS

Entradas: CHIRPS v2.0, NDVI de MODIS (MOD13Q1), DEM Copernicus GLO-30 y, para comparar, CR2MET v2.5.

En resumen

Decisiones

  • Mes sin datos ≠ mes seco
  • Marcar, no borrar
  • climatol solo cambia lo medido donde hay un salto

Controles

  • Mínimo 48 meses para rellenar
  • inht solo del exploratorio
  • Antes/después, mismos meses

Validación de los rásters

  • LOOCV y exclusión de 25–100 km
  • ANUSPLIN: mismo modelo, mismas celdas
  • Plausibilidad: 1,5 × el máximo medido