We introduce stochastic Interaction-Round-a-Face (IRF) models that are related to representations of the elliptic quantum group E τ,η (sl 2 ). For stochasic IRF models in a quadrant, we evaluate averages for a broad family of observables that can be viewed as higher analogs of q-moments of the height function for the stochastic (higher spin) six vertex models.In a certain limit, the stochastic IRF models degenerate to (1+1)d interacting particle systems that we call dynamic ASEP and SSEP; their jump rates depend on local values of the height function. For the step initial condition, we evaluate averages of observables for them as well, and use those to investigate one-point asymptotics of the dynamic SSEP.The construction and proofs are based on remarkable properties (branching and Pieri rules, Cauchy identities) of a (seemingly new) family of symmetric elliptic functions that arise as matrix elements in an infinite volume limit of the algebraic Bethe ansatz for E τ,η (sl 2 ).SYMMETRIC ELLIPTIC FUNCTIONS, IRF MODELS, AND DYNAMIC EXCLUSION PROCESSES 2 foundational works. This article is concerned with two of its very recent and closely related applications to:(a) theory of symmetric functions, with new families of symmetric functions being introduced and new summation identities for older ones being proved, see [7,35,31,32,15,16,26];(b) deriving new exact formulas for averages of observables in two-dimensional integrable lattice models and (1+1)-dimensional random growth models, and using those to study large scale and time asymptotics [18,15,16,4,1,8,14,9,22].All the above cited papers deal with the R-matrix for the (higher spin) six vertex model and its degenerations (in other words, with representations of the affine quantum group U q (ŝl 2 ) and their limits). Certain progress has been achieved for the U q (ŝl n ) case as well, with new Markov chains introduced via the corresponding R-matrices [30,17] and a duality functional for them provided in [29].The goal of this work is to climb higher in the hierarchy and to extend some of the recent progress from vertex models to the so-called Interaction-Round-a-Face (IRF) models, also known as face and solid-on-solid (SOS) models. The corresponding R-matrices satisfy a face version of the star-triangle relation also known as the dynamical Yang-Baxter equation.The IRF models were originally introduced by Baxter [5] as a tool to analyze the eight vertex model, but they quickly became a subject on their own, see, e. g., [6,36] and references therein. We will only be concerned with the sl 2 case, the basic instance of which is due to the original work [5] and is often called the eight vertex SOS model, and whose fused versions were introduced and extensively studied in [21,19,20].The IRF models were framed in a representation theoretic language by Felder [23], who introduced the concept of an elliptic quantum group. The simplest instance is the quantum elliptic group E τ,η (sl 2 ); the corresponding representation theory and the algebraic Bethe ansatz were developed by F...