{"id":"d04227f6-2934-4cd9-817c-aed52d4128ae","arxiv_id":"1908.03534","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"For the 2014 March 29 Moreton wave, the authors measure a 4 km/s downward chromospheric velocity and derive coronal shock Mach numbers of roughly 1.02 to 1.28.","lead":"This paper measures the plasma motion in the Sun's chromosphere during a Moreton wave from a solar flare, finding downward speeds up to 4 km/s and using that to estimate the strength of the coronal shock. The study combines ground-based H-alpha images with space-based EUV and X-ray data to derive Mach numbers from two independent methods.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Equation (4) is not derivable from Eqs. (2)-(3) as printed; it silently assumes cch/cco = 1/sqrt(a), and using the paper's stated cch = 10 km/s and cco = 185 km/s changes vi from 22 to ~13 km/s, dropping the H-alpha-based Mach numbers below the DEM-based range.","rationale":"The paper is a careful multi-instrument study of a well-observed Moreton wave. Its strongest quantitative payload is the weak-shock Mach numbers, and the authors present two routes: the H-alpha Doppler velocity (Sections 5.1-5.2) and the DEM compression ratio (Section 5.3). The propagation speeds (640-859 km/s) and the 4 km/s downward Doppler detection are plausible, and the DEM analysis is an independent diagnostic with real support. I considered whether the DEM-based headline range is the weaker link: it is an upper envelope (maximum compression ratios selected at times of strongest downward motion) and is restricted to 6.1 <= log T/K <= 6.4, where heating can masquerade as compression. That is a legitimate secondary concern. However, the paper itself flags the temperature restriction as a limitation, and a modest overestimate of X would still leave a weak shock (Mach roughly 1.1-1.3). The least secure premise is the analytic bridge Eq. (4), because it is the only place where an arbitrary density ratio and an implicit sound-speed relation convert the observed chromospheric velocity into a coronal shock property, and because the printed formula is not the direct consequence of Eqs. (2)-(3) that the text claims. Correcting this changes the H-alpha-based Mach numbers and removes the claimed agreement with the DEM route. This is a falsifiable, fixable problem rather than a fatal flaw, so the conditional verdict stands.","tokens_in":30746,"tokens_out":14132,"duration_ms":143153,"concrete_test":"Independently derive Eq. (4) from Eqs. (2) and (3) without assuming cch/cco = 1/sqrt(a), then recompute Table 2's vi column and Table 3's Mach numbers using the paper's stated cco = 185 km/s, cch = 10 km/s, and a = 100, with a = 25 and a = 400 as a sensitivity check. If the corrected H-alpha-based Mach numbers fall below the DEM-based range of 1.06-1.28 at the corresponding times, the claimed consistency between the two diagnostics is not supported and Sections 5.1-5.2 need revision or the consistency claim must be softened.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim rests on two diagnostics said to agree: the H-alpha Doppler velocity and the DEM compression ratio. The load-bearing step is the H-alpha-to-coronal-shock conversion in Sections 5.1-5.2. Equation (4), vi = ((1+sqrt(a))/2) vt, is stated as following from Eqs. (2) and (3), but direct division of Eq. (3) by Eq. (2) gives vi = (1 + a*cch/cco) vt / 2, not the printed form. The printed form holds only under the unstated condition cch/cco = 1/sqrt(a). The authors instead treat cch/cco as negligible in Eq. (5) and in the surrounding discussion, which is a different and inconsistent assumption. With the paper's own numbers (cch = 10 km/s for the chromosphere at 10^4 K, cco = 185 km/s for the corona at 2.5 MK, and a = 100), a*cch/cco is about 5.4, so vi is approximately 3.2*vt = 13 km/s for vt = 4 km/s, not 5.5*vt = 22 km/s. Recomputing Table 2 and Table 3 with this corrected relation lowers the H-alpha-based Mach numbers from up to 1.12 to roughly 1.07. The claimed quantitative agreement with the DEM-derived range (1.06-1.28) is therefore not established at the stated precision. The DEM numbers themselves are not directly invalidated by this algebra, but the two-diagnostic consistency argument at the center of the paper is weakened.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports a multiwavelength study of the Moreton wave associated with the X1.0 flare on 2014 March 29, using Flare Monitoring Telescope Hα wing observations, SDO/AIA EUV data, and Hinode/XRT. The authors measure the Moreton wave propagation speed (640–859 km/s) and, from Hα wing Doppler signals, infer a maximum downward chromospheric velocity of about 4 km/s at the wave front. They then use a weak-shock transmission model to convert this chromospheric velocity amplitude into an incident-shock Mach number in the corona, and independently derive Alfvén and fast-mode Mach numbers from AIA-based differential emission measure compression ratios. The two diagnostics yield overlapping Mach number ranges (Hα: up to ~1.12; DEM: 1.05–1.27), which the authors interpret as consistent evidence for a weak fast-mode shock.","tokens_in":31070,"tokens_out":10198,"duration_ms":92816,"significance":"If the results stand, the paper provides a valuable multi-diagnostic case study of a Moreton wave: it quantifies the chromospheric velocity amplitude with Hα wing imaging, derives shock Mach numbers from two independent datasets, and connects the chromospheric and coronal signatures of a global disturbance. The use of the FMT wing data to extract Doppler velocities and the combination with DEM analysis are genuine strengths. The paper also carefully discusses limitations such as line-of-sight projection and the restricted DEM temperature range. However, the derivation of the central transmission relation contains an algebraic error that must be corrected before the quantitative results can be accepted.","major_comments":[{"comment":"Equation (4) is not correctly derived from Eqs. (2) and (3). Dividing Eq. (3) by Eq. (2) gives vi = [1 + a(cch/cco)] vt / 2, not vi = [(1 + sqrt(a))/2] vt. The printed form holds only under the unstated assumption cch/cco = 1/sqrt(a). Using the paper's own values (cch = 10 km/s, cco = 185 km/s, a = 100), the coefficient is about 3.2, not 5.5, so the incident velocities vi in Table 2 and the velocity-transmittance discussion in §6.3 and Figure 14 are overestimated by a factor of about 1.7. Please correct Eq. (4), Table 2, Figure 14, and all related text, or explicitly state and justify the additional assumption on the sound-speed ratio.","section":"§5.1, Eq. (4)"},{"comment":"The treatment of the sound-speed ratio cch/cco is internally inconsistent. In §5.2 the term cch/cco is neglected as small in Eq. (5), but Eq. (4) as printed effectively requires cch/cco = 1/sqrt(a) ≈ 0.1, which is not negligible. With the stated chromospheric and coronal sound speeds (10 and 185 km/s), cch/cco ≈ 0.054, so neglecting it in Eq. (5) is reasonable, but the same ratio cannot also be set to 0.1 in Eq. (4). The manuscript should reconcile these two approximations and state which values of cch and cco are used in each relation.","section":"§5.1 vs. §5.2"},{"comment":"The Hα-based Mach numbers in Table 3 are computed from Eq. (5) using the observed vt, not from Eq. (4). Therefore the algebraic error in Eq. (4) does not directly invalidate the Hα Mach numbers. However, the strong claim of consistency with the DEM-derived Mach numbers rests entirely on the validity of Eq. (5), which assumes a one-dimensional, purely hydrodynamic weak shock with the transition region approximated as a contact discontinuity. The paper should explicitly state the regime of validity of this approximation and discuss how neglecting magnetic fields and the oblique nature of the shock could affect the derived Mach numbers, in order to support the cross-diagnostic comparison.","section":"§5.2, Eq. (5) and Table 3"}],"minor_comments":[{"comment":"There are several typographical errors, e.g., 'does exit' in the introduction should be 'does exist', and 'Filter trasmission' in Figure 7 should be 'Filter transmission'.","section":"General"},{"comment":"Equation (8) assumes that the line-of-sight depth l is the same ahead of and behind the shock, and that the proton number density is constant along the line of sight. The paper should explicitly note that this is an approximation and comment on how inhomogeneities or temperature-dependent filling factors might bias the compression ratio.","section":"§5.3, Eq. (8)"},{"comment":"The Doppler velocities are line-of-sight velocities, but the event is at N11, W32, so the radial (vertical) component is foreshortened. The paper does not discuss this projection effect or its impact on the maximum downward velocity of 4 km/s; a brief quantitative statement would be useful.","section":"§4.2 and Figure 8"},{"comment":"The Mach numbers derived from the DEM compression ratio are presented without error bars or uncertainties. Since the compression ratio X enters through Eq. (7) nonlinearly, it would clarify the robustness of the quoted ranges (1.06–1.28 and 1.05–1.27) if some estimate of the uncertainty from the DEM inversion and the choice of β were provided.","section":"§5.3 and Figure 13"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a solid observational study of a well-observed event, and the main scientific conclusions are likely of interest to the solar physics community. The principal issue is the algebraic error in Eq. (4) and the accompanying inconsistency in the treatment of cch/cco; these are fixable within the scope of the manuscript. The recommendation is major revision rather than reject because the central two-diagnostic comparison is not directly destroyed by the error, but the quantitative results for vi and the transmission discussion must be corrected."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: this is a careful single-event study with a new Doppler measurement for the 2014 March 29 Moreton wave, and the qualitative picture is convincing. But the quantitative link from the H-alpha velocity to the coronal Mach number has a derivation problem that undercuts the main consistency claim.\n\nWhat is genuinely new: the FMT two-wing H-alpha data give clean chromospheric Doppler signals—propagation speeds of 640–859 km/s along four paths and a maximum downward velocity of about 4 km/s. That is a useful addition, especially because the concurrent Long et al. study used only H-alpha core observations. The DEM analysis with AIA provides independent compression ratios and hence Alfven/fast-mode Mach numbers around 1.05–1.28. The multiwavelength time-distance comparison across chromosphere, transition region, and corona is well executed, and the co-spatiality argument is credible. The authors cite the relevant prior H-alpha Doppler work (Warmuth, Balasubramaniam, etc.) and do not overstate novelty.\n\nWhere it gets soft: Eq. (4) is not derivable from Eqs. (2)–(3) as printed. Dividing Eq. (3) by Eq. (2) gives vi = (1 + a*cch/cco)*vt/2, not (1+sqrt(a))*vt/2. The printed form silently requires cch/cco = 1/sqrt(a), which is inconsistent with neglecting cch/cco later in Eq. (5). Using the paper's own numbers (cch=10 km/s, cco=185 km/s, a=100) yields vi ~13 km/s rather than 22 km/s. The H-alpha-based Mach numbers consequently drop from ~1.12 to ~1.07. That still overlaps the DEM range at its lower end, but the two-diagnostic agreement is not as clean as the paper claims. This is fixable: show the derivation, state the assumption, and recompute the tables.\n\nThe other soft spots are minor: Table 2 has no uncertainty estimates for the Doppler velocities; the choice of maximum compression peaks for Figure 13 is a selection that could bias Mach numbers upward; and the 1D hydrodynamic weak-shock model ignores line-of-sight projection and magnetic fields—though the authors acknowledge both limitations. The restricted DEM temperature range is also acknowledged in the text.\n\nOverall, the observations and the qualitative scenario (weak fast-mode shock, chromospheric downward push, relaxation) hold up. The paper deserves serious peer review and a revision that fixes the algebra and adds uncertainties.\n\nRecommendation: send it to review, but with a referee who will check the transmission equations carefully.","headline":"A solid single-event Moreton wave study with a genuinely new H-alpha Doppler measurement, but the transmission-equation algebra in Eq. (4) has a real slip that weakens the claimed agreement with the DEM-based Mach numbers.","tokens_in":31682,"tokens_out":3905,"would_cite":false,"duration_ms":42663,"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":"The 2014 March 29 Moreton wave pushed chromospheric plasma downward at up to 4 km/s, implying a weak coronal shock with Alfvén Mach number 1.06–1.28 and fast-mode Mach number 1.05–1.27.","keywords":["solar flares","Moreton waves","H-alpha Doppler","MHD shock waves","chromosphere","coronal waves","differential emission measure","Mach number"],"falsifier":"Observe the same wave front with full H-$\\alpha$ line-profile spectroscopy: if the red-wing excess is a true Doppler shift, the line core shifts redward by a corresponding few km/s; if the profile only deepens without a core shift, the transmitted-velocity interpretation, and therefore the Mach numbers $M_A \\approx 1.06$–$1.28$ and $M_f \\approx 1.05$–$1.27$, would need revision. Independently, type II radio burst drift rates for this flare give a separate shock-speed estimate to compare with the derived Mach numbers.","tokens_in":30532,"feed_emoji":"☀️","tokens_out":8498,"duration_ms":81748,"temperature":0.7,"pith_summary":"This paper reconstructs the dynamics of the Moreton wave from the 2014 March 29 X1.0 flare using ground-based H-$\\alpha$ wing images. It finds that the wave front swept across the chromosphere at 640–859 km/s while the plasma at the front was pushed downward at up to 4 km/s, followed by an upward relaxation. Treating the transition region as a contact discontinuity and applying the weak-shock approximation, the authors convert this small chromospheric push into an incident coronal shock that is barely super-Alfvénic and barely super-fast-mode. They independently obtain Mach numbers of $M_A \\approx 1.06$–$1.28$ and $M_f \\approx 1.05$–$1.27$ from emission-measure compression ratios. A sympathetic reader would care because this connects a faint, difficult-to-observe chromospheric wave to the basic strength of coronal shocks and supports the standard picture of Moreton waves as footprints of coronal fast-mode shocks.","feed_headline":"Moreton wave exposes a weak coronal shock, Mach 1.05-1.27","feed_subtitle":"A 4 km/s downward push in H-alpha wings places the 2014 X-class shock just above the fast-mode speed.","key_machinery":"The load-bearing identities are the Doppler signal $DS = (I_r - I_b)/(I_r + I_b)$ from the H-$\\alpha$ $\\pm 0.8$ Å wing images and two shock-transmission formulas: $v_i = (1+\\sqrt{a})v_t/2$ for the incident coronal velocity amplitude in terms of the observed transmitted chromospheric amplitude, and the weak-shock relation $v_t = -\\frac{4}{\\gamma+1}(M_i^2-1)\\frac{c_{\\mathrm{ch}}}{1+c_{\\mathrm{ch}}/c_{\\mathrm{co}}}$ that converts that amplitude into the incident shock Mach number. The same machinery includes the oblique MHD shock jump relation, solved for the Alfvén Mach number as a function of the DEM-derived compression ratio and field inclination. Together these turn two observable quantities—a small red-wing excess at the wave front and a coronal density jump—into a diagnostic of coronal shock strength.","core_discovery":"The central claim is that the Moreton wave of 2014 March 29 was the chromospheric footprint of a weak, almost perpendicular fast-mode shock propagating in the corona. In H-$\\alpha$ wing data the leading edge appears as a red-wing absorption and blue-wing brightening, which the paper calibrates through a synthetic Doppler signal into a downward line-of-sight velocity reaching $-4$ km/s, with a delayed upward swing. Using the momentum- and energy-flux matching at the corona–chromosphere boundary and the weak-shock relation, a 1–4 km/s chromospheric amplitude translates into coronal incident velocities up to about 22 km/s and Mach numbers near 1.1; AIA differential-emission-measure maps give compression ratios that yield the same nearly sonic range ($M_A \\approx 1.06$–$1.28$, $M_f \\approx 1.05$–$1.27$). The paper therefore claims that a single large-scale MHD disturbance can explain the quasi-simultaneous response of the corona, transition region, and chromosphere, and that the weak shock approximation is sufficient to diagnose coronal shock strength from ground-based H-$\\alpha$ Doppler measurements.","pith_inferences":["A testable extension: compare these Mach numbers with type II radio-burst drift speeds for the same event; agreement would confirm the weak-shock mapping, while systematic disagreement would point to projection or non-Doppler contamination in the H-alpha signal.","The same Doppler-signal method applied to other Moreton waves observed with H-alpha wing telescopes could produce a statistical sample of coronal shock strengths and their dependence on flare energy.","Because the H-alpha wings at $\\pm 0.8$ Å sample a fixed wavelength, line-profile broadening or intensity changes can masquerade as Doppler shift; full-profile spectroscopy at the wave front would separate true mass motion from thermodynamic changes.","If the direction-dependence of the wave speed reflects the Alfvén-speed map, Moreton-wave kinematics themselves could be inverted to constrain coronal magnetic-field topology in the low corona."],"forward_implications":["If the central claim is right, Moreton-wave Doppler amplitudes can be used as a ground-based measure of coronal shock Mach number without EUV data.","The observed propagation speeds (640–859 km/s in H-alpha, faster along some coronal paths) imply that the shock strength and direction are shaped by the local Alfvén speed, favoring propagation into weak-field regions.","The near-simultaneous response in 304 Å, H-alpha, 211 Å, and X-ray data strengthens the standard picture where the chromospheric Moreton wave is the footprint of a globally expanding coronal fast-mode shock.","The small velocity transmittance (roughly one fifth to one quarter) quantifies why only weak chromospheric Doppler signatures are seen even when the coronal disturbance is large.","For this event the shock is weak ($M_A, M_f \\lesssim 1.3$), so the weak-shock approximation is internally consistent."],"supporting_citations":[{"why":"Established H-alpha wing observations as the standard way to detect Moreton waves and set this paper's observational method.","marker":"Moreton & Ramsey (1960)"},{"why":"Supplies the theoretical model this paper confirms: a coronal fast-mode shock whose footprint is the chromospheric Moreton wave.","marker":"Uchida (1968)"},{"why":"Provides the incident-to-transmitted velocity amplitude relation used to convert the H-alpha velocity into a coronal incident velocity.","marker":"Takahashi et al. (2015)"},{"why":"Gives the Rankine-Hugoniot and weak-shock basis for the Mach-number formula derived in the appendix.","marker":"Landau & Lifshitz (1987)"},{"why":"Supplies the oblique MHD shock jump relation that underlies the Alfvén and fast-mode Mach number calculation.","marker":"Priest (2000)"},{"why":"Presents the same algebraic solution of the oblique shock relation used here for the horizontal and perpendicular cases.","marker":"Vršnak et al. (2002b)"},{"why":"Provides the DEM inversion method used to derive compression ratios from AIA observations.","marker":"Cheung et al. (2015)"},{"why":"Supplies the XRT filter-ratio method used to fix the pre-shock coronal temperature and sound speed.","marker":"Narukage et al. (2011)"}],"fun_headline_variants":["Moreton wave shock nearly sonic, Mach 1.05-1.27","H-alpha wing data reveal Moreton wave's weak coronal shock","Downward 4 km/s at front: Moreton wave's chromospheric shock","Barely supersonic: Moreton wave shock from ground H-alpha","X-class flare's Moreton wave: shock just above fast-mode speed"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The assumption that carries the result is that the H-$\\alpha$ Doppler velocity at the wave front equals the transmitted velocity amplitude $v_t$ of a one-dimensional hydrodynamic shock crossing the transition region, with a density ratio $a=\\rho_{\\mathrm{ch}}/\\rho_{\\mathrm{co}} \\approx 100$ and a negligible sound-speed ratio; if projection or brightness effects contaminate the H-$\\alpha$ signal, the derived Mach numbers shift.","fun_headline_variants_meta":{"raw":{"variants":["Moreton wave shock nearly sonic, Mach 1.05-1.27","H-alpha wing data reveal Moreton wave's weak coronal shock","Downward 4 km/s at front: Moreton wave's chromospheric shock","Barely supersonic: Moreton wave shock from ground H-alpha","X-class flare's Moreton wave: shock just above fast-mode speed"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000673,"raw_usage":{"total_tokens":3147,"prompt_tokens":1111,"completion_tokens":2036,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":727,"completion_tokens_details":{"reasoning_tokens":1936}},"tokens_in":727,"tokens_out":2036,"duration_ms":18923,"temperature":1.0,"reasoning_tokens":1936,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:10:42.017510+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Observe the same wave front with full H-$\\alpha$ line-profile spectroscopy: if the red-wing excess is a true Doppler shift, the line core shifts redward by a corresponding few km/s; if the profile only deepens without a core shift, the transmitted-velocity interpretation, and therefore the Mach numbers $M_A \\approx 1.06$–$1.28$ and $M_f \\approx 1.05$–$1.27$, would need revision. Independently, type II radio burst drift rates for this flare give a separate shock-speed estimate to compare with the derived Mach numbers.","supporting_citations":[{"cited_title":"1968, , 4, 30","cited_arxiv_id":null,"evidence_quote":"Supplies the theoretical model this paper confirms: a coronal fast-mode shock whose footprint is the chromospheric Moreton wave."},{"cited_title":"2015, , 801","cited_arxiv_id":null,"evidence_quote":"Provides the incident-to-transmitted velocity amplitude relation used to convert the H-alpha velocity into a coronal incident velocity."},{"cited_title":"1987, Fluid Mechanics 2nd ed., (Pergamon Press)","cited_arxiv_id":null,"evidence_quote":"Gives the Rankine-Hugoniot and weak-shock basis for the Mach-number formula derived in the appendix."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the oblique MHD shock jump relation that underlies the Alfvén and fast-mode Mach number calculation."},{"cited_title":"2015, , 807, 143","cited_arxiv_id":null,"evidence_quote":"Provides the DEM inversion method used to derive compression ratios from AIA observations."},{"cited_title":"2011, , 269, 169","cited_arxiv_id":null,"evidence_quote":"Supplies the XRT filter-ratio method used to fix the pre-shock coronal temperature and sound speed."}],"review_version":1}