NWP verification against own analysis by using a Data Assimilation confidence mask
Authors: Casati Barbara , Fortin Vincent , Lespinas Franck , Khedhaouiria Dikraa 
Published: June 2021

&lt;p&gt;Numerical Model Prediction (NWP) verification against station measurements from a surface network is affected by sub-tile representativeness issues. Moreover, the station network is often not representative of the whole verification domain (e.g. usually coastal stations are predominant) and large unpopulated regions (such as oceans, Polar regions, deserts) are under-sampled. Verification against gridded analyses mitigate these issues, since they partially address the sub-tile representativeness, and sample homogeneously the verification domain. Moreover, gridded analyses merge station network measurements to radar and satellite retrieval estimates, in a physical coherent fashion, over the same NWP grid. Verification against own analysis, despite quite convenient, is however hampered by its dependence on the NWP background model, which renders the verification &amp;#8220;incestuous&amp;#8221;, further than being affected by the uncertainties introduced by retrieval algorithms and Data Assimilation (DA) procedures.&lt;/p&gt;&lt;p&gt;In this study we investigate the use of a gridded NWP own analysis for verification, by applying a mask to reduce the background model contribution. The mask weights the verification scores to account for the amounts of observations assimilated and their associated uncertainty, as estimated from DA. We illustrate the approach by using the Canadian Precipitation Analysis (CaPA), which assimilates station measurements, radar and satellite-based (IMERG) observations. The CaPA confidence (weighting) mask is dynamic and changes depending on the daily available (assimilated) observations, and on their corresponding DA error statistics; it is defined as&lt;/p&gt;&lt;p&gt;&amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160; &amp;#160;mask = 1 - var(A-O)/var(B-O)&lt;/p&gt;&lt;p&gt;where A=analysis, B=Background, O=observations. Where the analysis is identical to the background model, the weighting mask is zero.&lt;/p&gt;&lt;p&gt;We evaluate the Canadian Regional Deterministic Prediction System (RDPS), which is the NWP system used as background model for CaPA. As expected, the verification results obtained by using the weighting mask lay between the verification results obtained verifying against the analysis over the full domain, and the results obtained verifying against station measurements.&amp;#160;The effects of sub-tile representativeness are quantified by comparing verification results against station measurements to verification results against CaPA for the grid-points co-located with the stations. Finally, the comparison of the verification results against CaPA over the full domain versus the verification results against CaPA for the grid-points co-located with stations, estimates to which extent the station network is representative of the full domain.&lt;/p&gt;&lt;p&gt;The approach aims to propose a simple -yet effective- better practice for verification against own analysis.&lt;/p&gt;