Análisis satelital independiente · grace_mx v1.2.0

Almacenamiento de Agua Terrestre
en México, 2002–2026

Anomalías de almacenamiento de agua (TWSA) calculadas a partir de 251 meses de datos GRACE/GRACE-FO (JPL RL06.3 V4) para las 13 Regiones Hidrológicas Administrativas de México. Periodo: abril 2002 – enero 2026.

Reporte técnico completo Ver resultados
−259 km³
TWSA acumulada respecto
al baseline 2004–2009
12/13
RHA con tendencia negativa
significativa (S/R > 2)
−10.9 km³/año
Tasa de cambio integrada
(tendencia secular OLS)
251
meses de observación
abr 2002 – ene 2026
+3.1 mm/año
Tendencia RHA XII
Península de Yucatán

TWSA integrada — México (2002–2026)

Anomalía de almacenamiento de agua terrestre (TWSA) sumada sobre las 13 RHAs (1.96 × 10⁶ km²), expresada en km³ respecto al baseline 2004–2009. El valor de enero 2026 es −177 km³ — el más bajo del periodo de registro GRACE.

Anomalía de almacenamiento de agua en México 2002-2026
Anomalía de almacenamiento de agua terrestre (km³) integrada sobre las 13 regiones hidrológicas de México. La línea roja indica la tendencia secular. Fuente: GRACE/GRACE-FO JPL RL06.3 V4, procesado con grace_mx v1.2.0.

Fuente de datos y procesamiento

GRACE/GRACE-FO mide variaciones del campo gravitacional terrestre causadas por redistribución de masa de agua (subterránea, superficial y en suelo). La señal integrada de TWS no permite desagregar componentes individuales del ciclo hidrológico.

Trazabilidad: del satélite al dato procesado
🛰
Principio de medición

GRACE y GRACE-FO son pares de satélites en órbita polar separados ~220 km. Variaciones en la distancia inter-satelital (del orden de micrómetros) permiten inferir anomalías en el campo gravitacional local, que se relacionan con cambios en el almacenamiento de agua terrestre.

📡
Archivo fuente: JPL RL06.3 V4 Mascon CRI

El producto utilizado (GRCTellus.JPL.200204_202601.GLO.RL06.3M.MSCNv04CRI.nc) es el producto mascon de referencia más robusto y ampliamente citado en la literatura, descargado directamente del servidor PO.DAAC de la NASA (Watkins et al., 2015, JGR).

📐
Regiones oficiales de CONAGUA

Los límites geográficos corresponden a las 13 Regiones Hidrológicas Administrativas (RHA) definidas por la Comisión Nacional del Agua, en proyección Equal Area EPSG:6933.

🔢
Metodología estándar internacional

Agregación área-ponderada con gain factors CLM4.0, modelo OLS con armónicos anuales y semianuales, corrección AR(1). La misma metodología usada en Rodell et al. (2018, Nature) y Scanlon et al. (2012, WRR).

💻
Pipeline reproducible con trazabilidad completa

Cada salida registra el SHA1 del archivo GRACE, el hash del commit git y la versión del software. Cualquier persona puede descargar los datos de la NASA y reproducir estos resultados exactos.

Alcance y limitaciones

Qué mide este análisis
  • Almacenamiento total de agua terrestre (TWS): agua subterránea + humedad de suelo + agua superficial + nieve
  • Variación mensual respecto al baseline 2004–2009
  • Tendencia secular 2002–2026 por región hidrológica
  • Señal integrada a escala regional (~100,000–400,000 km²)
Qué no mide / limitaciones
  • No desagrega acuíferos individuales ni niveles piezométricos
  • No es precipitación, nivel de presas ni caudal de ríos
  • No permite inferir causas (extracción, clima, uso de suelo)
  • Resolución espacial mínima ~300 km — no válido para análisis locales
  • No es una predicción futura

Tendencias seculares por RHA

Tasas de cambio en TWSA (mm/año) estimadas mediante OLS con armónicos anuales/semianuales y corrección AR(1). Valores negativos indican reducción del almacenamiento respecto al baseline 2004–2009. Periodo: abril 2002 – enero 2026 (251 meses).

← Desliza para ver toda la tabla

RHA Región Tendencia mm/año ± 1σ S/R Magnitud relativa
XGolfo Centro−8.631.764.9×
IVBalsas−7.971.176.8×
XIFrontera Sur−7.342.463.0×
VIRío Bravo−7.211.325.5×
IINoroeste−6.801.076.4×
IIIPacífico Norte−6.431.504.3×
VIIILerma Santiago Pacífico−6.171.484.2×
VIICuencas Centrales del Norte−5.631.284.4×
XIIIAguas del Valle de México−5.201.164.5×
VPacífico Sur−5.120.995.2×
IXGolfo Norte−3.970.984.1×
IPenínsula de Baja California−2.140.375.8×
XIIPenínsula de Yucatán+3.111.861.7×
Pérdida severa (> 6 mm/año) Pérdida moderada (4–6 mm/año) Pérdida baja (< 4 mm/año) Ganancia

S/R = señal/ruido (|tendencia ÷ 1σ|). S/R > 2 indica tendencia estadísticamente significativa. Periodo: abril 2002 – enero 2026 (251 meses).

Metodología

Pipeline científico reproducible desarrollado en Python, siguiendo el estándar de la comunidad GRACE.

🛰
Datos fuente

GRACE/GRACE-FO JPL Mascon RL06.3 V4 con filtro CRI. Variable lwe_thickness (cm) convertida a mm.

⚖️
Gain Factors CLM4.0

Scale factors embebidos en el archivo JPL aplicados antes de la normalización de longitudes. Corrigen atenuación espacial.

📐
Agregación regional

Promedio área-ponderado de celdas mascon dentro de cada RHA. Pesos por intersección geométrica exacta.

📈
Tendencia secular

OLS con armónicos anuales y semianuales. Corrección AR(1) de Matalas-Jacobs para incertidumbre de tendencia.

📅
Baseline 2004–2009

Mismo periodo de referencia que JPL. Calculado dinámicamente desde los datos — nunca hardcodeado.

🔗
Trazabilidad

SHA1 del archivo GRACE, hash del commit git y versión del pipeline registrados en cada salida.

Fuente de datos Watkins, M.M., Wiese, D.N., Yuan, D.-N., Boening, C., Landerer, F.W. (2015). Improved methods for observing Earth's time variable mass distribution with GRACE using spherical cap mascons. JGR Solid Earth, 120, 2648–2671.

Datos descargados de NASA Physical Oceanography DAAC (PO.DAAC): podaac.jpl.nasa.gov

Reporte técnico

Incluye series de tiempo por región, estadísticas de tendencia, meses con z-score < −1.0, gain factors aplicados y metadata completa de reproducibilidad.

Ver reporte técnico