Assimilation of GPM-retrieved Ocean Surface Meteorology Data for Two Snowstorm Events during ICE-POP 2018
- 1University of Alabama in Huntsville, Huntsville, Alabama, USA
- 2NASA Marshall Space Flight Center, Huntsville, Alabama, USA
- 3Universities Space Research Association, Huntsville, Alabama, USA
- 4ENSCO, Inc./NASA SPoRT Center, Huntsville, Alabama, USA
- 5Kyungpook National University, Daegu, Republic of Korea
- 1University of Alabama in Huntsville, Huntsville, Alabama, USA
- 2NASA Marshall Space Flight Center, Huntsville, Alabama, USA
- 3Universities Space Research Association, Huntsville, Alabama, USA
- 4ENSCO, Inc./NASA SPoRT Center, Huntsville, Alabama, USA
- 5Kyungpook National University, Daegu, Republic of Korea
Abstract. As a component of the National Aeronautics and Space Administration (NASA) Weather Focus Area and GPM Ground Validation participation in the International Collaborative Experiments for PyeongChang 2018 Olympic and Paralympic Winter Games (ICE-POP 2018) field research and forecast demonstration programs, hourly ocean surface meteorology properties were retrieved from the Global Precipitation Measurement (GPM) microwave observations for January – March 2018. In this study, the retrieved ocean surface meteorological products – 2-m temperature, 2-m specific humidity, and 10-m wind speed were assimilated into a regional numerical weather prediction (NWP) framework to explore the application of these observations for two heavy snowfall events during the ICE-POP 2018: 27–28 February, and 7–8 March 2018. The Weather Research and Forecasting (WRF) model and the community Gridpoint Statistical Interpolation (GSI) were used to conduct high resolution simulations and data assimilation experiments. The results indicate that the data assimilation has a large influence on surface thermodynamic and wind fields in the model initial condition for both events. With cycled data assimilation, positive influence of the retrieved surface observation was found for the March case with improved quantitative precipitation forecast and reduced error in temperature forecast. A slightly smaller yet positive impact was also found in the forecast of the February case.
Xuanli Li et al.
Status: final response (author comments only)
-
RC1: 'Comment on gmd-2021-161', Anonymous Referee #1, 24 Oct 2021
Review ofAssimilation of GPM-retrieved Ocean Surface Meteorology Data for Two Snowstorm Events during ICE-POP 2018by Li et al.General commentsThis manuscript intends to investigate the impact of the assimilation of the Global Precipitation Measurement (GPM) retrieved ocean surface observations on two winter storms during 2018 using GSI and WRF model. The discussions primarily focus on the analysis impacts in the observation space, and the precipitation predictions for each case. This study has potentials in improving our understanding of how to optimally utilize the GPM retrievals. But the unclear motivation and lack of in-depth discussions make the paper reading frustrating. I would not recommend publication of this work until substantial revisions have been made. The detailed comments are listed below.Major comments/concerns:
- One major concern I have about this manuscript is the unclear goal. If this study is focusing on the assimilation of GPM retrievals as suggested by the title, then the introduction should provide more information about the GPM retrievals, such as how it was (or was not) previously applied/researched in the field of snowstorm predictions (or other similar fields), and what’s the novelty of this current study (not just the technical specifics in section 2.1). Also, the results section should contain more discussions on how the corrections from this dataset are physically benefiting the storm predictions, not just showing there is a difference and the model is improved. This is a case study after all. Please consider re-organize the manuscript/title to better represent your point and outstand your novelty.
- The results section spends large portions of contents on how the background is modified towards the observations for each snowstorm event. But the discussions on either A-B or O-B can just tell us that the DA is working properly since the verification is dependent. Any successful DA should make the analysis closer to observations than the background, this is from the mathematical nature of the minimization of the cost function. As indicated between L126-132, there should be plenty of other observations (e.g. D3R) available in addition to the retrievals. Why not use those independent observations to verify your analysis? It will objectively measure if the DA is really improving the analysis.
Minor but still important comments:- L220: The cost function is not just to measure the difference between the model and observations (otherwise it should be just yobs-Hxb). It is the sum of weighted (B and R) differences between the analysis estimation xa and model background xb (Jb), and between xa and observations y (Jo). The goal of DA is to find an optimal estimation xa that minimizes this J.
- L225: xguess is xb, please use just one subscription for consistency. This y term is also commonly known as innovation in recent DA literatures.
- L240: This comment may be trivial, but can you provide more details on how you perform the cycling? There are two start mode in WRF, restart and cold start. In the restart mode, WRF is able to continue the model integration without interruption by using both the analysis at the current time step and the tendencies from the previous time step. However, our current DA methods usually only updates the analysis, not the tendencies. I'm wondering if you saw any discontinuity issue if you are using this mode. On the other hand, if you are using cold start mode (restart=.false. by default, which only uses the analysis at the current step), it will interrupt previous integrations and perform differently with your control experiments even without DA. You have to let the control exp stop at the same time for consistency.
- L261: Current Fig. 5 may not be the best way to show those info. I don't think it's possible to retrieve the percentile info directly from the figure.
- L261-271: What do these statistics mean? Some of the o-b standard deviations are almost comparable or smaller than the corresponding observation errors (e.g. specific humidity). Are you suggesting your model background is very good already?
- L291-293: Maybe just show the RMSD of O-A and O-B to better present the point if you are not going to physically explain how those improvements can benefit the storm analysis or predictions?
- L299-301: I’m confused by the motivation of the discussion here. Differences should be there as long as you changed the initial analysis. Of cause the differences will spread out as the model tries to balance the changes. But without observations in those areas, how do we know if this “spread of information” is correct or not?
- L308-315: What is your localization length scale? Are these reducing differences in the area of your observations and increasing differences in the remote area reasonable? What’s the physical explanations or guesses?
- L324-329: How should we interpret this result from Fig. 8? Is it reasonable for a continuously cycled DA experiment to become more alike the NoDA experiment after several cycles? It appears to me that the impact of DA is fading after cycling. Is this true? Also, the increasing RMSD of wind speed at almost all levels from 21 UTC to 22 UTC (Fig. 8d) seems to be inconsistent with the reducing total RMSD in Fig. 8a. Why?
- Fig. 9 -14: Can you use something like ETS score to quantitively verify the precipitation predictions? Also, in your discussions with the precipitation patterns, can you physically relate them to your previous analysis impacts? E.g. how does a warmer temperature in the analysis within certain area result in more precipitation predictions.
- L371-380 and others: Why does the accuracy of this IMERG matter? How is this related to the focus of this paper? It reads to me that the discussions on the IMERG is out of no where. If you think this is a important part of your project and is somehow related to the focus of this paper, please add more information in the intro.
- L407-410: Why does the specific humidity increase while the innovations are mostly negative?
- L444-445: What makes the Feb case less significant than the Mar case? Any hypothesis? It appears to me that the O-B differences in Feb case is larger than the Mar case (Fig. 6 v.s. Fig. 12). Could it be related to my previous comment on the specific humidity?
-
AC2: 'Reply on RC1', Xuanli Li, 24 Nov 2021
The comment was uploaded in the form of a supplement: https://gmd.copernicus.org/preprints/gmd-2021-161/gmd-2021-161-AC2-supplement.pdf
-
CEC1: 'Comment on gmd-2021-161', Juan Antonio Añel, 25 Oct 2021
Dear authors,
After checking your manuscript, it has come to our attention that it does not comply with our Code and Data Policy.
https://www.geoscientific-model-development.net/policies/code_and_data_policy.htmlThe GPM-retrieved ocean surface meteorology data is stored in "https://www.nsstc.uah.edu/public/xuanli.li/GMD/". However, this is not a suitable repostory. Therefore, please, publish your data in one of the appropriate repositories (e.g., Zenodo).
In this way, before the Discussions stage is closed, you must reply to this comment with the link to the repository for the code and the corresponding DOI.Also, remember including in a potential reviewed version of your manuscript the modified 'Code and Data Availability' section with the DOI of the new repository.
Regards,
Juan A. Añel
Geosc. Mod. Dev. Executive Editor
- AC1: 'Reply on CEC1', Xuanli Li, 25 Oct 2021
-
RC2: 'Comment on gmd-2021-161', Anonymous Referee #2, 29 Oct 2021
The comment was uploaded in the form of a supplement: https://gmd.copernicus.org/preprints/gmd-2021-161/gmd-2021-161-RC2-supplement.pdf
-
AC3: 'Reply on RC2', Xuanli Li, 24 Nov 2021
The comment was uploaded in the form of a supplement: https://gmd.copernicus.org/preprints/gmd-2021-161/gmd-2021-161-AC3-supplement.pdf
-
AC3: 'Reply on RC2', Xuanli Li, 24 Nov 2021
Xuanli Li et al.
Xuanli Li et al.
Viewed
HTML | XML | Total | BibTeX | EndNote | |
---|---|---|---|---|---|
689 | 138 | 27 | 854 | 6 | 5 |
- HTML: 689
- PDF: 138
- XML: 27
- Total: 854
- BibTeX: 6
- EndNote: 5
Viewed (geographical distribution)
Country | # | Views | % |
---|
Total: | 0 |
HTML: | 0 |
PDF: | 0 |
XML: | 0 |
- 1