Using the ERA5 land/sea mask data in WRF

[ccalkins]

I have noticed that there is an issue with the land sea mask that you get from ECMWF is not binary.



This is also noted on ECMWF webpage:
Land-sea mask in ERA5 and ERA-Interim - Copernicus User Support Forum - ECMWF Confluence Wiki

Because the mask is not binary I believe it creates an issue with metgrid log
INFORM: Field LANDSEA.mask does not have a valid mask and will not be checked for missing values

which I believe causes issues with TSLB in the wrfinput files


Is there a standard with dealing with these issues? Such as modifying the SFC grib file and changing the values to be either 0/1 based or modifying the TSLB or some other variable in the met_em files to be the fill_value of 285? I thought it could possiblyy be a issue with the model level data but I've done it with Pressure levels from ERA5 and TSLB also shows up as 0K along areas of water. I should note that TSLB is 0K in the wrfout 00 hour file but in other files, the 0K issue has been fixed to be a single constant value over areas of water.

I have also provided a MFT link to some of my files from my run, that I left uncompleted. Let me know if other files would be needed that you also wish to examine.
https://mft.aer.com/download?domain...b8ec5d5b903b-104f5fa9fc1045ef8b22d2518bb0efa4

Chase

 

[Ming Chen]

Hi Chase,
NCAR RDA provides ERA5 dataset (NCAR RDA Dataset ds633.0). From this website you can find the link to "ERA5 atmospheric invariant fields", where you can download the datafile

e5.oper.invariant.128_172_lsm.ll025sc.1979010100_1979010100.grb

WPS can ungrib this file and yield correct landsea mask.

Would you please try and let me know if you have any issue? Thanks.

 

[MLandreau]

Hi,
I think I faced the same issue. If yes, the answer is in this post :

Issue with ERA5 SST fields

Hello all, I'm making a WRF run forced wit ERA5 reanalysis, and I noted that ERA5 SST fields look weird... First, the mask values over land is 0ºC, which when WPS will interpolate the fields in metgrid will give irrealistic surface temperature fields over coastal areas. obviously, this will in...

forum.mmm.ucar.edu

You need to change the metgrid.TBL file.
I join a modified version that works for me.
Mathieu

 

[Ming Chen]

Mathieu,
Would you please provide more information regarding the changes in METGRID.TBL? Specifically,
(1) It is based on which version of WPS?
(2) Where did you download the ERA5 reanalysis data? What is the resolution of the ERA5 data?
(3) Does this modified METGRID.TBL work for ERA5 data of various resolution?
Thank you.

 

[MLandreau]

Hi,
1- I faced this issue on version 4.4, I am still using this metgrid.TBL on version 4.5
2- I downloaded ERA5 data on Copernicus with a resolution of 31 km, with 37 pressure levels
3- I tried only with my ERA5 dataset

More informations :
With the default metgrid.TBL file, the landmask used in metgrid.exe seems to be the ERA5 landmask. SST was defined only at points considered as water in ERA5 (see figures a_LANDMASK and a_SST, the coastline resolution is not the same). This cause wrong values in TSLB but only in water points (which is not that much a problem). However, it also creates wrong SST values near the coastline (extrapolated from land points TSK). These wrong SST are kept all along the simulation if sst_update isn't used.
With the modified Metgrid.TBL, these problems are fixed.

Mathieu

 

[Ming Chen]

Thank you for the kind information !

 

[wallis]

This works for ERA5 (0.25deg), but what about IFS? The IFS landmask has the same problem with fractional values, and I don't think NOAA/RDA has a 0.1degree landmask available.

 

[Ming Chen]

I don't have a solution to IFS yet ....

 

[jhegarty]

Hi Ming,

I am working with Chase who posted the first message in this thread. I tried your solution using the file e5.oper.invariant.128_172_lsm.ll025sc.1979010100_1979010100.grb to create a constants file with the landsea mask. I also tried Mathieu's solution modifying the interpolation of SST in METGRID.TBL. Interestingly, when I looked at the TSK and T2 fields after 1-hour of simulation time on a 12 km grid they were almost identical for the two cases even though the SST fields at the coast were still very wrong in the landsea mask case. I have attached some files to illustrate this point. I am wondering how the TSK and T2 can be correct if SST is wrong and which solution is preferable. Thanks for your help. Jen

 

[Ming Chen]

Hi Jen,

Let's try the option: remove the following line from Vtable.ECMWF

34 | 1 | 0 | | SST | K | Sea-Surface Temperature |

Then rerun ungrib.exe. In this case you will have no SST. SKINTEMP will be used as SST later. In this way I expect that there is no unreasonable SST along the coastline.

Let me know if you have more questions.

 

[Josh]

I would like to add that while this issue seems insignificant, it could cause an obscure WRF runtime error that will crash the entire WRF model. That is why this issue should be corrected by using one of the above described methods.

I will provide an example of the WRF runtime error and my troubleshooting.

1) I was running a 1 year WRF simulation with ERA5 data. 1 domain with 16 processors. I am using the fractional LANDSEA field described above. Can WPS handle this field correctly (or handle it well enough where it doesn't really matter)?
Or should I use the "e5.oper.invariant.128_172_lsm.ll025sc.1979010100_1979010100.grb" file @Ming Chen mentioned?

2) The WRF model suddenly crashes on Jaunary 26th (about 1 month into the simulation) for no reason.

3) There are no CFL errors in any of the rsl files.

4) The main rsl files give no hint of what caused the error.

5) Some (but not all) of the rsl files specific to each processor have strange repeated output:

Flerchinger USEd in NEW version. Iterations= 10

6) This is the first clue.

7) I tried starting the WRF simulation just a few hours before the hour in January 26th when it fails. It still fails at exactly the same hour as when starting the run ~40 days earlier. From this, I know the problem is related to the input data at this specific time only (possibly more times later in the year).

8) I have sst_update = 1. If I change my namelist.input to use sst_update = 0, the run no longer fails at this time. This means sst is somehow related to the error.

9) I compare the SST and LANDSEA variables in the met_em file that gives the WRF error downstream, and the met_em file the hour before this file that does not give the WRF error downstream.
Surprisingly SST = SST and LANDSEA = LANDSEA between both files.

10) I change all the SST values that are 0 to 310.1 (a valid value) in the met_em file that gives the error. I rerun real.exe and wrf.exe, and WRF does not crash anymore. This indicates that SST is not really the issue but it indicates there is an issue with a different variable (maybe TSLB)?

11) The issue is not significant enough or doesn't exist frequently enough to cause all hours in the WRF simulation to fail. For my situation, it ran for over a month without issue, before a sudden issue. This likely indicates the issue could be related to missing data for a variable, and how that interacts with the LANDSEA mask and the WRF calculations related to sea surface temperature.

12) I have not tried any of the above methods to fix/get around this issue, but I plan to first try to remove SST from the Vtable so WRF will use SKINTEMP instead.

Josh