Modeling of geochemical processes Equilibrium modeling J. Faimon Modeling of geochemical processes Equilibrium modeling A set of functions, f, are defined that must be solved simultaneously to determine equilibrium for a given set of conditions. Many of these functions are derived from * mole-balance equations for each element or element valence state, exchange site, and surface site * mass-action equations for pure phases and solid solutions. * aqueous charge balance Additional functions are derived for * alkalinity * activity of water * gas-phase equilibria * ionic strength * surface-complexation equilibria. Each function is reduced to contain a minimum number of variables, such that the number of functions equals the number of variables. The program uses a modified Newton-Raphson method to solve the simultaneous nonlinear equations. Modeling of geochemical processes Equilibrium modeling In the following equations, lack of a subscript or the subscript "(aq)" will refer to entities in the aqueous phase, "(e)" refers to exchangers, "(g)" refers to gases, "(s)" refers to surfaces, "(ss)" refers to solid solutions, and "(p)" refers to phases. The unknowns for each aqueous species i are the activity, a[i], activity coefficient, g[i], molality, m[i], and moles in solution, n[i]. The activity of water, a[H2O], the ionic strength, m, and the mass of solvent water in an aqueous solution, W[aq]. The following relationships apply a[i] = g[i]m[i] n[i] = m[i] W[aq ] Modeling of geochemical processes Equilibrium modeling Components of base: * A[aq] -- water * A[i] -- dissolved species, base components * A[p] -- mineral phases in equilibrium * A[g] -- gases of known fugacity Another species: * A[j] -- dissolved species (secondary) * A[pj] -- secondary minerals, (they can be virtual) Formation of mineral, Aj, from the base components B = (A[aq], A[i], A[p], A[g]) Modeling of geochemical processes Equilibrium modeling Example Let us have the base B = (H[2]O, H^+, HCO[3]^-), Secondary species are CO[2(aq)] and CO[3]^2-. Build the secondary species from the base components: and Reaction between the secondary species can be expressed by linear combination of both equations (e.g. by subtracting of second equation from the first one): It gives: Modeling of geochemical processes Equilibrium modeling Mass-action equations In general, mass-action equations can be written as where K[i] is a temperature-dependent equilibrium constant, n[m,i] is the stoichiometric coefficient of master species m in species i and M[aq] is the total number of aqueous master species. The values of n[m,i] may be positive or negative. The total moles of an aqueous species i can be derived from the mass-action expression: Modeling of geochemical processes Equilibrium modeling Activity coefficients of aqueous species are defined with the Davies equation: or the extended Debye-Hückel equation: where z[i] is the ionic charge of aqueous species i, and A and B are constants dependent only on temperature. a[i]^o and b[i] are ion-specific parameters in the Debye-Hückel equation. The ionic strength of the aqueous solution is Modeling of geochemical processes Equilibrium modeling Example Calcite dissolution Modeling of geochemical processes Equilibrium modeling Mole-balance equation The total moles of an element in the system are the sum of the moles initially present in the * pure-phase * solid-solution assemblages, * aqueous phase, * exchange assemblage, * surface assemblage, * gas phase, and * diffuse layers of the surfaces. The following function is the general mole-balance equation: Modeling of geochemical processes Equilibrium modeling f[m] is zero when mole-balance is achieved T[m] is the total moles of the element N[p] is the number of phases SS is the number of solid solutions N[ss] is the number of components in solid solution ss N[aq] is the number of aqueous species E is the number of exchangers N[e] is the number of exchange species for exchange site e S is the number of surfaces K[s] is the number of surface types for surface s N[sk] is the number of surface species for surface type s[k ]N[g] is the number of gas-phase components Modeling of geochemical processes Equilibrium modeling Example * Distribution of carbonate species in closed carbonate system Components of base, B = [CaCO[3], H^+, H[2]O, HCO[3]^-] Secondary species, Ca^2+, OH^-, CO[3]^2-, H[2]CO[3], CO[2(g)]... OH^- = H[2]O -- H^+ CO[3]^2- = HCO[3]^- - H^+ H[2]CO[3] = HCO[3]^- + H^+ CO[2(g)] = HCO[3]^- - (H[2]O - H^+) Ca^2+ = CaCO[3] -- HCO[3]^- + H^+ Total mass of carbonate species, m[TOT] = 10^-3 mol Modeling of geochemical processes Equilibrium modeling Aqueous Charge Balance where f[z] is zero when charge balance has been achieved. If the diffuse-layer composition is explicitly calculated, a separate charge-balance equation is included for each surface and the sum of the terms in the parentheses will be zero when surface charge balance is achieved. If the diffuse-layer composition is not calculated, the second term inside the parentheses is zero. Modeling of geochemical processes Equilibrium modeling Example Charge balance of carbonate system Components of base, B = [Ca^2+, H^+, H[2]O, HCO[3]^-] Secondary species OH^-, CO[3]^2-, H[2]CO[3], CO[2(g)]... OH^- = H[2]O -- H^+ CO[3]^2- = HCO[3]^- - H^+ H[2]CO[3] = HCO[3]^- + H^+ CO[2(g)] = HCO[3]^- - (H[2]O - H^+) Charge balance