GRACE — Almacenamiento de Agua Terrestre

Regiones

hidrológicas analizadas

Periodo

253

meses · 2002-04 – 2026-03

Mayor pérdida

-8.6 mm/año (X)

tendencia secular más negativa

Mayor ganancia

+3.2 mm/año (XII)

tendencia secular más positiva

Meses sequía severa

268

z-score < −1.5 en todas las regiones

Validación de datos

✓ Unidades: cm (lwe_thickness attrs)

✓ Conversión ×10 aplicada (cm → mm)

✓ Gain factors CLM4.0 (scale_factor embebido)

✓ CRI-filtered (leakage costero corregido)

✓ SHA-1: dc826f343daa… (integridad del archivo)

✓ CRS áreas: EPSG:6933 (Equal Earth — igual-área)

Distribución geográfica — Tendencia secular TWSA

Haz clic en una región para ir a su análisis detallado.

⬇ Descargar mapa (PNG)

Tendencias seculares por región

ID	Nombre	Área (km²)	Tendencia	Incert.	Amp. anual (mm)	Amp. semi-anual (mm)	Meses	Sequía severa	Calidad
X	Golfo Centro	104,780	-8.58 mm/año	±1.73	153.1	21.0	253	20	~ caution
IV	Balsas	119,234	-7.98 mm/año	±1.15	114.5	19.0	253	21	~ caution
XI	Frontera Sur	101,219	-7.33 mm/año	±2.43	143.6	7.4	253	19	~ caution
VI	Río Bravo	379,582	-7.30 mm/año	±1.32	12.0	6.1	253	28	✓ safe
II	Noroeste	206,454	-6.84 mm/año	±1.05	9.3	7.7	253	24	✓ safe
III	Pacífico Norte	151,976	-6.49 mm/año	±1.47	54.3	22.7	253	23	~ caution
VIII	Lerma Santiago Pacífico	190,733	-6.23 mm/año	±1.46	125.8	28.9	253	17	~ caution
VII	Cuencas Centrales del Norte	202,459	-5.67 mm/año	±1.27	32.4	11.0	253	22	✓ safe
XIII	Aguas del Valle de México	16,423	-5.18 mm/año	±1.14	70.4	11.4	253	22	! context
V	Pacífico Sur	77,506	-5.11 mm/año	±0.98	107.7	12.9	253	23	! context
IX	Golfo Norte	127,120	-3.98 mm/año	±0.97	30.5	6.4	253	19	~ caution
I	Península de Baja California	147,630	-2.14 mm/año	±0.36	2.7	3.2	253	20	~ caution
XII	Península de Yucatán	138,697	+3.20 mm/año	±1.85	75.9	2.5	253	10	~ caution

Calidad GRACE según área de la región (Handbook JPL L3, §7.3): ✓ safe ≥ 200,000 km² — señal confiable; ~ caution 100,000–200,000 km² — usar con precaución; ! context < 100,000 km² — SNR bajo, posible contaminación por leakage. Resolución nativa del mascon JPL: ~300–350 km (~110,000 km² por celda).

Peores meses de estrés hídrico (z-score < −1.5)

