Uncoupled high-latitude wave models in COAMPS

J. Yu; R.A. Allard; T.J. Campbell; W.E. Rogers

REVIEW 3 major objections 4 minor 1 references

Uncoupled high-latitude wave models in COAMPS

T0 review · 3 major / 4 minor · reviewed 2026-08-05 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Standard uncoupled runs of SWAN and WAVEWATCH III in the COAMPS regional system, using a new empirical ice-thickness-dependent dissipation formula, produce high-latitude wave hindcasts whose skill is quantifiable, with four non-standard set

desk verdict Abstract describes a useful operational test of a thickness-dependent ice dissipation formula in COAMPS, but the full text is unreadable mojibake in this version, so judgment rests on the abstract alone. read the letter →

arxiv 2508.16755 v1 pith:YXAXZSXH submitted 2025-08-22 physics.ao-ph

W.E. Rogers , T.J. Campbell , J. Yu , R.A. Allard This is my paper · ORCID

classification physics.ao-ph

keywords sea-icewavedissipationhindcasthigh-latitudeoceanuncoupledmodelingCOAMPSSWANWAVEWATCHIIIicethickness

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This report tries to establish that a standard, uncoupled configuration of the wave models SWAN and WAVEWATCH III, run inside the COAMPS regional coupled modeling system, can produce credible high-latitude wave hindcasts in ice-covered seas, and that the choice of sea-ice dissipation parameterization measurably matters. It uses a new empirical dissipation formula that depends on ice thickness, not just concentration. Six cycling hindcast cases are run at 12-hour increments for 4–26 days each, across the Bering Strait, Gulf of Bothnia, Sea of Okhotsk, and Barents Sea. Verification combines visual field inspection, controlled experiments with four altered settings (ice dissipation off, an older thickness-independent formula, higher-resolution forcing, and no surface currents), and, in the two Barents Sea cases, spectral wave data from satellite and from instruments deployed on the ice. The paper's case is that the new settings are workable and that these sensitivity experiments quantify how each choice affects skill.

What carries the argument

The central machinery is the uncoupled wave-model workflow in COAMPS, hosting the phase-averaged spectral wave models SWAN (Simulating WAves Nearshore) and WAVEWATCH III, together with the new empirical dissipation formula that scales wave-energy loss with sea-ice thickness. This formula is what separates the standard configuration from the alternatives, so the experiments that switch it off or replace it with an older thickness-independent formula are testing the physical core of the report. The COAMPS workflow supplies archived global forcing (winds, ice concentration, ice thickness, surface currents) and handles setup, cycling, and post-processing in 12-hour increments.

What would settle it

Run the same two-model configuration and the same four sensitivity experiments in an independent ice-covered sea with its own in-ice wave records, either motion sensors or another satellite mission. If the thickness-dependent dissipation formula is no better than the older thickness-independent formula there, or if its advantage disappears when surface currents are included, the central claim about thickness dependence would be contradicted.

Watch

Extended reading notes

Core claim

The central claim is that the thickness-dependent ice dissipation formula, embedded in the standard uncoupled configuration, yields wave hindcasts whose skill can be quantified, and that each of the four non-standard settings changes that skill in a measurable way. The two Barents Sea demonstrations provide the main external check: satellite-retrieved wave spectra and in-ice motion sensors. For the other four cases, support comes from visual inspection of output fields and from comparing alternative settings against the standard run. The paper also compares against a global lower-resolution wave model, giving the regional runs a baseline. The claim is not that every high-latitude case is val

Load-bearing premise

The skill measured in the two Barents Sea cases, plus visual and model-versus-model checks on the other four, is enough to conclude that the thickness-dependent dissipation formula and the standard settings improve or correctly represent high-latitude wave-ice interactions more generally.

Editorial extensions

If this is right

