This study offers a valuable approach to mapping hazards and identifying sampling points in the lunar south polar region. One aspect I found intriguing is how the thermal models used (Paige and Oxford) rely on certain assumptions, such as uniform regolith properties and static thermal conductivity values. Given the dynamic nature of the lunar environment, with local variations in regolith layering and transient shadowing effects, I wonder how these factors might influence the accuracy of ice stability predictions. Were any steps taken to incorporate or account for these environmental complexities in the models, or could this be a focus for future refinements?