Región	Nombre	Mes	TWSA	Z-score
VIII	Lerma Santiago Pacífico	Sep 2023	-142.6 mm	-2.91 🔴
VIII	Lerma Santiago Pacífico	Aug 2023	-175.5 mm	-2.79 🔴
VIII	Lerma Santiago Pacífico	Oct 2023	-109.3 mm	-2.74 🔴
IV	Balsas	Oct 2023	-122.8 mm	-2.70 🔴
XIII	Aguas del Valle de México	Oct 2023	-123.1 mm	-2.68 🔴
V	Pacífico Sur	Oct 2023	-65.2 mm	-2.68 🔴
V	Pacífico Sur	Jul 2023	-177.8 mm	-2.62 🔴
XIII	Aguas del Valle de México	Jul 2023	-175.0 mm	-2.61 🔴
VIII	Lerma Santiago Pacífico	May 2024	-346.0 mm	-2.60 🔴
V	Pacífico Sur	Nov 2023	-70.3 mm	-2.58 🔴
XIII	Aguas del Valle de México	Sep 2023	-137.8 mm	-2.54 🔴
X	Golfo Centro	Nov 2023	-125.4 mm	-2.54 🔴
XIII	Aguas del Valle de México	Aug 2023	-146.9 mm	-2.53 🔴
V	Pacífico Sur	Apr 2024	-230.1 mm	-2.51 🔴
IX	Golfo Norte	Aug 2023	-110.4 mm	-2.48 🔴
V	Pacífico Sur	Dec 2023	-103.9 mm	-2.46 🔴
X	Golfo Centro	Oct 2023	-149.8 mm	-2.45 🔴
IV	Balsas	Jul 2023	-234.5 mm	-2.43 🔴
IV	Balsas	Apr 2024	-298.3 mm	-2.43 🔴
X	Golfo Centro	Jun 2024	-410.4 mm	-2.41 🔴

Z-score estandarizado por mes calendario. Muestra los 20 meses más severos.

Volumen de anomalía TWSA — México total (km³)

Metodología — Conversión TWSA a volumen
V [km³] = TWSA [mm] × A [km²] × 10⁻⁶
(1 mm de agua × 1 km² = 10⁻³ m × 10⁶ m² = 10³ m³ = 10⁻⁶ km³)

Las áreas se calculan a partir de los polígonos reprojectados a EPSG:6933 (WGS 84 / Equal Earth), proyección cilíndrica de igual área que conserva superficies con error < 0.1% en latitudes medias. El volumen total de México es la suma directa de las 13 RHA (regiones no traslapadas).

Equivalente metodológico: Khorrami et al. (2023, GRL) — convierte mascon TWS a km³/año multiplicando por el área de los tiles JPL.

Contribución por región hidrológica (km³)

ID	Nombre	Área (km²)	TWSA mín (km³)	TWSA máx (km³)	std (km³)
VI	Río Bravo	379,582	-77.845	24.375	23.135
VIII	Lerma Santiago Pacífico	190,733	-65.995	42.141	22.114
X	Golfo Centro	104,780	-43.003	38.501	15.172
II	Noroeste	206,454	-39.883	11.906	11.717
IV	Balsas	119,234	-39.809	25.539	13.136
VII	Cuencas Centrales del Norte	202,459	-38.687	24.634	11.970
III	Pacífico Norte	151,976	-38.239	28.692	12.035
XI	Frontera Sur	101,219	-37.485	33.068	14.363
XII	Península de Yucatán	138,697	-22.832	38.835	12.069
V	Pacífico Sur	77,506	-20.245	16.365	7.324
IX	Golfo Norte	127,120	-18.065	12.666	5.956
I	Península de Baja California	147,630	-9.504	4.091	2.786
XIII	Aguas del Valle de México	16,423	-3.776	2.771	1.254

Detalle por región

+ Toca cada región para ver sus estadísticas, mapa de sequía y serie temporal completa.

Golfo Centro

-8.58 mm/año ± 1.73 mm/año ~ caution

TWSA mínimo

-410.4 mm

TWSA máximo

367.4 mm

Z mínimo

-2.54 🔴

Sequía severa

20 meses

Sequía moderada

20 meses

Peor mes

Nov 2023

Z-score mensual (pasa el cursor sobre cada celda)

sequía severa moderada normal excedente

↑ Volver al mapa

Balsas

-7.98 mm/año ± 1.15 mm/año ~ caution

TWSA mínimo

-333.9 mm

TWSA máximo

214.2 mm

Z mínimo

-2.70 🔴

Sequía severa

21 meses

Sequía moderada

26 meses

Peor mes

Oct 2023

Z-score mensual (pasa el cursor sobre cada celda)

sequía severa moderada normal excedente

↑ Volver al mapa

