Journal cover Journal topic
Geoscientific Model Development An interactive open-access journal of the European Geosciences Union
https://doi.org/10.5194/gmd-2017-323
© Author(s) 2018. This work is distributed under
the Creative Commons Attribution 4.0 License.
Development and technical paper
16 Jan 2018
Review status
This discussion paper is a preprint. It is a manuscript under review for the journal Geoscientific Model Development (GMD).
A conservative reconstruction scheme for the interpolation of extensive quantities in the Lagrangian particle dispersion model FLEXPART
Sabine Hittmeir1,2, Anne Philipp2, and Petra Seibert3 1Faculty of Mathematics, University of Vienna, Vienna, Austria
2Department of Meteorology and Geophysics, University of Vienna, Vienna, Austria
3Institute of Meteorology, University of Natural Resources and Life Sciences, Vienna, Austria
Abstract. Lagrangian particle dispersion models require interpolation of all meteorological input variables to the position in space and time of computational particles. The widely used model FLEXPART uses linear interpolation for this purpose, implying that the discrete input fields contain point values. As this is not the case for precipitation (and other fluxes) which represent cell averages or integrals, a preprocessing scheme is applied which ensures the conservation of the integral quantity with the linear interpolation in FLEXPART, at least for the temporal dimension. However, this mass conservation is not ensured per grid cell, and the scheme thus has undesirable properties such as temporal smoothing of the precipitation rates. Therefore, a new reconstruction algorithm was developed, in two variants. It introduces additional supporting grid points in each time interval and is to be used with a piecewise linear interpolation to reconstruct the precipitation time series in FLEXPART. It fulfils the desired requirements by preserving the integral precipitation in each time interval, guaranteeing continuity at interval boundaries, and maintaining non-negativity. The function values of the reconstruction algorithm at the sub-grid and boundary grid points constitute the degrees of freedom which can be prescribed in various ways. With the requirements mentioned it was possible to derive a suitable piecewise-linear reconstruction. To improve the monotonicity behaviour, two versions of a filter were also developed that form a part of the final algorithm. Currently, the algorithm is meant primarily for the temporal dimension. It was shown to significantly improve the reconstruction of hourly precipitation time series from three-hourly input data. Preliminary considerations for the extension to additional dimensions are also included as well as suggestions for a range of possible applications beyond the case of precipitation in a Lagrangian particle model.

Citation: Hittmeir, S., Philipp, A., and Seibert, P.: A conservative reconstruction scheme for the interpolation of extensive quantities in the Lagrangian particle dispersion model FLEXPART, Geosci. Model Dev. Discuss., https://doi.org/10.5194/gmd-2017-323, in review, 2018.
Sabine Hittmeir et al.
Sabine Hittmeir et al.
Sabine Hittmeir et al.

Viewed

Total article views: 313 (including HTML, PDF, and XML)

HTML PDF XML Total Supplement BibTeX EndNote
246 44 23 313 6 5 20

Views and downloads (calculated since 16 Jan 2018)

Cumulative views and downloads (calculated since 16 Jan 2018)

Viewed (geographical distribution)

Total article views: 313 (including HTML, PDF, and XML)

Thereof 313 with geography defined and 0 with unknown origin.

Country # Views %
  • 1

Saved

Discussed

Latest update: 24 Feb 2018
Publications Copernicus
Download
Short summary
Model output of quantities such as precipitation usually represents integrals, for example sums over 3 h. It is not trivial to interpolate a time series of such integral values to instantaneous precipitation rates conserving the integral values. A piecewise-linear reconstruction is presented which fulfils the conservation, is non-negative and continuous at interval boundaries. It shall be used in the FLEXPART Lagrangian dispersion model, but has many other possible applications.
Model output of quantities such as precipitation usually represents integrals, for example sums...
Share