I commend the authors for addressing all my comments -- which was quite some work! I am satisfied with most of the revisions and I feel the results are very robust now.
I have a number of remaining minor/technical comments, and only a couple of bigger 'complains' from my side. While due to the overall amount of comments it almost sums to another major revision, it all should be nevertheless fast and straightforward to address.
#
# My two remaining complains:
1) The continuous color scale in most figures, would be much easier to interpret and guess the numbers by eye, if the color scale had discrete steps.
2) Too much focus and speculation about the MJO and related topics in the discussion, while I can't notice any MJO-like signal in any of the figures.
There are numerous specific comments related to 1) and 2) listed below:
#
# Specific comments:
- l.28: specify the study period.
--> same for the 1st paragraph of section 2.1, the analysis time period should be specified also there.
- l.29: remove "referred to as SO3 in"
- About the GWD scheme used in ERA5, this should be mentioned again in the 2.1 subsection and in your discussion regarding Fig.3.
- l.83-84: perhaps describe more specifics of RBS in your region of study in the previous sentence (like what resolution exactly it has in the UTLS height range), and here just state that RBS are different in the extratropics/polar regions.
- l.104: interpolated how? (refer to methods subsection)
- l.124-125: reference needed
- l.132-133: to make the sentence more concise, substitute "is intentional. This... by treating it as" --> serves as a comparison with
- l.167-168: Ern et al (2023) proved this, also Zagar et al. (2025) detail more about Kelvin waves and shear.
- l.226-227: "The lower bound is set..." --> sentence can be removed without loss of any info
- l.228: tropopause height, specify how it is defined, e.g. cold-point from model levels?
- l.278-279: cite Lux et al (2022) for this
- Figure 2: please improve the color scale, make separate colors every 1 or 2 J/Kg -- i.e. make it discrete, not continuous -- add another color beyond yellow.
- l.314: "Boreal spring 2019 to Austral summer 2020" is an extremely awkward way to state your analysis period. Simply JJA 2019 to MAM 2021 does the job, without confusing the reader.
- l.315: patterns, describe which ones.
- l.318-319: --> The reader should note that also some variance from equatorial waves, centered at the equator (by definition), will be inevitably present to some small degree.
- l.324-325: In ERA5 a lot is missing especially in the active monsoon regions // further away from the Equator, please note this in the text.
- l.334-338: Again, I feel seasonal monsoon convection is missing in this discussion.
- l.349-356: Perhaps mention that this will be looked in more detail in the next subsection in relation to OLR.
--> Figure 3 confirms your GW hotspots follow deep convective systems
- Figure 3:
--> Especially in panels a-b, I see the same problem as with Fig.2: please make the color scale discrete (in every panels would be best)
--> Longrange --> Longwave!
--> I think there might be some mask where positive values (red) are not completely transparent to the stippling, so it appears light gray.
Stippling within blue regions looks fine.
I assume the stippling is intended to be the same everywhere
- l.373: "the observations" --> more precision needed here -- Aeolus and ERA5 (HLOS) GW and their differece.
- l.375-377: need to mention the OLR patterns that coincide with your GW hotspots
- l.386-388: again your figure shows this is the case -- GW hotspots following low OLR regions, please discuss this and state the OLR values shown in the figure.
- l.397-399: Better to simply state it's a good estimator of cloud top temperature and thereby convection depth.
- l.400: "primary weakness"
--> not sure about calling it a primary weakness...
These convective GW's are not parameterized basically, from what I understand from your introduction?
This result from Fig. 3 highlights the need to parameterize convectively generated GWs better.
- l.410-411: Either be more specific about what metrics you refer to, or remove this sentence (no info lost if removed)
- Figure 4: Please improve color scale with a discrete separation of colors every 0.1.
As it is now, it's impossible to tell by eye where exactly the 2 ratio is, or whether it's 1.8 instead.
- l.423-424: the caption about OLR --> Can be shortened to "white and black contour lines represent 210 and 220 W/m2 OLR, respectively."
- l.438: Oscillations --> Oscillation (no plurals with MJO)
- l.447: dynamically active, you mean convectively?
- l.456: each line--> each row
- l.457: I don't see any white bins
- Figures 5-6-7: same as previous figures, please make a discrete color scale.
- l.489: "lesser convective areas" --> you mean non-convective?
Perhaps consider adding a contour line with high OLR indicating stable conditions into Fig. 6.
- l.506-512: I'd like to see a longer discussion comparing it with Fig. 3 in terms of general values, peaks and sign of the differences
--> Good agreement with ERA5 in Ep, but I see lots of (light red) color in Fig.6c.
--> Poor agreement with ERA5 in Fig.3 --> clear underestimation of Ek in convective regions (red colors), but what about the general blue color (although insignificant?) elsewhere.
-l. 535: "including orographic influences": there is no orography at 200 lon
- l.541: "RO temperature data" --> maybe simply state "RO measurements".
what is actually being assimilated are GNSS-RO bending angles (which contains the temperature information within).
- l.558: to me inconsistent implies lots of ups and downs.
--> perhaps decaying performance is more suitable here?
- l.564-565: Seasonal cycle or monsoon could be as well.
- l.560-573: Whole MJO discussion is very speculative and not backed by any result of yours. What one can see clearly in Figs. 2 and 3 is the seasonal cycle.
MJO has a timescale of 30-90 days, meaning full positive-negative phase within 3 months at most, moving form the Indian Ocean into the Pacific: sorry but I don't see anything on that timescale in these figures.
- l.574-583: this is actually an important result and should belong in an earlier result section and not here.
Expand the discussion in lines 230-235 with the infos from this paragraph.
You may refer to this in the discussion/conclusions in a summarized manner later.
- l.601-607: Newer papers from Zagar on Aeolus and equatorial waves (2021, and especially 2025) would help with this discussion.
The IFS has evolved quite a lot from 2004.
- l.612-620: too speculative, was there any remarkable MJO event during your analysis time period? I again suggest to focus the discussion more on seasonal monsoon convection, and existing results on Aeolus capturing equatorial waves (e.g. the previously mentione Ern and Zagar papers), and the QBO (Banyard).
Also, spending almost an entire paragraph on QBO-MJO modulation is going off-topic.
- l.632-633: No works that are more recent? It's surprising to me, that the newest reference here is 10y old and that there were no follow-ups (although I'm not an expert in this particular field)
- l.645: "increasingly challenged by observations" --> references please, or say "(see references above / previous paragraph)"
- l.652-653: refer to the section where you come up with these numbers
- l.657-665: regarding the hot pixels, wasnt this taken care of in the latest reprocessing that you use?
Even if you state later that the latest baseline improves this, the current phrasing implies it's still a major source of uncertainties -- which is not really the case right?
--> I overall don't like the fit of this sentence for the study's research focus, I suggest to remove even the whole paragraph -- it's mixing many things in a very vague way, without any way forward.
--> Oscillating perturbations misiterpreted as GW signals (Ratysnki etal 2023) can be mentioned in the Data/Methods section to give more nuance on the reliability of the dataset.
- l.676: include lower bound of the wavelengths that are band-pass filtered. Also, band-pass is the more common expression (instead of passband)
- l.687-693: please also mention Aeolus' swath size vs whole globe -- small scales are only sampled very locally every day --> how many equator crossings every 12h, I believe it's in the range of 14-16?
# References
#
Ern et al (2023), The quasi-biennial oscillation (QBO) and global-scale tropical waves in Aeolus wind observations, radiosonde data, and reanalyses , https://doi.org/10.5194/acp-23-9549-2023
Lux et al. (2022), Quality control and error assessment of the Aeolus L2B wind results from the Joint Aeolus Tropical Atlantic Campaign, https://doi.org/10.5194/amt-15-6467-2022
Zagar et al. (2025), ESA’s Aeolus mission reveals uncertainties in tropical wind and wave-driven circulations, http://dx.doi.org/10.1029/2025GL114832 |
The manuscript presents an analysis of gravity wave (GW) kinetic energy distributions, derived from new Aeolus satellite wind profiles, that shows great promise in pushing the needle forward in the construction of observational constraints of gravity waves and their impacts on upper troposphere/lower stratosphere circulation. A methodology is presented for deriving the kinetic energy associated with small-scale GWs in regions of deep convection in the tropics over a period spanning June 2019 to August 2022. Comparisons with ERA5 suggest that the reanalysis product underestimates GW-associated kinetic energy; conversely, GW-associated potential energy comparisons between ERA5 and temperature-profiles from an independent instrument (GNSS-RO) show much more consistency, suggesting that the use of kinetic energy highlights a distinct feature of the GW energy spectrum that is not typically assessed (and, incidentally, is not well represented in ERA5). The authors further speculate that this underestimate may reflect lack of assimilated direct wind observations, in contrast to temperatures, which are assimilated. All in all, the manuscript does a good job of presenting a new dataset with all necessary caveats, while also making a generally convincing case that this new data will be valuable. To this end, I recommend acceptance, pending that minor revisions be made to address the following concerns:
#1. Page 4: There is no description of the GW drag parameterization employed in ERA5. In particular, does the model have an explicit parameterization for non-orographic GW drag due to parameterized convection? If so, what is it and how has it been evaluated/performed in past assessments? This will be important in terms of interpreting the dearth of kinetic energy in the model, relative to the Aeolus-derived energy.
#2. Page 7: Presumably the definition of "background" based on "the arguments presented in Alexander et al. (2008b)" apply to past analysis of temperature, not wind, profiles, no? More generally, it would be good for the reader to have a better sense of the sensitivity of the profiles depicted in Figure 1a to choice of grid box averaging domain, the temporal period over which profiles are averaged (currently set to 7 days, etc.), etc. I imagine the authors have already done this sensitivity analysis, so they could consider showing in an appendix figure.
#3. Equation (1): This notation becomes slightly confusing/counterintuitive as the text moves on, since the meridional component often goes to zero due to the pointing vector retaining its approximate angle at ~100 degrees. In other words, V_HLOS would be more intuitively referred to as U_HLOS (or something similar) since, indeed, it primarily reflects the zonal component of the flow. Is there any particular reason why "v" is used instead of something more generic? I suggest changing.
#4. Figure 14, lines 354-355: The first sentence of this paragraph does not make sense to me. In particular, the bit referring to "ERA5 shows a considerable reduction" is vague. Reduction relative to what? Please clarify.
#5. Figure 5: The temporal resolution labeled on the y-axes of these hovmoller plots is too high/unnecessary as it crowds the figures. Please show only every other two or three months. Same comment applies to Figure 7.
#6. Figure 16, Discussion concluding Section 3.2: The discussion here seems weak and understates the disagreement between the Aeolus and ERA5 Ek temporal patterns. The second-to-last paragraph highlights the common features between Aeolus and ERA5, but I think the plots look very different. In particular, the hotspots coincident with low OLR are totally missing in ERA5 (Fig. 5b). The phrasing in the text, however, seems to suggest that the differences are only minor. Please rephrase.
#7. Section 4: Doesn't the ratio of Ek/Ep (shown for ERA-5 in Fig. 8a) suggest that these two quantities are extremely different and not meaningful to compare with each other? I appreciate that the authors want to move beyond traditional (conservative) analysis and attempt to do a bit more, but Figure 8a suggests that the two quantities are in much more disagreement than the discrepancy predicted by llinear wave theory (i.e., factor of 4, not factor of 2). My suggestion here is to introduce Figure 8a earlier as a way to more directly address the concerns with comparing potential and kinetic energy (within a self-consistent product like ERA-5).
#8. Last paragraph on page 19 (lines 467-470): How do you know it's the failure to assimilate the winds directly that's causing the poor representation of GW-associated EK? In principle, one might be able to capture these features using a convective non-orgraphic gravity wave drag parameterization within the ERA-5 model, no? In other words, the assimilation is one way to correct the problem, but an alternative approach is to tackle the model bias directly. However, without having more knowledge about the underlying GW drag parameterization in the model it's hard for the reader to know how many degrees of freedom are afforded to the modeler. Can the authors please comment on the role played here by model bias? And how this is/is not handled by the GW drag parameterization?
#9. Discussion: No mention is made of how these observations might be used to develop constraints on the momentum fluxes (which is what modelers seek most). Is that something that the author has considered? This is a challenging question, so I am not seeking any complete answers here; I am just wondering if the author can speculate in a sentence or two how to potentially bridge V_HLOS with the momentum fluxes.