TY - BOOK A1 - Van Leeuwen, Peter Jan A1 - Cheng, Yuan A1 - Reich, Sebastian T1 - Nonlinear data assimilation T3 - Frontiers in applied dynamical systems: reviews and tutorials ; 2 N2 - This book contains two review articles on nonlinear data assimilation that deal with closely related topics but were written and can be read independently. Both contributions focus on so-called particle filters. The first contribution by Jan van Leeuwen focuses on the potential of proposal densities. It discusses the issues with present-day particle filters and explorers new ideas for proposal densities to solve them, converging to particle filters that work well in systems of any dimension, closing the contribution with a high-dimensional example. The second contribution by Cheng and Reich discusses a unified framework for ensemble-transform particle filters. This allows one to bridge successful ensemble Kalman filters with fully nonlinear particle filters, and allows a proper introduction of localization in particle filters, which has been lacking up to now. Y1 - 2015 SN - 978-3-319-18346-6 SN - 978-3-319-18347-3 U6 - https://doi.org/10.1007/978-3-319-18347-3 PB - Springer CY - Cham ER - TY - GEN A1 - Acevedo, Walter A1 - Reich, Sebastian A1 - Cubasch, Ulrich T1 - Towards the assimilation of tree-ring-width records using ensemble Kalman filtering techniques T2 - Postprints der Universität Potsdam : Mathematisch Naturwissenschaftliche Reihe N2 - This paper investigates the applicability of the Vaganov–Shashkin–Lite (VSL) forward model for tree-ring-width chronologies as observation operator within a proxy data assimilation (DA) setting. Based on the principle of limiting factors, VSL combines temperature and moisture time series in a nonlinear fashion to obtain simulated TRW chronologies. When used as observation operator, this modelling approach implies three compounding, challenging features: (1) time averaging, (2) “switching recording” of 2 variables and (3) bounded response windows leading to “thresholded response”. We generate pseudo-TRW observations from a chaotic 2-scale dynamical system, used as a cartoon of the atmosphere-land system, and attempt to assimilate them via ensemble Kalman filtering techniques. Results within our simplified setting reveal that VSL’s nonlinearities may lead to considerable loss of assimilation skill, as compared to the utilization of a time-averaged (TA) linear observation operator. In order to understand this undesired effect, we embed VSL’s formulation into the framework of fuzzy logic (FL) theory, which thereby exposes multiple representations of the principle of limiting factors. DA experiments employing three alternative growth rate functions disclose a strong link between the lack of smoothness of the growth rate function and the loss of optimality in the estimate of the TA state. Accordingly, VSL’s performance as observation operator can be enhanced by resorting to smoother FL representations of the principle of limiting factors. This finding fosters new interpretations of tree-ring-growth limitation processes. T3 - Zweitveröffentlichungen der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe - 892 KW - proxy forward modeling KW - data assimilation KW - fuzzy logic KW - ensemble Kalman filter KW - paleoclimate reconstruction Y1 - 2020 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:517-opus4-436363 SN - 1866-8372 VL - 46 IS - 892 SP - 1909 EP - 1920 ER - TY - JOUR A1 - Escribano, Bruno A1 - Akhmatskaya, Elena A1 - Reich, Sebastian A1 - Azpiroz, Jon M. T1 - Multiple-time-stepping generalized hybrid Monte Carlo methods JF - Journal of computational physics N2 - Performance of the generalized shadow hybrid Monte Carlo (GSHMC) method [1], which proved to be superior in sampling efficiency over its predecessors [2-4], molecular dynamics and hybrid Monte Carlo, can be further improved by combining it with multi-time-stepping (MTS) and mollification of slow forces. We demonstrate that the comparatively simple modifications of the method not only lead to better performance of GSHMC itself but also allow for beating the best performed methods, which use the similar force splitting schemes. In addition we show that the same ideas can be successfully applied to the conventional generalized hybrid Monte Carlo method (GHMC). The resulting methods, MTS-GHMC and MTS-GSHMC, provide accurate reproduction of thermodynamic and dynamical properties, exact temperature control during simulation and computational robustness and efficiency. MTS-GHMC uses a generalized momentum update to achieve weak stochastic stabilization to the molecular dynamics (MD) integrator. MTS-GSHMC adds the use of a shadow (modified) Hamiltonian to filter the MD trajectories in the HMC scheme. We introduce a new shadow Hamiltonian formulation adapted to force-splitting methods. The use of such Hamiltonians improves the acceptance rate of trajectories and has a strong impact on the sampling efficiency of the method. Both methods were implemented in the open-source MD package ProtoMol and were tested on a water and a protein systems. Results were compared to those obtained using a Langevin Molly (LM) method [5] on the same systems. The test results demonstrate the superiority of the new methods over LM in terms of stability, accuracy and sampling efficiency. This suggests that putting the MTS approach in the framework of hybrid Monte Carlo and using the natural stochasticity offered by the generalized hybrid Monte Carlo lead to improving stability of MTS and allow for achieving larger step sizes in the simulation of complex systems. KW - Force splitting KW - Mollification KW - Generalized hybrid Monte Carlo KW - Molecular dynamics KW - Modified Hamiltonians Y1 - 2015 U6 - https://doi.org/10.1016/j.jcp.2014.08.052 SN - 0021-9991 SN - 1090-2716 VL - 280 SP - 1 EP - 20 PB - Elsevier CY - San Diego ER - TY - JOUR A1 - Aizinger, Vadym A1 - Korn, Peter A1 - Giorgetta, Marco A1 - Reich, Sebastian T1 - Large-scale turbulence modelling via alpha-regularisation for atmospheric simulations JF - Journal of turbulence N2 - We study the possibility of obtaining a computational turbulence model by means of non-dissipative regularisation of the compressible atmospheric equations for climate-type applications. We use an -regularisation (Lagrangian averaging) of the atmospheric equations. For the hydrostatic and compressible atmospheric equations discretised using a finite volume method on unstructured grids, deterministic and non-deterministic numerical experiments are conducted to compare the individual solutions and the statistics of the regularised equations to those of the original model. The impact of the regularisation parameter is investigated. Our results confirm the principal compatibility of -regularisation with atmospheric dynamics and encourage further investigations within atmospheric model including complex physical parametrisations. KW - hydrostatic atmosphere KW - non-dissipative regularisations KW - Lagrangian-averaged equations Y1 - 2015 U6 - https://doi.org/10.1080/14685248.2014.991443 SN - 1468-5248 VL - 16 IS - 4 SP - 367 EP - 391 PB - Routledge, Taylor & Francis Group CY - Abingdon ER -