the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Tomographic reconstruction algorithms for retrieving two-dimensional ice cloud microphysical parameters using along-track (sub)millimeter-wave radiometer observations
Abstract. The submillimeter-wave radiometer operating in the along-track scanning mode continuously collects brightness temperature (TB) data over a two-dimensional (2D) cloud cross-section as the platform moves forward. TB observations from multiple positions and viewing angles show great promise in better constraining the 2D cloud microphysical properties compared to single-angle observations. In this study, we develop two types of tomographic reconstruction algorithms to retrieve 2D ice water content (IWC) profiles using multi-angle TB observations. The one-dimensional (1D) tomographic algorithm performs 1D retrievals beam by beam using each TB observation at a specific position and angle to derive cloud properties along the propagation path. It then integrates the 1D retrieval results to construct 2D cloud distributions. The 2D tomographic algorithm directly constrains the 2D cloud microphysical properties using multi-angle scanning TB observations. Starting with an initial assumption, the algorithm iteratively refines the 2D cloud microphysical quantities by minimizing discrepancies between TB simulations and observations under prior constraints. Both tomographic algorithms are developed based on a hybrid of Bayesian Monte Carlo Integration (MCI) and Optimal Estimation Method (OEM). A simulation experiment is conducted to evaluate the performance of two tomographic reconstruction algorithms. The experiment demonstrates stable convergence of both tomographic methods, with the 2D tomographic algorithm exhibiting superior performance. The experiment results highlight the significant advantage of using multi-angle observations to constrain 2D cloud structure. Compared to nadir-only retrievals, the tomographic technique provides a detailed reconstruction of ice clouds’ inner structure with high spatial resolution. Also, the technique significantly improves retrieval accuracy by correcting systematic biases and reducing the derivation of retrieval errors. Furthermore, the tomography technique effectively increases detection sensitivity for small ice cloud particles.
- Preprint
(4524 KB) - Metadata XML
- BibTeX
- EndNote
Status: final response (author comments only)
-
RC1: 'Comment on amt-2024-188', Anonymous Referee #1, 21 Dec 2024
General comments
The manuscript provides a novel sub-millimeter ice cloud retrieval algorithm that uses radiometer along-track scanning simulation data and a hybrid Bayesian MCI and OEM optimization approach to reconstruct two-dimensional ice cloud scenes. The manuscript is well-written and rich in information. I believe the manuscript deserves to be published in AMT. I only have a few minor issues that need to be explained by the authors.
Specific comments
- Line 79. Is the current CoSSIR not capable of along-track scanning?
- Line 99. Since CloudSat is not scanning, meaning that an orbit of 2C-ICE is striped, how to understand the 3D scene here?
- Line 101. How to understand the IBA technique here, whether it is by averaging 1D simulated observations within the field of view or by extracting multiple atmospheric profiles at different locations on the path to synthesize a scene?
- Figure 5. What spatial extent of cloud parameters are reconstructed in one along-track scan, is it the entire scan? How do you handle the front and back scans where there is a crossover?
- Line 190. If the BMCI performs successfully, does the OEM not execute, why the minimum database case is 25?
- Line 197. Does the local Gaussian a priori constrain mean the a prior IWC and covariance matrix used for OEM?
- Line 197. That is, even if the BMCI fails, at the start of the OEM, it still has to find 25 cases in the database to update the a posteriori PDF. I'm a little confused here.
- Line 237. How to understand auto-adjustable BMCI? Is it any different than the BMCI in 1D?
- Line 265. Computational resources and speed are really an issue with this algorithm, roughly how much computational resources and computation time are needed at the moment?
Citation: https://doi.org/10.5194/amt-2024-188-RC1 -
AC1: 'Reply on RC1', Yuli Liu, 15 Jan 2025
The comment was uploaded in the form of a supplement: https://amt.copernicus.org/preprints/amt-2024-188/amt-2024-188-AC1-supplement.pdf
-
RC2: 'Comment on amt-2024-188', Anonymous Referee #2, 23 Dec 2024
This is a well written manuscript that illustrates the benefits of performing 2D tomographic retrievals as a way of improving on historical 1D retrievals. While the results between 1D-nadir, 1D-tomographic and 2D-tomographic retrievals could have been anticipated, the value of this paper is that it clearly outlines viable approaches for the 1-D and 2-D tomographic inversions. The treatment of uncertainties and correlated uncertainties, while perhaps not the only way to deal with the problem, is sound and pushes our treatment of correlations length in OEM retrievals beyond our current methods. I think this paper is definitely worth publishing since it not only explores a new technique, but also opens up the possibility of occasionally flying cross-track scanning satellites sideways in order to obtain better microphysical insight as suggested here.
Â
Specific issues the authors should address
Â
Line ~110 – The results presented here depend to some extent on the details of the particles that were assumed to make up the liquid and ice clouds, and how supercooled water was handled. To make the results more reproducible, I would encourage the authors to make their PSD assumptions explicit. Because the profiles are artificially simple, the authors should also add a caveat to the conclusions. It may be that changes in PSD properties or the inclusion of super-cooled water in the cloud may be another reason to perform 2D tomography but that has not been demonstrated.
Â
Line ~140 – the suggestion is that Tb are computed for each CloudSat profile at multiple angles needs a little more explanation. Because slant path computations cut through multiple horizontally adjacent profiles, this leads to uneven layering of the slant path profiles. Exactly how the authors handled this was not clear from the description.Â
Â
Line ~150 – If Tbs are constructed from single profiles but simply for different view angles, then the vertical correlations are assumed in the prior data. Some discussion here about how this is ultimately handled in the retrieval would be appropriate here.
Â
Line 220 – Simply averaging the same voxels retrieved from each view angle in the 1-D scheme seems to artificially simplistic. Given that the authors have the goodness of fit from the OE, should this not be a weighted average?
Citation: https://doi.org/10.5194/amt-2024-188-RC2 -
AC2: 'Reply on RC2', Yuli Liu, 15 Jan 2025
The comment was uploaded in the form of a supplement: https://amt.copernicus.org/preprints/amt-2024-188/amt-2024-188-AC2-supplement.pdf
-
AC2: 'Reply on RC2', Yuli Liu, 15 Jan 2025
Viewed
HTML | XML | Total | BibTeX | EndNote | |
---|---|---|---|---|---|
182 | 39 | 7 | 228 | 8 | 8 |
- HTML: 182
- PDF: 39
- XML: 7
- Total: 228
- BibTeX: 8
- EndNote: 8
Viewed (geographical distribution)
Country | # | Views | % |
---|
Total: | 0 |
HTML: | 0 |
PDF: | 0 |
XML: | 0 |
- 1