Scheduled Downtime
On Friday 21 April 2023 @ 5pm MT, this website will be down for maintenance and expected to return online the morning of 24 April 2023 at the latest

SST Bias Underneath Sea Ice

CalvinC

New member
Hello MPAS-A Forum! This is my first time posting here so apologies if I forget to include some crucial information.

I am attempting to use MPAS-A (uniform 60km resolution mesh) to replay some Arctic cold air intrusions and warm air outbreaks, and I am investigating how long MPAS-A stays faithful to the ERA5 reanalysis conditions. I noticed that both the 2 Meter Temperature and the Skin Temperature are initialized at reasonable values, and then become very cold very quickly. I investigated the sea surface temperature, and found that, despite using initial conditions and update files with MPAS, the SST beneath the sea ice is not being initialized correctly or being updated.

I have attached two plots. The first is an Arctic area weighted average of skin temperature and SST of an MPAS-A run initialized on January 20th, 2013. It shows that the skin temperature is initialized at a reasonable value and then becomes much colder in MPAS than ERA5, and it shows that the MPAS SST starts well below the ERA5 initial conditions. The second plot is a map of SSTs at the first model timestep (January 20th, 2013, 00h), in ERA5, and MPAS-A (uniform 60km resolution), as well as their anomaly. I have plotted the 15% sea ice concentration contour. This shows that the MPAS-A cold bias only occurs underneath the sea-ice.

I would appreciate any tips on why this is happening! Is there a way to have a skin temperature update file to periodically push the MPAS-A skin temp back to ERA5 values, like how we update SST and sea ice concentrations? On maybe a more concerning note, why are the SSTs beneath the sea-ice edge not being initialized at the ERA5 values? Have I set something incorrectly in my namelists?

Namelist and streams files have been attached as well.
I look forward to hearing from you all!
 

Attachments

  • area_avg.pdf
    17.4 KB · Views: 3
  • era5_MPAS_comp_2013-01-20T00.png
    era5_MPAS_comp_2013-01-20T00.png
    144.7 KB · Views: 4
  • namelist.atmosphere.txt
    2 KB · Views: 2
  • streams.atmosphere.txt
    4.7 KB · Views: 2
Hi,

I apologize for answering so late. There are many posts pouring in every day, and somehow this post was skipped and I just found that it was never answered.

Thanks a lot for the detailed description of the SST issue. In your 2nd plot, the SST anomaly is the difference between MPAS and ERA at the first time step, is this correct?

Would you please plot SST difference between MPAS initial condition and ERA5? It seems that SST in MPAS is much lower at the right beginning. I would like to make sure whether it is derived from the process of initialization.

Also, please let me me which version of MPAS you are using. Thanks.
 
Last edited:
Hi Ming Chen!

No worries about your reply time, I appreciate the help!

It is correct that the plot I sent before is from the first time step. I have made a new plot (attached here) of the ERA5 SST from the time of the start of the run, and the SST from the MPAS init file (x1.163842.init.nc). I used the ERA5 file to generate my initial conditions, so it seems that something is going wrong with the generation of initial conditions. Do you have any tips on how I could diagnose this issue?

I am using MPAS-v8.3.1 on the Narval HPC, based out of Montreal.

Thank you so much!
Calvin
 

Attachments

  • initial_conditions.png
    initial_conditions.png
    144.9 KB · Views: 2
Hi Calvin,

Thank you for the plots of the initial SST. I do 't have an immediate explanation why SST in MPAS initialization is so severely deviated from ERA5. I will look into this issue and get back to you.

Since your case is a region MPAS case, would you please send me your regional domain information (i.e. points.txt you used to create the regional mesh).

I suppose that you didn't run grid_rotate to relocate any refinement. Please let me know if I am wrong.

Thanks.
 
Hi Ming,

I should clarify, this is not a regional MPAS case, I am simulating the entire globe but I am only plotting the Arctic as this is my area of interest (and where I noticed the SST discrepancy). I did not relocate any refinement because this is a uniform 60km grid.

Other details: I used the era5_to_int.py tool to convert my era5 data to initial conditions, and I am using the convert_mpas tool to interpolate my MPAS data to lat/lon coordinates for plotting purposes.

If there are any other files (logs, initial conditions, etc) that you need from me, please let me know and I can provide them!

Thank you for looking in to this.
Calvin
 
Calvin,

Thank you for the clarification. I will look into this issue and get back to you once I know for sure what is wrong. It may take some time partly because I exceeded the threshold of HPC usage and cannot submit job at present. Thank you for your patience.

Ming
 
Top