Journal cover Journal topic
Nonlinear Processes in Geophysics An interactive open-access journal of the European Geosciences Union
Journal topic

Journal metrics

IF value: 1.558
IF 5-year value: 1.475
IF 5-year
CiteScore value: 2.8
SNIP value: 0.921
IPP value: 1.56
SJR value: 0.571
Scimago H <br class='widget-line-break'>index value: 55
Scimago H
h5-index value: 22
© Author(s) 2020. This work is distributed under
the Creative Commons Attribution 4.0 License.
© Author(s) 2020. This work is distributed under
the Creative Commons Attribution 4.0 License.

  19 Jun 2020

19 Jun 2020

Review status
This preprint is currently under review for the journal NPG.

Ensemble-based statistical interpolation with Gaussian anamorphosis for the spatial analysis of precipitation

Cristian Lussana, Thomas N. Nipen, Ivar A. Seierstad, and Christoffer A. Elo Cristian Lussana et al.
  • Norwegian Meteorological Institute, Oslo, Norway

Abstract. Hourly precipitation over a region is often simultaneously simulated by numerical models and observed by multiple data sources. An accurate precipitation representation based on all available information is a valuable result for numerous applications and a critical aspect of climate. Inverse problem theory offers an ideal framework for the combination of observations with a numerical model background. In particular, we have considered a modified ensemble optimal interpolation scheme, that takes into account deficiencies of the background. An additional source of uncertainty for the ensemble background has been included. A data transformation based on Gaussian anamorphosis has been used to optimally exploit the potential of the spatial analysis, given that precipitation is approximated with a gamma distribution and the spatial analysis requires normally distributed variables. For each point, the spatial analysis returns the shape and rate parameters of its gamma distribution. The Ensemble-based Statistical Interpolation scheme with Gaussian AnamorPhosis (EnSI-GAP) is implemented in a way that the covariance matrices are locally stationary and the background error covariance matrix undergoes a localization process. Concepts and methods that are usually found in data assimilation are here applied to spatial analysis, where they have been adapted in an original way to represent precipitation at finer spatial scales than those resolved by the background, at least where the observational network is dense enough. The EnSI-GAP setup requires the specification of a restricted number of parameters and specifically the explicit values of the error variances are not needed, since they are inferred from the available data. The examples of applications presented provide a better understanding of the characteristics of EnSI-GAP. The data sources considered are those typically used at national meteorological services, such as local area models, weather radars and in-situ observations. For this last data source, measurements from both traditional and opportunistic sensors have been considered.

Cristian Lussana et al.

Interactive discussion

Status: final response (author comments only)
Status: final response (author comments only)
AC: Author comment | RC: Referee comment | SC: Short comment | EC: Editor comment
[Login for Authors/Editors] [Subscribe to comment alert] Printer-friendly Version - Printer-friendly version Supplement - Supplement

Cristian Lussana et al.

Cristian Lussana et al.


Total article views: 330 (including HTML, PDF, and XML)
HTML PDF XML Total BibTeX EndNote
229 76 25 330 26 27
  • HTML: 229
  • PDF: 76
  • XML: 25
  • Total: 330
  • BibTeX: 26
  • EndNote: 27
Views and downloads (calculated since 19 Jun 2020)
Cumulative views and downloads (calculated since 19 Jun 2020)

Viewed (geographical distribution)

Total article views: 308 (including HTML, PDF, and XML) Thereof 308 with geography defined and 0 with unknown origin.
Country # Views %
  • 1



No saved metrics found.


No discussed metrics found.
Latest update: 30 Sep 2020
Publications Copernicus
Short summary
An unprecedented amount of rainfall data is nowadays available, such as ensemble model output, weather radar estimates and in-situ observations from networks of both traditional and opportunistic sensors. Nevertheless, the exact amount of precipitation, to some extent, eludes our knowledge. The objective of our study is the precipitation reconstruction through the combination of numerical model output with observations from multiple data sources.
An unprecedented amount of rainfall data is nowadays available, such as ensemble model output,...