This question involves two related simulations I am reviewing (not creating) that both use MODFLOW-NWT and MT3DMS (i.e., NOT MT3D-USGS). In the results of the first simulation (steady-state flow), cells in the upper two layers are saturated and have finite concentration values, as calculated by MT3DMS. These concentration results are then used as the initial concentrations for the second simulation. In the second simulation, the water table is lowered below the upper two layers (also steady-state flow). These cells are now above the water table and are "inactive," but the results show non-zero concentrations. The unsaturated zone is not being simulated. I understand that water can flow through such "dry" cells in NWT. However, because the original MT3DMS is used I would have thought that the concentration in the dry cells would be zero, i.e., MT3DMS would mark these dry cells as inactive. But when viewing the results in GMS, cells that are clearly marked as dry cells have finite concentration values. I would think that these dry cells would have zero concentration even if MT3D-USGS was being used. If anyone can help me understand this I would appreciate it. Thanks.