Evaluation of the High Altitude Lidar Observatory Methane Retrievals During the Summer 2019 ACT-America Campaign
- 1NASA Langley Research Center, Hampton, VA, USA
- 2Science Systems and Applications, Inc., Hampton, VA, USA
- 3Department of Meteorology and Atmospheric Science, and Earth and Environmental Systems Institute, The Pennsylvania State University, University Park, PA, USA
- 1NASA Langley Research Center, Hampton, VA, USA
- 2Science Systems and Applications, Inc., Hampton, VA, USA
- 3Department of Meteorology and Atmospheric Science, and Earth and Environmental Systems Institute, The Pennsylvania State University, University Park, PA, USA
Abstract. The NASA Langley Research Center High Altitude Lidar Observatory (HALO) is a multi-function and modular lidar developed to address the observational needs of NASA’s weather, climate, carbon cycle, and atmospheric composition focus areas. HALO measures atmospheric H2O mixing ratios, CH4 mole fractions, and aerosol/cloud optical properties using the Differential Absorption Lidar (DIAL) and High Spectral Resolution Lidar (HSRL) techniques, respectively. In 2019 HALO participated in the NASA Atmospheric Carbon and Transport – America campaign on board the NASA C-130 to compliment a suite of greenhouse gas in-situ sensors and provide, for the first time, simultaneous measurements of column CH4 and aerosol/cloud profiles. HALO operated in 18 of 19 science flights where the DIAL and Integrated Path Differential Absorption lidar (IPDA) techniques at 1645 nm were used for column and multi-layer measurements of CH4 mole fractions, the HSRL and backscatter techniques at 532 and 1064 nm, respectively, for retrievals of aerosol backscatter, extinction, depolarization, and mixing layer heights. In this paper we present HALO’s measurement theory for the retrievals of column and multi-layer XCH4, retrieval accuracy and precision including methods for bias correction, and a comprehensive total column XCH4 validation comparison to in-situ observations. Comparisons of HALO XCH4 to in-situ derived XCH4, collected during spiral ascents and descents, indicates mean difference of 2.54 ppb and standard deviation of the differences of 16.66 ppb when employing 15 s along track averaging (< 3 km). A high correlation coefficient of R = 0.9058 was observed for the 11 in-situ spiral comparisons. Column XCH4 measured by HALO over regional scales covered by the ACT-America campaign are compared against in-situ CH4 measurements carried out within the planetary boundary layer (PBL) from both the C-130 and B200 aircraft. Favorable correlation between the in-situ point measurements within the PBL and the remote column measurements from HALO elucidates the sensitivity of a column integrating lidar to CH4 variability within the PBL, where surface fluxes dominate the signal. Novel capabilities for CH4 profiling in regions of clear air using the DIAL technique are presented and validated for the first time. Additionally, profiling of CH4 is used to apportion the PBL absorption from the total column and is compared to previously reported IPDA cloud slicing techniques that estimate PBL columns using strong echoes from fair weather cumulus. The analysis presented here points towards HALO’s ability to retrieve accurate and precise CH4 columns with the prospects for future multi-layer profiling in support of future suborbital campaigns.
Rory A. Barton-Grimley et al.
Status: final response (author comments only)
-
RC1: 'Comment on amt-2022-106', Anonymous Referee #1, 20 Apr 2022
Review of Evaluation of the high altitude lidar observatory methane retrievals during the summer 2019 ACT-America campaign by Barton-Grimley et al.
General comments.
First I would like to mention that this paper show impressive work and results with respect to XCH4 estimate using lidar IPDA and even DIAL methods. The paper is well written and the figures are clear and well detailed which make easy the understanding of the measurements.
Although previous XCH4 airborne IPDA measurements with CHARM-F is mentioned some discussion of performances (precision, resolution, biaises) with respect to CHARM-F and maybe future validation of MERLIN space mission is missing. The paper is quite long and one may think that the HSRL measurements, as much as a detailed geophysical analysis of XCH4 data should be kept for a second paper.
As I said the paper seems to me a little long, however, in the same time, a fundamental part is missing: the consideration and correction of statistical biases in DAOD and DIAL estimates. Such omission can lead to misunderstanding of differences between in situ and lidar measurements. Therefore, some of the results in this paper should be calculate again and corrected and this is the reason why I indicated "major revision". I recommend to make the following correction before consideration for publication.
Specific comments.
- Instrument:
- L162. Note that main difference with CHARM-F is the OFF wavelength 1645.86 nm -> 1645.37 nm on the other side of the CH4 multiplet. Recent spectroscopic data on H2O absorption lines [Delahaye JQRST 2019] show that minimising H2O impact on CH4 DAOD requires to use the OFF-line at 1645.86 nm. This should be maybe indicated or at least taken into account if the authors plan to contribute to a future validation of MERLIN CH4 space lidar mission
- Opical depth biais correction:
- L334 and Figure 6c. Why the biais correction depends on the gain Figure 6c ? I don’t see a potential explanation in all is described L334-342.
Moreover, as the main difference of biais correction is shown for the LOLE channel (with the lowest SNR) one may wonder if the necessary correction of averaged optical depth with low SNR have been taken into account in the signal processing (Tellier et al. AMT, 2018) ?
The reviewer notes that nothing is said about statistical biaises in IPDA/DIAL in theoretical paragraph 2.2 and further in the paper which is not acceptable.
If no statistical biais is considered, the a posteriori biais correction used by the authors in paragraph 3.1.1 is clearly SNR dependent and this should be indicated even if it is not obvious in the DAOD estimates.
- L406. Precision performances of HALO should be compared and discussed with respect to Amediek et al. AO 2017 paper and CHAM-F results.
- Figure 9 and L422. XCH4 noise statistic decreases less than the square root law. We can find a similar result in Amediek et al. 2017. The reason that is suggested by the authors is « harsh operating condition in the C-130 » and « high vibrational environment » which is fully possible and may entail a « degraded laser frequency stability » … what about optical misalignment issue? did the authors make some vibrational tests of their system?
- Regional scale observations
- Figure 13. and L530-535. The unexpected result of larger IPDA XCH4 than PBL in situ CH4 is very unusual. This shows that in situ data measured both in the free troposphere and in the PBL may not be sufficient to make a validation of space-based measurement such as MERLIN. Co-located airborne measurement and maybe XCH4 profiling with ground-based lidar should be used to explain such enhancement of CH4 in the free troposphere.
- I think that the second part L550 to L600 is not necessary in this paper (although really interesting!). Also, HSRL measurements seem not to be so essential in this paper as the authors proved that 1.645 µm backscatter is sufficient to give the vertical structure of the atmosphere and even, I guess, the height of the PBL.
- Advanced CH4 products - Atmospheric profiling
L 646- 663. and Figure 17
- L 646. The estimate of the DAOD profile is confusing. Did the author slice average the backscattered profiles to 350 m and 15 s first before using equation 6 ? It does not look this way given the variations of DAOD profile in Figure 17b… Signal processing should be clarified here.
- A decreasing DAOD is of course not expected and I agree that this may be a manifestation of low SNR. I have then the same question as for IPDA measurement: did the authors make an estimate of the statistical biais on the DAOD (and thus a correction) due to the non-linearity of equation 6? this question is linked to the question just above giving that an averaging enables to increases the SNR and then to decrease such biais… once again please read Tellier et al. AMT 2018 but this issue was also mentioned in early DIAL measurements with high precision such as for CO2 (Gibert et al. JTECH 2008)
As for IPDA, a correction of DAOD with statistical biais is a basis for modern DIAL measurements and the authors should includes and discuss in details the impact on SNR on their measurements. This is to my mind essential.
However the authors should be aware that the correction of DAOD with SNR linked biais might not be sufficient to remove entirely the decrease of DAOD seen in Figure 17b. At low SNR, especially for ON line signal the impact of Pb removal in Equation 1 may entail other issue due to the electronic baseline and linearity of the detection.
- L 658. A linear regression on the DAOD that is not weighted by error bars on each DAOD point is biased for the reason mentioned above and non linearity of equation 6. Gibert et al. AO 2006 used such likelihood estimate to make accurate XCO2 measurement in the PBL. In Figure 17b the DAOD will then not impact so much a likelihood calculated slope coefficient and I expect that there will be a better agreement with IPDA and in situ DAOD.
In conclusion the difference is not at this point explained by spectroscopy as the authors wrote (this sentence should be removed) but clearly by the non consideration of statistical biaises in their estimates.
- L710 - 720. Of course what is mentioned above should be considered in all this paragraph, i.e. the different statistical biaises should be corrected before the comparison of PBL XCH4 using the cloud slicing method and the DIAL profiles.
As already said before, spectroscopy induced error should be mentioned, if necessary, only in a second step.
Technical corrections
- L390. please remove one « is » in the sentence.
- L696. Please add an error bar for each retrieval IPDA ground, cloud, PBL
- 5.1 paragraph. As there is no 5.2 paragraph I guess that this title should be removed
- Conclusion must be re-written in agreement with statistical bias corrected results.
-
AC1: 'Reply on RC1', Rory Barton-Grimley, 28 Jun 2022
The comment was uploaded in the form of a supplement: https://amt.copernicus.org/preprints/amt-2022-106/amt-2022-106-AC1-supplement.pdf
-
RC2: 'Comment on amt-2022-106', Anonymous Referee #2, 05 May 2022
This paper describes the first results from a lidar system deployed onboard a research aircraft measuring atmospheric methane. The paper is well-written, and fits well within the scope of AMT. However, a few issues listed below should be addressed before the paper can be recommended for publication.
General comments:
Note that the term “mole fraction” is recommended rather than “mixing ratio”, see e.g. https://www.empa.ch/web/s503/gaw_glossary#recommendations. I suggest simply replacing throughout the text. Also I would recommend to consistently use ppb for dry air mole fractions of CH4. Using both ppb and ppm (e.g. Fig. 10 (b) ) is confusing to the reader.
Comparison to in-situ measurements: The deployment of the different aircraft sampling different altitude regimes really has potential, as indicated e.g. by Figs. 12 and 13 and the associated discussion. I suggest a simple combination of the in-situ measurements within the free troposphere from the C-130 aircraft, the boundary layer in-situ measurements from the B200 aircraft, and the estimate of the boundary layer height derived from the HSRL measurements onboard the C-130, to calculate a partial column XCH4 based on in-situ observations that can directly be compared to HALO XCH4. The assumption is that CH4 is well mixed within the PBL and also within the free troposphere. Any advection of air masses with enhanced CH4 above the PBL would clearly stick out as differences between HALO XCH4 and the aircraft derived XCH4.
Dry air mole fraction - impact from H2O: In the in-situ measurement community there is much discussion on drying/conditioning samples before measurement vs. correcting based on simultaneous H2O measurement within the exact same sample. As the authors describe, MERRA humidity is used in the retrieval of XCH4 (the column average dry air mole fraction). The uncertainty in XCH4 introduced by this choice should be assessed, e.g. by comparing MERRA water vapor to that of the in-situ observations.
Specific comments:
L271 “over samples” -> “oversamples”
L280 “Altitude is used in lieu of MSL for all figures” this is not clear. May be “Altitude above MSL is used in lieu of Altitude above ground level”?
L280: “post- flight reanalysis” may be drop “post-flight”? I guess reanalysis products are available only for past periods, i.e. after the flights, anyway.
L331: “The superscript will be dropped for simplicity.” Which superscript?
L362: the matrix T should contain the elements that the beta-vector elements are multiplier with, i.e. 0th, 1st 2nd and 3rdorder terms as formulated in the Eq. on line 361. May be simply write down the first and last row of the matrix, and the few elements
Fig. 6: please adjust color selection for the different gains so that color blind people can read the figure. To me HOLE and LOLE are identical, LOLE is very slightly different.
Fig. 6 caption: please explain DEM (I know what it is, but it should be mentioned once)
Fig. 11: the Y-intercept is not clear. It should be negative, given the slope is larger than one, and the regression line crosses the 1:1 line at around 1900 ppb.
L522: “PA region” – to make this clear to non-US readers (AMT it is a European journal) may be add a label to Figs 11 (a) and (c)
L571: I don’t see any cross-hatched area in Fig. 14, may be I am misunderstanding something
L612 “spatial” -> “spatially” or drop
L659: the dial DAOD estimated at SSE shown in the inset of Fig. 17 (b) (magenta symbol) is around 0.2875, not at 0.9243 as given in the text. Please clarify.
L661: “un-bias-corrected” may be use non-bias-corrected
L736 “PBL fluxes” use PBL mole fractions or concentrations
-
AC2: 'Reply on RC2', Rory Barton-Grimley, 28 Jun 2022
The comment was uploaded in the form of a supplement: https://amt.copernicus.org/preprints/amt-2022-106/amt-2022-106-AC2-supplement.pdf
-
AC2: 'Reply on RC2', Rory Barton-Grimley, 28 Jun 2022
Rory A. Barton-Grimley et al.
Rory A. Barton-Grimley et al.
Viewed
HTML | XML | Total | BibTeX | EndNote | |
---|---|---|---|---|---|
341 | 122 | 12 | 475 | 6 | 5 |
- HTML: 341
- PDF: 122
- XML: 12
- Total: 475
- BibTeX: 6
- EndNote: 5
Viewed (geographical distribution)
Country | # | Views | % |
---|
Total: | 0 |
HTML: | 0 |
PDF: | 0 |
XML: | 0 |
- 1