We consider weak solutions of second order nonlinear elliptic systems of divergence type under subquadratic growth conditions. Via the method of $${\mathcal{A}}$$ -harmonic approximation we give a characterization of regular points up to the boundary which extends known results from the quadratic and superquadratic case. The proof yields directly the optimal higher regularity on the regular set.