A fully 3D finite-element model used to simulate densi- ty-driven coupled flow and contaminant transport in saturated and unsaturated zones. FEMWATER allows modeling of salinity intrusion and other density-dependent contaminants. Complex stratigraphy can be developed in GMS and directly represented in the model. Solutions can be displayed using realistic 3D plots and animation sequences.