Frontera Sur

-7.33 mm/año ± 2.43 mm/año ~ caution

TWSA mínimo

-370.3 mm

TWSA máximo

326.7 mm

Z mínimo

-2.16 🔴

Sequía severa

19 meses

Sequía moderada

21 meses

Peor mes

Sep 2019

Z-score mensual (pasa el cursor sobre cada celda)

sequía severa moderada normal excedente

↑ Volver al mapa

Río Bravo

-7.30 mm/año ± 1.32 mm/año ✓ safe

TWSA mínimo

-205.1 mm

TWSA máximo

64.2 mm

Z mínimo

-2.36 🔴

Sequía severa

28 meses

Sequía moderada

24 meses

Peor mes

May 2025

Z-score mensual (pasa el cursor sobre cada celda)

sequía severa moderada normal excedente

↑ Volver al mapa

Noroeste

-6.84 mm/año ± 1.05 mm/año ✓ safe

TWSA mínimo

-193.2 mm

TWSA máximo

57.7 mm

Z mínimo

-2.36 🔴

Sequía severa

24 meses

Sequía moderada

28 meses

Peor mes

Jan 2025

Z-score mensual (pasa el cursor sobre cada celda)

sequía severa moderada normal excedente

↑ Volver al mapa

III

Pacífico Norte

-6.49 mm/año ± 1.47 mm/año ~ caution

TWSA mínimo

-251.6 mm

TWSA máximo

188.8 mm

Z mínimo

-2.35 🔴

Sequía severa

23 meses

Sequía moderada

21 meses

Peor mes

Feb 2025

Z-score mensual (pasa el cursor sobre cada celda)

sequía severa moderada normal excedente

↑ Volver al mapa

VIII

Lerma Santiago Pacífico

-6.23 mm/año ± 1.46 mm/año ~ caution

TWSA mínimo

-346.0 mm

TWSA máximo

220.9 mm

Z mínimo

-2.91 🔴

Sequía severa

17 meses

Sequía moderada

27 meses

Peor mes

Sep 2023

Z-score mensual (pasa el cursor sobre cada celda)

sequía severa moderada normal excedente

↑ Volver al mapa

VII

Cuencas Centrales del Norte

-5.67 mm/año ± 1.27 mm/año ✓ safe

TWSA mínimo

-191.1 mm

TWSA máximo

121.7 mm

Z mínimo

-2.26 🔴

Sequía severa

22 meses

Sequía moderada

29 meses

Peor mes

May 2025

Z-score mensual (pasa el cursor sobre cada celda)

sequía severa moderada normal excedente

↑ Volver al mapa

XIII

Aguas del Valle de México

-5.18 mm/año ± 1.14 mm/año ! context

⚠ Señal de baja confiabilidad — Área 16,423 km² < 100,000 km² (umbral mínimo recomendado por el Handbook JPL L3 §7.3). La resolución nativa del mascon es ~300–350 km (~110,000 km² por celda): esta región es submascon, con SNR bajo y posible contaminación por leakage. Interpretar tendencias y anomalías con precaución.

TWSA mínimo

-229.9 mm

TWSA máximo

168.7 mm

Z mínimo

-2.68 🔴

Sequía severa

22 meses

Sequía moderada

16 meses

Peor mes

Oct 2023

Z-score mensual (pasa el cursor sobre cada celda)

sequía severa moderada normal excedente

↑ Volver al mapa

Pacífico Sur

-5.11 mm/año ± 0.98 mm/año ! context

⚠ Señal de baja confiabilidad — Área 77,506 km² < 100,000 km² (umbral mínimo recomendado por el Handbook JPL L3 §7.3). La resolución nativa del mascon es ~300–350 km (~110,000 km² por celda): esta región es submascon, con SNR bajo y posible contaminación por leakage. Interpretar tendencias y anomalías con precaución.

TWSA mínimo

-261.2 mm

TWSA máximo

211.2 mm

Z mínimo

-2.68 🔴

