Knowledge of the solubility of quartz over a broad spectrum of aqueous fluid compositions and T–P conditions is essential to our understanding of water–rock interaction in the Earth’s crust. We propose an equation to compute the molality of aqueous silica, mSiO2(aq), mol·(kg H 2 O) −1 , in equilibrium with quartz and water–salt–CO 2 fluids, as follows:logmSiO2=A(T)+B(T)·log18.0152VH2O∗+2logxH2OHere A(T) and B(T) are polynomials from Manning’s (GCA 58 (1994), 4831) equation for quartz solubility in pure water, and xH2O and VH2O∗ stand for the mole fraction and effective partial molar volume of H 2 O in the fluid, respectively. The value of VH2O∗ is computed from the relation Vmix=xH2OVH2O∗+∑xsVs, where V mix is the molar volume of the fluid mixture (in cm 3 mol −1 ), and x s and V s denote the mole fraction and the intrinsic volume of the solute, s, respectively. Values of V mix may be obtained from experimental data on the fluid mixture or from a reliable equation of state for the mixture. Adoption of the V s values V NaCl =30.8cm 3 mol −1 and VCO2=29.9cm 3 mol −1 permits satisfactory prediction of quartz solubility both in binary and ternary aqueous systems. In lieu of experimental data V s can be estimated from pure substance properties: the intrinsic volumes of molten salts yield V s for the electrolyte components, whereas the excluded volumes of gas species in Redlich–Kwong–Soave-type equations of state yield V s for the volatiles.The accuracy of our density model is only slightly inferior to the empirical regressions that experimentalists have used to interpolate their measurements of quartz solubility. The strength of our model lies in its ability to predict trends in quartz solubility in fluid mixtures over an extremely wide range of T–P–x s conditions relevant to the Earth’s crust, including conditions hitherto unexplored experimentally. This success is attributable to our model having only one adjustable parameter per solute.