Estimating parameters of the nonlinear cloud and rain equation from a large-eddy simulation

Abstract Predator–prey dynamics have been suggested as simplified models of stratocumulus clouds, with rain acting as a predator of the clouds. We describe a mathematical and computational framework for estimating the parameters of a simplified model from a large eddy simulation (LES). In our method, we extract cycles of cloud growth and decay from the LES and then search for parameters of the simplified model that lead to similar cycles. We implement our method via Markov chain Monte Carlo. Required error models are constructed based on variations of the LES cloud cycles. This computational framework allows us to test the robustness of our overall approach and various assumptions, which is essential for the simplified model to be useful. Our main conclusion is that it is indeed possible to calibrate a predator–prey model so that it becomes a reliable, robust, but simplified representation of selected aspects of a LES. In the future, such models may then be used as a quantitative tool for investigating important questions in cloud microphysics.

[1]  A. Chorin,et al.  Stochastic Tools in Mathematics and Science , 2005 .

[2]  David J. Spiegelhalter,et al.  Introducing Markov chain Monte Carlo , 1995 .

[3]  Christopher Jones,et al.  A coherent structure approach for parameter estimation in Lagrangian Data Assimilation , 2017, 1706.04834.

[4]  C. Fox,et al.  A general purpose sampling algorithm for continuous distributions (the t-walk) , 2010 .

[5]  Marc Bocquet,et al.  Data Assimilation: Methods, Algorithms, and Applications , 2016 .

[6]  R. Wood,et al.  Spatial variability of liquid water path in marine low cloud : The importance of mesoscale cellular convection , 2006 .

[7]  A. Sokal Monte Carlo Methods in Statistical Mechanics: Foundations and New Algorithms , 1997 .

[8]  G. Feingold,et al.  On the reversibility of transitions between closed and open cellular convection , 2015 .

[9]  G. Feingold,et al.  Stratocumulus to Cumulus Transition by Drizzle , 2017 .

[10]  G. Feingold,et al.  Adaptive behavior of marine cellular clouds , 2013, Scientific Reports.

[11]  Matthias Morzfeld,et al.  Feature-based data assimilation in geophysics , 2017 .

[12]  Daniel Foreman-Mackey,et al.  emcee: The MCMC Hammer , 2012, 1202.3665.

[13]  Jonathan R Goodman,et al.  Ensemble samplers with affine invariance , 2010 .

[14]  G. Feingold,et al.  Aerosol–cloud–precipitation system as a predator-prey problem , 2011, Proceedings of the National Academy of Sciences.

[15]  Ulli Wolff Monte Carlo errors with less errors , 2004 .

[16]  E. Agee Observations from Space and Thermal Convection: A Historical Perspective. , 1984 .

[17]  Colin J. Cotter,et al.  Probabilistic Forecasting and Bayesian Data Assimilation , 2015 .

[18]  G. Feingold,et al.  A model of coupled oscillators applied to the aerosol–cloud–precipitation system , 2013 .

[19]  Ilan Koren,et al.  Exploring the nonlinear cloud and rain equation. , 2017, Chaos.

[20]  A. P. Siebesma,et al.  Climate goals and computing the future of clouds , 2017 .

[21]  G. Feingold,et al.  Network approach to patterns in stratocumulus clouds , 2017, Proceedings of the National Academy of Sciences.