the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Impact of aeolian sediments on the ice-covered lake surface temperature of Zonag Lake, Tibetan Plateau
Abstract. Aeolian sediments on ice-covered lakes can modify surface energy conditions and potentially influence the thermal environment of lake ice. However, the quantitative relationship between aeolian sediments and lake ice surface temperature remains insufficiently understood. This study investigates the relationship between aeolian sediments and ice-covered lake surface temperature (LST) in Zonag Lake on the Tibetan Plateau using multi-temporal Landsat 8 OLI imagery and MOD11A1 LST data from 2013 to 2018. The Normalized Difference Sandy Land Index (NDSLI) was applied to identify the distribution of aeolian sediments on lake ice. To quantify the sediment impact on LST, spatial regression analyses were conducted using ordinary least squares (OLS), geographically weighted regression (GWR) and multi-scale geographically weighted regression (MGWR). The results show that NDSLI effectively captures the spatial distribution and migration patterns of aeolian sediments on lake ice. Sediment-covered areas (NDSLI < 0) consistently exhibit higher temperature than uncovered areas (NDSLI > 0), with observed temperature differences ranging from 0.83 to 2.36 °C. Finer-scale NDSLI classification indicates a progressive increase in surface temperature with increasing sediment accumulation. Spatial regression results reveal a significant positive relationship between sediment and LST, with MGWR providing the best model performance and explanatory capability among the tested models. This study provides quantitative evidence that aeolian sediment deposition can modify the thermal conditions of lake ice surfaces, offering new insights into cryosphere–aeolian interactions in arid-cold regions and improving the understanding of thermal variability in ice-covered lakes.
- Preprint
(1727 KB) - Metadata XML
- BibTeX
- EndNote
Status: final response (author comments only)
-
CC1: 'Comment on egusphere-2026-2129', Lei Zheng, 24 Jul 2026
-
AC1: 'Reply on CC1', Guangyin Hu, 18 Aug 2026
The comment was uploaded in the form of a supplement: https://egusphere.copernicus.org/preprints/2026/egusphere-2026-2129/egusphere-2026-2129-AC1-supplement.pdf
-
AC1: 'Reply on CC1', Guangyin Hu, 18 Aug 2026
-
CC2: 'Comment on egusphere-2026-2129', Ke Zhang, 01 Aug 2026
This is a well-conducted and timely study that provides quantitative evidence for the thermal perturbation of ice-covered lake surfaces by aeolian sediments on the Tibetan Plateau. The application of spatial regression models, particularly MGWR, to quantify sediment-induced thermal anomalies is appropriate and helps address an important knowledge gap regarding cryosphere–aeolian interactions. The results indicate that sediment-covered areas exhibit surface temperature increases of 0.83–2.36°C, highlighting the potential role of aeolian sediments as an important non-climatic driver of lake ice thermal heterogeneity. The manuscript is generally well organized, and the methodology appears robust. I have a few comments that I believe will further strengthen the manuscript.
Specific comments
1. The Landsat 8 and MODIS datasets are not always acquired on the same date. Given that aeolian sediment streaks can migrate rapidly under strong wind conditions, this temporal mismatch may introduce uncertainty into the spatial correlation analysis. The authors should discuss the likely migration rate or persistence of these sediment streaks and evaluate how such temporal offsets may affect the reported relationships. If available, daily meteorological observations from nearby weather stations could also be used to demonstrate that background thermal conditions remained relatively stable during the interval between the two satellite acquisitions, thereby helping to constrain this source of uncertainty.
2. The 2011 Zonag Lake outburst and the subsequent reorganization of the regional lake system, which resulted in extensive exposure of former lakebeds, represent a major environmental change in the study area. Although this event is briefly introduced in the study area description, its potential implications are not fully explored in the subsequent analysis. Since the newly exposed lakebed likely became the dominant source of aeolian sediments, it would be informative to examine whether the observed sediment–temperature relationship became stronger or more evident after 2011 compared with earlier periods. At a minimum, the authors should discuss this possibility in Section 5, as it would provide valuable environmental context for interpreting the results.
3. Line 77: A punctuation mark is needed between "conditions" and "Monitoring"because these are two independent sentences. The sentence should read:"...regional aeolian conditions. Monitoring data revealed..."
4. The first paragraph of Section 5.2 would fit better in the Introduction. It clearly outlines the scientific significance of the study and identifies the key research gap, providing a stronger rationale for the work at an earlier stage of the manuscript. Moving this paragraph to the Introduction would improve the overall flow and better frame the objectives of the study.
Citation: https://doi.org/10.5194/egusphere-2026-2129-CC2 -
AC2: 'Reply on CC2', Guangyin Hu, 18 Aug 2026
The comment was uploaded in the form of a supplement: https://egusphere.copernicus.org/preprints/2026/egusphere-2026-2129/egusphere-2026-2129-AC2-supplement.pdf
-
AC2: 'Reply on CC2', Guangyin Hu, 18 Aug 2026
-
CC3: 'Comment on egusphere-2026-2129', Ke Zhang, 01 Aug 2026
This is a well-conducted and timely study that provides quantitative evidence for the thermal perturbation of ice-covered lake surfaces by aeolian sediments on the Tibetan Plateau. The application of spatial regression models, particularly MGWR, to quantify sediment-induced thermal anomalies is appropriate and helps address an important knowledge gap regarding cryosphere–aeolian interactions. The results indicate that sediment-covered areas exhibit surface temperature increases of 0.83–2.36°C, highlighting the potential role of aeolian sediments as an important non-climatic driver of lake ice thermal heterogeneity. The manuscript is generally well organized, and the methodology appears robust. I have a few comments that I believe will further strengthen the manuscript.
Specific comments
1. The Landsat 8 and MODIS datasets are not always acquired on the same date. Given that aeolian sediment streaks can migrate rapidly under strong wind conditions, this temporal mismatch may introduce uncertainty into the spatial correlation analysis. The authors should discuss the likely migration rate or persistence of these sediment streaks and evaluate how such temporal offsets may affect the reported relationships. If available, daily meteorological observations from nearby weather stations could also be used to demonstrate that background thermal conditions remained relatively stable during the interval between the two satellite acquisitions, thereby helping to constrain this source of uncertainty.
2. The 2011 Zonag Lake outburst and the subsequent reorganization of the regional lake system, which resulted in extensive exposure of former lakebeds, represent a major environmental change in the study area. Although this event is briefly introduced in the study area description, its potential implications are not fully explored in the subsequent analysis. Since the newly exposed lakebed likely became the dominant source of aeolian sediments, it would be informative to examine whether the observed sediment–temperature relationship became stronger or more evident after 2011 compared with earlier periods. At a minimum, the authors should discuss this possibility in Section 5, as it would provide valuable environmental context for interpreting the results.
3. Line 77: A punctuation mark is needed between "conditions" and "Monitoring"because these are two independent sentences. The sentence should read:"...regional aeolian conditions. Monitoring data revealed..."
4. The first paragraph of Section 5.2 would fit better in the Introduction. It clearly outlines the scientific significance of the study and identifies the key research gap, providing a stronger rationale for the work at an earlier stage of the manuscript. Moving this paragraph to the Introduction would improve the overall flow and better frame the objectives of the study.
Citation: https://doi.org/10.5194/egusphere-2026-2129-CC3 -
AC3: 'Reply on CC3 (The same with the "Reply on CC2")', Guangyin Hu, 18 Aug 2026
The comment was uploaded in the form of a supplement: https://egusphere.copernicus.org/preprints/2026/egusphere-2026-2129/egusphere-2026-2129-AC3-supplement.pdf
-
AC3: 'Reply on CC3 (The same with the "Reply on CC2")', Guangyin Hu, 18 Aug 2026
-
RC1: 'Comment on egusphere-2026-2129', Anonymous Referee #1, 02 Sep 2026
Overall comment: The topic is interesting and the spatial patterns shown in the paper are worth investigating. My main concern is that the current analysis shows a strong spatial relationship between sediment distribution and LST, but this is sometimes presented as proof that sediments are causing the temperature differences. I think the paper will be much stronger if the authors explain the Landsat-MODIS matching clearly, provide more information about the satellite processing and regression models, address the uncertainty in the data, and make the conclusions more careful.
Major Comments
- I am concerned about how the Landsat and MODIS images were matched. Some of the images are several days apart. For example, the Landsat image from December 26, 2017 is matched with MODIS LST from December 20, 2017. There is a six-day difference. Since the paper shows that sediments can move across the lake surface with wind, how can we be sure that the sediment distribution six days later represents the conditions when the MODIS temperature was measured? Weather conditions could also change during those days. The authors should show the time difference for every Landsat-MODIS pair and explain why these differences are acceptable. If possible, the main analysis should use same-day observations.
- This part needs more explanation. Landsat NDSLI is at 30 m while MODIS LST is at 1 km. The paper says that MODIS LST was resampled to “higher resolution (1 km × 1 km),” but MOD11A1 is already 1km. How exactly were the two datasets combined? Was the 30 m NDSLI aggregated to 1 km? If so, was the mean NDSLI used, the percentage of sediment-covered pixels, or something else? If MODIS was resampled to 30 m, this would not create new temperature and could result in many Landsat pixels sharing the same MODIS temperature. The authors need to clearly explain the complete spatial matching and resampling process.
- The paper says MOD11A1 V6 was used, but important processing information is missing. Was daytime or nighttime LST used? What MODIS quality flags were applied? How were cloudy pixels removed? How were mixed shoreline pixels handled? What was the acquisition time? This is especially important because some of the reported temperature differences are less than 1 °C. The authors should also discuss the uncertainty of MODIS LST over snow- and ice-covered surfaces.
- The authors use NDSLI < 0 to identify sediments and NDSLI < -0.15 for stronger sediment accumulation. However, frozen lake surfaces can have different spectral conditions because of snow, bare ice, cracks, wet ice, snow depth, and surface roughness. How can we be sure that negative NDSLI values are caused mainly by sediments? The visual comparison in Figure 2 is useful, but I do not think it is enough to say that the sediment boundaries were “accurately extracted.” The authors should provide additional validation or test how the results change when different NDSLI thresholds are used or soften the language
- The text refers to Gao et al. (2025), including when discussing the previous validation of NDSLI. However, the reference list gives the related paper as Gao et al. (2026). Please check and make the citation year consistent.
- The authors report a 6.28 °C difference by comparing the maximum temperature in a sediment-covered area with the minimum temperature in a sediment-free area. I do not think this is the best way to estimate the sediment effect because these are two extreme values from different locations. The mean or median temperature difference between sediment and non-sediment areas is more useful. Confidence intervals or another measure of uncertainty would also help.
- Figure 7 shows that LST decreases as NDSLI increases, meaning that the relationship between NDSLI and LST is negative. Lower NDSLI means more sediment, so the relationship between sediment amount and LST is positive. These are not the same statement. The paper sometimes describes the relationship simply as a positive relationship. The authors should make this distinction clear throughout the paper.
- The paper mainly presents R², adjusted R², AICc, bandwidth, and sigma. However, the actual NDSLI regression coefficients are not shown. If the purpose of the model is to understand how NDSLI relates to LST, the coefficient is very important. Please report the coefficient values and their uncertainty or significance. For GWR and MGWR, maps showing the local coefficients would also be useful.
- The paper explains that MGWR is useful because different explanatory variables can operate at different spatial scales. However, the model appears to have only one main predictor, NDSLI. The authors should explain more clearly what MGWR adds in this particular case compared with GWR. The paper should also provide the kernel type, whether the bandwidth is fixed or adaptive, how bandwidth was selected, the distance method, and other important model settings.
- The manuscript repeatedly says that MGWR consistently performs better than GWR. However, this is not true for March 31, 2018. GWR has an adjusted R² of 0.442 and AICc of 355.202, while MGWR has an adjusted R² of 0.438 and AICc of 357.377. In this case, GWR performs slightly better. The authors should change “consistently outperformed” to something like “generally outperformed.”
- The manuscript repeatedly says that MGWR consistently performs better than GWR but that is not the case. Please change the language
- The authors do not explain how the spatial weights were created. Was queen contiguity, rook contiguity, distance, or k-nearest neighbours used? Were the weights standardized? How many permutations were used to determine significance? These details are necessary to reproduce the analysis. Also, the equations for Moran’s I should be checked carefully because Equation 2 does not appear to be written correctly.
- NDSLI appears to be the only main explanatory variable in the regression. However, lake ice temperature can also be influenced by air temperature, snow cover, radiation. Etc. The conclusion currently says that the results “robustly confirm” I think this should be changed to something more careful, such as saying that sediment coverage is strongly associated with the spatial variation in LST.
- The Methods say that images were selected when the lake was in a stable freezing period. However, the Results later say that one image represents a non-freezing condition and another represents snowfall conditions. If these images were intentionally included as comparison or control cases, please explain this clearly in the Methods.
Minor Comments
The Results move from 4.1 to 4.3, then 4.4, 4.5, and later back to 4.4. This is confusing.
The caption to Figure 8 reads OSL where OLS is intended.
L26: “Sediments accumulation” should be “sediment accumulation”. There are other places in the document with the same issue
L32: “Researches” should be “Research.”. Other parts of the document have this issue as well
L42: “Beside” should be “Besides.”
L47: “The” is missing. Change to: “At the Headwater Region of the Yellow River”
L50: Change “In Mars” to “On Mars.”
L31: “The” is missing in two places.
Change to: “...between aeolian sediments and the ice-covered period of the lake.”L77: A period is missing after “conditions.”
L189: Remove “was.” Change to: “NDSLI < 0 accounted for 46.1% of the total lake surface.”
L231: “These evidences indicate that the aeolian sediments accumulation significantly changes the thermal distribution of ice-covered lake surface...” Check this sentence
L249: “Exhibit” should be “exhibits,” and “a” is missing before “bimodal pattern.”
L267: “Only one date showed relatively good fitting result (R² = 0.567) is December 31, 2013, whereas other dates exhibited poor fitting results...” Check this sentence
L347: “Researches” should be “research.”
L408: “The authors declare that they have no conflict of interest..”
There is an extra period at the end.
There are some date inconsistencies
Table 1: Landsat 2017-11-08 / MODIS 2017-11-07
Figure 3(d): 2017-11-07
Figure 4(e): 2017-11-20
Section 4.4, Line 208: November 20, 2017The dates are not consistent. Please check whether November 20 is a typo or a different observation date that was not included in Table 1.
This is important because this date is linked to the reported Moran’s I value of 0.8327.
Citation: https://doi.org/10.5194/egusphere-2026-2129-RC1 -
RC2: 'Comment on egusphere-2026-2129', Anonymous Referee #2, 16 Oct 2026
1 Initial view
The article “Impact of aeolian sediments on the ice-covered lake surface temperature of Zonag Lake, Tibetan Plateau” describes the relationship between aeolian sediments and ice-covered lake surface temperature. The main goal is to examine the sediment impact on the lake surface temperature using spatial regression. The article demonstrates that the chosen index (NDSLI) captures the spatial distribution and migration patterns of aeolian sediments on the lake ice. The paper is well structured and very interesting to read. The proposed methodology is sensible, the spatial patterns of sediments on sea ice and their effects on surface albedo and its spatial variability is scientifically interesting and motivates future perspectives.
2 Conceptual overview
The researchers stated that MGWR provides best model performance among two other methods: OLS and GWR. The study contributes to the knowledge on the freeze-thaw cycle, supporting earlier findings on the positive linkage between sediments and surface temperature. The spatial pattern that aeolian sediment forms and drives the albedo in different ways, when comparing two observational datasets is itself very interesting and is the main focus of the discussion.
This paper addresses a relevant scientific question that fully aligns with the extent of TC interests. It offers a wide range of comparable methods, searches for the best fitting method to observe sediment spatial patterns and their impact on radiative properties of the surface. A previously unexplored regional relationships between the spatial patterns of aeolian sediments and surface temperature of the Zonag Lake are presented. The authors reach a clear conclusion, suggest possible future directions and consider their limitations. The scientific method of this study is valid, but could be clarified with greater care. The result is fully sufficient and supports the main conclusion. The interpretation is coherent but it is missing possible influence of physical drivers of sedimentation and sediment transport. The description of data retrieval and processing is unfortunately insufficient, but with greater detail it should be reproducible. The researchers follow on recent literature, contributing new measures to the field, which is directly stated in the title. The abstract fully grasps the articles’ structure, which is sufficient and clear as well. The language can partially be improved in structure and coherency. On the other hand, the mathematical equations, used symbols and abbreviations are well explained and properly used in the text. The general amount of description, provided visualizations and references in the text is appropriate.
3 Major issues
Starting in 3.1 on Data Sources, the data choice and matching should be explained in further detail:
- I am concerned about the quality of the LST product on sediment covered lake ice. Sediments will have a very different thermal emissivity than clean lake ice (which can also be seen in the differences in albedo) which might confuse the LST algorithm. Can you rule out, that sediments create a positive bias in the LST product which would also create a spurious correlation with NDSLI? Please discuss this in the manuscript and give references for the quality of MOD11A1 on lake ice and, ideally, sediment covered lake ice.
- LST was obtained from two datasets, MOD11A1 V6 with temporal resolution of 1 day and spatial resolution of 1km. OLI has spatial resolution of 30 meters, there is no clear information on the temporal resolution which should be 16 days. The first paragraph mentioned that the MODIS LST was resampled to higher resolution (1km x 1km), but the data already comes with 1km spatial resolution. Researchers compare this dataset directly to a higher resolution OLI (30m resolution) to provide detailed identification of surface features. A clarification on down-scaling and smoothing out the finer features could be provided.
- The temporal matching between MODIS and OLI is not persuasive, even if no other is accessible. Five out of seven are one day apart, which should be acceptable with careful physical reasoning. However, the remaining two are several days apart e.g. 2017-12-26 (OLI) and 2017-12-20 (MODIS) with a six-day difference. Given the nature of dust and sand particles, they can be easily moved by the wind masses within hours. This concerns already sedimented dust, which can be uplifted again by stronger wind gusts and transported far from the inspected location. The short-life of the particles does not support the long-term effect on the surface temperature. How can one be sure that the LST/NDSLI on both days come from the same spatial pattern of the sediment? I would advise a short description concerning possible matches and the temporal scales of physical drivers of sedimentation.
- Crucial information on data processing is missing as well. What are the errors that come already with the datasets? Especially, given the big difference in spatial resolution, if OLI is down-scaled, isn’t the result smoothed out? MODIS v6 provides LST from the nighttime and daytime, which one was used?
In the first part of the results, is there an uncertainty of the NDSLI retrieval from the images? The retrieval of NDSLI is explained in 3.2.1 and the formula is presented. The mask of NDSLI 0 would be then reasonable, but why does one put another threshold at NDSLI -0.15? Is the final threshold at 0 or -0.15?
Is there any seasonality to the wind directions? Does it explain the patterns? It is said that the analysis does not provide continuous time-series analysis. The spatial patterns can change on day-to-day basis, which depends on the daily wind direction. The only information is in the section 2 on the dominant wind direction which is southwest-west.
On Figure 3 we observe spatial distribution of resampled LST based on MODIS data. The text mentions that the uniform temperature pattern is visible on days without sediment accumulation, possibly referring to figure 2, which is based on OLI data (of different temporal resolution). The dates are different, even if it’s just one day of difference, the connection may be weak. It would be generally nice to have those two figures next to each other, so one can directly compare the corresponding days having the satellite images as in fig. 2 next to the spatial distribution of LST from fig. 3
The section numbering is mixed up, there are two sections 4.4 (one leading to and one following section 4.5)
In section 4.5, the authors mention the effect of sediment accumulation on LST. The visualization is really clear, the result is reasonable and explainable.
Part 4.4.5 (‘Spatial distribution characteristics of Local R² in the MGWR’), including Figure 9 should be explained in more detail. For each location and date there is only one NDSLI and one LST value, correct? Where do the R² values actually come from?
Part 4.4 (“Regression analysis of NDSLI and LST”) is good, but there is explanation missing on the actual regression values. OLS is poorly explained (you give the model but e.g. do not mention what is minimized), the F-statistic is not mentioned in the text. Only fitting is described, based on the and the adjusted . What does the * mean next to p-values? Same issue is with table about GWR and MGWR, the results are nicely presented in a clear table but the description could be extended.
4 Minor issues
L1, L7 - The researchers refer to Lake Surface Temperature as LST. However, the LST retrieved from MOD11A1 is usually Land Surface Temperature. This should be clarified in detail, whether it is just likewise naming. Is the MOD11A1 LST fully able to see the lake surface if it is specified for land observation?
L50 - "On Mars" and not "In Mars".
L52 - the sentence "Aeolian sediments in its high-altitude..." is grammatically incorrect. It lacks a main verb (predicate) to complete the thought.
L57 - what scientific questions? None were directly asked in the text before. It could be "knowledge gaps" that are mentioned.
L77 - point (.) missing, after "regional aeolian conditions"
L93 - Sentence starting with "Secondly, ..." and the following are stylistically inconsistent. Please choose either passive voice or active/imperative voice when describing step-by-step processes.
L103 - Gao et al. 2025 should be Gao et al. 2026 as in the references.
L124 - citation is debatable. The equation (2) from L125 in this paper and equation (9) from Guo et al. 2020, page 4 are the same equations, explained very similarly. I understand that the authors wanted to provide a relevant and recent source, but perhaps more detailed information is missing and maybe a citation of the original paper: doi.org/10.2307/2332142 would be a good addition. Information is missing on the interpretation of the Moran’s I and Local Moran’s indices and the range of the possible values, i.e. -1 to 1. ... generally, the methodology could be explained in more detail
L167-8 - incorrect sentence, starting with "The observation data shows" which is traditionally plural, "The observational data show" sounds more natural. "under the driven of the wind" should be changed to "driven by the wind" or "under the drive of the wind".
L173 - While the data is arranged chronologically, the seven dates are spaced far apart from one another, therefore it is not the best example of a "time-series" which might confuse some readers.
L175 - "calculated the remote sensing images" Images cannot be calculated. They can be processed, or the spectral indices can be calculated, such as NDSLI from the images.
L189 - remove "was"
L208 - November 20, 2017 mentioned in the text does not seem to be analyzed in figures 2, 3. Since this date is not showing any ice cover (right?) it might be best to include it as a test, but not in the following analysis.L210 - April 1, 2018 is not the date on figure 4, but on 2.
L256 - A "demonstrate a clear association" is debatable, a change of words is advised. The relationship between NDSLI and LST is somehow visible, but the interpretation (p. 13) seems very loose.
L262 - study used... and L265 - we constructed... - is should be either passive voice or active voice, mixing both is not a good practiceOn Figure 2 the dataset is OLI, but fig. 2a does not match the date from table 1 (2013-12-31 and 2013-12-30).
Figure 4 is very small.
The label of Figures 5 and 6 are rather unclear.
In Figure 7 some dates seem to come from one dataset, some from the other? If both datasets are used, the dates could be removed or there should be a common label for comparison
In Figure 7 the legend in the lower-right corner is added, even though each subfigure has its own.
Figure 8 - the colour choice is questionable. Especially the little difference between GWR and MGWRis not well visible.
Citation: https://doi.org/10.5194/egusphere-2026-2129-RC2
Viewed
| HTML | XML | Total | BibTeX | EndNote | |
|---|---|---|---|---|---|
| 297 | 162 | 43 | 502 | 29 | 27 |
- HTML: 297
- PDF: 162
- XML: 43
- Total: 502
- BibTeX: 29
- EndNote: 27
Viewed (geographical distribution)
| Country | # | Views | % |
|---|
| Total: | 0 |
| HTML: | 0 |
| PDF: | 0 |
| XML: | 0 |
- 1
This manuscript addresses a timely and important topic at the intersection of aeolian processes and ice-covered lake dynamics. The study’s focus on the Tibetan Plateau is particularly relevant given the region’s sensitivity to climate change and the significant sediment sources. I would like to offer a few comments on the current version.
1. The introduction provides a broad review of light-absorbing particles and snow/ice processes, but it remains largely descriptive. The manuscript would benefit from a sharper articulation of the specific knowledge gap this study addresses. What is currently known versus unknown about the thermal effects of aeolian sediments on lake ice, and how the present work fills that gap. Strengthening this framing would better justify the novelty and significance of the research.
2. Several methodological choices require clearer justification. It would be helpful to explain why MOD11A1 was selected over other available surface temperature products, and why Landsat 8 OLI was chosen for sediment mapping given that higher-resolution alternatives (e.g., Sentinel-2) exist. The thresholds for NDSLI are introduced without quantitative validation, providing accuracy metrics based on at least manual interpretation would strengthen the sediment extraction results.
3. The advantages over conventional GWR are not clearly demonstrated. Is the multiscale GWR framework truly necessary? Moreover, the broadband albedo is used as a key argument for the proposed radiative mechanism. If this analysis is essential, the corresponding data, methods, and results should be moved to the main body.