Abstract: We present a method for surface approximation that yields closed surfaces within geometric guarantees. Through refinement we construct a piecewise linear function on a 3D triangulation within the tolerance volume, whose isosurface is an intersection-free surface triangle mesh. Further, the algorithm provides strict guarantees over global approximation error of the simplification. Each iteration of the algorithm can be parallelized.