The purpose of this paper is to propose the quasi-linear theory of tsunami run-up and run-down on a beach with complex bottom topography. We begin with the one-dimensional nonlinear shallow-water wave equations, which we consider over a beach of complex geometry that can be modeled by a piecewise continuous function, along with several natural initial and boundary conditions. The primary obstacle in solving this problem is the moving boundary associated with the shoreline motion. To avoid this difficulty, we replace the moving boundary with a stationary boundary by applying a transformation to the spatial variable of the computational domain. A characteristic feature of any tsunami problem is the smallness of the parameter ε=η0/h0, where η0...