Sequía severa

23 meses

Sequía moderada

19 meses

Peor mes

Oct 2023

Z-score mensual (pasa el cursor sobre cada celda)

sequía severa moderada normal excedente

↑ Volver al mapa

Golfo Norte

-3.98 mm/año ± 0.97 mm/año ~ caution

TWSA mínimo

-142.1 mm

TWSA máximo

99.6 mm

Z mínimo

-2.48 🔴

Sequía severa

19 meses

Sequía moderada

27 meses

Peor mes

Aug 2023

Z-score mensual (pasa el cursor sobre cada celda)

sequía severa moderada normal excedente

↑ Volver al mapa

Península de Baja California

-2.14 mm/año ± 0.36 mm/año ~ caution

TWSA mínimo

-64.4 mm

TWSA máximo

27.7 mm

Z mínimo

-2.39 🔴

Sequía severa

20 meses

Sequía moderada

37 meses

Peor mes

Feb 2025

Z-score mensual (pasa el cursor sobre cada celda)

sequía severa moderada normal excedente

$TWSA I$

↑ Volver al mapa

XII

Península de Yucatán

+3.20 mm/año ± 1.85 mm/año ~ caution

TWSA mínimo

-164.6 mm

TWSA máximo

280.0 mm

Z mínimo

-2.17 🔴

Sequía severa

10 meses

Sequía moderada

36 meses

Peor mes

Oct 2009

Z-score mensual (pasa el cursor sobre cada celda)

sequía severa moderada normal excedente

↑ Volver al mapa

Metodología y reproducibilidad

1. Datos de entrada

Archivo GRACEGRCTellus.JPL.200204_202603.GLO.RL06.3M.MSCNv04CRI.nc

SHA-1 (integridad)dc826f343daa…

ProductoJPL GRACE/GRACE-FO RL06.3 V4 · Mascon CRI-filtered

Variable leídalwe_thickness — Liquid Water Equivalent Thickness

Unidades originalescm (confirmado en atributo units del NetCDF)

Conversión a mm× 10 (1 cm = 10 mm)

Periodo temporal2002-04-01 → 2026-03-01 (253 meses; 0 meses sin datos)

Resolución nativa~300–350 km (mascon de ~3° de soporte; la resolución efectiva varía con la latitud); muestreado en grilla 0.5°

2. Anomalía de almacenamiento (TWSA)

Referencia interna JPL El archivo ya es anomalía: JPL removió la media de 2004.000 to 2009.999 (atributo time_mean_removed). No se lee TWS absoluto.

Baseline del pipeline 2004-01-01 / 2009-12-01 — Se calcula la media regional en este periodo y se sustrae de toda la serie. Si coincide con el periodo JPL, la sustracción es ≈ 0 (neutro). Si difiere, re-centra la anomalía. Incluso cuando coinciden, el paso garantiza consistencia a nivel de agregación regional, ya que el promedio espacial ponderado puede introducir desviaciones leves respecto a cero.

Fórmula TWSA_usuario(t) = TWSA_JPL(t) − mean(TWSA_JPL, baseline)
No es una anomalía de TWS absoluto, sino un re-centrado de la anomalía JPL a un periodo de referencia diferente.

Nota sobre tws_mm vs twsa_mm La columna tws_mm en los CSV es TWSA_JPL: la agregación directa de lwe_thickness × 10, cuya media temporal ha sido removida por JPL para el periodo 2004.0–2009.999. twsa_mm es el re-centrado al baseline especificado por el usuario. Ambas columnas son anomalías (no almacenamiento absoluto); la diferencia es únicamente el periodo de referencia.

Componentes físicos incluidos TWSA representa la anomalía del almacenamiento terrestre de agua integrado verticalmente, incluyendo: aguas subterráneas, humedad del suelo, agua superficial (ríos, lagos, humedales), nieve y agua en el dosel vegetal. GRACE no separa estas componentes individualmente; interpolar la señal como "déficit de agua subterránea" exclusivamente constituye una sobreinterpretación.