An operational regional wave-hindcast configuration for ice-covered high latitudes—using archived forcing rather than a full coupled run—is plausible for Barents-like conditions, with documented sensitivity to four settings.
The ice-thickness-dependent dissipation formula is a concrete candidate for other world regions where wave forecasts must cross marginal ice zones of varying thickness.
Because disabling sea-ice dissipation and replacing the formula with an older one both change skill, ice-dissipation parameterization should be treated as a first-order choice, not a detail, in high-latitude wave modeling.
The measured response to omitting surface currents and to raising forcing resolution tells users which components of the forcing chain are worth investing in.
The comparison against a global lower-resolution model frames the value of regional resolution at high latitudes.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

If the thickness-dependent formula is doing the work, similar gains should show up in other ice-covered seas with contrasting thickness regimes, such as the Weddell Sea or the Canadian Arctic; a direct test there would separate the formula's physical merit from Barents-specific conditions.
The verification design leaves room for a stronger test: the satellite spectra and in-ice motion data used here as targets could be redivided—part for calibration, part for validation—to show the empirical formula is not being tuned to the same observations that score it.
Because the four sensitivity experiments change one setting at a time against archived global forcing, the report implies that forcing provenance matters for operational forecasts; publishing forcing fields alongside wave products would let users judge when the configuration is out of its tested envelope.

Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the letter

Two things you should know. First, the abstract describes a useful, bounded operational study: six hindcast cases of SWAN and WW3 in COAMPS, with a new ice-thickness-dependent dissipation formula toggled against older or thinner settings, and verification against CFOSAT/SWIM and ice-mounted sensors for two Barents Sea cases. Second, the body text we received is unreadable mojibake, with a header from a different astro-ph paper, so everything below rests on the abstract. That is not the authors' fault, but it caps what I can endorse.

What the paper does well, if the abstract is honest: it is not a wild claim. It documents a practical configuration where the wave model runs uncoupled inside COAMPS and tests four settings that matter operationally. The external comparisons for the Barents cases are real evidence and the right kind of evidence for a skill claim. The comparison against a global lower-resolution model is also a sensible baseline.

Soft spots, in proportion: the 'quantified skill' is only fully external for two of six cases; the other four are visual inspection and model-versus-model, which is fine for a demonstration but weak for the general claim about settings. The abstract does not report the skill numbers, error bars, or sample sizes. The larger risk is the one in the stress-test note: the new formula is empirical, but the abstract does not say how its coefficients were obtained. If they were tuned to the same Barents observations used in validation, the apparent advantage of thickness-dependent dissipation is partly in-sample. That is a serious, referee-level question, not a reason to reject out of hand. The right fix is for the authors to state calibration data and report a blind or holdout evaluation, even a simple one.

Bottom line: if the full text is readable and the equations/figures match the abstract, this is a competence contribution to operational wave forecasting in the marginal ice zone. It will be useful to people running COAMPS or similar systems, less so to someone looking for new physics. Cite it only after checking the calibration question. I would send it to peer review, with a request to address the calibration/validation overlap and to provide actual numbers. For a reading group, it is a maybe: good as a case study of how operational workflows are validated, but not a breakthrough.

Referee Report

3 major / 4 minor

Summary. The manuscript describes six demonstration hindcasts of two uncoupled wave models, SWAN and WW3, run within the Navy's COAMPS system for high-latitude regions with sea ice. The standard configuration uses archived global forcing (wind, sea ice concentration and thickness, surface currents) and a new empirical, ice-thickness-dependent dissipation formula. Four cases are assessed by visual inspection and by comparing model runs with alternative settings; only the two Barents Sea cases are also compared with external observations (CFOSAT/SWIM satellite spectra and Norwegian Meteorological Institute ice-mounted motion sensors). Experiments vary four settings: disabling ice dissipation, using an older thickness-independent dissipation formula, using higher-resolution forcing, and omitting surface currents. The abstract claims that the impact of these settings on model skill is quantified and that skill is also compared with a global wave model.

