Problems in hydrology frequently have moving fronts and dynamic driving mechanisms such as wells. Since the location of important features changes during a simulation, accurate modeling requires uniformly fine resolution or the ability to change resolution during the simulation. We will describe an algorithm for refinement and unrefinement of tetrahedral/triangular meshes that has been implemented in the adaptive hydrology (ADH) code. The codes including the refinement/unrefinement algorithms are implemented in parallel to accommodate problems with large run time and memory requirements. In this paper, we describe the parallel, adaptive grid algorithm used in ADH and show the resulting grids from some example problems.