3. Factores de ganancia (Gain Factors CLM4.0)

AplicadosSí — scale_factor embebido (CLM4.0, JPL RL06.3 V4)

Propósito Recuperan la amplitud de la señal TWS atenuada por el filtrado espacial de la solución mascon. Factor espacialmente variable (un valor por pixel de 0.5°), derivado del modelo CLM4.0 (Community Land Model).

Fórmula lwe_corr(x,y,t) = lwe(x,y,t) × scale_factor(x,y) — aplicado pixel a pixel antes del promedio regional.

CRI (Coastline Resolution Improvement) El archivo GRCTellus.JPL.200204_202603.GLO.RL06.3M.MSCNv04CRI.nc incluye el filtro CRI (CRIv04) que reduce el leakage entre tierra y océano en regiones costeras. Esto hace que las RHA costeras (Noroeste, Frontera Sur, Pacífico Norte/Sur, Golfo Norte/Centro, Yucatán) sean más confiables que con productos sin CRI. Ref: Wiese et al. (2016, Water Resources Research).

Advertencia CLM4.0 no incluye agua subterránea ni actividad humana. Esta limitación puede afectar tanto la magnitud como la distribución espacial de las tendencias seculares, así como la comparabilidad con estudios que usan otros modelos de gain factors (p. ej. WGHM, VIC). El impacto es mayor en regiones con extracción intensa de agua subterránea (p. ej. Valle de México, Río Bravo), donde el factor CLM4.0 podría sub-corregir la señal (Handbook JPL L3 §7.3; Landerer & Swenson 2012).

Incertidumbre NO se aplican gain factors a la variable uncertainty: la incertidumbre formal del mascon ya está en el espacio observacional (Watkins et al. 2015; Wiese et al. 2016).

4. Áreas de las regiones hidrológicas

Origen Calculadas en el pipeline a partir de los polígonos del GPKG, no leídas de un atributo pre-existente del shapefile.

CRS de cálculo EPSG:6933 — WGS 84 / Equal Earth. Proyección cilíndrica de igual área (Šavrič et al. 2019). Error de superficie < 0.1% en latitudes medias. Fórmula: área_km² = geometry.area / 1e6 después de reprojectar.

Promedio regional Promedio ponderado por área de las celdas GRACE que intersectan cada región. Peso de celda = fracción del área de la celda dentro del polígono. Los pesos se precalculan y almacenan en caché (.npz).

Supuesto en propagación de incertidumbre La incertidumbre regional se propaga como raíz del promedio ponderado de varianzas, asumiendo independencia entre celdas mascon. En realidad, los errores GRACE exhiben correlación espacial, por lo que la incertidumbre formal podría estar ligeramente subestimada para regiones grandes.

ID	Nombre	Área EPSG:6933 (km²)	Celdas GRACE	Soporte efectivo (km²)	Leakage index	Cobertura baseline	Calidad
I	Península de Baja California	147,630	107	147,630	1.00	100.0%	~ caution
II	Noroeste	206,454	108	206,454	1.00	100.0%	✓ safe
III	Pacífico Norte	151,976	89	151,976	1.00	100.0%	~ caution
IV	Balsas	119,234	67	119,234	1.00	100.0%	~ caution
IX	Golfo Norte	127,120	70	127,120	1.00	100.0%	~ caution
V	Pacífico Sur	77,506	50	77,506	1.00	100.0%	! context
VI	Río Bravo	379,582	189	379,582	1.00	100.0%	✓ safe
VII	Cuencas Centrales del Norte	202,459	111	202,459	1.00	100.0%	✓ safe
VIII	Lerma Santiago Pacífico	190,733	104	190,733	1.00	100.0%	~ caution
X	Golfo Centro	104,780	59	104,780	1.00	100.0%	~ caution
XI	Frontera Sur	101,219	52	101,219	1.00	100.0%	~ caution
XII	Península de Yucatán	138,697	83	138,697	1.00	100.0%	~ caution
XIII	Aguas del Valle de México	16,423	13	16,423	1.00	100.0%	! context

