{"id":"c35cc4ae-0ad6-4bac-93b9-677b53b36476","arxiv_id":"1908.07789","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":3.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Acoustic-gravity waves in a stratified isothermal atmosphere can be represented as a combination of two harmonic oscillations at the acoustic and gravitational eigenfrequencies, with explicit real-variable solutions.","lead":"The paper derives real-valued solutions for acoustic-gravity waves in an isothermal atmosphere, showing that each wave is a sum of two oscillations at distinct eigenfrequencies. It offers a two-frequency interpretation that may explain why satellite-observed wave trains look like beats, though the comparison is only qualitative.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central mathematical derivation is sound, but Section 6's claim that satellite observations are explained by two eigenfrequencies is not quantitatively tested: a single-wavenumber, two-line fit to the observed wave trains is needed before the explanation can carry weight.","rationale":"The reader's weakest assumption and my independent read coincide: Section 6 is the least secure part of the central claim. In good faith, the theoretical derivation deserves credit: the two-real-oscillator representation follows from the initial-value problem, Eq. (10) is algebraically equivalent to the standard isothermal AGW dispersion relation, and the solution form is self-consistent (initial acceleration is zero, matching p'=ρ'=0 at t=0). The known typo in Eq. (12) (the second trigonometric factor in the ξz sum should be sin(k·r), not cos) is a mechanical error that does not change the structure of the argument. The observational explanation, however, is asserted rather than demonstrated. Since real AGW observations involve broadband spectra and atmospheric variability, the qualitative beat resemblance is insufficient to single out the two-frequency mechanism. A concrete spectral fit would settle this. Because this is the same concern as the reader's weakest assumption, I do not change the CONDITIONAL verdict.","tokens_in":6105,"tokens_out":10926,"duration_ms":96388,"concrete_test":"Analyze the Fedorenko et al. (2015) satellite wave-train events with kxH≈0.5 and kzH→0. Compute the power spectrum of Vx(t) and Vz(t), then fit two spectral lines with frequencies ω1(kx,kz) and ω2(kx,kz) from Eq. (10) for a common (kx,kz), with free amplitudes V0x, V0z and phases φx, φz. Compare the fitted line frequencies and the predicted phase/amplitude relations from Eq. (13) against the observed time series, and test the two-line model against a broadband alternative (e.g., by residual analysis or an F-test). If the two-line fit is rejected or the predicted beat period/phase relation is absent, the Section 6 explanation fails; if it survives, the observational claim is supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim bundles a mathematical result (Eqs. (10)-(13): linear isothermal AGW disturbances at fixed (kx,kz) evolve as the sum of two real eigen-oscillations) with an observational explanation (Section 6). The mathematical part is internally consistent and reproduces the standard AGW dispersion relation, so my concern is not with the algebra. The load-bearing weak spot is the observational interpretation. Section 6 selects observed wave trains with kxH≈0.5, kzH→0 and asserts that they 'may indicate' the two-frequency regime, citing morphological similarity to beats. But Eq. (13) predicts specific quantitative signatures - the beat period set by ω1-ω2, the amplitude modulation envelope, and the phase relationship between Vx and Vz - and none of these are compared with the satellite data. Real observed AGWs are broadband, dissipative, wind-affected, and temperature-dependent; a superposition of many modes can also produce beat-like time series. Without a quantitative two-spectral-line fit at the frequencies of Eq. (10) and a check of the predicted polarization/phase relations, the explanation is an unsupported qualitative claim. This is load-bearing because the paper's main result is explicitly framed as explaining observations that classical theory cannot.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript develops a two-frequency description of acoustic-gravity waves in an isothermal atmosphere. Starting from the linearized hydrodynamic equations, the authors derive coupled second-order equations for the horizontal and vertical displacement components, impose initial conditions of zero displacement and prescribed velocity, and obtain real-time solutions expressed as sums of two eigen-oscillations with the frequencies given by Eq. (10). The paper further presents velocity solutions in Eq. (13), examines special limiting cases, and claims that the two-frequency regime explains beat-like wave trains and phase properties observed by satellites (Section 6).","tokens_in":6349,"tokens_out":6274,"duration_ms":59846,"significance":"The paper's mathematical core is a self-contained, parameter-free derivation. It recovers the standard acoustic-gravity dispersion relation and offers real-variable, physically interpretable solutions that reduce to single-frequency oscillations in appropriate limits. The main weakness is the observational section, which asserts explanatory power without a quantitative comparison. If the observational analysis were made rigorous, the framework could be a useful interpretive tool for AGW data.","major_comments":[{"comment":"The second spatial term in the expression for ξz is printed as cos(k·r) twice; consistency with Eq. (8) and with the derivation of Eq. (13) requires the second term to be sin(k·r). This should be corrected.","section":"Section 4, Eq. (12)"},{"comment":"The claim that the two-frequency regime explains satellite observations is not supported by a quantitative comparison. The paper cites morphological similarity and the qualitative properties kxH≈0.5, kzH→0, but it does not compute the predicted beat period (ω1−ω2), the modulation envelope, or the phase relationships between Vx and Vz for the observed events, nor does it show that a single-wavenumber superposition is adequate. Without such a test, the statement that these observations 'can be explained' remains an unsubstantiated assertion. I recommend adding a quantitative fit or rephrasing the claim as a qualitative suggestion.","section":"Section 6, Eqs. (10) and (13)"},{"comment":"The reduction of Eqs. (13) to the compact beat form (14) is presented without derivation, and the notation is confusing because Vx,z denotes both the velocity component and its amplitude envelope. Please provide the algebra and introduce separate symbols for the envelope and phase.","section":"Section 5, Eq. (14)"}],"minor_comments":[{"comment":"The paper promises spectral characteristics of the perturbed velocity, but no spectra are computed or plotted; either add spectral analysis or remove the claim.","section":"Abstract and Section 7"},{"comment":"The cited observed properties come mainly from the authors' previous papers; the manuscript should summarize the measurement procedure and error bars so that the qualitative comparison can be evaluated.","section":"Section 6"},{"comment":"In the text, 'the dimensionless velocity components Vx/V0x and Vx/V0x' appears to contain a typo; the second quantity should likely be Vz/V0z.","section":"Section 5 and Figure captions"},{"comment":"Several references are abbreviated or contain inconsistencies (e.g., 'Prist' appears for what is likely Priest); the reference list should be standardized.","section":"References"},{"comment":"At x=0 the expression Vx,z = V0(t) cos(ω+t) cos(ω−t) is ambiguous because Vx(t) in Eq. (14) is nonnegative, whereas cos(ω−t) changes sign; please clarify whether V0(t) is meant to include an absolute value.","section":"Section 5, Eq. (14)"}],"recommendation":"major_revision","confidential_remarks":"The manuscript's observational support relies heavily on the authors' own prior papers; the editor may wish to ensure that the cited datasets are accessible and independently corroborated before publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know. The math is fine, and the real-variable two-frequency solution is a convenient repackaging of the standard AGW result, not a new wave mode. The observational section is the weak part: it claims to explain satellite wave trains, but never compares the predicted beat period, envelope, or polarization against the data.\n\nWhat the paper does well: it derives the coupled displacement equations for an isothermal atmosphere, solves them in real variables with initial conditions, and obtains explicit superpositions of the two eigenfrequencies. The special cases in Section 5 are genuinely nice—for example, at kz = 0 and kxH = 1/2, the solution reduces to the evanescent modes from the authors' earlier paper, with clean div V conditions. The beat behavior in Eq. (14) is also transparent. For teaching or for quick data interpretation, having a real-variable expression is handy.\n\nThe soft spots, in order of importance. First, Section 6 is not quantitative. The authors show that observed wave trains with kxH ≈ 0.5, kzH → 0 look beat-like, and that observed in-phase velocity components contradict Hines's π/2 phase shift. Then they say these 'can be explained' by the two-frequency approach. But they do not fit Eq. (13) to the data. There is no comparison of the observed beat period to ω1 − ω2, no check of the predicted amplitude envelope, no test of the Vx/Vz phase relation. Broadband, dissipative, wind-affected AGWs can produce beat-like time series from many modes, so the single-wavenumber, two-frequency story is a hypothesis, not an explanation. The paper should say that clearly, or better, add a fit.\n\nSecond, the novelty is overstated. Eq. (10) is the standard two-branch AGW dispersion relation; the two roots are the acoustic and gravity branches known since Hines (1960). The authors cite Hines and Yeh & Liu, so this is not a hidden issue, but the abstract and Section 7 frame the two-frequency representation as a new result. It's a real-variable repackaging, useful but not new physics.\n\nThird, there is a typo in Eq. (12): the second spatial term in ξz is written as cos(k·r) but should be sin(k·r) to match Eq. (8). Mechanical, but with a heavily typed formula it matters.\n\nWho this is for: AGW specialists who want a self-contained real-variable derivation and a clear picture of the two-frequency superposition. It deserves a serious referee: the algebra is sound and the tutorial value is real, but the observational claims need to be either scaled back or tested. I'd send it out, with a strong request to fix the typo, temper the novelty language, and make Section 6 explicitly a qualitative suggestion.","headline":"A clean real-variable derivation of the standard AGW two-branch solution, useful as a tutorial, but the satellite 'explanation' is a qualitative guess that needs a quantitative test.","tokens_in":6895,"tokens_out":3020,"would_cite":false,"duration_ms":29055,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Atmospheric acoustic-gravity waves beat at two frequencies at once","keywords":["acoustic-gravity waves","isothermal atmosphere","two coupled oscillators","two eigenfrequencies","beat phenomena","real-variable solutions","satellite observations","thermosphere"],"falsifier":"Take a satellite or lidar time series of an acoustic-gravity wave packet with $k_xH\\approx0.5$ and $k_zH\\to0$ and fit Eq. (13) allowing only amplitude, phase, and the two frequencies; if the fitted envelope period disagrees with $\\omega_{-}=\\omega_{0}[(1+2\\varepsilon)^{1/2}-(1-2\\varepsilon)^{1/2}]/2$ by more than the observational uncertainty, or if the predicted in-phase relation between vertical velocity and density fails, the single-mode two-frequency explanation is ruled out for that event.","tokens_in":1950,"feed_emoji":"🌊","tokens_out":4225,"duration_ms":134526,"temperature":0.7,"pith_summary":"The paper reworks linear acoustic-gravity wave theory in an isothermal atmosphere, treating a wave packet as a system of two coupled oscillators rather than a single monochromatic wave. It derives a dispersion relation with two eigenfrequencies and constructs the general solution for displacement as a sum of the two corresponding oscillations. The resulting real-valued formulas for perturbed velocity display beats when the two frequencies are close, which the paper argues explains observed wave trains in satellite data that standard single-frequency theory cannot account for. A sympathetic reader would care because this offers a parameter-free way to see beat-like acoustic-gravity signatures as intrinsic to the linear wave equation rather than as nonlinear or external effects.","feed_headline":"Atmospheric waves beat at two frequencies at once","feed_subtitle":"Real solutions produce beat patterns that match satellite wave trains.","key_machinery":"The machinery is a two-oscillator description of the coupled horizontal and vertical displacement equations. Writing displacements as $\\sin(\\omega t)$ times plane waves in $x,z$ reduces the dynamics to a $4\\times4$ algebraic system in the coefficients $a,b,c,d$; nontrivial solvability gives the dispersion relation (10) with the two eigenfrequencies. Each eigenfrequency yields a polarization relation linking vertical coefficients to horizontal ones, and the general solution is the sum of the two independent eigen-solutions. Initial conditions on perturbed velocity select the amplitudes, producing real, not complex, formulas for $V_x$ and $V_z$ whose interference generates beats.","core_discovery":"The central claim is that a small-amplitude acoustic-gravity disturbance in an isothermal stratified atmosphere is not a single mode but a superposition of two eigen-oscillations. The frequencies are the roots of Eq. (10), $\\omega_{1,2}^{2}=\\frac{1}{2}(\\omega_{01}^{2}+\\omega_{02}^{2}\\pm\\sqrt{D})$ with $D=(\\omega_{01}^{2}-\\omega_{02}^{2})^{2}+4(\\omega_{03}^{4}+\\omega_{04}^{4})$, where $\\omega_{01}^{2}=k_{x}^{2}c_{s}^{2}$, $\\omega_{02}^{2}=(k_{z}^{2}+1/(4H^{2}))c_{s}^{2}$, $\\omega_{03}^{2}=\\varepsilon k_{x}c_{s}^{2}/H$, $\\omega_{04}^{2}=k_{x}k_{z}c_{s}^{2}$, and $\\varepsilon=1/\\gamma-1/2$. Starting from the coupled second-order displacement equations, the paper shows that the general real solution is the sum of the two eigen-solutions, with coefficients fixed by initial conditions on velocity, giving closed-form expressions (13) for $V_x$ and $V_z$. For equal initial amplitudes and phases, these factor into $\\cos(\\omega_{+}t)\\cos(\\omega_{-}t)$ beats, where $\\omega_{\\pm}=\\omega_{0}[(1+2\\varepsilon)^{1/2}\\pm(1-2\\varepsilon)^{1/2}]/2$. The paper claims these beats match the morphology of satellite-observed acoustic-gravity wave trains, including typical $k_xH\\approx0.5$, $k_zH\\to0$ scaling, near-in-phase $V_x$ and $V_z$, and in-phase vertical velocity with density.","pith_inferences":["A natural extension beyond the paper is to ask whether the same two-frequency structure survives in slowly varying non-isothermal backgrounds; one testable prediction is that beat envelopes should persist along ray paths with locally defined scale height.","In the editor's reading, the paper's comparison with satellite data is qualitative; a quantitative test would be to fit Eq. (13) to individual wave-train time series and check whether the envelope period and the in-phase polarization relation hold within observational uncertainty.","The two-frequency decomposition may also apply to other stratified media, including stellar or solar atmospheric layers, where beat-like velocity signals could be misidentified as separate modes rather than as two eigenfrequencies of a single disturbance."],"forward_implications":["If the two-frequency picture is correct, beat-like intensity modulations in satellite observations of acoustic-gravity waves follow from linear theory alone, without invoking nonlinear interactions or external modulation.","The two-frequency solution contains the standard single-frequency regime as a limit: for $k_zH\\approx k_xH\\gg1$ the general solution reduces to a high-frequency acoustic oscillation, and at the evanescent point $k_z=0$, $k_x=1/(2H)$ it recovers previously studied evanescent acoustic-gravity modes.","At horizontal location $x=0$ the velocity components are exactly $\\cos(\\omega_{+}t)\\cos(\\omega_{-}t)$, a pure beat, while at $x=\\pi H/2$ the same solution gives two harmonic oscillations with slightly different frequencies, so the degree of beating is spatially dependent.","The real-variable form provides direct predictions for phase relations between $V_x$ and $V_z$ and between $V_z$ and density; the paper claims these can account for observed in-phase behaviour near $k_z\\to0$ where classical theory predicts a $\\pi/2$ phase shift.","Because the construction is parameter-free given wavenumbers, scale height, adiabatic index, and initial amplitudes, Eq. (13) can be compared directly with time-series data rather than through dispersion fits alone."],"supporting_citations":[{"why":"Establishes the classical single-frequency acoustic-gravity dispersion and polarization relations that the two-frequency result is compared against.","marker":"Hines, 1960"},{"why":"Provides the coupled second-order displacement equations used as the starting point for the two-frequency derivation.","marker":"Tolstoy, 1963"},{"why":"Supplies the linear-oscillation method by which the general solution is written as a sum of two independent eigen-solutions.","marker":"Landau and Lifshitz, 1969"},{"why":"Defines the evanescent acoustic-gravity modes whose single-frequency limit is recovered from the two-frequency solutions.","marker":"Cheremnykh et al., 2019"},{"why":"Supplies the satellite wave-train observations with $k_xH\\approx0.5$, $k_zH\\to0$ that the paper interprets as two-frequency beats.","marker":"Fedorenko et al., 2015"},{"why":"Documents the observed near-in-phase velocity components and vertical-velocity/density phasing that the paper claims the two-frequency approach explains.","marker":"Fedorenko, 2013"}],"fun_headline_variants":["Two-tone atmosphere: gravity waves beat like a chord","Acoustic-gravity waves are two-frequency beats","Atmospheric waves show dual-frequency beating","Beat patterns explain satellite wave trains","Two eigenfrequencies tune atmospheric wave beats"],"cache_read_input_tokens":9088,"weakest_assumption_plain":"The comparisons with satellite data assume that an observed wave train is well represented by a single plane-wave component at one fixed wavenumber pair $(k_x,k_z)$ that is a superposition of the two eigenfrequencies, even though real atmospheric acoustic-gravity waves are broadband and affected by winds, dissipation, and temperature variability.","fun_headline_variants_meta":{"raw":{"variants":["Two-tone atmosphere: gravity waves beat like a chord","Acoustic-gravity waves are two-frequency beats","Atmospheric waves show dual-frequency beating","Beat patterns explain satellite wave trains","Two eigenfrequencies tune atmospheric wave beats"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000205,"raw_usage":{"total_tokens":1405,"prompt_tokens":967,"completion_tokens":438,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":583,"completion_tokens_details":{"reasoning_tokens":380}},"tokens_in":583,"tokens_out":438,"duration_ms":51416,"temperature":1.0,"reasoning_tokens":380,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:55:42.685230+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take a satellite or lidar time series of an acoustic-gravity wave packet with $k_xH\\approx0.5$ and $k_zH\\to0$ and fit Eq. (13) allowing only amplitude, phase, and the two frequencies; if the fitted envelope period disagrees with $\\omega_{-}=\\omega_{0}[(1+2\\varepsilon)^{1/2}-(1-2\\varepsilon)^{1/2}]/2$ by more than the observational uncertainty, or if the predicted in-phase relation between vertical velocity and density fails, the single-mode two-frequency explanation is ruled out for that event.","supporting_citations":[{"cited_title":"Internal atmospheric gravity waves at ionospheric heights //Canadian J.Phys., 1960, V.38, P.1441-1481","cited_arxiv_id":null,"evidence_quote":"Establishes the classical single-frequency acoustic-gravity dispersion and polarization relations that the two-frequency result is compared against."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the coupled second-order displacement equations used as the starting point for the two-frequency derivation."},{"cited_title":"Evanescent acoustic -- gravity modes in the isothermal atmosphere: systematization and applications to the Earth and solar atmospheres","cited_arxiv_id":null,"evidence_quote":"Defines the evanescent acoustic-gravity modes whose single-frequency limit is recovered from the two-frequency solutions."},{"cited_title":"Geophys., 33, 101-108, 2015","cited_arxiv_id":null,"evidence_quote":"Supplies the satellite wave-train observations with $k_xH\\approx0.5$, $k_zH\\to0$ that the paper interprets as two-frequency beats."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Documents the observed near-in-phase velocity components and vertical-velocity/density phasing that the paper claims the two-frequency approach explains."}],"review_version":1}