<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "https://jats.nlm.nih.gov/nlm-dtd/publishing/3.0/journalpublishing3.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" article-type="research-article" specific-use="SMUR" dtd-version="3.0" xml:lang="en">
<front>
<journal-meta>
<journal-id journal-id-type="publisher">EGUsphere</journal-id>
<journal-title-group>
<journal-title>EGUsphere</journal-title>
<abbrev-journal-title abbrev-type="publisher">EGUsphere</abbrev-journal-title>
<abbrev-journal-title abbrev-type="nlm-ta">EGUsphere</abbrev-journal-title>
</journal-title-group>
<issn pub-type="epub"></issn>
<publisher><publisher-name>Copernicus Publications</publisher-name>
<publisher-loc>Göttingen, Germany</publisher-loc>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.5194/egusphere-2026-1817</article-id>
<title-group>
<article-title>A robustness diagnostic framework for NMIP ensembles: Application to NMIP2 soil N&lt;sub&gt;2&lt;/sub&gt;O emission estimates</article-title>
</title-group>
<contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Inatomi</surname>
<given-names>Motoko</given-names>
<ext-link>https://orcid.org/0009-0005-3592-9797</ext-link>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
</contrib>
</contrib-group><aff id="aff1">
<label>1</label>
<addr-line>Institute for Agro-Environmental Science, NARO, Tsukuba, 305-8604, Japan</addr-line>
</aff>
<pub-date pub-type="epub">
<day>11</day>
<month>05</month>
<year>2026</year>
</pub-date>
<volume>2026</volume>
<fpage>1</fpage>
<lpage>17</lpage>
<permissions>
<copyright-statement>Copyright: &#x000a9; 2026 Motoko Inatomi</copyright-statement>
<copyright-year>2026</copyright-year>
<license license-type="open-access">
<license-p>This work is licensed under the Creative Commons Attribution 4.0 International License. To view a copy of this licence, visit <ext-link ext-link-type="uri"  xlink:href="https://creativecommons.org/licenses/by/4.0/">https://creativecommons.org/licenses/by/4.0/</ext-link></license-p>
</license>
</permissions>
<self-uri xlink:href="https://egusphere.copernicus.org/preprints/2026/egusphere-2026-1817/">This article is available from https://egusphere.copernicus.org/preprints/2026/egusphere-2026-1817/</self-uri>
<self-uri xlink:href="https://egusphere.copernicus.org/preprints/2026/egusphere-2026-1817/egusphere-2026-1817.pdf">The full text article is available as a PDF file from https://egusphere.copernicus.org/preprints/2026/egusphere-2026-1817/egusphere-2026-1817.pdf</self-uri>
<abstract>
<p>Multi-model ensembles, such as the global N&lt;sub&gt;2&lt;/sub&gt;O Model Intercomparison Project (NMIP), are essential for quantifying terrestrial nitrous oxide (N&lt;sub&gt;2&lt;/sub&gt;O) fluxes; however, interpreting where ensemble results are reliable and where they are not remains challenging. Here, we propose a robustness diagnostic framework (RDF) that classifies each grid cell into four categories: Robust-increase, Robust-decrease, Divergent and Uncertain, based on the three metrics of ensemble-mean change, model agreement on the sign of change and inter-model standard deviation. By explicitly separating directional disagreement (Divergent) from quantitative disagreement (Uncertain), the RDF reveals the nature of model uncertainty rather than its magnitude alone. We applied the RDF to soil N&lt;sub&gt;2&lt;/sub&gt;O emission estimates from eight terrestrial biosphere models participating in NMIP phase 2 (NMIP2), comparing pre-industrial (1850s) and contemporary (2010s) periods. Of 56,852 valid grid cells, 40.8 % were classified as Robust-increase, 2.6 % as Robust-decrease, 36.1 % as Divergent and 20.5 % as Uncertain. Stratification by land-use type revealed that the nature of uncertainty differed fundamentally: cropland-dominated regions were dominated by Uncertain (72.6 %), indicating an agreement on the direction of N&lt;sub&gt;2&lt;/sub&gt;O increase but a large quantitative spread, whereas forest-dominated regions were dominated by Divergent (50.0 %), indicating disagreement on the direction of change itself. Pasture-dominated regions exhibited the highest robustness (59.1 % Robust-increase). The inter-model spread correlated strongly with nitrogen input intensity (Spearman &amp;rho; = 0.75), and the Divergent to Uncertain transition followed a gradient of cropland fraction. These contrasting patterns implied that different model improvement strategies were needed: observational benchmarking of emission magnitudes for croplands, and improved process understanding of the competition between carbon dioxide (CO&lt;sub&gt;2&lt;/sub&gt;) fertilization and warming effects in forests. The proposed framework is general and applicable to any multi-model ensemble of biogeochemical change.</p>
</abstract>
<counts><page-count count="17"/></counts>
</article-meta>
</front>
<body/>
<back>
</back>
</article>