Leakage index = soporte efectivo / área región (~1.0 = cobertura completa; <0.9 = posible subestimación de señal). Cobertura baseline = % de meses del periodo de referencia con datos válidos.

5. Conversión a volumen (km³)

Fórmula V [km³] = TWSA [mm] × A [km²] × 10⁻⁶
Derivación: 1 mm × 1 km² = 10⁻³ m × 10⁶ m² = 10³ m³ = 10⁻⁶ km³

Supuesto espacial La conversión asume distribución espacialmente uniforme de la anomalía dentro de cada región. Esto es una aproximación dado que la resolución nativa del mascon JPL es ~300–350 km y las regiones son heterogéneas internamente.

Total México Suma directa de twsa_km3 de las 13 RHA (regiones no traslapadas). Esta agregación asume independencia espacial entre regiones; sin embargo, la solución mascon JPL exhibe correlación espacial y leakage residual que puede introducir un sesgo menor en la suma nacional. Equivalente metodológico: Khorrami et al. (2023, GRL) — suma de mascon tiles JPL convertidos a km³/año.

Columna en CSV twsa_km3 en TWSA_baseline_*.csv

6. Tendencia secular

Modelo Mínimos cuadrados ordinarios (OLS) con modelo estándar GRACE:
TWSA(t) = a + b·t + c·sin(2πt) + d·cos(2πt) + e·sin(4πt) + f·cos(4πt)
donde t = años decimales. El ajuste simultáneo del ciclo anual y semi-anual elimina su influencia sobre el coeficiente secular b (Rodell et al. 2018 Nature; Scanlon et al. 2018 Nature Comm.).

Incertidumbre 1σ de la matriz de covarianza OLS, corregida por autocorrelación temporal AR(1): σ_b_corr = σ_b × √VIF, donde VIF = (1+φ)/(1−φ) y φ es el coeficiente AR(1) de los residuos (Tourian et al. 2018, JGR). Solo se aplica cuando hay ≥ 6 pares consecutivos.

Mínimo de datos 24 meses válidos para estimar tendencia.

Descomposición estacional El mismo ajuste OLS entrega la amplitud del ciclo anual A_anual = √(c²+d²) y semi-anual A_semianual = √(e²+f²). La amplitud anual refleja la recarga estacional (precipitación verano − evapotranspiración). Regiones con A_anual alta (~100–150 mm: Balsas, Lerma, Golfo Centro) tienen ciclo estacional dominante; las áridas (Noroeste, Baja California, Río Bravo: < 15 mm) tienen ciclo débil. Estas amplitudes se reportan en la tabla de tendencias.

ID	Nombre	Tendencia (mm/año)	Incertidumbre 1σ	N meses
I	Península de Baja California	-2.14	±0.36	253
II	Noroeste	-6.84	±1.05	253
III	Pacífico Norte	-6.49	±1.47	253
IV	Balsas	-7.98	±1.15	253
IX	Golfo Norte	-3.98	±0.97	253
V	Pacífico Sur	-5.11	±0.98	253
VI	Río Bravo	-7.30	±1.32	253
VII	Cuencas Centrales del Norte	-5.67	±1.27	253
VIII	Lerma Santiago Pacífico	-6.23	±1.46	253
X	Golfo Centro	-8.58	±1.73	253
XI	Frontera Sur	-7.33	±2.43	253
XII	Península de Yucatán	+3.20	±1.85	253
XIII	Aguas del Valle de México	-5.18	±1.14	253

7. Z-score mensual (estrés hídrico)

Fórmula z(t) = (TWSA(t) − μ_mes) / σ_mes
donde μ_mes y σ_mes son la media y desviación estándar de todos los valores del mismo mes calendario (enero, febrero, …) para cada región. Elimina la estacionalidad y expresa cuán anómalo es cada mes en relación a ese mes históricamente.

