This thesis analyzes partial differential equations related to the coupled surface and subsurface flows and develops efficient high order discontinuous Galerkin (DG) methods to solve them numerically. Specifically, the ...