{"id":"69597f5d-0af6-4a27-b72a-e38bae21275a","arxiv_id":"1908.03764","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Zonal spectral MSE variance budgets show that radiation builds variance at scales over 1,000 km while advection and surface fluxes damp it, consistently across convection-permitting models and observations.","lead":"This paper compares how radiation, surface heat exchange, and air motion create and destroy large-scale moisture patterns in the tropics, using the same spectral budget on computer simulations and on observed/reanalysis data. It finds that idealized storm-resolving models reproduce the sign and scale of the observed process balances, supporting their use as a simplified laboratory for tropical variability.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Observed advection variance rate is computed as a residual of a non-closing budget in ERA; the claimed damping by advection across wavelengths may be an artifact of budget residuals rather than a real atmospheric process.","rationale":"The reader's weakest assumption exactly identifies the residual-based advection in the observed budget, which is also acknowledged by the paper in Section 3.3. This is the most load-bearing concern because the advection term is one of the three main similarity claims (radiation increases variance at long wavelengths, advection damps across wavelengths, surface fluxes damp at intermediate scales). The other two terms (radiation and surface fluxes) are directly computed from flux products, though they still depend on the MSE field from ERA. The advection residual is the only term that could be entirely an artifact of budget non-closure in ERA. The paper's own statement that 'the fine variability of its variance rate may not be resolved' underscores this limitation. The proposed direct computation would settle the concern by providing an independent estimate of the advection variance rate. If the direct estimate confirms damping across all wavelengths, the central claim stands; if not, the analogy is weakened. The reader gave CONDITIONAL, which remains appropriate pending this test. We agree with the reader's assessment and find no additional load-bearing concern beyond the residual-advection issue.","tokens_in":13533,"tokens_out":5126,"duration_ms":54836,"concrete_test":"Recompute the advection variance rate in ERA5 by explicitly calculating the horizontal and vertical divergence of MSE flux using 3D wind, temperature, and humidity fields, rather than treating total advection as the residual of Eq. 7; then compare the sign and scale-selectivity with the model curves in Figure 3d. If the direct computation differs substantially, the observed advection damping is an artifact of the residual method, and the central claim would need revision.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim that MSE tendencies share similar signs and scale-selectivity in models and observations rests heavily on the advection term, which damps variance at all wavelengths in both ERA and CERES (Figure 3d). However, Section 3.3 states that total advection is calculated as a residual of Eq. 7, and that ERA does not close the MSE budget. Thus the observed advection variance rate absorbs all budget imbalances, analysis increments, and unrepresented processes. If the directly computed advection term differs in sign or spectral shape, the claimed analogy between idealized RCE and the real world would lose one of its three main pillars. The paper itself acknowledges that direct 3D advection is needed, but the abstract and conclusions present advection damping as a robust observed fact. This is load-bearing because the qualitative agreement in Figure 3d could be an artifact of the residual calculation rather than a real physical process.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper introduces a zonal spectral budget for column-integrated transient moist static energy (MSE), following the authors' earlier work, and applies it to four datasets: ERA5 reanalysis, CERES satellite radiation, a non-rotating long-channel RCE simulation (LC), and a near-global rotating aquaplanet simulation (NG). For each dataset, the authors compute the rate of MSE variance injection/removal by longwave radiation, shortwave radiation, surface enthalpy fluxes, and advection, as a function of zonal wavelength. They find that radiation (especially longwave) injects variance at wavelengths above about 1000 km, advection removes variance at all resolved wavelengths, and surface fluxes mostly remove variance between 1000 and 10,000 km. On this basis, they argue that idealized RCE simulations provide a valid analogy for the real tropical atmosphere, despite differences in the spectral shape (model spectra show a peak, observations do not). The spectral budget derivation is clean, and the authors check that the left-hand side of their budget equation is small relative to the dominant terms.","tokens_in":13651,"tokens_out":7098,"duration_ms":69612,"significance":"If the results hold, the paper provides a compact multi-scale framework for comparing process-level MSE variance budgets across models and observations, and supports the use of idealized RCE for tropical variability. The main strengths are: the diagnostic is parameter-free (no fitted constants); the spectral budget derivation is straightforward and checked (LHS = 0.1–5.1% of the longwave variance rate); the observed MSE spectrum is corroborated by two independent datasets (ERA and CERES); and the model sensitivity experiments (UNI-RAD, UNI-SEF) provide independent physical context. The central weakness is that the observed advection term is a residual of an unbalanced reanalysis budget, which the paper acknowledges but does not fully resolve.","major_comments":[{"comment":"The observed advection variance rate is not directly computed; it is derived as the residual of Equation 7 after subtracting the radiation and surface flux terms. Because the ERA reanalysis does not close the moist static energy budget, as the paper itself notes in the same paragraph, this residual includes analysis increments and unrepresented processes. The sign and spectral shape of the advection term in Figure 3d may therefore be contaminated, and the conclusion that advection damps MSE variance across all wavelengths in observations is not directly supported. Please compute the advection tendency explicitly from three-dimensional ERA fields (horizontal and vertical advection of column MSE) and compare it with the residual, or alternatively state in the abstract and conclusions that the observed advection damping is inferred from a residual of an unbalanced budget.","section":"§3.3, Figure 3d"},{"comment":"The paper's headline claim that 'the results affirm the validity of the RCE simulations as an analogy to the real world' rests in part on the observed advection term (Figure 3d), which is subject to the residual uncertainty described above. Since the paper itself acknowledges in §3.3 that 'explicitly calculating the horizontal and vertical components of MSE advection from three-dimensional data will be needed to clarify its scale-selectivity in observations,' the abstract and conclusions should be tempered to reflect this unresolved uncertainty, or the direct calculation should be performed.","section":"Abstract and Section 4"}],"minor_comments":[{"comment":"The sentence 'The new reanalysis dataset has a better hydrological cycle and sea surface temperatures in the Tropics and is calibrated for climate applications' is vague; please specify what is meant by 'better' and 'calibrated,' or provide a reference for the improvements.","section":"§2.1"},{"comment":"The NG simulation is described as 'convection-permitting,' but its saved output resolution is 156.25 km. Please clarify whether the model is integrated at convection-permitting resolution with output saved at coarser resolution, or whether the model itself uses this coarse grid. This distinction affects the interpretation of the short-wavelength part of the NG spectra.","section":"§2.3"},{"comment":"In the caption, 'The UNI-SEF rates of variance injection (dashed lines) have been divided by a factor of 5 because the denominator of equation 7... is smaller for non-aggregated simulations' is unclear: dividing by an arbitrary factor makes the curves comparable visually but obscures the actual amplitude difference. Please explain the scaling more fully and consider plotting the true rates on a separate axis or with a different normalization.","section":"Figure 3 caption"},{"comment":"The time-averaging operation used to go from Equation (6) to Equation (7) is not explicitly defined. Please define the average, for example as \\overline{(\\cdot)} = (1/t_H) \\int_{t_0}^{t_0+t_H} (\\cdot) \\, dt, so that the notation is unambiguous.","section":"Equation (7)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is well-suited to GRL in length and scope, and the code/data availability is commendable. The main risk is the observed advection residual issue, which the authors themselves flag but do not resolve; I recommend requiring either a direct calculation of advection from ERA winds or a clearly qualified statement in the abstract and conclusions. The paper should not be rejected, as the framework and model-observation comparison are valuable and the central derivation is sound."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is the first scale-by-scale comparison of MSE variance budget terms between idealized convection-permitting RCE and observations/reanalysis, and it is a useful diagnostic contribution rather than a breakthrough. The spectral budget derivation is clean, and the sensitivity experiments (UNI-RAD, UNI-SEF) provide clean controls. The paper does what it sets out to do, and the authors are honest about the main caveats.\n\nWhat it does best: it shows that radiation consistently injects MSE variance at wavelengths above ~1000 km and surface enthalpy fluxes damp variance between ~1000 and 10000 km in both models and observations. The left-hand side of Eq. 7 is checked to be under a few percent of the longwave variance rate, and the code and data are public. The comparison with CERES and ERA using the same zonal spectral diagnostics is a genuine new result.\n\nThe soft spot is real and load-bearing: advection in the observations is computed as a residual of Eq. 7, and ERA does not close the MSE budget. So the claim that advection damps variance across wavelengths in the abstract and Figure 3d may be contaminated by budget residuals, analysis increments, and unrepresented processes. The stress-test note is right about this. The authors do flag it in Section 3.3 and say explicit three-dimensional advection is needed, but the abstract and conclusions state it as a robust finding, which overstates the evidence. Also absent: any confidence intervals or uncertainty quantification on the spectral rates, and the analysis is zonal-only, which excludes meridional transport and makes the LC comparison harder to interpret.\n\nThose caveats do not sink the central qualitative result. The radiation and surface-flux results are computed from direct observations, and the model/observation agreement there is credible. The paper is a useful methodological template for global storm-resolving model intercomparison, and it should go to review. My recommendation: send it to peer review, ask for a direct computation of 3D advection in a reanalysis that closes the budget, or at least a quantitative estimate of the residual contamination, plus uncertainty bounds on the spectral rates. Without the advection pillar the conclusion loses one of its three supports, but the other two are solid enough to make this a worthwhile contribution for tropical convection specialists and modeling centers.","headline":"Useful first scale-resolved comparison of MSE variance budget terms between idealized RCE and observations; the observed advection term rests on a non-closing budget residual, so that pillar is softer than the abstract implies.","tokens_in":14211,"tokens_out":1908,"would_cite":true,"duration_ms":20110,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Idealized cloud-resolving simulations reproduce observed tropical moisture variance across scales","keywords":["convective self-aggregation","moist static energy","spectral budget","radiative-convective equilibrium","tropical variability","convection-permitting models","scale selectivity","reanalysis"],"falsifier":"Explicitly compute the three-dimensional horizontal and vertical advection of moist static energy from ERA5 fields and evaluate its spectral variance rates; if the explicit advection rates are not negative at all wavelengths, the residual-based claim of uniform advection damping would be falsified. A second check would compare the ERA5 spectral budget against an independent reanalysis or a fully instrument-based product.","tokens_in":1765,"feed_emoji":"🌧️","tokens_out":4000,"duration_ms":69034,"temperature":0.7,"pith_summary":"This paper argues that idealized, convection-permitting simulations of radiative-convective equilibrium can serve as valid analogies for real tropical moisture variability. By applying a zonal spectral budget to column moist static energy, the authors compare how radiation, advection, and surface enthalpy fluxes create or destroy variance across wavelengths in idealized simulations, a reanalysis, and satellite observations. They find the same scale-selective pattern in all datasets: radiation increases variance at wavelengths above 1,000 km, advection damps variance across wavelengths, and surface fluxes mostly reduce variance between 1,000 and 10,000 km. If true, this supports the use of idealized self-aggregation experiments to understand real tropical rainfall organization and to evaluate global cloud-resolving models.","feed_headline":"Idealized cloud runs mirror observed tropical moisture variability","feed_subtitle":"In reanalysis and satellite data, radiation builds moisture contrast at long scales while advection and surface fluxes smooth it.","key_machinery":"The central object is the zonal power spectrum $\\phi_H$ of transient column moist static energy, together with the spectral variance budget of equation 7, which divides each tendency's co-spectrum $\\Re(\\hat H'^*\\hat{\\dot H}'_i)$ by the time-mean spectrum $\\phi_H$ to obtain per-wavelength variance rates with units of inverse time. Because the budget normalizes by the spectrum, it does not explicitly depend on the time-mean zonal structure, which lets the authors make direct analogies between observations and zonally symmetric radiative-convective equilibrium.","core_discovery":"The paper establishes that when column-integrated moist static energy is decomposed into zonal scales, idealized radiative-convective equilibrium simulations run with convection-permitting models produce the same scale-selective moist static energy variance tendencies as observed in the tropics. Longwave radiation injects variance at all wavelengths, shortwave radiation injects variance at long wavelengths, advection removes variance across scales, and surface enthalpy fluxes mostly remove variance between roughly 1,000 and 10,000 km. The observed spectra and the control simulations agree at long wavelengths, while simulations in which radiation is homogenized lose variance at wavelengths above 1,000 km and disagree with observations by more than an order of magnitude. The paper therefore affirms the radiative-convective equilibrium analogy and identifies the long-channel configuration as an inexpensive, reduced-size framework for studying the processes that maintain convective aggregation.","pith_inferences":["The paper's advection result in observations rests on a residual calculation, so an explicit three-dimensional diagnosis of horizontal and vertical advection could revise the claim that advection damps variance at all wavelengths.","If surface-enthalpy-flux damping is driven by near-surface enthalpy disequilibrium, then ocean coupling strength may control both the degree of aggregation and the presence or absence of a spectral peak in observed moist static energy.","Applying the same spectral budget to lower-tropospheric water vapor, or to output from global cloud-resolving models, would test whether the moist static energy result extends to variables that are more directly tied to precipitation extremes.","The absence of a distinct spectral peak in observations, despite positive radiation variance injection, suggests that external forcing and lateral mixing may mask self-aggregation in spectral space; this could be tested by filtering observed fields into intraseasonal versus higher-frequency components."],"forward_implications":["If the radiative-convective equilibrium analogy holds, convective self-aggregation likely plays a role in generating observed tropical moisture variability, not just in idealized models.","Positive longwave variance injection at scales above about 1,000 km supports the idea that radiative-convective feedbacks are key to producing realistic moisture variability from homogeneous boundary conditions.","Stronger surface-enthalpy-flux damping in the reanalysis, relative to fixed-sea-surface-temperature models, is consistent with interactive ocean coupling or meridional sea surface temperature gradients damping self-aggregation patterns.","The long-channel configuration could serve as a relatively inexpensive framework for studying convective-aggregation processes across different climates.","The spectral framework generalizes to three-dimensional tracer variance budgets and to limited-area or localized analyses through discrete cosine or wavelet transforms."],"supporting_citations":[{"why":"Provides the long-channel control simulation and the long-channel geometry used as an idealized, reduced-size configuration.","marker":"Wing & Cronin, 2016"},{"why":"Provides the near-global rotating simulation whose control and sensitivity runs are compared against observations.","marker":"Khairoutdinov & Emanuel, 2018"},{"why":"Supplies the spectral variance budget derivation and the uniform-radiation and uniform-surface-flux sensitivity experiments.","marker":"Beucler & Cronin, 2018"},{"why":"Establishes the moist static energy variance budget approach for diagnosing convective self-aggregation in convection-permitting models.","marker":"Wing & Emanuel, 2014"},{"why":"Foundational moist static energy variance framework that this paper generalizes to a per-wavelength spectral budget.","marker":"Bretherton et al., 2005"},{"why":"Earlier comparison of MSE advection and radiation in realistic convective-scale simulations, extended here to larger domains and observations.","marker":"Holloway, 2017"},{"why":"Documents the ERA5 reanalysis dataset that provides one of the two observational estimates of the MSE spectrum and budget.","marker":"Hersbach & H., 2016"},{"why":"Introduces the CERES satellite dataset that provides the independent observational radiative flux data used in the spectral budget.","marker":"Wielicki et al., 1996"}],"fun_headline_variants":["Model clouds echo real tropical moisture patterns","Idealized runs mirror observed moisture variance","Simulated storms match real moisture swings","Cloud models mimic tropics' moisture variability","RCE simulations reproduce observed MSE spectra"],"cache_read_input_tokens":16384,"weakest_assumption_plain":"In the observational data, the advection term is not measured directly but inferred as whatever remains after accounting for radiation and surface fluxes, and the reanalysis does not fully close the moist static energy budget, so the conclusion that advection damps variance at all wavelengths in observations could be contaminated by these budget residuals.","fun_headline_variants_meta":{"raw":{"variants":["Model clouds echo real tropical moisture patterns","Idealized runs mirror observed moisture variance","Simulated storms match real moisture swings","Cloud models mimic tropics' moisture variability","RCE simulations reproduce observed MSE spectra"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000858,"raw_usage":{"total_tokens":3699,"prompt_tokens":896,"completion_tokens":2803,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":512,"completion_tokens_details":{"reasoning_tokens":2741}},"tokens_in":512,"tokens_out":2803,"duration_ms":23136,"temperature":1.0,"reasoning_tokens":2741,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:02:06.562824+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Explicitly compute the three-dimensional horizontal and vertical advection of moist static energy from ERA5 fields and evaluate its spectral variance rates; if the explicit advection rates are not negative at all wavelengths, the residual-based claim of uniform advection damping would be falsified. A second check would compare the ERA5 spectral budget against an independent reanalysis or a fully instrument-based product.","supporting_citations":[{"cited_title":"\\ Emanuel, K","cited_arxiv_id":null,"evidence_quote":"Provides the near-global rotating simulation whose control and sensitivity runs are compared against observations."},{"cited_title":"\\ Emanuel, K a","cited_arxiv_id":null,"evidence_quote":"Establishes the moist static energy variance budget approach for diagnosing convective self-aggregation in convection-permitting models."},{"cited_title":", Blossey, P N","cited_arxiv_id":null,"evidence_quote":"Foundational moist static energy variance framework that this paper generalizes to a per-wavelength spectral budget."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Documents the ERA5 reanalysis dataset that provides one of the two observational estimates of the MSE spectrum and budget."},{"cited_title":", Barkstrom, B R","cited_arxiv_id":null,"evidence_quote":"Introduces the CERES satellite dataset that provides the independent observational radiative flux data used in the spectral budget."}],"review_version":1}