Interpretación El z-score es un indicador de anomalía relativa de almacenamiento, no un déficit hídrico directo ni un índice de sequía operacional. Valores bajos indican almacenamiento inusualmente bajo en relación al historial del mismo mes calendario, lo que puede correlacionarse con condiciones de sequía pero no equivale a balance hídrico.

Umbrales z < −1.5 → almacenamiento severo por debajo de lo normal | −1.5 ≤ z < −1.0 → almacenamiento moderadamente por debajo de lo normal | z > +1.5 → almacenamiento excepcionalmente alto

Salida twsa_zscore en TWSA_baseline_*.csv; drought_months.csv con todos los meses z < −1.0

8. Clasificación de calidad por región

Criterio Área de la región vs resolución nativa del mascon JPL (~300–350 km ≈ 110,000 km²).

Niveles ✓ safe ≥ 200,000 km² — señal confiable (≥ 2 mascons completos)
~ caution 100,000–200,000 km² — usar con precaución (1–2 mascons)
! context < 100,000 km² — SNR bajo, posible leakage de señal adyacente

Referencia Handbook JPL L3 §7.3 (2021); Watkins et al. (2015, GRL); Wiese et al. (2016, GRACE-FO).

9. Software y reproducibilidad

grace_mx versiónv1.2.0 · git: 487a0da

Fecha de procesamiento2026-05-19T20:50:22Z

Comando python3 scripts/run_rha13.py

Dependencias clave xarray, geopandas, numpy, pandas, matplotlib, pyproj, shapely

Correcciones pre-aplicadas en L3 Geocentro (grado-1, Swenson et al. 2008), C₂₀ (SLR, Cheng et al. 2011), GIA (modelo secular), atmósfera (ECMWF), océano (AOD1B), destriping + filtrado gaussiano 300 km, corrección elipsoidal (V4).

SHA-1 del archivo GRACE dc826f343daab0cc81c4c9eafe33f4796090b4dc

10. Referencias

Watkins, M.M. et al. (2015). Improved methods for observing Earth's time variable mass distribution with GRACE. JGR Solid Earth. doi:10.1002/2014JB011547
Wiese, D.N. et al. (2016). Quantifying and reducing leakage errors in the JPL RL05M GRACE mascon solution. Water Resources Research. doi:10.1002/2016WR019344
Landerer, F.W. & Swenson, S.C. (2012). Accuracy of scaled GRACE terrestrial water storage estimates. Water Resources Research. doi:10.1029/2011WR011453
Rodell, M. et al. (2018). Emerging trends in global freshwater availability. Nature. doi:10.1038/s41586-018-0123-1
Scanlon, B.R. et al. (2018). Global models underestimate large decadal declining and rising water storage trends. Nature Communications. doi:10.1038/s41467-018-04475-5
Khorrami, M. et al. (2023). Groundwater Volume Loss in Mexico City. Geophysical Research Letters. doi:10.1029/2022GL101962
Castellazzi, P. et al. (2016). Groundwater depletion in Central Mexico. Water Resources Research. doi:10.1002/2015WR018211
Tourian, M.J. et al. (2018). A quantile function approach to discharge estimation from satellite altimetry. JGR. doi:10.1029/2018WR023166
Šavrič, B. et al. (2019). The Equal Earth map projection. International Journal of Geographical Information Science. doi:10.1080/13658816.2018.1504949
JPL GRACE Tellus (2021). GRACE-FO Level-3 Data Handbook. NASA Jet Propulsion Laboratory. GRACE D-103133.

GRACE / GRACE-FO — Anomalía de Almacenamiento de Agua Terrestre (TWSA)

Distribución geográfica — Tendencia secular TWSA

Tendencias seculares por región

Peores meses de estrés hídrico (z-score < −1.5)

Volumen de anomalía TWSA — México total (km³)

Contribución por región hidrológica (km³)

Detalle por región

Metodología y reproducibilidad