A multispecies numerical code was developed to simulate flow and mass transport with kinetic adsorption in variable-density flow systems The two-dimensional code simulated the transport of bromide (Br−)  a nonreactive tracer  and lithium (Li+)  a reactive tracer  in a large-scale tracer test performed in a sand-and-gravel aquifer at Cape Cod  Massachusetts