Significant petrogenetic processes governing the geochemical evolution of magma bodies include magma Recharge (including formation of ‘quenched inclusions’ or enclaves), heating and concomitant partial melting of country rock with possible ‘contamination’ of the evolving magma body (Assimilation), and formation and separation of cumulates by F ractional Crystallization (RAFC). Although the importance of modeling such open-system magma chambers subject to energy conservation has been demonstrated, the effects of concurrent removal of magma by eruption and/or variable assimilation (involving imperfect extraction of anatectic melt from wall rock) have not been considered. In this study, we extend the EC-RAFC model to include the effects of Eruption and variable amounts of assimilation, Ax. This model, called EC-E 0RAcFC, tracks the compositions (trace elements and isotopes), temperatures, and masses of magma body liquid (melt), eruptive magma, cumulates and enclaves within a composite magmatic system undergoing simultaneous eruption, recharge, assimilation and fractional crystallization. The model is formulated as a set of 4 þ t þ i þ s coupled nonlinear differential equations, where the number of trace elements, radiogenic and stable isotope ratios modeled are t, i and s, respectively. Solution of the ECE 0RAcFC equations provides values for the average temperature of wall rock (Ta), mass of melt within the magma body (Mm), masses of cumulates (Mct), enclaves (Men) and wall rock (M o a) and the masses of anatectic melt generated (M a ) and assimilated (cM a ). In addition, t trace element concentrations and i þ s isotopic ratios in melt and eruptive magma (Cm, em, dm), cumulates (Cct, em, dm), enclaves (Cen, er , d o r ) and anatectic melt (Ca, ea, d o a) as a function of magma temperature (Tm) are also computed. Input parameters include the (user-defined) equilibration temperature (Teq), a factor describing the efficiency of addition of anatectic melt (c) from country rock to host magma, the initial temperature and composition of pristine host melt (Tom, C o m, em, d o m), recharge melt (T o r , C o r , er , dor ) and wall rock (T o a, C o a, ea, d o a), distribution coefficients (Dm, Dr, Da) and their temperature dependences (DHm, DHr, DHa), latent heats of transition (melting or crystallization) for wall rock (Dha), pristine magma (Dhm) and recharge magma (Dhr) as well as the isobaric specific heat capacity of assimilant (Cp,a), pristine (Cp,m) and recharge (Cp,r) melts. The magma recharge mass and eruptive magma mass functions,Mr(Tm) andMe(Tm), respectively, are specified a priori. Mr(Tm) and Me(Tm) are modeled as either continuous or episodic (step-like) processes. Melt productivity functions, which prescribe the relationship between melt mass fraction and temperature, are defined for end-member bulk compositions characterizing the local geologic site. EC-E 0RAcFC has potential for addressing fundamental questions in igneous petrology such as: What are intrusive to extrusive ratios (I/E) for particular magmatic systems, and how does this factor relate to rates of crustal growth? How does I/E vary temporally at single, long-lived magmatic centers? What system characteristics are most profoundly influenced by eruption? What is the quantitative relationship between recharge and assimilation? In cases where the extraction efficiency can be shown to be less than unity, what geologic criteria are important and