Significance. If the claims are fully supported, the paper would provide a practically useful, operational quantification of how ice dissipation, forcing resolution, and surface currents affect uncoupled wave-model skill in ice-covered high-latitude seas. The use of two independent observational datasets, cycling runs of 4--26 days, and systematic alternative configurations are strengths. However, the evidence as described supports quantitative skill claims for only two of six cases, and the calibration provenance of the new empirical dissipation formula is not disclosed. The broader generalization of the configuration's skill and of the settings' sensitivity therefore remains an open question.

major comments (3)

[Abstract] The abstract states that six demonstration cases are presented and that 'the impact of these settings on model skill is quantified by comparison to the observations.' But only the two Barents Sea cases have external observations; the other four are verified by visual inspection and model-versus-model comparisons. Visual inspection and run-to-run differences do not quantify skill against independent truth. Please either report quantitative skill metrics for all cases where observations exist, or explicitly restrict the skill-quantification claim to the Barents Sea cases and present the other four as qualitative demonstrations.
[Abstract / new empirical formula] The dissipation formula is called 'new empirical' and is central to the standard configuration, but the abstract does not say how its coefficients were obtained or on which data they were calibrated. If any coefficients were tuned to the same CFOSAT/SWIM and motion-sensor data used for the Barents Sea verification, the reported superiority over the no-dissipation and thickness-independent settings would be partly in-sample. Please report the formula's coefficients, identify the calibration data, and clearly separate calibration from validation. Adding out-of-sample quantitative metrics for at least one non-Barents case would materially strengthen the claim of transferability.
[Full text] The submitted body text is not readable: it appears as corrupted/mojibake content and includes an extraneous line 'arXiv:2508.16754v1 [astro-ph.EP] 22 Aug 2025'. Without readable methods, equations, tables, and results, it is impossible to check the implementation, the dissipation formula, the verification statistics, or the numerical setup. This is a blocking issue independent of the scientific content. Please resubmit a legible manuscript, and ensure the correct arXiv identifier is attached.

minor comments (4)

[Abstract] The claim that the settings' impact 'is quantified' would be easier to assess if the abstract included at least one numerical skill measure (e.g., bias, root-mean-square error, scatter index) and the sample sizes or time periods used for the Barents Sea comparisons.
[Abstract / Introduction] Please define 'uncoupled' clearly: the wave models are run within COAMPS but do not feed back to the atmosphere/ocean; the current abstract does not make this explicit.
[Abstract / Results] The choice of model (SWAN vs WW3) by region should be justified; if the two models are not expected to behave identically, the aggregation of 'six cases' into a single configuration claim may be misleading.
[Full text] The page headers/footers and the embedded 'arXiv:2508.16754v1 [astro-ph.EP]' line appear to belong to a different submission or are pipeline artifacts; they should be removed or corrected.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity demonstrable: validation is against external observations and alternative settings, not against the formula's own inputs.

full rationale

The central claim is an operational demonstration: six uncoupled COAMPS/SWAN/WW3 hindcast cases, with the standard configuration using a new empirical ice-thickness-dependent dissipation formula, are compared against external observations (CFOSAT/SWIM satellite spectra and Norwegian Meteorological Institute ice-mounted motion sensors) and against alternative settings. The abstract does not assert that the formula is derived from first principles or that the verification data are predicted from the formula by construction. The sensitivity experiments (disabling ice dissipation, using an older thickness-independent formula, higher-resolution forcing, omitting surface currents) are independent perturbations whose skill differences are measured against the same external observations; they are not restatements of the fitted formula. No equation or definition in the available abstract makes a predicted quantity equal to an input by construction, and no load-bearing self-citation chain is visible. The only plausible circularity risk—that the 'new empirical formula' may have been calibrated on the same Barents Sea observations used for verification—is a calibration/validation overlap concern that the abstract neither confirms nor rules out; under the hard rules, speculation about unstated fitting data is not sufficient to claim circularity. The full text is largely unreadable due to encoding corruption, but the abstract's stated design is self-contained against external benchmarks, so the appropriate finding is no significant circularity (score 0).

Assumptions & free parameters 2 free parameters · 3 assumptions · 0 invented entities

