A 41-year (1979–2019) passive microwave derived lake ice phenology data record of the Northern Hemisphere
- 1Jiangsu Provincial Key Laboratory of Geographic Information Science and Technology, Key Laboratory for Land Satellite Remote Sensing Applications of Ministry of Natural Resources, School of Geography and Ocean Science, Nanjing University, Nanjing, China
- 2Department of Geography and Environmental Management, University of Waterloo, Ontario, Canada
- 3H2O Geomatics Inc., Waterloo, Ontario, Canada
- 1Jiangsu Provincial Key Laboratory of Geographic Information Science and Technology, Key Laboratory for Land Satellite Remote Sensing Applications of Ministry of Natural Resources, School of Geography and Ocean Science, Nanjing University, Nanjing, China
- 2Department of Geography and Environmental Management, University of Waterloo, Ontario, Canada
- 3H2O Geomatics Inc., Waterloo, Ontario, Canada
Abstract. Seasonal ice cover is one of the important attributes of lakes at middle and high latitude regions. The annual freeze-up and break-up dates and the durations of ice cover (i.e., lake ice phenology) are sensitive to the weather and climate, and hence can be used as an indicator of climate variability and change. The Calibrated Enhanced Resolution Brightness Temperature (CETB) dataset available from the National Snow and Ice Data Center (NSIDC) provides an alternate source of passive microwave brightness temperature (TB) measurements for the determination of lake ice phenology on a 3.125 km grid. This study used Scanning Multi-channel Microwave Radiometer (SMMR), Special Sensor Microwave Image (SSM/I) and Special Sensor Microwave Imager/Sounder (SSMIS) data from the CETB dataset to extract the ice phenology for 56 lakes across the Northern Hemisphere from 1979 to 2019. According to the differences in TB between lake ice and open water, a threshold algorithm based on the moving t test method was applied to determine the lake ice status for grids located at least 6.25 km away from the lake shore, and the ice phenology dates for each lake were then extracted. When ice phenology could be extracted from more than one satellite over overlapping periods, results from the satellite offering the largest number of observations were prioritized. The lake ice phenology results showed strong agreement with an existing product derived from Advanced Microwave Scanning Radiometer for EOS (AMSR-E) and Advanced Microwave Scanning Radiometer 2 (AMSR-2) data (2002 to 2015), with mean absolute errors of ice dates ranging from 2 to 4 days. Compared to near-shore in-situ observations, the lake ice results, while different in terms of spatial coverage, still showed overall consistencies. The produced lake ice record also displayed significant consistencies when compared to a historical record of annual maximum ice cover of the Laurentian Great Lakes of North America. From 1979 to 2019, the average complete freezing duration and ice cover duration for lakes forming a complete ice cover on an annual basis were 153 and 161 days, respectively. The lake ice phenology dataset – a new climate data record (CDR) – will provide valuable information to the user community about the changing ice cover of lakes in the last four decades. The dataset is available at https://www.pangaea.de/tok/c8fc0eab3d30777fc38979ad514217b6b7e86a65 (Cai et al., 2021).
Yu Cai et al.
Status: final response (author comments only)
-
CC1: 'Comment on essd-2021-435', Yubao Qiu, 26 Feb 2022
Lake ice is good indicator for the regional changes of climate or enviromental or even the water and energy exchanges. Author would like to extract the big lake (as we know there are much more records also pulished as dataset for the large lake or even small lakes) , the method and imput data are not so much difference with the previous develoment, so encourage the author articulate the value that this added to the dataset.
Others questions:
1) Why not engage the AMSR2, as this is really a good record for the lake observations?
2)The input brightness temperature wasnot calibrated at the same basis, so how do you use the threshold (interuptive changes detector), which may induce the errors to your decision of dates.
3) the optical image would be much better for the validation to the ice cover, while compared with Du etc.(2017) result, this is also from the passive microwave remote sensing, this introduce the error to it.
4) The data opened include other parameters of lake ice, or looks only the start and end dates, while the passive microwave remote sensing (for those larger lakes) can contribute starting of freezing dates, and end of freezing dates also.Highly appreicated to include the last development in this area, or may be compared with existing records, it would be giving more value to this work. Thanks.
-
AC1: 'Reply on CC1', Yu Cai, 10 May 2022
Thank you for the comments.
To our knowledge, there are three published lake ice phenology datasets based on passive microwave data. Du et al. (2017) obtained lake ice phenology of 71 lakes from 2003 to 2015 based on AMSR-E/2 data; Su et al. (2021) obtained lake ice phenology of 22 lakes from 1979 to 2018 based on SMMR and SSM/I-SSMIS data at a grid spacing of 25 km; Wang et al. (2021) obtained lake ice phenology of 409 lakes from 1978 to 2020 and 344 lakes from 2002 to 2020 based on SMMR (1978-1987, 6.25 km), SSM/I (1987-2002, 6.25km), AMSR-E (2002-2011, 6.25 km), MWRI (2011-2012, 50×30 km) and AMSR2 (2012-2020, 22×14 km).
Therefore, our dataset has a longer time period than Du's product and contains more lakes than Su's product. In comparison with Wang's product, which used only one central pixel from multi-source passive microwave data at different grid spacings for each lake, and used visual interpretation to obtain lake ice phenology, our dataset considered all the pixels 6.25 km away from the lake shore, and the lake ice phenology were extracted by an automatic algorithm based on the calibrated enhanced resolution passive microwave (CETB) dataset. In this case, we think that Wang's product provides precious records for small lakes, but the larger the lake, the worse the representation of the single pixel. In the contrast, our dataset is more suitable for the analysis of spatiotemporal changes in large lakes.
References:
Du, J., Kimball, J. S., Duguay, C., Kim, Y. and Watts, J. D.: Satellite microwave assessment of Northern Hemisphere lake ice phenology from 2002 to 2015, Cryosphere, 11(1), 47–63, doi:10.5194/tc-11-47-2017, 2017.
Su, L., Che, T. and Dai, L.: Variation in ice phenology of large lakes over the northern hemisphere based on passive microwave remote sensing data, Remote Sens., 13(7), doi:10.3390/rs13071389, 2021.
Wang, X., Qiu, Y., Zhang, Y., Lemmetyinen, J., Cheng, B., Liang, W., and Leppäranta, M.: A lake ice phenology dataset for the Northern Hemisphere based on passive microwave remote sensing, 00, 1–19, https://doi.org/10.1080/20964471.2021.1992916, 2021.
The followings are point-to-point response to the comments.
- SMMR, SSM/I-SSMIS data from CETB dataset provide continuous brightness temperature records with consistent grid resolution, the obtained lake ice phenology results are more comparable and suitable for time series analysis. Moreover, a lake ice phenology product based on AMSR-E/2 data has already been published (Du et al., 2017). Therefore, if users prefer to use the lake ice phenology from AMSR2 data, they can replace the results after 2012 with the results from Du's product. Furthermore, we compared the lake ice phenology from the two products in the manuscript (see Section 3.3.2), providing users with a reference for consistency and bias in using the two datasets.
- The brightness temperature of SSM/I and SSMIS data from the CETB dataset were calibrated, while the SMMR data were not cross-calibrated with the SSM/I-SSMIS data. However, we used unique threshold for each pixel from each satellite during the extraction of lake ice phenology (see Line 173). Therefore, whether or not the brightness temperature data were calibrated would not affect the extraction of lake ice phenology.
- We did not compare the lake ice phenology results with optical data because there are currently no lake ice phenology datasets in the northern hemisphere based on optical data that directly provide calendar date results. Instead, in addition to the AMSR-E/2 product, we also compared the lake ice phenology results with GLRIPD records and GLERL ice cover records. We believe that in-situ observations and ice charts can better represent the actual situation of lake ice status, and the comparisons with these records are more meaningful. Example of day-to-day comparison of ice coverage were also presented in the section of comparisons with GLERL ice cover records.
- We provide freeze-up start date, freeze-up end date, break-up start date, break-up end date, complete freezing duration, and ice cover duration in the dataset. For seven intermittently ice-covered lakes, annual maximum ice cover are also provided. We have made detailed description in the Data Availability section (see Lines 413-416), you can also open the link to look over the dataset (https://doi.org/10.1594/PANGAEA.937904).
-
AC1: 'Reply on CC1', Yu Cai, 10 May 2022
-
RC1: 'Comment on essd-2021-435', Tamlin Pavelsky, 10 Mar 2022
Review
A 41-year (1979-2019) passive microwave derived lake ice phenology data record of the Northern Hemisphere
By Cai, Duguay, and Ke
For Earth System Science Data
By Tamlin Pavelsky, University of North Carolina
Summary
This manuscript presents a new dataset of lake ice phenology for several dozen large lakes in the Northern Hemisphere based on passive microwave data. This dataset, based on a new passive microwave dataset beginning in 1979, represents the most comprehensive passive microwave lake ice phenology dataset to date. The dataset is compared against in situ observations for several lakes, including the Laurentian great lakes, and another passive microwave dataset from the AMSR-E and AMSR-2 instruments and generally found to be consistent with these other datasets, though with greater consistency for larger lakes at higher latitudes.
Overall Review
Fundamentally, the dataset presented in this paper is likely to be useful to the scientific community. I recommend some alterations to the paper and to the dataset itself, which I believe will improve it. However, I believe it should likely be published after consideration of reviewer comments. The strengths of the dataset are strongly related to the strengths of passive microwave remote sensing in general. There is a long record of global data, it is not impeded by clouds or other atmospheric effects, and there is a considerable literature suggesting that it can be used to detect ice status in large lakes. The dataset presented here is quite clearly the most comprehensive passive microwave lake ice phenology dataset, and it is likely to be of use to researchers interested in ice phenology patterns, including as a point of comparison for datasets collected using other methods (e.g. ground-based surveys, optical remote sensing, active-source radar).
There are 2 primary weaknesses that I believe should be addressed before the dataset and paper are finally published:
- The authors refer multiple times to uncertainty in ice phenology for various lakes, but there is no uncertainty field included in the dataset itself for any of the data fields. It is probably impossible to include all sources of uncertainty, but that should not stop the authors from include those sources they can quantify. I have some experience in this area with optical datasets, and what we have generally included are gaps between viable observations (see, for example, Pavelsky and Smith, 2004 and Zhang et al., 2021). These gaps are also clearly present in the passive microwave data (though for different reasons) and are discussed throughout the paper. The authors should at least be able to represent this source of uncertainty. However, they may also be able to represent other sources related to lake size. For example, there are several lakes that are represented by only one passive microwave pixel. In these cases, ice flagging is binary and likely to be less accurate. In any case, I would like to suggest that the authors quantify uncertainty on all dates as fully as possible and include those estimates in the dataset. Of course this also entails including a description in the paper of what sources of uncertainty are included in the flags.
- There are a number of places in the paper where statements are made that may be true but which are not supported by any data or analysis in the paper. These include:
- Line 308: “The main reason for the difference between lake ice dates from in-situ observations and passive microwave is their different observation ranges. In-situ records rely on observations of lake ice status visible from lake shores by human observers, while passive microwave satellites record TB from the entire lake surface (here within the pre-defined buffer).” While this seems reasonable to me, there is no evidence in the paper that it is true, and no other work is cited.
- Line 336: “Therefore, AMSR-E/2 data can capture more information near the lake shore than SMMR, SSM/I & SSMIS data, which led to the directional differences between the lake ice phenology dates extracted from the two datasets.” Same as above.
- Line 355: “The correlations between the ice cover for Huron and Ontario were lower than that of the other three lakes, which is also because the ice cover of these two lakes were usually small, and the ice first forms near shore which may not be covered in the set buffer (Figure 5).” This assessment potentially conflicts with the assessment shown in the next paragraph indicating that problems with the MTT algorithm my result in lower accuracies for Lake Ontario.
- Line 435: “The differences between the in-situ observations and the lake ice phenology extracted in this study were mainly due to the different fields of view of human observers versus satellite instruments.” Again, while this is plausible, no evidence of this explanation is explicitly shown in the paper.
I would strongly recommend that the authors either qualify these statements, provide evidence for them, or remove them. I would tend to hope for one of the two former options, as long as more concrete evidence suggests that they are correct.
Other than these two substantial areas for potential improvement, all other comments (listed below) are minor.
Specific Points
Line 14: I would write “. . . in middle and high latitude regions.”
Line 17: An alternate to what?
Line 28: I would write “consistency” rather than “consistencies” in both cases in this line.
Lines 39-40: I would recommend citing a few papers here, as there are a number of papers that have looked at this. For example, Smejkalova et al., 2016; Sharma et al., 2016.
Line 42: Might also consider citing Knoll et al., 2019 here.
Line 50: I might change the wording here, as a reader could be confused into thinking that all 865 site records begin in 1443, when most of them begin much latter.
Line 72: I might, again, cite the Smejkalova et al. 2016 paper here.
Line 101: I would probably delete “Satellite” before Nimbus-7.
Line 104: I would replace “data of the” with “data from the”
Line 108: “is against” should be “is compared against” or similar.
Figure 1: perhaps it’s unnecessary to do, but I believe there are a number of Patagonian lakes that might be close to the necessary size threshold for inclusion in this study. If the authors have not looked at these lakes to see if they are viable, I would recommend doing so given the paucity of similar records in the southern hemisphere.
Line 141: Typo at the end of the line. Should be “For comparison with lake ice phenology. . .”
Line 144: I would write “on 7 lakes” instead of “of 7 lakes.”
Line 149: I would write “from the CETB dataset.”
Line 155: I would write “from the National. . .”
Line 184: I would write “For the remaining pixels. . .”
Line 185: Why were thresholds of 5% and 95% chosen? Is there any sensitivity compared to say, 80-90% and 10-20%?
Line 194: This choice to omit ice covered periods of <30 days is significant. Would you then say that your estimates of ice duration likely underestimate total ice duration, as they ignore any intermittent ice cover occurring during the breakup or freezeup seasons? If so, I would explicitly mention this.
Line 228: I would write “For years with multiple lake ice records. . .”
Lines 246-249: A minor point, but given the size of these lakes, I’m not sure there’s any need to include the fourth significant digit in latitude. All are nearly 0.5 degree (or more) in N-S extent.
Line 270: I would recommend making sure you stay in consistently the present or past tense for this sentence and the following one.
Line 276: I would write “periodically missing data” instead of “periodical data missing.”
Line 301: I would write “compared” rather than “used to compare”
Line 302: I would write “complete ice cover” rather than “completely ice covered.”
Lines 313-316: while I agree with the statement in this sentence to a degree, I would argue that the agreement with GLRIPD records is not particularly strong in many cases. As such, I’m not sure I agree that the analysis presented here provides strong evidence for this statement. Rather, I would say that remotely sensed observations can complement in situ measurements.
Lines 335-336: Make sure to keep verb tenses the same in this sentence.
Line 429: The conclusion that lakes at low latitudes and/or small areas tend to have larger uncertainties would be much more robust if there were a more consistent uncertainty quantification, as mentioned above.
References
Knoll, L. B., Sharma, S., Denfeld, B. A., Flaim, G., Hori, Y., Magnuson, J. J., ... & Weyhenmeyer, G. A. (2019). Consequences of lake and river ice loss on cultural ecosystem services. Limnology and Oceanography Letters, 4(5), 119-131.
Pavelsky, T. M., & Smith, L. C. (2004). Spatial and temporal patterns in Arctic river ice breakup observed with MODIS and AVHRR time series. Remote sensing of environment, 93(3), 328-338.
Sharma, S., Magnuson, J. J., Batt, R. D., Winslow, L. A., Korhonen, J., & Aono, Y. (2016). Direct observations of ice seasonality reveal changes in climate over the past 320–570 years. Scientific Reports, 6(1), 1-11.
Šmejkalová, T., Edwards, M. E., & Dash, J. (2016). Arctic lakes show strong decadal trend in earlier spring ice-out. Scientific reports, 6(1), 1-8.
Zhang, S., Pavelsky, T. M., Arp, C. D., & Yang, X. (2021). Remote sensing of lake ice phenology in Alaska. Environmental Research Letters, 16(6), 064007.
- AC2: 'Reply on RC1', Yu Cai, 10 May 2022
-
RC2: 'Comment on essd-2021-435', J. Uusikivi, 12 Apr 2022
Overall review
The authors have used an automated method to extract ice phenology data from passive microwave data. The data set presented here and explained in the article is generally very interesting and will be useful for research community. The data set is likely the longest and most comprehensive ice phenology data set from satellite-based observations for that large number of lakes. This data covers multiple climatological areas and lake sizes and is therefore well worth publication. Data set is usable in the present format.
Comments:
I would like to present 2 recommendations to improve the usability of the data and the manuscript.
- Data set does not include any sort of error estimates for dates, duration, or maximum ice cover area. In the manuscript is long discussion on the errors and their possible sources, but these should be quantified in the data, or at least in the manuscript. It is very difficult to compare this data set to other similar data sets without this information. In the manuscript one major target for this data is climate research, it is difficult to draw conclusion if error marginals are unknown. To use this data to complement data gaps of in situ archives of ice phenology, more precise definition of the errors and their sources compared to the GRLIPD and GLERL ice cover data sets should be included.
- Using 37 GHz H-polarized data has some limitations in distinguishing ice and open water. Signal can be strongly affected by open water surface roughness from wind (for example, K.-K. Kang et al.: Estimating ice phenology on large northern lakes from AMSR-E; doi:10.5194/tc-6-235-2012). This problem and its implications to the data is not discussed in the manuscript at all, and it is not covered in any of the references provided. By discussing this matter or providing references that discuss this, will make this data much more reliable and usable.
I also have some minor comments to consider:
- on line 176: “When the lake is water covered, the TB for land-contaminated pixels will be higher than that of a pure pixel, while when the lake is ice covered, the TB will be lower than that of pure pixel.” Last 2 words: Is it pure pixel of ice/water/land?
- on line 271: “When the lake area was large enough, the gradual freeze-up or break-up within the pixel can be ignored, but for small lakes, it may lead to certain deviations in the lake ice phenology results.” What are the certain deviations?
- on line 312: “Overall, if the overlapping time between the two dataset was longer, the lake ice dates could show a higher consistency.” How or Why that could be the case?
- on line 353:” This is because a buffer of 6.25 km was used to exclude pixels near the lake shore, which happens to be the place where lake ice forms first.” If this is the only explanation in the difference between GLERL data and this data set, one would expect the difference to gradually wannish as one nears 100% ice coverage. This is not the case in all the lakes in all the years. Why is that?
- AC3: 'Reply on RC2', Yu Cai, 10 May 2022
Yu Cai et al.
Data sets
Lake ice phenology in the Northern Hemisphere extracted from SMMR, SSM/I and SSMIS data from 1979 to 2020 Cai, Yu; Duguay, Claude R.; Ke, Chang-Qing https://www.pangaea.de/tok/c8fc0eab3d30777fc38979ad514217b6b7e86a65
Yu Cai et al.
Viewed
HTML | XML | Total | BibTeX | EndNote | |
---|---|---|---|---|---|
324 | 112 | 20 | 456 | 6 | 10 |
- HTML: 324
- PDF: 112
- XML: 20
- Total: 456
- BibTeX: 6
- EndNote: 10
Viewed (geographical distribution)
Country | # | Views | % |
---|
Total: | 0 |
HTML: | 0 |
PDF: | 0 |
XML: | 0 |
- 1