Topoflow is a powerful and spatially distributed hydrological model with a user-friendly click interface. Its main purpose is to simulate many different physical processes in a watershed, and its goal is to accurately predict how various hydrological variables will evolve in response to climate forces.

PHREEQC is a computer software (Parkhurst and appelo, 1999) for the calculation of various low-temperature hydrogeochemical reactions. PHREEQC, based on the ion linked water model, can (1) calculate the species and saturation index of matter formation; (2) simulate geochemical inversion process; (3) calculate batch reaction and one-dimensional migration reaction. In addition, PHREEQC coupled with multi-component solute transport model can generate PHAST, a three-dimensional reaction transport simulator for simulating groundwater flow system. PHREEQC is also useful in the investigation of unpolluted and polluted environment, laboratory, field and regional investigation. PHREEQC can investigate various natural and human influenced environmental geochemical reactions, including acid mine water discharge, radioactive waste isolation, pollutant migration, nutrient enrichment, aquifer natural and artificial restoration, aquifer reserve recovery, drinking water treatment, indoor experiment and regional aquifer system. PHREEQC has two graphical user interfaces available. Phreeqcl is a screen driven interface that provides a descriptive explanation for each data item (Charlton and Parkhurst, 2002). PHREEQC used in windows is a row editing object interface. Its additional input features enable it to generate graphic simulation results (post, 2004). 　

The n3dade model is used to evaluate solute transport in three-dimensional nonequilibrium porous media.

HydroGeoSphere model is jointly developed by water loo University, Laval University and hydrogeologic company. It can simulate the flow and solute transport process of groundwater and surface water in three dimensions. The model has two coupling modes for groundwater and surface water: the common node mode assumes that the instantaneous consistency of water head is maintained between the surface and the groundwater body, and the water exchange between the surface and groundwater is not calculated in the coupling process of the model; the dual node mode does not assume the continuity of the water head between the two systems, but describes the surface and ground by Darcy flow relationship The flow exchange between the lower water bodies.