🏔️ Snow Cover & Energy Flux Analysis

Northern Pakistan (2015-2025) | ERA5-Land Data Analysis

23 Districts 123 Months Spatial Clustering Detected

📊 Key Statistics (2020-2024)

Moran's I (Snowfall)

0.42
p = 0.003 ⭐ Significant

Districts Analyzed

23
Northern Pakistan

Snow Depth Loss

-9,241%
Nagar & Gilgit (Extreme)

Correlation (Snow-Energy)

r = 0.89
Snowfall ↔ Snowmelt

Summary Statistics (2020-2024)

Variable Mean Std Dev Min Max Unit
Snow Cover -1.80 0.91 -3.93 -0.10 %
Snowfall -1.00 0.54 -2.09 -0.10 m × 10⁴
Snowmelt -0.66 0.47 -1.84 0.07 m × 10⁴
Latent Heat Flux -67,700 28,760 -117,299 5,062 J/m²

📈 Interactive Visualizations

District Rankings: Snowmelt Change

Temporal Trends (2016-2025)

Correlation Matrix

🔍 Key Findings

🌍 Spatial Clustering

  • Snowfall: I = 0.42 (p < 0.01)
  • Snowmelt: I = 0.21 (p < 0.05)
  • Pattern: Not random distribution

❄️ Extreme Changes

  • Nagar: -9,241% snow depth
  • Gilgit: -9,241% snow depth
  • Diamir: -2,241% snow depth

⚡ Energy-Snow Coupling

  • Snowfall-Snowmelt: r = 0.89
  • Latent-Solar: r = -0.89
  • Model R²: 0.25

📊 Period Comparison

  • 2015-2019: Western high snowfall
  • 2020-2024: Eastern shift
  • Change: Pattern reversal

🚨 Priority Districts

  • Gilgit: Extreme loss + hotspot
  • Hattian Bala: -1.84 × 10⁴ snowmelt
  • Muzaffarabad: High vulnerability

💧 Water Resources

  • 80%: Reduced snowmelt
  • Impact: Indus River system
  • Risk: Agricultural water stress

LISA Hotspots (p < 0.05)

District Cluster Type p-value Interpretation
Lower Dir (Temergara) High-High 0.015 High snowmelt, high neighbors
Upper Dir High-High 0.028 High snowmelt, high neighbors
Diamir (At Chilas) High-Low 0.022 Isolated high snowmelt
Muzaffarabad Low-Low 0.003 Low snowmelt, low neighbors
Neelum Low-Low 0.023 Low snowmelt, low neighbors