TY - GEN A1 - Boer, H. Marike T. A1 - Stötzel, Claudia A1 - Röblitz, Susanna A1 - Deuflhard, Peter A1 - Veerkamp, Roel F. A1 - Woelders, Henri T1 - A simple mathematical model of the bovine estrous cycle: follicle development and endocrine interactions N2 - Bovine fertility is the subject of extensive research in animal sciences, especially because fertility of dairy cows has declined during the last decades. The regulation of estrus is controlled by the complex interplay of various organs and hormones. Mathematical modeling of the bovine estrous cycle could help in understanding the dynamics of this complex biological system. In this paper we present a mathematical model of the bovine estrous cycle that includes the processes of follicle and corpus luteum development and the key hormones that interact to control these processes. Focus in this paper is on development of the model, but also some simulation results are presented, showing that a set of equations and parameters is obtained that describes the system consistent with empirical knowledge. Even though the majority of the mechanisms that are included are only known qualitatively as stimulatory or inhibitory effects, the model surprisingly well features quantitative observations made in reality. This model of the bovine estrous cycle could be used as a basis for more elaborate models with the ability to study effects of external manipulations and genetic differences. T3 - ZIB-Report - 10-06 KW - cow KW - reproduction KW - hormone patterns KW - ordinary differential equations KW - systems biology Y1 - 2010 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-11682 SN - 1438-0064 ER - TY - GEN A1 - Röblitz, Susanna A1 - Stötzel, Claudia A1 - Deuflhard, Peter A1 - Jones, Hannah M. A1 - Azulay, David-Olivier A1 - van der Graaf, Piet A1 - Martin, Steven T1 - A mathematical model of the human menstrual cycle for the administration of GnRH analogues N2 - This study presents a differential equation model for the feedback mechanisms between Gonadotropin-releasing Hormone (GnRH), Follicle-Stimulating Hormone (FSH), Luteinizing Hormone (LH), development of follicles and corpus luteum, and the production of estradiol (E2), progesterone (P4), inhibin A (IhA), and inhibin B (IhB) during the female menstrual cycle. In contrast to other models, this model does not involve delay differential equations and is based on deterministic modelling of the GnRH pulse pattern, which allows for faster simulation times and efficient parameter identification. These steps were essential to tackle the task of developing a mathematical model for the administration of GnRH analogues. The focus of this paper is on model development for GnRH receptor binding and the integration of a pharmacokinetic/pharmacodynamic model for the GnRH agonist Nafarelin and the GnRH antagonist Cetrorelix into the menstrual cycle model. The final mathematical model describes the hormone profiles (LH, FSH, P4, E2) throughout the menstrual cycle in 12 healthy women. Moreover, it correctly predicts the changes in the cycle following single and multiple dose administration of Nafarelin or Cetrorelix at different stages in the cycle. T3 - ZIB-Report - 11-16 KW - human menstrual cycle KW - systems biology KW - mathematical modelling KW - GnRH analogues KW - differential equations Y1 - 2011 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-12732 ER - TY - GEN A1 - Stötzel, Claudia A1 - Plöntzke, Julia A1 - Röblitz, Susanna T1 - Advances in modelling of the bovine estrous cycle: Administration of PGF2alpha N2 - Our model of the bovine estrous cycle is a set of ordinary differential equations which generates hormone profiles of successive estrous cycles with several follicular waves per cycle. It describes the growth and decay of the follicles and the corpus luteum, as well as the change of the key substances over time. In this work we describe recent improvements of this model, including the introduction of new components, and elimination of time delays. We validate our model by showing that the simulations agree with observations from synchronization studies and with measured progesterone data after a single dose administration of synthetic prostaglandin F2alpha. T3 - ZIB-Report - 11-17 KW - cow KW - reproduction KW - hormone patterns KW - differential equations KW - systems biology Y1 - 2011 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-12740 ER - TY - GEN A1 - Plöntzke, Julia A1 - Berg, Mascha A1 - Olany, Amiya A1 - Leonhard-Marek, Sabine A1 - Müller, Kerstin-Elisabeth A1 - Röblitz, Susanna T1 - Modeling potassium balance in dairy cows N2 - Potassium is fundamental for cell functioning including signal transduction, acid-base- and water-metabolism. Since diet of dairy cows is generally rich in potassium, hypokalemia was not in the focus of research for long time. Furthermore, hypokalemia was not frequently diagnosed because blood potassium content is difficult to measure. In recent years, measurement methods have been improved. Nowadays hypokalemia is increasingly diagnosed in cows with disorders such as abomasal displacement, ketosis or down cow syndrome, calling for intensified research on this topic. In this report we describe the development of a basic mechanistic, dynamic model of potassium balance based on ordinary differential and algebraic equations. Parameter values are obtained from data of a clinical trial in which potassium balance and the influence of therapeutic intervention in glucose and electrolyte metabolism on potassium balance in non-lactating dairy cows were studied. The model is formulated at a high abstraction level and includes information and hypotheses from literature. This work represents a first step towards the understanding and design of effective prophylactic feed additives and treatment strategies. T3 - ZIB-Report - 13-09 KW - systems biology KW - differential equations KW - metabolism Y1 - 2013 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-17749 SN - 1438-0064 ER - TY - THES A1 - Moualeu-Ngangue, Dany Pascal T1 - A Mathematical Tuberculosis Model in Cameroon N2 - This thesis firstly presents a nonlinear extended deterministic model for the transmission dynamics of tuberculosis, based on realistic assumptions and data collected from the WHO. This model enables a comprehensive qualitative analysis of various aspects in the outbreak and control of tuberculosis in Sub-Saharan Africa countries and successfully reproduces the epidemiology of tuberculosis in Cameroon for the past (from 1994-2010). Some particular properties of the model and its solution have been presented using the comparison theorem applied to the theory of differential equations. The existence and the stability of a disease free equilibrium has been discussed using the Perron-Frobenius theorem and Metzler stable matrices. Furthermore, we computed the basic reproduction number, i.e. the number of cases that one case generates on average over the course of its infectious period. Rigorous qualitative analysis of the model reveals that, in contrast to the model without reinfections, the full model with reinfection exhibits the phenomenon of backward bifurcation, where a stable disease-free equilibrium coexists with a stable endemic equilibrium when a certain threshold quantity, known as the basic reproduction ratio (R0), is less than unity. The global stability of the disease-free equilibrium has been discussed using the concepts of Lyapunov stability and bifurcation theory. With the help of a sensitivity analysis using data of Cameroon, we identified the relevant parameters which play a key role for the transmission and the control of the disease. This was possible applying sophisticated numerical methods (POEM) developed at ZIB. Using advanced approaches for optimal control considering the costs for chemoprophylaxis, treatment and educational campaigns should provide a framework for designing realistic cost effective strategies with different intervention methods. The forward-backward sweep method has been used to solve the numerical optimal control problem. The numerical result of the optimal control problem reveals that combined effort in education and chemoprophylaxis may lead to a reduction of 80\% in the number of infected people in 10 years. The mathematical and numerical approaches developed in this thesis could be similarly applied in many other Sub-Saharan countries where TB is a public health problem. Y1 - 2013 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:kobv:188-fudissthesis000000095327-7 UR - http://www.diss.fu-berlin.de/diss/receive/FUDISS_thesis_000000095327 ER - TY - GEN A1 - Moualeu-Ngangue, Dany Pascal A1 - Röblitz, Susanna A1 - Ehrig, Rainald A1 - Deuflhard, Peter T1 - Parameter Identification for a Tuberculosis Model in Cameroon N2 - A deterministic model of tuberculosis in sub-Saharan Africa in general and Cameroon in particular including lack of access to the treatment and weak diagnose capacity is designed and analyzed with respect to its transmission dynamics. The model includes both frequency- and density-dependent transmissions. It is shown that the model is mathematically well-posed and epidemiologically reasonable. Solutions are non-negative and bounded whenever the initial values are non-negative. A sensitivity analysis of model parameters is performed and most sensitive parameters of the model are identified using a state-of-the-art Gauss-Newton Method. In particular, parameters representing the proportion of individuals having access to medical facilities have a large impact on the dynamics of the disease. It has been shown that an increase of these parameter values over the time can significantly reduce the disease burden in the population within the next 15 years. T3 - ZIB-Report - 13-72 KW - Tuberculosis, parameter identification, sensitivity analysis, Gauss-Newton method, Nonlinear dynamical systems. Y1 - 2013 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-43087 SN - 1438-0064 ER - TY - GEN A1 - Moualeu-Ngangue, Dany Pascal A1 - Weiser, Martin A1 - Ehrig, Rainald A1 - Deuflhard, Peter T1 - Optimal control for a tuberculosis model with undetected cases in Cameroon N2 - This paper considers the optimal control of tuberculosis through education, diagnosis campaign and chemoprophylaxis of latently infected. A mathematical model which includes important components such as undiagnosed infectious, diagnosed infectious, latently infected and lost-sight infectious is formulated. The model combines a frequency dependent and a density dependent force of infection for TB transmission. Through optimal control theory and numerical simulations, a cost-effective balance of two different intervention methods is obtained. Seeking to minimize the amount of money the government spends when tuberculosis remain endemic in the Cameroonian population, Pontryagin's maximum principle is used to characterize the optimal control. The optimality system is derived and solved numerically using the forward-backward sweep method (FBSM). Results provide a framework for designing cost-effective strategies for diseases with multiple intervention methods. It comes out that combining chemoprophylaxis and education, the burden of TB can be reduced by 80 % in 10 years T3 - ZIB-Report - 13-73 KW - Tuberculosis KW - Optimal control KW - Nonlinear dynamical systems Y1 - 2013 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-43142 SN - 1438-0064 ER - TY - GEN A1 - Meisel, Kristin A1 - Röblitz, Susanna T1 - Mathematische Modellierung der Knochenremodellierung N2 - Knochenremodellierung beinhaltet den Auf- und Abbau der Knochenmasse durch die verschiedenen Knochenzellen und findet fast überall, auch am Skelett Erwachsener, statt. Für die Erneuerung der Knochensubstanz sind die Osteoblasten zuständig. Sie ersetzen exakt die Menge der Knochenmasse, welche zuvor durch Osteoklasten abgebaut wurde. Störungen dieses, durch viele Faktoren beeinflussten Prozesses führen zu pathologischen Veränderungen, beispielsweise zu Osteoporose oder Arthritis. Auf der Grundlage analysierter publizierter Modelle wurde ein Modell entwickelt, welches den Einfluss einiger dieser Faktoren realistisch abbildet. Auf diese Weise kann die Wirkung von basalem PTH, des RANKL-OPG-Systems ( RANKL: Receptor activator of nuclear factor kappa-B ligand, OPG: Osteoprotegerin) und des Estradiols auf den Knochenstoffwechsel am vorgestellten System untersucht werden. Außerdem wurde durch Estradiolmangel hervorgerufene Osteoporose und der kurative Effekt von synthetischen Medikamenten wie Estradiol oder intermittierend verabreichtem PTH modelliert. Mit der Parameterschätzung anhand des Gauß-Newton-Verfahrens wird des Weiteren eine Methode vorgestellt, die es ermöglicht, mathematische Modelle bestmöglich durch Variation der Parameterwerte an experimentelle Daten anzupassen. Das am Zuse-Institut Berlin entwickelte Softwarepaket POEM wendet diesen Algorithmus an und wird ebenfalls erläutert. T3 - ZIB-Report - 13-41 Y1 - 2013 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-42010 SN - 1438-0064 ER - TY - JOUR A1 - Plöntzke, Julia A1 - Berg, Mascha A1 - Ehrig, Rainald A1 - Leonhard-Marek, Sabine A1 - Müller, Kerstin-Elisabeth A1 - Röblitz, Susanna T1 - Model Based Exploration of Hypokalemia in Dairy Cows JF - Scientific Reports N2 - Hypokalemia, characterized by too low serum potassium levels, is a severe mineral disorder which can be life threatening. It is increasingly diagnosed in veterinarian healthcare and a topic of ongoing research. In this paper, we explore the different originating conditions of hypokalemia: reduced potassium intake, increased excretion, acid-base disturbances, or increased insulin, by using a dynamic mathematical model for potassium balance in non-lactating and lactating cows. Simulations are compared with literature. The results give insights into the network dynamics and point to scenarios on which experimental effort should be focused. Application of mathematical models can assist in experimental planning as well as the reduction, refinement and replacement of animal experiments. Y1 - 2022 U6 - https://doi.org/10.1038/s41598-022-22596-0 VL - 12, 19781 ER - TY - JOUR A1 - Schwieger, Robert A1 - Siebert, Heike A1 - Röblitz, Susanna T1 - Correspondence of Trap Spaces in Different Models of Bioregulatory Networks JF - SIAM Journal on Applied Dynamical Systems N2 - Mathematical models for bioregulatory networks can be based on different formalisms, depending on the quality of available data and the research question to be answered. Discrete boolean models can be constructed based on qualitative data, which are frequently available. On the other hand, continuous models in terms of ordinary differential equations (ODEs) can incorporate time-series data and give more detailed insight into the dynamics of the underlying system. A few years ago, a method based on multivariate polynomial interpolation and Hill functions has been developed for an automatic conversion of boolean models to systems of ordinary differential equations. This method is frequently used by modellers in systems biology today, but there are only a few results available about the conservation of mathematical structures and properties across the formalisms. Here, we consider subsets of the phase space where some components stay fixed, called trap spaces, and demonstrate how boolean trap spaces can be linked to invariant sets in the continuous state space. This knowledge is of practical relevance since finding trap spaces in the boolean setting, which is relatively easy, allows for the construction of reduced ODE models. Y1 - 2018 U6 - https://doi.org/https://doi.org/10.1137/17M1155302 VL - 17 IS - 2 SP - 1742 EP - 1765 ER - TY - JOUR A1 - Lange, Alexander A1 - Schwieger, Robert A1 - Plöntzke, Julia A1 - Schäfer, Stefan A1 - Röblitz, Susanna T1 - Follicular competition in cows: the selection of dominant follicles as a synergistic effect JF - Journal of Mathematical Biology N2 - The reproductive cycle of mono-ovulatory species such as cows or humans is known to show two or more waves of follicular growth and decline between two successive ovulations. Within each wave, there is one dominant follicle escorted by subordinate follicles of varying number. Under the surge of the luteinizing hormone a growing dominant follicle ovulates. Rarely the number of ovulating follicles exceeds one. In the biological literature, the change of hormonal concentrations and individually varying numbers of follicular receptors are made responsible for the selection of exactly one dominant follicle, yet a clear cause has not been identified. In this paper, we suggest a synergistic explanation based on competition, formulated by a parsimoniously defined system of ordinary differential equations (ODEs) that quantifies the time evolution of multiple follicles and their competitive interaction during one wave. Not discriminating between follicles, growth and decline are given by fixed rates. Competition is introduced via a growth-suppressing term, equally supported by all follicles. We prove that the number of dominant follicles is determined exclusively by the ratio of follicular growth and competition. This number turns out to be independent of the number of subordinate follicles. The asymptotic behavior of the corresponding dynamical system is investigated rigorously, where we demonstrate that the ω-limit set only contains fixed points. When also including follicular decline, our ODEs perfectly resemble ultrasound data of bovine follicles. Implications for the involved but not explicitly modeled hormones are discussed. Y1 - 2018 U6 - https://doi.org/https://doi.org/10.1007/s00285-018-1284-0 ER - TY - CHAP A1 - Mancini, Toni A1 - Mari, Federico A1 - Massini, Annalisa A1 - Melatti, Igor A1 - Salvo, Ivano A1 - Sinisi, Stefano A1 - Tronci, Enrico A1 - Ehrig, Rainald A1 - Röblitz, Susanna A1 - Leeners, Brigitte T1 - Computing Personalised Treatments through In Silico Clinical Trials. A Case Study on Downregulation in Assisted Reproduction T2 - Experimental Evaluation of Algorithms for Solving Problems with Combinatorial Explosion, Workshop of the Federated Logic Conference 2018 (FLoC 2018), Oxford N2 - In Silico Clinical Trials (ISCT), i.e., clinical experimental campaigns carried out by means of computer simulations, hold the promise to decrease time and cost for the safety and efficacy assessment of pharmacological treatments, reduce the need for animal and human testing, and enable precision medicine. In this paper we present a case study aiming at quantifying, by means of a multi-arm ISCT supervised by intelligent search, the potential impact of precision medicine approaches on a real pharmacological treatment, namely the downregulation phase of a complex clinical protocol for assisted reproduction. Y1 - 2018 U6 - https://doi.org/10.29007/g864 SP - 1 EP - 16 ER - TY - JOUR A1 - Moewis, Philippe A1 - Boeth, Heide A1 - Heller, Markus O. A1 - Yntema, Caroline A1 - Jung, Tobias A1 - Doyscher, Ralf A1 - Ehrig, Rainald A1 - Zhong, Yanlin A1 - Taylor, William R. T1 - Towards understanding knee joint laxity: Errors in non-invasive assessment of joint rotation can be corrected JF - Medical Engineering and Physics Y1 - 2014 U6 - https://doi.org/doi.org/10.1016/j.medengphy.2014.03.017 VL - 36 IS - 7 SP - 889 EP - 895 ER - TY - GEN A1 - Winkelmann, Stefanie T1 - Markov Control with Rare State Observation: Average Optimality N2 - This paper investigates the criterion of long-term average costs for a Markov decision process (MDP) which is not permanently observable. Each observation of the process produces a fixed amount of \textit{information costs} which enter the considered performance criterion and preclude from arbitrarily frequent state testing. Choosing the \textit{rare} observation times is part of the control procedure. In contrast to the theory of partially observable Markov decision processes, we consider an arbitrary continuous-time Markov process on a finite state space without further restrictions on the dynamics or the type of interaction. Based on the original Markov control theory, we redefine the control model and the average cost criterion for the setting of information costs. We analyze the constant of average costs for the case of ergodic dynamics and present an optimality equation which characterizes the optimal choice of control actions and observation times. For this purpose, we construct an equivalent freely observable MDP and translate the well-known results from the original theory to the new setting. T3 - ZIB-Report - 16-59 KW - Markov decision process KW - partial observability KW - average optimality KW - information costs Y1 - 2016 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-60981 SN - 1438-0064 ER - TY - GEN A1 - Klebanov, Ilja A1 - Sikorski, Alexander A1 - Schütte, Christof A1 - Röblitz, Susanna T1 - Empirical Bayes Methods, Reference Priors, Cross Entropy and the EM Algorithm N2 - When estimating a probability density within the empirical Bayes framework, the non-parametric maximum likelihood estimate (NPMLE) usually tends to overfit the data. This issue is usually taken care of by regularization - a penalization term is subtracted from the marginal log-likelihood before the maximization step, so that the estimate favors smooth solutions, resulting in the so-called maximum penalized likelihood estimation (MPLE). The majority of penalizations currently in use are rather arbitrary brute-force solutions, which lack invariance under transformation of the parameters(reparametrization) and measurements. This contradicts the principle that, if the underlying model has several equivalent formulations, the methods of inductive inference should lead to consistent results. Motivated by this principle and using an information-theoretic point of view, we suggest an entropy-based penalization term that guarantees this kind of invariance. The resulting density estimate can be seen as a generalization of reference priors. Using the reference prior as a hyperprior, on the other hand, is argued to be a poor choice for regularization. We also present an insightful connection between the NPMLE, the cross entropy and the principle of minimum discrimination information suggesting another method of inference that contains the doubly-smoothed maximum likelihood estimation as a special case. T3 - ZIB-Report - 16-56 KW - parameter estimation KW - Bayesian inference KW - Bayesian hierarchical modeling KW - hyperparameter KW - hyperprior KW - EM algorithm KW - NPMLE KW - MPLE KW - DS-MLE KW - principle of maximum entropy KW - cross entropy KW - minimum discrimination information KW - reference prior KW - Jeffreys prior Y1 - 2016 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-61230 SN - 1438-0064 ER - TY - GEN A1 - Klebanov, Ilja A1 - Sikorski, Alexander A1 - Schütte, Christof A1 - Röblitz, Susanna T1 - Empirical Bayes Methods for Prior Estimation in Systems Medicine N2 - One of the main goals of mathematical modelling in systems medicine related to medical applications is to obtain patient-specific parameterizations and model predictions. In clinical practice, however, the number of available measurements for single patients is usually limited due to time and cost restrictions. This hampers the process of making patient-specific predictions about the outcome of a treatment. On the other hand, data are often available for many patients, in particular if extensive clinical studies have been performed. Therefore, before applying Bayes’ rule separately to the data of each patient (which is typically performed using a non-informative prior), it is meaningful to use empirical Bayes methods in order to construct an informative prior from all available data. We compare the performance of four priors - a non-informative prior and priors chosen by nonparametric maximum likelihood estimation (NPMLE), by maximum penalized lilelihood estimation (MPLE) and by doubly-smoothed maximum likelihood estimation (DS-MLE) - by applying them to a low-dimensional parameter estimation problem in a toy model as well as to a high-dimensional ODE model of the human menstrual cycle, which represents a typical example from systems biology modelling. T3 - ZIB-Report - 16-57 KW - Parameter estimation KW - Bayesian inference KW - Bayesian hierarchical modelling KW - NPMLE KW - MPLE KW - DS-MLE KW - EM algorithm KW - Jeffreys prior KW - reference prior KW - hyperparameter KW - hyperprior KW - principle of maximum entropy Y1 - 2016 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:0297-zib-61307 SN - 1438-0064 ER - TY - CHAP A1 - Lange, Alexander A1 - Plöntzke, Julia A1 - Schäfer, Stefan A1 - Röblitz, Susanna T1 - Follicular maturation in cows: mathematical models and data T2 - 10. European Conference on Mathematical and Theoretical Biology Y1 - 2016 ER - TY - GEN A1 - Omari, Mohamed A1 - Plöntzke, Julia A1 - Berg, Mascha A1 - Röblitz, Susanna T1 - Modelling of glucose-insulin metabolism and its effect on the estrous cycle in bovine T2 - Annual Meeting of the Population Approach Group in Europe Y1 - 2016 ER - TY - JOUR A1 - Winkelmann, Stefanie T1 - Markov Control with Rare State Observation: Average Optimality JF - Markov Processes and Related Fields N2 - This paper investigates the criterion of long-term average costs for a Markov decision process (MDP) which is not permanently observable. Each observation of the process produces a fixed amount of information costs which enter the considered performance criterion and preclude from arbitrarily frequent state testing. Choosing the rare observation times is part of the control procedure. In contrast to the theory of partially observable Markov decision processes, we consider an arbitrary continuous-time Markov process on a finite state space without further restrictions on the dynamics or the type of interaction. Based on the original Markov control theory, we redefine the control model and the average cost criterion for the setting of information costs. We analyze the constant of average costs for the case of ergodic dynamics and present an optimality equation which characterizes the optimal choice of control actions and observation times. For this purpose, we construct an equivalent freely observable MDP and translate the well-known results from the original theory to the new setting. Y1 - 2017 VL - 23 SP - 1 EP - 34 ER - TY - JOUR A1 - Winkelmann, Stefanie A1 - Schütte, Christof T1 - The Spatiotemporal Master Equation: Approximation of Reaction-Diffusion Dynamics via Markov State Modeling JF - Journal of Chemical Physics N2 - Accurate modeling and numerical simulation of reaction kinetics is a topic of steady interest.We consider the spatiotemporal chemical master equation (ST-CME) as a model for stochastic reaction-diffusion systems that exhibit properties of metastability. The space of motion is decomposed into metastable compartments and diffusive motion is approximated by jumps between these compartments. Treating these jumps as first-order reactions, simulation of the resulting stochastic system is possible by the Gillespie method. We present the theory of Markov state models (MSM) as a theoretical foundation of this intuitive approach. By means of Markov state modeling, both the number and shape of compartments and the transition rates between them can be determined. We consider the ST-CME for two reaction-diffusion systems and compare it to more detailed models. Moreover, a rigorous formal justification of the ST-CME by Galerkin projection methods is presented. Y1 - 2016 U6 - https://doi.org/10.1063/1.4971163 VL - 145 IS - 21 ER -