diff --git a/experiment/meteo_france_data/scm_models_data/abstract_study.py b/experiment/meteo_france_data/scm_models_data/abstract_study.py
index d1180ec8316e7c07de295129c90128f17546b8ca..173cfa44dada9acfb85d0b91db49ffc2db7448f3 100644
--- a/experiment/meteo_france_data/scm_models_data/abstract_study.py
+++ b/experiment/meteo_france_data/scm_models_data/abstract_study.py
@@ -1,5 +1,8 @@
 import datetime
-from matplotlib.patches import Polygon
+
+from matplotlib.lines import Line2D
+from matplotlib.markers import MarkerStyle
+from matplotlib.patches import Polygon, Patch
 import io
 import os
 import os.path as op
@@ -396,6 +399,19 @@ class AbstractStudy(object):
         if axis_off:
             plt.axis('off')
 
+        # Add legend for the marker
+        if massif_name_to_marker_style is not None:
+            labels = ['\mathcal{M}_{\mu_1}', '\mathcal{M}_{\sigma_1}', '\mathcal{M}_{\mu_1, \sigma_1}']
+            markers = ["s", "^", "D"]
+            legend_elements = [
+                Line2D([0], [0], marker=marker, color='w', label='${}$'.format(label),
+                           markerfacecolor='w', markeredgecolor='k', markersize=8)
+                for label, marker in zip(labels, markers)
+            ]
+            ax.legend(handles=legend_elements, bbox_to_anchor=(0.01, 0.03), loc='lower left')
+            ax.annotate("Filled symbol = significant trend ", xy=(0.05, 0.015), xycoords='axes fraction',
+                         fontsize=7)
+
         if show:
             plt.show()
 
diff --git a/experiment/paper_past_snow_loads/result_trends_and_return_levels/main_result_trends_and_return_levels.py b/experiment/paper_past_snow_loads/result_trends_and_return_levels/main_result_trends_and_return_levels.py
index 4a8827d1bd1abddb3101b4dd6365a6b9017bce67..f86d97aac1770497c5eb409030435ed2565d9bc6 100644
--- a/experiment/paper_past_snow_loads/result_trends_and_return_levels/main_result_trends_and_return_levels.py
+++ b/experiment/paper_past_snow_loads/result_trends_and_return_levels/main_result_trends_and_return_levels.py
@@ -58,9 +58,9 @@ def major_result():
 
 
 if __name__ == '__main__':
-    major_result()
+    # major_result()
     # intermediate_result(paper_altitudes)
-    # minor_result(altitude=600)
+    minor_result(altitude=600)
     # intermediate_result(altitudes=[1500, 1800], massif_names=['Chartreuse'],
     #                     uncertainty_methods=[ConfidenceIntervalMethodFromExtremes.ci_mle,
     #                                          ConfidenceIntervalMethodFromExtremes.ci_bayes],