Skip to content

Fix an issue in snow initialization#1438

Open
smileMchen wants to merge 1 commit intoMPAS-Dev:release-v8.4.0from
smileMchen:initsnow
Open

Fix an issue in snow initialization#1438
smileMchen wants to merge 1 commit intoMPAS-Dev:release-v8.4.0from
smileMchen:initsnow

Conversation

@smileMchen
Copy link
Copy Markdown

@smileMchen smileMchen commented Apr 1, 2026

This PR fixes the bug in snow initialization over Arctic Ocean and other areas covered by seaice.

Snow water equivalent (as well as snow depth and snow cover) in the initial conditions of MPAS is currently set to zero over sea ice areas. This leads to a significant warm bias of skintemp during model integration. When MPAS is initialized in winter, the simulated conditions over polar regions become physically unrealistic.

This PR corrects snow initialization, making it physically reasonable over Arctic Ocean and other areas covered by seaice in winter.

@smileMchen
Copy link
Copy Markdown
Author

smileMchen commented Apr 1, 2026

image

@smileMchen
Copy link
Copy Markdown
Author

smileMchen commented Apr 1, 2026

The upper two panels show initial snow before and after the bug fix ( initial time is 2024-02-01_00)

@mgduda mgduda changed the base branch from develop to release-v8.4.0 April 2, 2026 16:05
@mgduda mgduda self-requested a review April 2, 2026 19:59
This commit fixes an issue in snow initialization over the Arctic Ocean and
other areas covered by seaice.

Snow water equivalent (as well as snow depth and snow cover) in the initial
conditions of MPAS was set to zero over sea ice areas. This leads to a
significant warm bias of skintemp during model integration. When MPAS is
initialized in winter, the simulated conditions over polar regions become
physically unrealistic.

This commit corrects snow initialization, making it physically reasonable over
Arctic Ocean and other areas covered by seaice in winter.
@mgduda mgduda self-requested a review April 2, 2026 23:04
@mgduda
Copy link
Copy Markdown
Contributor

mgduda commented Apr 2, 2026

@smileMchen When testing the changes in this PR with a global 60-km mesh using ERA5 data valid on 2025-01-01, it looks like snow over land points is getting removed. In the screenshot below, the top panel is the snowc field from the current release-v8.4.0 branch, and the bottom panel is the snowc field after locally merging the changes from this PR.
snowc_pr1438
Is this expected?

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Projects

None yet

Development

Successfully merging this pull request may close these issues.

2 participants