Snow Cover Dynamics with MODIS and DeepSeek

// Define the region of interest (North Morocco).
var geometry = ee.Geometry.Rectangle([-6, 32, -2, 35]);

// Load the MODIS MOD10A1 snow cover dataset.
var snowCover = ee.ImageCollection('MODIS/061/MOD10A1')
  .filterDate('2020-01-01', '2023-12-31')
  .filterBounds(geometry);

// Filter for December, January, and February.
var winterSnow = snowCover.filter(ee.Filter.calendarRange(12, 2, 'month'));

print('MODIS Snow Cover Collection:', winterSnow);

// Calculate the mean NDSI snow cover for each year.
var annualSnow = ee.List.sequence(2020, 2023).map(function(year) {
  var startDate = ee.Date.fromYMD(year, 12, 1);
  var endDate = ee.Date.fromYMD(ee.Number(year).add(1), 2, 28);
  var annualImage = winterSnow.filterDate(startDate, endDate)
    .mean()
    .clip(geometry);
  return annualImage.set('year', year);
});

annualSnow = ee.ImageCollection.fromImages(annualSnow);

// Define visualization parameters.
var snowVisParams = {
  min: 0.0,
  max: 100.0,
  palette: ['00008B', '0000FF', 'ADD8E6', 'FFFFFF']
};

// Display the annual mean NDSI snow cover as a map.
annualSnow.sort('year').iterate(function(image) {
  var year = image.get('year');
  Map.addLayer(image, snowVisParams, 'Mean Snow Cover ' + year);
});

// Export the mean snow cover image to Google Drive.
Export.image.toDrive({
  image: annualSnow.first(),
  description: 'Mean_Snow_Cover_2020',
  folder: 'GEE_Exports',
  scale: 500,
  region: geometry,
  maxPixels: 1e13
});