Because the body text was unreadable, this ledger is inferred from the abstract. The dominant ledger item is the new thickness-dependent ice dissipation formula: an empirical parameterization whose coefficients are fitted to data by construction, and whose fitting set could not be checked against the validation set from the abstract alone. The remaining items are standard domain assumptions of uncoupled regional wave modeling: external forcing quality, model physics validity in ice, and the absence of feedbacks. No new physical entities (particles, forces, dimensions) are introduced.

free parameters (2)

Coefficients of the new thickness-dependent sea-ice dissipation formula
The abstract calls it 'a new empirical formula for dissipation of wave energy by sea ice that is dependent on ice thickness.' Empirical dissipation formulas carry tunable constants; the constants and their fitting data are not given in the abstract, and the full text is unreadable.
Coefficients of the older thickness-independent dissipation formula
Used as the comparison baseline ('an older formula... which does not depend on ice thickness'). Its constants come from prior literature and are inputs this paper does not control, but they still parameterize the comparison.

assumptions (3)

domain assumption Archived global model output (wind vectors, sea ice concentration and thickness, surface current vectors) is accurate enough to drive regional wave hindcasts at the scales studied.
Abstract: 'forcing comes from archived global model output.' Any bias in these fields transfers directly into the wave model and is not independently corrected.
domain assumption SWAN and WW3 represent wave physics adequately in high-latitude, ice-covered seas at the resolutions used.
The models are run uncoupled below their standard validation envelope; any physics mismatch in the ice-covered regime is absorbed as model error rather than corrected or quantified.
domain assumption Neglect of two-way coupling (waves to ice, waves to currents) does not dominate the skill differences being measured.
The title and abstract emphasize 'uncoupled' operation; the error contribution of omitting feedback is not quantified in the abstract.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Uncoupled high-latitude wave models in COAMPS." pith.science (2026). https://pith.science/paper/YXAXZSXH

@misc{pith2026250816755,
  author       = {Pith},
  title        = {Pith review of: Uncoupled high-latitude wave models in COAMPS},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/YXAXZSXH}},
  note         = {Machine review of arXiv:2508.16755}
}

read the original abstract

This report describes six demonstration cases of two numerical ocean wave models in high latitude regions where waves interact with sea ice. The two wave models are SWAN (Simulating WAves Nearshore) and WW3 (WAVEWATCH III), run in uncoupled mode within the Navy's coupled regional modeling system, COAMPS(R) (Coupled Ocean Atmosphere Mesoscale Prediction System). The COAMPS software handles a large majority of the tasks associated with the setup, running, and post-processing of the wave models. All six cases are cycling runs with 12-hour increments, each providing a continuous hindcast of four to 26 days duration. SWAN is applied in a Bering Strait case and two Gulf of Bothnia cases. WW3 is applied in a Sea of Okhotsk case and two Barents Sea cases. Verification is performed by visual inspection of model output fields, and by comparing model runs with alternative settings. In the standard configuration, forcing comes from archived global model output, including information on surface wind vectors, sea ice concentration and thickness, and surface current vectors, and the wave model uses a new empirical formula for dissipation of wave energy by sea ice that is dependent on ice thickness. The Barents Sea cases are compared to spectral wave data from satellite (SWIM instrument on CFOSAT) and from motion sensors deployed on the ice by the Norwegian Meteorological Institute. Experiments are performed with non-standard settings, 1) disabling the dissipation by sea ice, 2) using an older formula for dissipation by sea ice which does not depend on ice thickness, 3) using higher resolution wind and sea ice concentration forcing fields, and 4) omitting surface currents. The impact of these settings on model skill is quantified by comparison to the observations. The skill is also compared to that from a global (thus, lower resolution) ocean wave model.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

1 extracted references · 1 canonical work pages

[1]

�� ...

work page Pith review arXiv 2025

Pith tools

Reviewed August 5, 2026 · model on record in the stance chip above.

[1] [1]

�� ...

work page Pith review arXiv 2025