the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
RICHARD 1.0 – Routine for the Isolation of Chemical Hotspots in Atmospheric Research Data
Abstract. Here, we introduce version 1.0 of the RICHARD algorithm, a Routine for the Isolation of Chemical Hotspots in Atmospheric Research Data, e.g. in satellite measurement datasets (level-2 and above) or atmospheric chemistry model output. The overall goal of the algorithm is to identify "hotspot" areas in which local enhancements of an atmospheric constituent can only occur due to strong local emissions. To detect hotspot areas, we use a mathematical method that combines spatiotemporal proxy data for calibration purposes and a selection algorithm in a novel way. For each input file family (e.g. the near surface mixing ratio or the tropospheric partial column of atmospheric constituents as a function of longitude and latitude) we define a structure quotient by which the algorithm decides - based on a threshold value - whether at a particular iteration step the hotspot area criteria are met and if they are kept for further calculations or not. The python based command line tool RICHARD 1.0 comes with a set of implemented features like an automated generator for user-defined patterns and an analysis tool to determine the optimal threshold value for a given dataset.
For testing purposes of RICHARD 1.0, we use simulations of the atmosphere and chemistry modeling framework ICON-ART, a joint development of the German Weather Service and Max-Planck-Institute for Meteorology in Hamburg. We comprehensively explore different aspects of RICHARD using ICON-ART model output datasets. We provide an analysis of the decision making process coded in RICHARD, and provide a detailed look at the competing effects of emissions and advection. Here, we also consider the direction and speed of the wind that affect the advection of prescribed (and thus known) emissions in the model and look at the resulting tracer mixing ratios as to evaluate the sensitivity of the algorithm and its ability to identify objectively hotspots of strong emissions, based on the self-determined threshold values.
The results show that RICHARD can identify frequently (or continuously) emitting localised sources as hotspots. Furthermore, the algorithm is able to distinguish between an actual emission source and other circumstances that lead to enhanced tracer concentrations, e.g. as caused by wind conditions and associated transport processes. In addition, the in ICON-ART prescribed emission source strengths are detected and quantified regardless of overlying transport features with only a small error of about 5 %. This increases significantly the accuracy of determined source strengths compared to other methods that we have explored.
RICHARD 1.0 is a novel comprehensive tool for the identification and quantification of emission hotspots and uses a novel workflow that includes spatiotemporal proxy data as well as a selection algorithm. Here, we present a model-based proof of concept that is already fully transferable to applications using satellite measurement data.
This preprint has been withdrawn.
-
Withdrawal notice
This preprint has been withdrawn.
-
Preprint
(1055 KB)
Interactive discussion
Status: closed
-
RC1: 'Comment on gmd-2023-91', Anonymous Referee #1, 17 Jul 2023
-
AC1: 'Reply on RC1', Christian Scharun, 30 Sep 2023
Dear referees,
I want to express my sincere gratitude for your time and effort in reviewing this paper. Your valuable feedback has been instrumental in helping me gain insights into possible areas for improvements, and I genuinely appreciate your dedication to enhancing the quality of my paper and its presentation. While I understand that the presentation of my work is rather fundamental and “bare” (and thus did not meet your expectations), the methodology itself delivers novel aspects and insights that certainly require further testing and refinement in “real world cases”. However, this is something that I will not be able to deliver, due to a change in career path. The code and documentation will still be available to interested parties and I am optimistic that revisions and adjustments (as suggested by you and worked on by a successor) will result in a more comprehensive and compelling manuscript that will advance the current idealized proof-of-concept case study to a practical tool directly applicable to real data.
All the best
Christian ScharunCitation: https://doi.org/10.5194/gmd-2023-91-AC1
-
AC1: 'Reply on RC1', Christian Scharun, 30 Sep 2023
-
RC2: 'Comment on gmd-2023-91', Anonymous Referee #2, 01 Sep 2023
The paper introduces the RICHARD 1.0 algorithm written in Python that aims to identify hot spots in satellite measurements and atmospheric transport simulations. While such a tool would be of scientific importance and would contribute to the scope of GMD, the scientific quality, reproducibility and presentation quality of the manuscript are poor . Therefore, I recommend that the manuscript should be rejected.
Some general comments:
The introduction is not well written. The cited literature is incomplete and it remains unclear how the cited examples motivate the development of the RICHARD algorithm. Furthermore, the introduction has several wrong statements. For example, Dumont Le Brazidec et al. (2023) do not analyze CO2 fields from the TROPOMI instrument but synthetic CO2M observations. TROPOMI does not have a CO2 product. While Dumont Le Brazidec et al. (2023) introduce a new loss function, this is not done to address the issue of low signal-to-noise (L41f), but to define the extent of the true CO2 plume in training dataset. Furthermore, instrument noise and background fields are added to the simulated CO2 plumes to generate the synthetic CO2M observations (L42f).
The description of the RICHARD algorithm (Section 2.1) is unclear. The algorithm starts with two patterns (P1 and P2), but it is not explained how they are initialized from the dataset. Based on Section 2.2, it seems to be related to a list of sources provided to the program. The patterns are used to define the structure quotient. An iterative step is mentioned to "drop the ones", but it is never explained how this is implemented and it is unclear what "the ones" refers to.
Since the goal and method used by the algorithm is unclear, I feel not able to review the results section.
Citation: https://doi.org/10.5194/gmd-2023-91-RC2
Interactive discussion
Status: closed
-
RC1: 'Comment on gmd-2023-91', Anonymous Referee #1, 17 Jul 2023
-
AC1: 'Reply on RC1', Christian Scharun, 30 Sep 2023
Dear referees,
I want to express my sincere gratitude for your time and effort in reviewing this paper. Your valuable feedback has been instrumental in helping me gain insights into possible areas for improvements, and I genuinely appreciate your dedication to enhancing the quality of my paper and its presentation. While I understand that the presentation of my work is rather fundamental and “bare” (and thus did not meet your expectations), the methodology itself delivers novel aspects and insights that certainly require further testing and refinement in “real world cases”. However, this is something that I will not be able to deliver, due to a change in career path. The code and documentation will still be available to interested parties and I am optimistic that revisions and adjustments (as suggested by you and worked on by a successor) will result in a more comprehensive and compelling manuscript that will advance the current idealized proof-of-concept case study to a practical tool directly applicable to real data.
All the best
Christian ScharunCitation: https://doi.org/10.5194/gmd-2023-91-AC1
-
AC1: 'Reply on RC1', Christian Scharun, 30 Sep 2023
-
RC2: 'Comment on gmd-2023-91', Anonymous Referee #2, 01 Sep 2023
The paper introduces the RICHARD 1.0 algorithm written in Python that aims to identify hot spots in satellite measurements and atmospheric transport simulations. While such a tool would be of scientific importance and would contribute to the scope of GMD, the scientific quality, reproducibility and presentation quality of the manuscript are poor . Therefore, I recommend that the manuscript should be rejected.
Some general comments:
The introduction is not well written. The cited literature is incomplete and it remains unclear how the cited examples motivate the development of the RICHARD algorithm. Furthermore, the introduction has several wrong statements. For example, Dumont Le Brazidec et al. (2023) do not analyze CO2 fields from the TROPOMI instrument but synthetic CO2M observations. TROPOMI does not have a CO2 product. While Dumont Le Brazidec et al. (2023) introduce a new loss function, this is not done to address the issue of low signal-to-noise (L41f), but to define the extent of the true CO2 plume in training dataset. Furthermore, instrument noise and background fields are added to the simulated CO2 plumes to generate the synthetic CO2M observations (L42f).
The description of the RICHARD algorithm (Section 2.1) is unclear. The algorithm starts with two patterns (P1 and P2), but it is not explained how they are initialized from the dataset. Based on Section 2.2, it seems to be related to a list of sources provided to the program. The patterns are used to define the structure quotient. An iterative step is mentioned to "drop the ones", but it is never explained how this is implemented and it is unclear what "the ones" refers to.
Since the goal and method used by the algorithm is unclear, I feel not able to review the results section.
Citation: https://doi.org/10.5194/gmd-2023-91-RC2
Viewed
HTML | XML | Total | BibTeX | EndNote | |
---|---|---|---|---|---|
430 | 116 | 38 | 584 | 34 | 35 |
- HTML: 430
- PDF: 116
- XML: 38
- Total: 584
- BibTeX: 34
- EndNote: 35
Viewed (geographical distribution)
Country | # | Views | % |
---|
Total: | 0 |
HTML: | 0 |
PDF: | 0 |
XML: | 0 |
- 1