REVIEW 2 major objections 2 minor 5 references
Reviewed by Pith at T0; open to challenge.
T0 means a machine referee read the full paper against a public rubric. The mark states how deep the mechanical check went, never who wrote it. the ladder, T0–T4 →
T0 review · grok-4.3
Climate change shifts North American heatwaves toward longer, persistence-dominated events in continental interiors through changes in low-frequency temperature modes.
2026-06-28 11:55 UTC pith:4EQ3GZQJ
load-bearing objection The fixed historical threshold likely drives much of the reported persistence increase mechanically, so the dynamical link to low-frequency modes needs stronger checks. the 2 major comments →
Multiscale Dynamics of Heatwave Persistence and Intensity Under Climate Change
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Results show a clear shift toward persistence-dominated heatwave regimes in the continental interior. By the late century, increases in seasonal heatwave days are accompanied by much longer events, with regional HWMD reaching about 26.66 days/event and HWD about 69 days, together with stronger above-threshold intensity, with HWI reaching about 6.88 K. Dynamical diagnostics indicate a redistribution of dominant activity toward lower-frequency levels and weaker effective damping in interior regions, while coastal and maritime regions show smaller changes. Heatwave-relevant low-frequency modes remain active during long events and align with persistence and intensity hotspots, supporting a proce
What carries the argument
multiresolution dynamic mode decomposition (mrDMD) of the daily mean temperature field, connected to heatwave metrics through heatwave-conditioned mode participation ratios and spatial alignment analyses.
Load-bearing premise
The fixed historical 90th percentile threshold, minimum-duration rule, and mrDMD mode participation on this particular regional model output capture the true causal drivers of persistence changes without major influence from model biases or threshold sensitivity.
What would settle it
Observational records or independent model ensembles for the late 21st century showing no substantial rise in interior heatwave duration or no spatial alignment between low-frequency modes and persistence hotspots would falsify the central claim.
If this is right
- Seasonal heatwave days increase alongside much longer individual events in continental interiors.
- Above-threshold intensity strengthens, reaching regional values around 6.88 K by late century.
- Dominant dynamical activity redistributes toward lower-frequency levels with weaker effective damping inland.
- Coastal and maritime regions exhibit comparatively small shifts in the same metrics.
- Low-frequency modes stay active during extended events and overlap with persistence and intensity hotspots.
Where Pith is reading between the lines
- Adaptation planning for prolonged heat events may need to focus more on interior continental zones than coastal ones.
- The workflow could be reapplied to other variables such as humidity or precipitation to check whether similar multiscale shifts appear.
- If low-frequency modes prove robust across models, targeted observational campaigns could test whether those modes strengthen in reality as projected.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper develops an event-dynamical workflow that identifies warm-season heatwaves over southern Canada using a fixed 2001-2010 90th-percentile threshold and minimum-duration criterion on CORDEX-NAM12 output, computes standard metrics (HWF, HWN, HWMD, HWD, HWI, HWM), and links them to multiresolution dynamic mode decomposition (mrDMD) via heatwave-conditioned mode participation ratios and spatial alignment analyses. It reports a shift to persistence-dominated regimes in the continental interior under SSP5-8.5, with late-century values reaching HWMD ≈ 26.66 days/event, HWD ≈ 69 days, and HWI ≈ 6.88 K, accompanied by redistribution toward lower-frequency mrDMD modes and weaker damping.
Significance. If the reported dynamical linkages are robust, the work offers a mechanistic bridge between multiscale temperature variability and regional heatwave amplification, moving beyond frequency-only projections; the mrDMD-based participation and alignment diagnostics represent a constructive approach for process-informed attribution in regional climate modeling.
major comments (2)
- [Abstract and Methods] Abstract and Methods (threshold definition): The central persistence claims (HWMD reaching 26.66 days/event and HWD ≈ 69 days) rest on events defined by a fixed 2001-2010 90th-percentile threshold applied to SSP5-8.5 runs. Because mean warming increases the frequency and duration of exceedances even if the variability spectrum is unchanged, the subsequent mrDMD conditioning on those days risks coupling the reported shift to lower-frequency modes and weaker damping directly to the threshold choice rather than to independent dynamical evolution. A sensitivity analysis with time-varying or relative thresholds is required to establish that the persistence amplification is not an artifact of the fixed-threshold definition.
- [Results] Results (numerical projections): The reported values (HWMD 26.66 days/event, HWI 6.88 K) are presented without accompanying uncertainty ranges, ensemble spread, or baseline comparisons against the reference period; this weakens the ability to judge whether the interior-coastal contrast exceeds internal variability or model bias.
minor comments (2)
- [Methods] The 15-day moving window for the percentile threshold is mentioned but its sensitivity to window length is not quantified.
- [Figures] Figure captions should explicitly state the reference period and threshold used for each panel to avoid ambiguity when comparing periods.
Simulated Author's Rebuttal
We thank the referee for the constructive and detailed comments. We address each major comment point by point below, indicating where revisions will be incorporated.
read point-by-point responses
-
Referee: [Abstract and Methods] Abstract and Methods (threshold definition): The central persistence claims (HWMD reaching 26.66 days/event and HWD ≈ 69 days) rest on events defined by a fixed 2001-2010 90th-percentile threshold applied to SSP5-8.5 runs. Because mean warming increases the frequency and duration of exceedances even if the variability spectrum is unchanged, the subsequent mrDMD conditioning on those days risks coupling the reported shift to lower-frequency modes and weaker damping directly to the threshold choice rather than to independent dynamical evolution. A sensitivity analysis with time-varying or relative thresholds is required to establish that the persistence amplification is not an artifact of the fixed-threshold definition.
Authors: We agree that a fixed historical threshold can increase exceedance frequency and duration under mean warming, potentially affecting the conditioning step. The mrDMD decomposition itself is performed on the full temperature field without reference to the threshold, and mode participation ratios are computed conditionally afterward. This design aims to separate dynamical changes from threshold effects. To directly address the concern, we will add a sensitivity analysis using a time-varying (period-specific) 90th-percentile threshold in the revised manuscript and report the resulting changes in HWMD, HWD, and mrDMD participation ratios. revision: yes
-
Referee: [Results] Results (numerical projections): The reported values (HWMD 26.66 days/event, HWI 6.88 K) are presented without accompanying uncertainty ranges, ensemble spread, or baseline comparisons against the reference period; this weakens the ability to judge whether the interior-coastal contrast exceeds internal variability or model bias.
Authors: We acknowledge that the reported point values lack explicit uncertainty quantification. The analysis is based on a single CORDEX-NAM12 simulation (CRCM5 driven by CanESM5), precluding ensemble spread. We will incorporate direct baseline comparisons of all metrics against the 2001-2010 reference period and explicitly discuss the single-model limitation when interpreting the interior-coastal contrasts. However, ensemble-based uncertainty ranges cannot be added without additional simulations. revision: partial
- Provision of ensemble spread or multi-model uncertainty ranges, as the study employs only one regional climate model simulation.
Circularity Check
No significant circularity detected
full rationale
The paper defines heatwave events via a fixed external historical 90th-percentile threshold (2001-2010) applied to model output, computes standard event metrics, performs mrDMD on the full daily temperature field, and then derives conditioned participation ratios and spatial alignments. None of these steps reduce a claimed result (e.g., HWMD ~26.66 days/event or mode redistribution) to its inputs by construction, nor invoke self-citation load-bearing uniqueness theorems or ansatzes. The derivation chain remains self-contained against the model data and independent decompositions.
Axiom & Free-Parameter Ledger
free parameters (2)
- 90th percentile threshold
- minimum-duration criterion
axioms (1)
- domain assumption CORDEX-NAM12 regional simulations under SSP5-8.5 provide a reliable representation of future temperature variability at the scales relevant to heatwave dynamics
read the original abstract
Climate change is expected to increase heatwave risk, but exceedance frequency alone cannot explain why some regions show stronger amplification in event persistence. This study develops an integrated event-dynamical workflow to diagnose changes in warm-season heatwaves and link them to coherent, multiscale structures of temperature variability. Heatwaves are identified over southern Canada using a fixed historical 90th percentile threshold (2001-2010 reference, 15-day moving window) and a minimum-duration criterion. Events are summarized using frequency (HWF, HWN), persistence (HWMD, HWD), and intensity (HWI, HWM) metrics. The daily mean temperature field is analyzed using multiresolution dynamic mode decomposition (mrDMD). Event and dynamical perspectives are connected through heatwave-conditioned mode participation ratios and spatial alignment analyses between mode-energy footprints and gridded heatwave metrics using hotspot overlap and Spearman rank association. The workflow is applied to CORDEX-NAM12 regional simulations (CRCM5 downscaling of CanESM5) under SSP5-8.5 for 2016-2025, 2051-2060, and 2091-2100. Results show a clear shift toward persistence-dominated heatwave regimes in the continental interior. By the late century, increases in seasonal heatwave days are accompanied by much longer events, with regional HWMD reaching about 26.66 days/event and HWD about 69 days, together with stronger above-threshold intensity, with HWI reaching about 6.88 K. Dynamical diagnostics indicate a redistribution of dominant activity toward lower-frequency levels and weaker effective damping in interior regions, while coastal and maritime regions show smaller changes. Heatwave-relevant low-frequency modes remain active during long events and align with persistence and intensity hotspots, supporting a process-informed interpretation of regional heatwave amplification under climate change.
Reference graph
Works this paper leans on
-
[1]
Adeyeri, O. E., Zhou, W., Ndehedehe, C. E., Ishola, K. A., Laux, P., Akinsanola, A. A., Dieng, M. D. B., & Wang, X. (2025). Global Heatwaves Dynamics Under Climate Change Scenarios : Multidimensional Drivers and Cascading Impacts. Earth’s Future , 13(6), e2025EF006486. https://doi.org/10.1029/2025EF006486 Ann, M., Behrens, J., & Sillmann, J. (2026). Dynam...
-
[2]
https://doi.org/10.1038/s41612-024-00779-y CanESM5 hourly dataset . (2026a, March 15). https://pavics.ouranos.ca/twitcher/ows/proxy/thredds/ncss/grid/datasets/simulations/RCM - CMIP6/CORDEX/NAM-12/1hr/NAM-12_CanESM5_ssp585_r1i1p2f1_OURANOS_CRCM5_v1- r1_1hr_201501010000-210012312300.ncml/dataset.html CanESM5 hourly dataset . (2026b, March 15). https://pavi...
-
[3]
https://doi.org /10.1038/s41467-025- 56067-7 Luo, M., Wu, S., Lau, G. N.-C., Pei, T., Liu, Z., Wang, X., Ning, G., Chan, T. O., Yang, Y., & Zhang, W. (2024). Anthropogenic forcing has increased the risk of longer -traveling and slower-moving large contiguous heatwaves. Science Advances, 10(13), eadl1598. https://doi.org/10.1126/sciadv.adl1598 Mankovich, N...
-
[4]
https://doi.org/10.1038/s41467-020-16970-7 Ranatunga, T., Tong, S. T. Y., & Yang, Y. J. (2016). An approach to measure parameter sensitivity in watershed hydrological modelling. Hydrological Sciences Journal , 1 –17. https://doi.org/10.1080/02626667.2016.1174335 33 Ren, Y., Liu, J., Zhang, T., Shalamzar i, M. J., Arshad, A., Liu, T., Willems, P., Gao, H.,...
-
[5]
https://doi.org/10.3390/rs15051369 Sato, T., Hino, H., & Kusaka, H. (2024). Separating urban heat island circulation and convective cells through dynamic mode decomposition. Atmospheric Science Letters , 25(12), e1279. https://doi.org/10.1002/asl.1279 Schmid, P. J. (2010). Dynamic mode decomposition o f numerical and experimental data. Journal of Fluid Me...
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.