NOAA-SWPC / GSMWAM-IPE

Coupled NEMS app for WAM-IPE
5 stars 11 forks source link

Tidal analysis on the recent WAM-IPE annual output with new eddy mixing parameterization #80

Open ZhuxiaoLi66 opened 1 year ago

ZhuxiaoLi66 commented 1 year ago

The IDL tidal analysis package from Svetlana has been revised and applied to do the tidal analysis on the output of the recently WAM-IPE annual run (year 2022) with our new eddy mixing parameterization (skeddy0=70, skeddy_semiann=-10) which already applied to the current operation. I will post the analysis results here and compare the results of Rashid (2008) and Svetlana's recent analysis on FV3WAM. The comparison with SABER data also is aiming to do in this summer.

ZhuxiaoLi66 commented 1 year ago

The following plots show the T amplitudes of DW1, and DE3 at different height levels, also attached the corresponding results in Rashid's 2008 paper. The DW1 at 100km of recent WAM run looks good except the fall peak is much larger than Rashid's result, also the SABER data. will double check DE3 output. Also attached Svetlana's analysis for recent FV3WAM at the bottom for comparison.

T_tides_DW1_2022_WAM_para5_cb18 T_tides_DE3_2022_WAM_para5 Rashid_2008_DW1_DE3

FV3WAM_T_tides_DW1

ZhuxiaoLi66 commented 1 year ago

T amplitude of SW2 is shown below, also attached Rashid's results for comparison. The two tides maps are quite similar in SW2. Again, the WAM runs seems has strong tides than obs data SABER. T_tides_SW2_2022_WAM_para5 Rashid_2008_SW2_T_DE3_U

ZhuxiaoLi66 commented 1 year ago

Here is the FV3WAM tides for SW2. FV3WAM_T_tides_SW2

ZhuxiaoLi66 commented 1 year ago

will do more analysis and double check the results.

ZhuxiaoLi66 commented 1 year ago

The following is the DW1 amplitude of U & V for the 2022 WAM output. U_tides_DW1_2022_WAM_para5_50 V_tides_DW1_2022_WAM_para5

ZhuxiaoLi66 commented 1 year ago

The following is the SW2 tides of U & V in 2022 for the new para WAM output. U_tides_SW2_2022_WAM_para5 V_tides_SW2_2022_WAM_para5

ZhuxiaoLi66 commented 1 year ago

TDE3_FV3_WAMGSM_INLINE

ZhuxiaoLi66 commented 1 year ago

With Svetlana's great help to do the Temperature DE3 analysis on the same 2022 output dataset as mine, we confirmed that my application of the tidal package is correct since the results look almost the same as follows, the different tidal features in the WAM-IPE 2022 free run output from Rashis's results should be caused by the model performance itself.

  1. Svetlana's results TDE3_WAMOP_Svetlana
  2. Zhuxiao's results T_tides_DE3_2022_WAM_para5_cb27
ZhuxiaoLi66 commented 1 year ago

Planning to do the tidal analysis of the operational WAM-IPE output for 2022 and 2021 respectively.

ZhuxiaoLi66 commented 1 year ago

The operational output of gsm10 for 2022 has been downloaded and processed, and a similar tidal analysis has been performed on the dataset. The peak amplitude of T DE3 in operational output is a little bit weak than in the free run, the annual pattern seems not improved much compared with the free run. T_tides_DE3_2022_WAM_opr_cb27

ZhuxiaoLi66 commented 1 year ago

T_DE3_116km_comparison

timfullerrowell commented 1 year ago

Thanks Zhuxiao,

This all looks good. It's encouraging Svetlana's and your results are consistent.

The DE3 from the operational WAM looks closer to previous Rashid/SABER with the single annual peak in northern summer months. Be nice to do the same "smoothing" as Rashid did to see if the amplitudes are similar.

The multi-day variability is also interesting. It would be good to correlate with a tropospheric weather/tropical convection "index", if there is one? Tim

On Wed, Jul 12, 2023 at 12:37 PM ZhuxiaoLi @.***> wrote:

[image: T_DE3_116km_comparison] https://user-images.githubusercontent.com/22546571/253067374-0e08b393-25c0-430b-b3ad-15e9aebdb5bc.png

— Reply to this email directly, view it on GitHub https://github.com/NOAA-SWPC/GSMWAM-IPE/issues/80#issuecomment-1633027162, or unsubscribe https://github.com/notifications/unsubscribe-auth/AH5BFOHE2EILQYMGKLHIZMLXP3VF5ANCNFSM6AAAAAAYECQGGY . You are receiving this because you were assigned.Message ID: @.***>

ZhuxiaoLi66 commented 1 year ago

Thanks for the reminder, Tim. Going to do the smoothing on the dataset. As for the 'multi-day variability', do you mean the multi-peaks with short periods in the T DE3 plot for the operational run?

timfullerrowell commented 1 year ago

Yes, exactly.

On Wed, Jul 12, 2023 at 2:29 PM ZhuxiaoLi @.***> wrote:

Thanks for the reminder, Tim. Going to do the smoothing on the dataset. As for the 'multi-day variability', do you mean the multi-peaks with short periods in the T DE3 plot for the operational run?

— Reply to this email directly, view it on GitHub https://github.com/NOAA-SWPC/GSMWAM-IPE/issues/80#issuecomment-1633170251, or unsubscribe https://github.com/notifications/unsubscribe-auth/AH5BFODYIWOP5TRW3GQETXLXP4CI3ANCNFSM6AAAAAAYECQGGY . You are receiving this because you were assigned.Message ID: @.***>

ZhuxiaoLi66 commented 1 year ago

tested the 30-day U_30d_DE3_95km_comparison running mean for our U DE3 tides at 95km.

