Here, u ∈ V := H 1 0 (Ω) is the state and q ∈ Q := BV (Ω) is the control, where BV (Ω) denotes the space of functions of bounded variation (BV) on the interval Ω := (0, 1). The operator A is elliptic and α is a positive real number. The two finite element schemes that will be analyzed are identical in regard to the discretization of state and adjoint state, but they differ in the treatment of the control. In the variational discretization the control is not discretized, while in the second scheme the control is discretized by piecewise constant functions.The significance of the above control problem is given by the use of the BVseminorm q ′ M(Ω) in the objective. This favors piecewise constant controls with only a limited number of jumps, which makes this problem type interesting in many practical applications. The precise functional analytic setting will be provided in the next section.Optimal control problems with BV-controls defined in one space dimension are strongly related to control problems with measures as controls. Both BV optimal control problems and optimal control problems with measures have attracted significant research interest in the recent past, see, e.g. [8,13,14,17,22,23] for the former and [11,12,15,16,29,30] for the latter.Error estimates for PDE-constrained optimal control problems involving measures have been presented in [11,30,31,34,35]. For error estimates of further sparsity promoting optimal control problems with PDEs see for example [19,30]. The literature on error estimates for optimal control problems with controls in BV is rather limited. We are only aware of [14,18]. Error estimates and numerical analysis for inverse problems involving BV-functions are studied in [5,6]. Related discussion of ODE-constrained control problems involving discontinuous functions and their numerical analysis can be found in, e.g., [1,2,10,21,25,26,37,38].The main difficulty in deriving error estimates for the above problem is given by the fact that it lacks certain coercivity properties that are usually employed to obtain error estimates for the controls, for instance by suitably testing the first order necessary optimality conditions. Hence, only error estimates for the state and the adjoint state can be proven in a rather direct manner; these are, however, suboptimal. To obtain an error estimate for the control and also to improve the error estimates for state and adjoint state, we make use of a structural assumption on the Lagrange multiplierΦ arising from the convex subdifferential of the term q ′ M(Ω) . Specifically, we assume thatΦ, which is a C 2 function inΩ, has only finitely many global extreme points and that it exhibits quadratic growth near those points (i.e.,Φ ′′ = 0 near those points; see Assumption 4.4 and Assumption 4.5). Since the jump set of the optimal control is contained in the set of global extreme points ofΦ, see Corollary 1, this assumption implies that the optimal control admits only finitely many jumps, which is a rather typical situation in practice. In addition, it ensures ...