Execution-grounded evaluation reveals hidden failures in language-model calculations for environmental science
Authors: Maohao Ran, Chendong Ma, Yanting Zhang, Dailing Jiang, Yusen Huang, Meng Gao, Jun Song
Organizations: Department of Geography, Hong Kong Baptist University, Hong Kong, China. · Hong Kong University of Science and Technology, Hong Kong, China. · The Chinese University of Hong Kong (Shenzhen), Shenzhen, China.
Large language models are increasingly used for quantitative work in the environmental sciences, yet existing evaluations score only final answers, leaving calculation process unobserved. Here we introduce AtmosCoder-Bench, an execution-grounded benchmark that makes the calculation process visible. Built through a transferable semi-automated pipeline (436 problems, 3,910 variants, 7,029 graded quantities), every problem is validated to be unambiguous and human-solvable, with uniquely verifiable answers. We find that (i) multiple-choice formats inflate measured accuracy by at least 12 percentage points; (ii) many failures arise not from missing knowledge but from models failing to apply known formulas and constraints consistently throughout multi-step computation; and (iii) even frontier models remain weak when task-specific conditions invalidate familiar methods, often reverting to canonical solution patterns rather than adapting methods to the relevant physical regime, leaving expert oversight essential.
Figures & tables
Category
Problems
Follow-up
Numeric
Paraphrase
queries
variants
variants
Atmospheric dynamics
116
47
520
580
Atmospheric thermodynamics
89
59
395
445
Atmospheric chemistry
51
76
160
255
Air quality
37
8
130
185
Boundary layer
32
16
145
160
Table 1: Composition of AtmosCoder-Bench by category. 436 base problems from 13 established textbooks spanning introductory to advanced level (undergraduate to graduate), each labelled with one of ten atmospheric–environmental categories. Follow-up queries counts the additional sub-questions a problem poses beyond its first (a problem asking for parts (1)–(3) contributes two); every sub-question is graded independently. Each perturbable problem carries five parameter-perturbed numeric variants (346 problems; contamination resistance) and every problem five paraphrase variants (phrasing robustness; answer unchanged).
Figure 1: The AtmosCoder-Bench admission pipeline. Textbook items pass a sequence of gates—extraction for self-contained numeric answerability, an executability gate, cross-vendor independent-derivation agreement (textbook-anchored or consensus-frozen), an AST + perturbation faithfulness audit, a well-posedness audit and triple-gated de-scaffolding—before certified, contamination-resistant variants are generated. Items failing a gate are removed (right), so the admitted set carries verifiable, execution-grounded ground truth.
Model
Option mode
Code mode (grading tolerance)
Δ
raw
chance-corr.
±5%
±10%
±20%
Gemini-3.1-Pro (r)
94.5%
92.7%
60.0%
61.0%
61.9%
+30.8
DeepSeek-R1 (r)
88.5% ∗
84.7%
53.0%
54.6%
55.7%
+29.0
GPT-5.5
80.1%
73.5%
49.6%
50.6%
53.0%
+20.5
DeepSeek-V4-flash
88.1%
84.1%
43.3%
44.0%
45.8%
+38.3
DeepSeek-V3
63.3% ∗
51.1%
27.2%
28.1%
29.6%
+21.5
Table 2: Multiple-choice inflation: option mode versus computed answers on identical problems. Δ is chance-corrected option accuracy minus computed accuracy at the most permissive ±20% tolerance; “(r)” denotes a reasoning configuration. ∗ Option accuracy from the source benchmark [ 33 ] .
Figure 2: Overall model performance under execution-grounded evaluation. a , Code-mode accuracy (mean ± s.d. over three runs) for the 16 leaderboard configurations, coloured by reasoning setting (labels ending in (R) denote reasoning enabled), together with the two domain-adapted ClimateGPT models shown separately as an out-of-distribution reference. b , Distribution of output tokens per problem under the two protocols, for the six configurations run under both, with the median direct-to-code ratio at the right. c , Accuracy against token cost; arrows join the non-reasoning and reasoning settings of the same backbone, dashed contours mark equal tokens per correct answer, and † marks the two configurations whose token counts are lower bounds because the endpoint returns only a summary of the reasoning.
Figure 3: Performance variation across atmospheric domains and failure modes. a , Accuracy across atmospheric subfields and difficulty levels. b , Per-trap outcomes for each configuration across the six trap mechanism families. Red denotes failures returning the exact shortcut value predicted during trap construction, whereas grey denotes other failures; pale grey marks a solved trap.
Figure 4: Failure modes on the Dogpatch Skunk Works Gaussian-plume task. No configuration solves it, and the failures differ in kind: strong models misreport the parameterisation used, intermediate models misclassify regimes or alter formulas, and weaker models show unstable sign, unit, and arithmetic errors.
SciCode is the standard measure of the scientific-coding ability of language models: research-level problems that demand both frontier scientific theory and its implementation as working numerical code. It is a component of the Artificial Analysis Intelligence Index and a standing evaluation in government and national-laboratory suites. Yet its scores have recently plateaued: the strongest 2026 models cluster tightly around 60% subproblem accuracy, and a successor model ties its predecessor. We trace this stagnation to defects in the benchmark itself. A per-problem, domain-expert audit of all 65 test problems uncovers 263 defects; 192 of them, spread across 91% of the main problems, cause correct, instruction-following solutions to be wrongly rejected---through non-reproducible gold answers, over-tight tolerances, or self-contradictory specifications. Critically, 78% of these score-suppressing defects require specialized physics or mathematics knowledge to detect, not mere clerical proofreading. We corrected every confirmable defect to produce SciCode-Verified. The corrections add only the specifications a well-posed problem requires, repair grading, and tighten the tests that were too lenient; every change is recorded with its justification and independently re-checked by a second domain expert. We re-evaluate twelve frontier model snapshots on the corrected benchmark and find a substantial recovery: subproblem accuracy rises from 45--60% to 84--98%, and main-problem accuracy from 9--27% to 69--92%. State-of-the-art models are far more proficient in scientific coding than SciCode has suggested---the bottleneck was not model capability, but the quality of the evaluation instrument. We release SciCode-Verified with its complete audit trail as the corrected public standard.
Sihan Hu, Lyuhan Huang, Youjin Deng +1
Hefei National Laboratory College of Mechanical and Electrical Engineering University of Science and Technology of China · Hefei National Laboratory Institute of Theoretical Physics University of Science and Technology of China Chinese Academy of Sciences Hefei 230026, China · Chinese Academy of Sciences Beijing 100190, China
Large language models (LLMs) often achieve strong performance on reasoning benchmarks, but final-answer accuracy alone does not show whether they faithfully execute the procedure specified in a prompt. We introduce a controlled diagnostic benchmark for procedural execution, where models are given a step-wise arithmetic procedure and two numeric inputs, and must return the final computed value. Complexity is varied through procedure length and look-back dependencies over intermediate variables. Average first-answer accuracy drops from 63% on 5-step procedures to 20% on 95-step procedures. Generation-level analysis shows that failures often involve missing answers, premature answers, self-correction after an initial error and under-executed traces. These findings suggest that apparent reasoning ability can mask substantial weaknesses in faithful long-horizon procedural execution.
Language-model benchmarks collapse two distinct measurement questions into a single accuracy score: whether a response reached an evaluable state, and whether its answer was judged correct. We introduce a two-layer evaluation framework that separates scorer-independent execution evidence, including termination, answer exposure, parseability, and completion length, from scorer-dependent correctness. Across 2,550 outputs from five fixed Qwen and DeepSeek configurations on MATH and ARC-Challenge, matched 2,048-token limits produce sharply different execution mixtures: 49 of 450 Qwen MATH outputs terminate without a final answer, compared with 5 of 300 DeepSeek MATH outputs and none of the 750 ARC outputs. Among the same 300 DeepSeek MATH question-model pairs, no missing-final length termination is observed at 8,192 tokens. A coverage-audited targeted verification study further shows that candidate-selection and aggregation policies can substantially alter comparative accuracy estimates. These results demonstrate that accuracy conflates execution case mix with verification policy. Evaluations of test-time methods should therefore report pre-intervention execution states, verification coverage, and scorer provenance alongside accuracy.
Zongyou Yang, Yinghan Hou
Dyson School of Design Engineering Imperial College London London, United Kingdom · Department of Electrical and Electronic Engineering Imperial College London London, United Kingdom