A parallel domain decomposition method for coupling of surface and groundwater flows