Indeed, the novel aspects of this study are somewhat hidden due to a heavier focus on the methods and results. We have added a short section titled "7.3 Significance of results" (page 17 line 32) to summarize the novel contributions of this paper. In this section we aim to highlight that we are not the first to use a LOWESS interpolation on MODIS time series data, however, this is the first study to apply it to the MODIS snow cover product. In addition the LOWESS allows us to detect the start and end dates of snow, which is an important contribution to our understanding of snow cover in British Columbia. We suggest that to extend this work further into the past, one could use lower resolution remote sensing data, or in-situ observations, and we include the provided references of (McClung, 2013, Barton 2017). These methods could likely be implemented in other regions, however the dates of the hydrological year will likely be different, and so could the optimal NDSI threshold. Finally, we highlight that the 500 m rasters that we have produced for this study could be of interest to other fields of science, and that our results could likely improve seasonal forecasting of snow. Here is the new section in its entirety: (7.