I saw this old mathoverflow question coming up again recently. Arguments in favor of Riemann integration stress its geometric transparency, its simplicity and no-nonsense character.
Q: Why then do we not dispense with all the business about partitions and sub- vs. supersums (sometimes called "lower" and "upper Darboux sums") and just define the integral $\int_a^b f$ of a bounded (piecewise) continuous function $f:[a,b] \to \mathbb{R}$ as the limit $$\lim_{n \to \infty}\frac{b-a}{n}\sum_{j=1}^n f\left(\frac{(n-j)a+jb}{n}\right)$$ or the limit $$\lim_{n \to \infty}\frac{1}{n}\sum_{j=\lceil na \rceil}^{\lfloor nb \rfloor} f\left(\frac{j}{n}\right)$$ or something similar?
Proving $\int_a^b f + \int_a^b g = \int_a^b (f+g)$ is much simpler with these definitions as it is for the Riemann integral. The proof of $\int_a^bf + \int_b^c f=\int_a^c f$ may have become more difficult with my first definition, but not with the second one. Proving the fundamental theorem of calculus and the change of variables theorem is also as easy with these definitions as it is for the Riemann integral. A basic version of Fubini's theorem is accessible with (a trivial generalization of) these definitions. These alternative definitions are equally transparent from a geometric point of view.
EDIT: of course this mode of integration is 'reliable' only for a limited class of functions, just as is the case for the Riemann integral. Staying true to the pragmatic spirit of the Riemann integral, I'm content if this modified version works fine for bounded piecewise continuous functions.