the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
A hybrid-grid global model for the estimation of atmospheric weighted mean temperature considering time-varying lapse rate in GNSS precipitable water vapor retrieval
Abstract. The atmospheric weighted mean temperature (Tm) is a key parameter in global navigation satellite system (GNSS) water vapor retrieval and can convert the zenith wet delay (ZWD) into precipitable water vapor (PWV). However, there are some shortcomings in the existing Tm models, such as the detailed time-varying lapse rate not being considered. Additionally, the spatiotemporal characteristics of Tm need to be further refined. Therefore, we developed a new global high-precision and high-spatiotemporal-resolution Tm model considering time-varying lapse rate using the latest European Centre for Medium-Range Weather Forecasts ReAnalysis 5 (ERA5) atmospheric reanalysis data. Firstly, a global multidimensional Tm lapse rate model (NGGTm-H model) was developed using the sliding window algorithm. Secondly, the daily variation characteristics of Tm and its relationships with geographical situation were investigated. Finally, a hybrid-grid global Tm model considering time-varying lapse rate (NGGTm model) was developed. To verify the effectiveness of the proposed model, the NGGTm model was compared with the Bevis and GPT3 models using the Tm data recorded at 378 radiosonde stations in 2017 and the surface grid Tm data calculated from the ERA5 reanalysis data. The results show that taking the surface grid Tm data of ERA5 as reference values, the average root mean square error (RMSE) value predicted by the NGGTm model was 2.84 K, which was higher with 0.50 K, 0.18 K and 0.06 K than those of the Bevis, GPT3-5 and GPT3-1 models, respectively. Meanwhile, taking the Tm data from the radiosonde stations as the reference values, the mean bias and RMSE of the NGGTm model were 0.10 K and 3.30 K, respectively, which exhibit the best accuracy and stability among the Bevis, GPT3-5 and GPT3-1 models.
- Preprint
(2068 KB) - Metadata XML
- BibTeX
- EndNote
Status: final response (author comments only)
-
RC1: 'Comment on gmd-2024-21', Anonymous Referee #1, 05 May 2024
The manuscript entitled "A hybrid-grid global model for the estimation of atmospheric weighted mean temperature considering time-varying lapse rate in GNSS precipitable water vapor retrieval" aims to improve the water vapor retrieval through a new global high-precision and high-spatiotemporal-resolution Tm (the atmospheric weighted mean temperature) model by considering time-varying lapse rate using the latest European Centre for Medium-Range Weather Forecasts ReAnalysis 5 (ERA5) atmospheric reanalysis data. The methods seem to work well and some impressive results were obtained. The manuscript is well organized, clearly illustrated as well as easy to read. I have no critical comments but recommend this manuscript for publication after some revisions as following suggestions.
- Please correct the grammar errors, typo, or missing words in the For example, in line 109, "we aim to global Tm model that takes into account..." and in line 298, “surface-level…”.
- In section 3.3, the NGGTm-H model is validated by radiosonde. I suggest describing the vertical resolutions or the altitude of the record at the 378 radiosonde stations.
- In line 256-260, do the authors try to explain the reason why the positive biases are smaller than the absolute value of the negative biases? If yes, please explain the reason in more detail.
- I suggest using the same color bar ranges in (a)-(d) of figure 5 and 6 to emphasize the result and avoid misunderstanding to readers, especially in figure 6.
Citation: https://doi.org/10.5194/gmd-2024-21-RC1 -
AC1: 'Reply on RC1', Jihong Zhang, 14 May 2024
Dear Referee,
Thank you very much for reviewing our manuscript and providing your valuable comments and suggestions. We have answered your questions one by one in order to improve our research results based on your feedback.
- we have correct the grammar errors, typo, or missing words in line 109, “we aim to develop a global Tm model that takes into account time-varying lapse rate and high-precision capabilities” and in line 298, “The surface gridded Tm data with a temporal resolution of 1 hour derived from the ERA5 reanalysis data in 2017 were selected as reference values”. In addition, we have carefully read the entire text and corrected the remaining errors.
- In section 3.3, we have added descriptions of the altitude of the record at the 378 radiosonde stations, “The altitude of radiosonde stations ranges from 0 to 4500 m, mostly within 2000 m”.
- In line 256-260, we have tried to explain the reason why the positive biases are smaller than the absolute value of the negative biases. Due to our limited research abilities, the reasons we analyzed are as follows: The vertical correction values of Tm obtained using the NGGTm-H model were slightly larger in land areas but smaller in marine areas than the reference values. However, a small number of radiosonde stations distributed in marine areas were susceptible to the influence of marine climate, resulting in the vertical correction values of the model was apparently smaller than the reference values. Therefore, the positive biases were smaller than the absolute value of the negative biases. We have added the reasons in the manuscript. In future work, we will further investigate the reasons and develop solutions.
- We have changed the color bar ranges in figures 5 and 6 to be consistent. Indeed, it will mislead readers if the color bar ranges are inconsistent.
Thank you for your review comments and suggestions, we will take them seriously and revise our manuscript according to your guidance. We promise to try our best to improve and further enhance our research results and quality in the new version. Thank you again for your support and evaluation of our research work.
Sincerely,
Dr. Jihong Zhang
jhzhang@stu.kust.edu.cn
Citation: https://doi.org/10.5194/gmd-2024-21-AC1 -
RC2: 'Reply on AC1', Anonymous Referee #1, 21 May 2024
Dear authors,
Thank you for your reply. Your explanation and corrections for the article are convincing. Hence I think the article is worth to be published.
Citation: https://doi.org/10.5194/gmd-2024-21-RC2
-
CEC1: 'Comment on gmd-2024-21 - No Compliance with GMD's policy', Juan Antonio Añel, 11 May 2024
Dear authors,
Unfortunately, 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.htmlYou have not published your code in a permanent repository that complies with our requirements before the submission of you manuscript. Also, your manuscript must contain such information. This is mandatory according to our policy. Therefore, please, publish your code in one of the appropriate repositories, and reply to this comment with the relevant information (link and DOI) as soon as possible, as it should be available before the Discussions stage. Also, please, include the relevant primary input/output data.
Please, note that the code must have a license. If you do not include a license, the code continues to be your property and can not be used by others, despite any statement on being free to use. Therefore, when uploading the model's code to the repository, you could want to choose a free software/open-source (FLOSS) license.
In this way, if you do not fix this problem, we will have to reject your manuscript for publication in our journal. I should note that, given this lack of compliance with our policy, your manuscript should not have been accepted in Discussions. Therefore, the current situation with your manuscript is irregular.
Juan A. Añel
Geosci. Model Dev. Executive Editor
Citation: https://doi.org/10.5194/gmd-2024-21-CEC1 -
AC2: 'Reply on CEC1', Jihong Zhang, 14 May 2024
Dear Editor,
Thank you for your recognition and reply. We have read the “Code and Data Policy” carefully. After our discussion, we have published the code and data related to our work in ZENODO (http://doi.org/10.5281/zenodo.11188772).
In addition, we have updated the “Code and dada availability” in our manuscript. It is shown below:
Code and data availability. The ERA5 reanalysis data used in this paper can be freely accessed at (http://cds.climate.copernicus.eu/cdsapp#!/search?text=ERA5). The radiosonde data can be accessed at (http://weather.uwyo.edu/upperair/sounding.html). All of the data generated during the current study and the code are available on ZENODO (http://doi.org/10.5281/zenodo.11188772).
Title: A hybrid-grid global model for the estimation of atmospheric weighted mean temperature considering time-varying lapse rate in GNSS precipitable water vapor retrieval
Author(s): Shaofeng Xie et al.
MS No.: gmd-2024-21
MS type: Development and technical paper
Thank you for your consideration. If there is anything else that needs to be modified or noted, please provide guidance and suggestions, and we will further improve it. I look forward to hearing from you.
Sincerely,
Dr. Jihong Zhang
jhzhang@stu.kust.edu.cn
Citation: https://doi.org/10.5194/gmd-2024-21-AC2 -
CEC2: 'Reply on AC2', Juan Antonio Añel, 18 May 2024
Dear authors,
Thanks for your reply. Unfortunately, I do not think it is a satisfactory solution for the issues I raised previously. First, what you have published are a few codes for the NGGT model that you propose. These are coded in M language. You should indicate what kind of interpreter you use to run the code (e.g. GNU Octave, Mathematica,...) and its version. Second. you have upload part of your code as binary files. Binary files can not be read unless the interpreter/software used to produce them is available. Therefore, they are not an optimal solution. It would be better if you can provide them as ASCII files. However, if it is not feasible, at least you should indicate what software is necessary to read them. It should be better if it is non-privative software.
Regards,
Juan A. Añel
Geosci. Model Dev. Executive Editor
Citation: https://doi.org/10.5194/gmd-2024-21-CEC2 -
AC3: 'Reply on CEC2', Jihong Zhang, 20 May 2024
Dear Editor,
First, thank you for giving us another opportunity and reply very much. We have carefully read the comments you provided. We apologize that our solution did not meet your satisfaction. The codes and datas were re-published on ZENODO. What we have published are a few core codes of our proposed NGGT model. The other codes lack universality and are not very meaningful to readers, so we have only published some key codes. These codes should be run using Matlab with 2017 or higher versions. Additionally, we have provided the datas as ASCII files. The codes and datas related to our work were published on ZENODO (http://doi.org/10.5281/zenodo.11218449).
In addition, the “Code and dada availability” was updated in our manuscript. It is shown below:
Code and data availability. The ERA5 reanalysis data used in this paper can be freely accessed at (http://cds.climate.copernicus.eu/cdsapp#!/search?text=ERA5). The radiosonde data can be accessed at (http://weather.uwyo.edu/upperair/sounding.html). All of the data generated during the current study and the code are available on ZENODO (http://doi.org/10.5281/zenodo.11218449).
Title: A hybrid-grid global model for the estimation of atmospheric weighted mean temperature considering time-varying lapse rate in GNSS precipitable water vapor retrieval
Author(s): Shaofeng Xie et al.
MS No.: gmd-2024-21
MS type: Development and technical paper
Thank you for your consideration. If there is anything else that needs to be modified or noted, please provide guidance and suggestions, and we will further improve it. I look forward to hearing from you.
Sincerely,
Dr. Jihong Zhang
jhzhang@stu.kust.edu.cn
Citation: https://doi.org/10.5194/gmd-2024-21-AC3 -
CEC3: 'Reply on AC3', Juan Antonio Añel, 20 May 2024
Dear authors,
Again, thanks for your quick reply. Regarding our request for the code, it is not an issue of whether it is "meaningful" (your view) or not but of having stored all the code you have used. Therefore, please publish the additional code independently of your perception of its usefulness.
Also, you still need to address the issue regarding the software version you use. It is not an issue of readers or reviewers running your code, but it is necessary to know precisely the version number of the software you are using to ensure the replicability of the results. Also, it is relevant in case of future detection of bugs in such software that could impact your results. "2017" or a newer version is not a version number. If you have used Matlab, you should identify its version, which should be listed in its release history (https://en.wikipedia.org/wiki/MATLAB#Release_history).
We appreciate your attention to these matters. Could you please address the issues mentioned and provide us with the necessary information? Your cooperation is crucial for the successful replication of your research.
Juan A. Añel
Geosci. Model Dev. Executive Editor
Citation: https://doi.org/10.5194/gmd-2024-21-CEC3 -
AC4: 'Reply on CEC3', Jihong Zhang, 24 May 2024
Thanks for giving us another opportunity and reply very much. We apologize again that our solution did not meet your satisfaction. All the additional codes have been re-published on ZENODO (http://doi.org/10.5281/zenodo.11255397). The version informations of Matlab we are using are as follows. Version: 9.2. Release name: R2017a. Number: 37. Year: 2017. Release data: March 9, 2017.
In addition, the “Code and dada availability” has been updated in our manuscript. It is shown below:
Code and data availability. The ERA5 reanalysis data used in this paper can be freely accessed at (http://cds.climate.copernicus.eu/cdsapp#!/search?text=ERA5). The radiosonde data can be accessed at (http://weather.uwyo.edu/upperair/sounding.html). All of the data generated during the current study and the code are available on ZENODO (http://doi.org/10.5281/zenodo.11255397).
Title: A hybrid-grid global model for the estimation of atmospheric weighted mean temperature considering time-varying lapse rate in GNSS precipitable water vapor retrieval
Author(s): Shaofeng Xie et al.
MS No.: gmd-2024-21
MS type: Development and technical paper
Again, thanks for your consideration. If there is anything else that needs to be modified or noted, please provide guidance and suggestions, and we will further improve it. I look forward to hearing from you.
Sincerely,
Jihong Zhang
jhzhang@stu.kust.edu.cn
Citation: https://doi.org/10.5194/gmd-2024-21-AC4 -
CEC4: 'Reply on AC4', Juan Antonio Añel, 24 May 2024
Dear authors,
Thanks for providing the information. We can consider now this issue solved. However, if you can, please, update the Zenodo repository with true code instead of .mat binary files. For example, I do not use Matlab, and therefore can not "load" and check the code inside such files.
Also, do not forget to provide clear and specific information about the reanalysis and radiosonde data files that you have used.
Regards,
Juan A. Añel
Geosci. Model Dev. Executive Editor
Citation: https://doi.org/10.5194/gmd-2024-21-CEC4 -
AC5: 'Reply on CEC4', Jihong Zhang, 27 May 2024
Dear Editor,
Thanks for giving us another opportunity and reply very much. We have updated the ZENODO repository with true code instead of .mat binary files (http://doi.org/10.5281/zenodo.11258940).
The specific information about the reanalysis and radiosonde data files that we have used can be obtaind at (https://cds.climate.copernicus.eu/cdsapp#!/dataset/reanalysis-era5-pressure-levels?tab=overview) and (http://weather.uwyo.edu/upperair/sounding.html).
In addition, the “Code and dada availability” has been updated in our manuscript. It is shown below:
Code and data availability. The ERA5 reanalysis data used in this paper can be freely accessed at (https://cds.climate.copernicus.eu/cdsapp#!/dataset/reanalysis-era5-pressure-levels?tab=overview). The radiosonde data can be accessed at (http://weather.uwyo.edu/upperair/sounding.html). All of the data generated during the current study and the code are available on ZENODO (http://doi.org/10.5281/zenodo.11258940).
Title: A hybrid-grid global model for the estimation of atmospheric weighted mean temperature considering time-varying lapse rate in GNSS precipitable water vapor retrieval
Author(s): Shaofeng Xie et al.
MS No.: gmd-2024-21
MS type: Development and technical paper
Again, thanks for your consideration. If there is anything else that needs to be modified or noted, please provide guidance and suggestions, and we will further improve it. I look forward to hearing from you.
Sincerely,
Jihong Zhang
jhzhang@stu.kust.edu.cn
Citation: https://doi.org/10.5194/gmd-2024-21-AC5
-
AC5: 'Reply on CEC4', Jihong Zhang, 27 May 2024
-
CEC4: 'Reply on AC4', Juan Antonio Añel, 24 May 2024
-
AC4: 'Reply on CEC3', Jihong Zhang, 24 May 2024
-
CEC3: 'Reply on AC3', Juan Antonio Añel, 20 May 2024
-
AC3: 'Reply on CEC2', Jihong Zhang, 20 May 2024
-
CEC2: 'Reply on AC2', Juan Antonio Añel, 18 May 2024
-
AC2: 'Reply on CEC1', Jihong Zhang, 14 May 2024
-
RC3: 'Comment on gmd-2024-21', Anonymous Referee #2, 07 Jun 2024
6 June 2024
Comments on “A hybrid-grid global model for the estimation of atmospheric weighted mean temperature considering time-varying lapse rate in GNSS precipitable water vapor retrieval” submitted to Geoscientific Model Development by Shaofeng Xie, Jihong Zhang, Liangke Huang, Fade Chen, Yongfeng Wu, Yijie Wang, Lilong Liu
This submission tries to model the spatial and temporal variation of Tm, which the authors claim to be a key parameter of GNSS PWV inversion. The authors claim they have developed the NGGTm-H model and the NGGTm model. They also compared the performance of their models with that of other three models: Bevis, GPT3-1, GPT3-5 and they claimed that their model has a better performance. Actually many Tm models (more than the Bevis and GPT3) can be found. I don’t know why they didn’t use. The result showed the so-called NGGTm model seems to have a fiddling improvement compared to Bevis and GPT3-1 and GPT3-5 models. However, the actual impact of such a fiddling improvement on the improvement of the GNSS PWV accuracy is unknown, which is the most desirable result the readers want to see. Moreover, the authors frequently stress the necessity of real-time. I know the NGGTm-H model and the NGGTm model are built using multiple-year ERA5 data (2012 to 2016). I am puzzled how the authors can achieve real-time if multiple years’ data are needed. In my viewpoint, the scientific and application contribution of this submission is fiddling or unknown. It has been known that Mr. Yao, Y. B. has conducted a lot of research on the Tm modeling. The authors have cited a large number of his papers too.
Some specific comments:
L15, What does “NGGTm-H” stand for?
L39-40, “Microwave radiometers and satellite remote sensing, which rely on infrared band detection, offer high detection accuracies.” This is logically misleading. Microwave radiometers operate in the microwave region of the electromagnetic spectrum, not the infrared.
L45, “a high spatiotemporal resolution” this is not completely right. Globally speaking, the number of GNSS stations is still low.
L54, “However, the variation law of ZWD influenced mainly by water…” What is the variation law of ZWD? If there is such a law, then the ZWD variation can be well known.
L58, “The accuracy of GNSS tropospheric water vapor retrievals can be significantly improved by using high-precision 𝑇𝑚 data.” How is the accuracy of PWV affected if different Tm values are used? Please provide numerical examples.
L63, “it is necessary to build a real-time and high-precision 𝑇𝑚.” I concur that a high-precision Tm is needed. Is it really necessary to build a real-time Tm model?
L64, “Existing 𝑇𝑚models can be divided into two categories: meteorological parameter models and nonmeteorological parameter models.” What are the representative models for the meteorological parameter models and what are the representative models for the nonmeteorological parameter models? Do you have reference papers?
L85-89, “The non-85 meteorological parameter 𝑇𝑚 model (named the Emardson model) was developed to take the annual cycle variation into account by using radiosonde data collected in Europe over many years, which was capable of meeting the requirement for GNSS water vapor detection (Emardson & Derks, 2000). Therefore, the model has been widely used in real-time GNSS meteorology research.” The Tm model was developed based on the radiosonde data collected in Europe. Does it really meet the accuracy requirement for GNSS water vapor detection out of Europe?
L98-100, “This new model can estimate the 𝑇𝑚 value at any location by simply inputting the station location and the annual product day, which have been applied to real-time GNSS PWV inversion studies worldwide.” Are there any reference papers? What is the accuracy of Tm used in GNSS PWV inversion?
L104-106, “Although the GPT3 model is currently the most representative empirical model with a high precision on the global scale, GPT3 model dose not take into account elevation correction or detailed 𝑇𝑚 lapse rate.” The current GPT3 can get a high precision on global scale, though it didn’t consider the Tm lapse rate. Justify why you need to do this work.
L109, “we aim to global 𝑇𝑚model that takes into account time-varying…” This sentence is grammatically erroneous.
Eqs. (5) and (6), why is the variable KP used? What is the rational of using KP, not Tm?
The paragraph (around L185) discussed the result of Fig. 2. However it is very hard to understand. You should cite the Fig. 2(a)… Fig. 2(f) in the discussion.
L211, “we focused on optimizing the model coefficients solely for these cycles to improve the calculation efficiency when developing the 𝑇𝑚 lapse rate model.” Why is the calculation efficiency so critical? What is the normal calculation efficiency? Is the current calculation not fast or efficient enough?
L219, “Note that, the sliding window algorithm has been used in the previous study, which exhibits a superior performance” What is the difference between the sliding window algorithm in the previous study and that in this submission?
L224, “by using the data of 9 gridded points in each window,” Explain what the 9 gridded points in each window? What is the window? It is better to have figure illustration.
Eq. (9), what is the # of windows for different grids 0.5°x0.5°, 1°x1° and 2°x2°?
L241, “Finally, a global real-time and high-precision 𝑇𝑚 lapse rate model was developed and” How can you get the gamma 𝛾 (the lapse rate of 𝑇𝑚) in real-time? How can you get the 𝑇𝑚𝐺 value at the height of the gridded points from the reanalysis data in real-time? In what applications, the real-time Tm is really needed?
L267, “…daily cycle amplitude, and semidaily cycle amplitude at all grid points using the least-squares adjustment using surface-level gridded 𝑇𝑚 data calculated from all the ERA5 reanalysis data” In L210, you wrote “Since the daily variation in the lapse rate of 𝑇𝑚can be overshadowed by the annual and semiannual variations.” Thus you didn’t model the daily cycle amplitude, or semidaily cycle amplitude at all grid points in the Section 3.2. I can’t understand why you bring up the daily cycle amplitude, and semidaily cycle amplitude in this Section 4.1.
In Eq. (9) of Section 3.2, you estimated the annual cycle amplitude and semiannual cycle amplitude of the lapse rate gamma 𝛾 of 𝑇𝑚. However in the Figure 4 of Section 4.1, you showed the annual cycle amplitude and semiannual cycle amplitude of Tm. I know the lapse rate gamma 𝛾 of 𝑇𝑚 is closely related to Tm. But they are not the same thing. You need to state clearly in the submission what you want to study: lapse rate or Tm.
L276, “In summary, 𝑇𝑚 not only undergoes significant annual and semiannual variations but also experiences significant daily and semidiurnal variation.” Again, it is quite perplexing to me that you stated that significant daily and semidiurnal variations here but you didn’t study it in Section 3.2.
L282, “Since the significant variations in the horizontal direction of 𝑇𝑚 compared to lapse rate, the estimation of 𝑇𝑚 at the gridded points did not use the sliding window algorithm.” It is hard to understand. Rephrase it.
This submission has an unusually high frequency of self-citation (e.g. Huang, L. K., Yao, Y. B.)
Citation: https://doi.org/10.5194/gmd-2024-21-RC3
Viewed
HTML | XML | Total | BibTeX | EndNote | |
---|---|---|---|---|---|
283 | 55 | 51 | 389 | 24 | 21 |
- HTML: 283
- PDF: 55
- XML: 51
- Total: 389
- BibTeX: 24
- EndNote: 21
Viewed (geographical distribution)
Country | # | Views | % |
---|
Total: | 0 |
HTML: | 0 |
PDF: | 0 |
XML: | 0 |
- 1