Snow and Ice Monitoring from Space: How NDSI Maps the Cryosphere

Why Snow Mapping Matters

Snow cover affects everything from water supply forecasting to climate modeling. Roughly one-sixth of the world's population depends on snowmelt for freshwater. Accurate snow extent mapping from satellites enables:

• Snowmelt runoff prediction — How much water will reservoirs receive in spring?
• Flood early warning — Rapid snowmelt combined with rain causes devastating floods
• Climate monitoring — Snow cover duration is a sensitive indicator of warming trends
• Ski industry and tourism — Snow depth and extent affect seasonal economies
• Transportation planning — Mountain road closures depend on snowfall and melt patterns

The Spectral Trick Behind NDSI

Snow is white — it reflects nearly all visible light. But in the shortwave infrared (SWIR), snow is surprisingly dark. This combination is unique among natural surfaces:

Surface	Visible (Green)	SWIR	NDSI
Fresh snow	Very high (~0.9)	Very low (~0.05)	+0.8 to +0.9
Old/dirty snow	High (~0.6)	Low (~0.15)	+0.5 to +0.7
Clouds	High (~0.7)	High (~0.5)	-0.1 to +0.2
Bare rock	Moderate (~0.3)	Moderate (~0.3)	-0.1 to +0.1
Vegetation	Moderate (~0.1)	Low (~0.15)	-0.3 to +0.1
Water	Low (~0.05)	Very low (~0.01)	+0.4 to +0.8 ⚠️

The formula:

NDSI = (Green - SWIR) / (Green + SWIR)

For Sentinel-2: (B3 - B11) / (B3 + B11)

The Cloud Discrimination Advantage

This is NDSI's most valuable property. Simple brightness-based snow detection fails because clouds are also white. But clouds are bright across all wavelengths, while snow drops dramatically in SWIR.

• Snow: Green ≈ 0.9, SWIR ≈ 0.05 → NDSI ≈ +0.9
• Cloud: Green ≈ 0.7, SWIR ≈ 0.5 → NDSI ≈ +0.2

With a threshold of 0.4, snow and clouds are clearly separated. This is the same principle used in NASA's MODIS snow products (MOD10A1) and is the standard approach for global snow mapping.

Practical Workflow

1. Select Imagery

Search for Sentinel-2 imagery over mountainous or high-latitude areas during winter or spring. Look for:

• Scenes with partial snow cover (fully snow-covered scenes are less interesting)
• Minimal cloud contamination
• Good sun illumination (low sun angles in winter can cause extensive shadows)

2. Apply NDSI

Switch to NDSI visualization. Snow-covered areas will appear bright (high NDSI), while snow-free areas appear dark.

3. Identify the Snowline

The transition zone where NDSI drops below the threshold (~0.4) marks the snowline — the elevation where snow cover begins. This boundary is a critical parameter for hydrological models.

4. Compare Dates

Load imagery from different dates to track:

• Snowline retreat during spring melt (snowline moves uphill)
• Snow accumulation during winter storms (snowline moves downhill)
• Year-to-year differences at the same date (climate indicators)

Common Pitfalls

Water Bodies Look Like Snow

Water has high NDSI values because it reflects some green light but absorbs SWIR almost completely. Without correction, lakes and rivers will be classified as "snow." This is the same spectral property exploited by NDWI for water detection.

Solution: Apply an additional NIR brightness threshold. Snow reflects strongly in NIR, while water absorbs it. Requiring NIR reflectance > 0.11 (the MODIS standard) effectively masks water.

Shadows Reduce NDSI

Mountain shadows reduce all reflectance values, pushing shadowed snow below the NDSI threshold even though snow is present.

Solution: Lower the threshold in shadowed areas, or use a shadow mask derived from a DEM and the sun angle. This is one of the most challenging aspects of accurate snow mapping in mountainous terrain.

Dirty or Aged Snow

