A first-order analytical solution for the transport of reactive solutes in physically and chemically heterogeneous porous media is derived and discussed. The solution relies on the assumption of chemical activity described by the local linear equilibrium assumption postulating the existence of a (spatially variable) retardation factor. Retardation factors and log permeabilities modeling heterogeneities are described statistically by random space functions with assigned correlation structure. Correlated as well as uncorrelated physical and chemical heterogeneities are studied. The analytical expressions derived reduce to Dagan's classic solution for the case of nonreactive solute transport and obey asymptotic limits already known from the literature. ¿ American Geophysical Union 1993