A mass balance model for mercury based on the fugacity concept is applied to Lake Superior, Lake Michigan, Onondaga Lake and Little Rock Lake to evaluate model performance, analyze cycling of three mercury species groups (elemental, divalent and methyl mercury), and identify important processes that determine the source-to-concentration relationship of the three mercury species groups in these lakes. This model application to four disparate ecosystems is an extension of previous applications of fugacity-based models describing mercury cycling. The model performs satisfactorily following site-specific parameterization, and provides an estimate of minimum rates of species interconversion that compare well with literature. Volatilization and sediment burial are the main processes removing mercury from the lakes, and uncertainty analyses indicate that air-water exchange of elemental mercury and water-sediment exchange of divalent mercury attached to particles are influential in governing mercury concentrations in water. Any new model application or field campaign to quantify mercury cycling in a lake should consider these processes as important.