As snow ages, dust, soot, and melt-freeze cycles reduce its visible reflectance while SWIR absorption remains strong. This lowers NDSI values, potentially below the detection threshold.

Solution: Use a lower threshold (0.3 instead of 0.4) for aged snow detection. The tradeoff is increased false positives from cloud edges and bright soil.

Glacier Ice vs Snow

Glacier ice has lower NDSI than fresh snow because ice absorbs more visible light (especially blue). Bare glacier ice may produce NDSI values of 0.3-0.5 — right at the threshold boundary.

Solution: For glacier mapping, consider combining NDSI with NIR/SWIR band ratios that better distinguish ice from snow.

Applications Beyond Snow Cover

Glacier Monitoring

Track glacier retreat by comparing NDSI maps from the same month across multiple years using difference analysis. The shrinking of persistent snow/ice areas reveals glacier loss.

Permafrost Indicators

Snow cover duration affects ground temperature and permafrost stability. Areas where snow cover is decreasing (shorter snow season) may experience permafrost degradation.

Albedo and Climate Feedback

Snow has high albedo (reflects sunlight back to space). When snow cover decreases, the darker underlying surface absorbs more solar energy, causing further warming. NDSI time-series help quantify this feedback loop.

Snow Cover Trends: What Satellites Reveal

Long-term satellite records (MODIS since 2000, Landsat since 1984) provide an unambiguous picture of cryosphere change. Key findings relevant for interpreting NDSI time series:

Region	Snow/Ice Trend (2000–2024)	Snowmelt Timing Shift	Climate Implication
Western US mountain snowpack	−10 to −25% area loss	1–3 weeks earlier spring melt	Colorado River system water security
European Alps glaciers	−2 to −3% area loss per decade	2–4 weeks shorter snow season	Tourism, hydropower impacts
Tibetan Plateau (high elevation)	Modest gain at highest elevations	Complex; elevation-dependent	Partially offset by increased precipitation
Arctic sea ice extent	−13% per decade (September minimum)	Seasonal-only: near ice-free summers by 2050?	Albedo feedback, shipping routes
Siberian snow cover duration	−5 to −10 days per decade	Earlier spring onset	Permafrost warming
Hindu Kush–Himalaya snowpack	−6% per decade (spring snow extent)	1–2 weeks earlier peak melt	Water security for 240 million people

NDSI accuracy context: MODIS MOD10A1 snow product achieves 90–97% overall accuracy under clear skies. Accuracy drops to 60–80% in dense forest (canopy obscures snow), and to near-zero under cloud cover (clouds must be masked before NDSI analysis). For the Alps and similar topographically complex terrain, shadow masking is the largest remaining source of error — improperly shadowed snow can reduce mountain NDSI accuracy by 5–15%.

Practical snowline detection accuracy: With Sentinel-2 at 20m resolution, snowline elevation can be mapped with ±50–100m vertical accuracy (corresponding to 100–200m horizontal distance on a typical alpine slope). This precision is sufficient for hydrological modeling but not for detailed glacier boundary mapping (which requires < 5m horizontal accuracy, achievable with commercial imagery or aerial LiDAR).

Combining NDSI with SAR

SAR imagery responds to snow differently than optical sensors:

• Dry snow is largely transparent to C-band SAR — the radar "sees through" it to the ground
• Wet snow strongly absorbs SAR energy, appearing very dark

This means SAR can detect wet snow (melt conditions) even under cloud cover. Combining NDSI (total snow extent) with SAR (wet snow detection) provides a more complete picture of snow conditions than either sensor alone.

Try It in Off-Nadir Delta

Off-Nadir Delta includes NDSI as a built-in visualization for Sentinel-2 data:

1. Search for Sentinel-2 imagery over a mountain range in winter or spring
2. Apply NDSI visualization — snow appears bright, everything else dark
3. Compare with true-color to verify the snow boundary
4. Load multiple dates to track snowline changes through the melt season
5. Overlay with Sentinel-1 SAR to identify wet snow zones