diff --git a/python/s2snow/snow_annual_map.py b/python/s2snow/snow_annual_map.py
index d868dbde1d1c62015d1a97b0c414e7b781eb09fc..bd93fcb68b172ddb9085f9a671f16f0065869b5d 100644
--- a/python/s2snow/snow_annual_map.py
+++ b/python/s2snow/snow_annual_map.py
@@ -102,7 +102,7 @@ def compute_SCD(binary_snow_mask_list, multitemp_cloud_vrt, input_dates_filename
     expression = "+".join(["im1b" + str(i) for i in band_index])
     logging.debug("expression: {}".format(expression))
     bandMathApp = band_math([gapfilled_timeserie],
-                            snow_coverage_duration,
+                            snow_coverage_duration  + "?&gdal:co:COMPRESS=DEFLATE",
                             expression,
                             ram,
                             otb.ImagePixelType_uint16)