ZhuxiaoLi66 commented 1 year ago

Aiming to investigate the tides performance of the operational WAM-IPE and its interannual variability, I downloaded so far two years (20210715-20230712) gsm10 3D output of the operational from HPSS, processed and 'suture' some small data gaps on the disk, did the tidal analysis with 30days smoothing on Jul.2021-Jul.2022, Jul.2022-Jul.2023 respectively to do some comparison. the results of the two years show consistency between the two years which is encouraging. the interannual variation exists but is not huge. On the other hand, the seasonal peaks of DE3 of T (Sep. & Oct.) & U (Jan. & Feb.) looks quite different. the results are shown as follows.

30d_T_DE3_202107_202207 30d_T_DE3_202207_202307 30d_U_DE3_202107_202207 U_DE3_20220715_por

ZhuxiaoLi66 commented 1 year ago

For SW2 tides at 100km compared with Rashid's former results.
30d_SW2_T_100km_comparison

ZhuxiaoLi66 commented 1 year ago

The comparison of T DE3 results. 30d_T_DE3_202207_202307_comparison

ZhuxiaoLi66 commented 1 year ago

T_DE3_ori_para5_opr_2022

ZhuxiaoLi66 commented 1 year ago

T_DE3_ori_para5_opr_U_DE3_para5

ZhuxiaoLi66 commented 1 year ago

With Astrid's great help, we have got the tidal analysis plot from Chihoko in LASP on the 2021 SABER data. SABER_amp_phase_90km_2021_Chihoko

ZhuxiaoLi66 commented 1 year ago

The following is my tidal analysis results on the 2022 WAM-IPE free run output with the new eddy mixing parameterization (para5) which is applied in current operational.
U_tide_amp_phase_90km_2022_para5_30d T_tide_amp_phase_90km_2022_para5

ZhuxiaoLi66 commented 1 year ago

This panel shows the similar results on our 2021 operational results, only a few months from July of 2021. just for a comparison. U_tide_90km_amp_phase_20210715_1231_opr_10d T_tide_90km_amp_phase_20210715_1231_opr_10d

ZhuxiaoLi66 commented 1 year ago

I think we just have got some encouraging results in the comparison of the SABER tidal analysis from Chihoko and our operational & para5 output in amp & phase of T tides. The following plots are in the almost same color-bar to compare. We can see the shape and magnitude of our results are very close to Chihoko's results. Astrid, I believe that the plot that Chihoko sent should be Temperature tides instead of U_winds. will send an email later to both of you. image T_tides_opr_202107_202207 T_tides_202207_202307 T_tides_2022_para5

ZhuxiaoLi66 commented 1 year ago

We will ask the plot for U_tides later from Chihoko, so far, it seems that there is some problems in SW1 tide phase of operational U_wind, not in free new para5 run though. post here and will do the double check later. U_tides_opr_202107_202207

U_tides_opr_202207_202307

U_tides_2022_para5

astridmaute commented 1 year ago

SABER measures temperatur, theref will be no neutral winds. The temperature amplitude comparison looks great.

ZhuxiaoLi66 commented 1 year ago

Thanks, Astrid. The following is the amplitude comparison between recent operational & para5 runs and Rashid's 2008 paper results at different height levels from above 90km. We can see temperature tides looks good, the U_wind_DE3 is kind of off.

image

image

image

image

ZhuxiaoLi66 commented 1 year ago

ICON_MIGHTI_DW1_yamazaki

ZhuxiaoLi66 commented 1 year ago

the following plots are the comparisons of amplitudes & phases among the SABER data, operational WAM-IPE output, and new-parameters WAM-IPE free runs output. The WAM-IPE results are very consistent with SABER observations. compari_amp_90km_4P compari_phase_90km_4P

twfang commented 1 year ago

@ZhuxiaoLi66 Thanks, but it is very difficult to see these plots. Could you please upload a better version?

ZhuxiaoLi66 commented 1 year ago

@twfang you mean the plot sizes are different?

astridmaute commented 1 year ago

the plots are too small to read the numbers, even when downloading and increasing size due to low resolutions the numbers are not legible.

ZhuxiaoLi66 commented 1 year ago

Got it, thanks, Astrid. Maybe these resulted from the cutting paste and the two times snip. will try some other method later for a clear version. As for the color bars, are plots are the same for the corresponding panels for comparison. By the way, I did this roughly at first just to let you all get a basic idea of the results, when I get data files from Chihoko, I will make all the panels in one plot for serious usage later.

ZhuxiaoLi66 commented 1 year ago

the following plots are the comparisons of amplitudes & phases among the SABER data, operational WAM-IPE output, and new-parameters WAM-IPE free runs output. The WAM-IPE results are very consistent with SABER observations.

compari_amp_90km_4P compari_phase_90km_4P

ZhuxiaoLi66 commented 1 year ago

The following is the horizontal wind tides based on our WAM-IPE operational output compared with the MIGHTI data tidal analysis by Chihoko.

WAM_MiGHTI_UV_202207_202307_DW1 WAM_MIGHTI_DW1_UV_202107_202207

WAM_MIGHTI_UV_SW2_202107_202207

WAM_MIGHTI_DE3_202107_202207

WAM_MIGHTI_UV_DE3_202207_202307

ZhuxiaoLi66 commented 1 year ago

zonal_wind_tides_MIGHTI_WAM-IPE-opr.pdf

ZhuxiaoLi66 commented 1 year ago

The plots below descript the more tidal analysis results of horizontal wind tides based on our WAM-IPE operational output compared with the MIGHTI data by Chihoko.

MIGHTI_WAM_UV_DE3_100km_2021 MIGHTI_WAM_UV_DW1_100km_2021 MIGHTI_WAM_UV_SW2_100km_2021