SST Bias Underneath Sea Ice
- Thread starter CalvinC
- Start date
-
- Tags
- sst_update
Hi Ming,
Thank you for your continued aid on this issue. I don't seem to be able to upload either the grib2 or the WPS intermediate files to this post as they are too large for the forum to handle. Is there another way I can send you these files? I can email a zip file to you if you would like? Let me know what you would like me to do.
I have tried both the GFS and GFS GDAS files (d084001 and d083003 respectively) with config_nfglevels set to both 27 and 34, and have received the same error each time. I have attached the most recent log files (both out and err) generated by my attempts, this is for GFS_GDAS files with config_nfglevels set to 27, and config_nvertlevels set to 55.
Let me know what you would like for me to do about getting you the grib2 and intermediate files that I am using. The zipped intermediate files are ~800 MB in size.
Thank you,
Calvin
Thank you for your continued aid on this issue. I don't seem to be able to upload either the grib2 or the WPS intermediate files to this post as they are too large for the forum to handle. Is there another way I can send you these files? I can email a zip file to you if you would like? Let me know what you would like me to do.
I have tried both the GFS and GFS GDAS files (d084001 and d083003 respectively) with config_nfglevels set to both 27 and 34, and have received the same error each time. I have attached the most recent log files (both out and err) generated by my attempts, this is for GFS_GDAS files with config_nfglevels set to 27, and config_nvertlevels set to 55.
Let me know what you would like for me to do about getting you the grib2 and intermediate files that I am using. The zipped intermediate files are ~800 MB in size.
Thank you,
Calvin
Hi Calvin,
I finally managed to finish two tests of global MPAS run using the same mesh as yours. I use GFS as input since it has snow information over the Arctic. Please see the plots I attach here. With the modified code I sent to you, skintemp is differenrt than that from the original code. I haven't compared with GFS data yet.
I finally managed to finish two tests of global MPAS run using the same mesh as yours. I use GFS as input since it has snow information over the Arctic. Please see the plots I attach here. With the modified code I sent to you, skintemp is differenrt than that from the original code. I haven't compared with GFS data yet.
Last edited:
I am not sure why your GFS data doesn't work. I run with gfs.0p25.2024020100.f000.grib2, which can be downloaded from NCAR RDA (d084001).
For this data, config_nfglevels = 34
Can you try to repeat what I did, just to make sure all works as expected on your side. Please let me know if you have any issues during this process.
For this data, config_nfglevels = 34
Can you try to repeat what I did, just to make sure all works as expected on your side. Please let me know if you have any issues during this process.
Hi Ming Chen!
I apologize for the delay in my response, however, I was able to get the GFS data to work, and I am able to get results similar to yours. The trick was to set the vertical levels of MPAS to only extend to 28000 meters, rather than 30000.
However, our issue seems to be further from being fixed with the new snow height. I performed a series of simulations changing the snow height. The results can be seen in the attached plot. The black line represents the ERA5 data, the orange line represents an MPAS run with ERA5 default initial conditions, the blue line is an MPAS run with the snow height over sea ice fixed at 1 meter, the green line is an MPAS run with the snow height over sea ice fixed at 5 meters, and the two olive lines represent MPAS runs using GFS initial conditions. The solid olive line is using the default MPAS code, and the dashed olive line is using your fix for the initial conditions processing.
As you can see, the larger the snow height, the colder the skintemp, and MPAS is already biased cold as compared to ERA5. So while it is probably an issue that MPAS is not capturing the snow height correctly from ERA5 initial conditions, the issue of the cold bias remains and becomes worse when we fix the snow issue.
Additionally, as you can see from the plot attached here and from the SST anomaly plot I attached on the first post, the sea surface temperatures in MPAS under sea ice are well below freezing, which is unphysical. I have found some attempts in the model code to prevent the SST under sea ice from going below 271 K, but evidently these are not working correctly.
All of this is to say, there might be a larger issue here that needs addressing, which is how MPAS initializes the sea surface temperature underneath sea ice points. If that were to be repaired, this might fix the skin temperature issue, as I believe that the skintemp above sea ice is sensitive to the temperature of the ice below it, which is being set to unphysically cold levels via the sea surface temperature initialization issue.
Let me know what you think of this. If you have any thoughts on how to fix the sea surface temperature under sea ice issue, I would love to hear about it. Thank you for your continued help on this issue!
Calvin
EDIT: There seems to be an issue with file uploads on the forum currently. I will send you an email with the plot and post it here when this issue is resolved.
I apologize for the delay in my response, however, I was able to get the GFS data to work, and I am able to get results similar to yours. The trick was to set the vertical levels of MPAS to only extend to 28000 meters, rather than 30000.
However, our issue seems to be further from being fixed with the new snow height. I performed a series of simulations changing the snow height. The results can be seen in the attached plot. The black line represents the ERA5 data, the orange line represents an MPAS run with ERA5 default initial conditions, the blue line is an MPAS run with the snow height over sea ice fixed at 1 meter, the green line is an MPAS run with the snow height over sea ice fixed at 5 meters, and the two olive lines represent MPAS runs using GFS initial conditions. The solid olive line is using the default MPAS code, and the dashed olive line is using your fix for the initial conditions processing.
As you can see, the larger the snow height, the colder the skintemp, and MPAS is already biased cold as compared to ERA5. So while it is probably an issue that MPAS is not capturing the snow height correctly from ERA5 initial conditions, the issue of the cold bias remains and becomes worse when we fix the snow issue.
Additionally, as you can see from the plot attached here and from the SST anomaly plot I attached on the first post, the sea surface temperatures in MPAS under sea ice are well below freezing, which is unphysical. I have found some attempts in the model code to prevent the SST under sea ice from going below 271 K, but evidently these are not working correctly.
All of this is to say, there might be a larger issue here that needs addressing, which is how MPAS initializes the sea surface temperature underneath sea ice points. If that were to be repaired, this might fix the skin temperature issue, as I believe that the skintemp above sea ice is sensitive to the temperature of the ice below it, which is being set to unphysically cold levels via the sea surface temperature initialization issue.
Let me know what you think of this. If you have any thoughts on how to fix the sea surface temperature under sea ice issue, I would love to hear about it. Thank you for your continued help on this issue!
Calvin
EDIT: There seems to be an issue with file uploads on the forum currently. I will send you an email with the plot and post it here when this issue is resolved.
Hi Calvin,
I did a wintertime test case and got the same results as you did, that is, with physically reasonable initial snow condition, skintemp becomes colder than that without snow in the initial condition. I haven't figured out what is the reason behind. I know that Noah LSM has a very simple snow cap module, which may be responsible for the results we have seen. However, more detailed analysis needs to be done before we know for sure what is going on.
I am extremely busy these days because I need to prepare for the release of WRFv4.8 and have little time to explore this issue. If convenient, can you (1) compare MPAS results (driven by GFS input, with and without snow in its initial condition) with GFS skintemp? Let's see which one yields smaller biases, (2) run the same case buy using NoahMP instead of Noah. Let's see how NoahMP performs.
Ming
I did a wintertime test case and got the same results as you did, that is, with physically reasonable initial snow condition, skintemp becomes colder than that without snow in the initial condition. I haven't figured out what is the reason behind. I know that Noah LSM has a very simple snow cap module, which may be responsible for the results we have seen. However, more detailed analysis needs to be done before we know for sure what is going on.
I am extremely busy these days because I need to prepare for the release of WRFv4.8 and have little time to explore this issue. If convenient, can you (1) compare MPAS results (driven by GFS input, with and without snow in its initial condition) with GFS skintemp? Let's see which one yields smaller biases, (2) run the same case buy using NoahMP instead of Noah. Let's see how NoahMP performs.
Ming