One can motivate the covariant differentiation using only vector calculus. It works for an oversimplified case though (but since the OP doesn't accept either the definition via Ehresmann connection nor the vector bundle definition, may be it's justified.)
Consider new coordinates $y^i$ on $\mathbb{R}^n$ (e.g. spherical). We require orthonormality on that coordinate system, $$\mathrm{d}s^2=(u_i)^2\mathrm{d}y^i.$$
One has $\mathbf{x}=x^i(y)e_i$ and defines $$e'_j=\frac{\partial\mathbf{x}(y) }{\partial y^j}.$$ Then the metric is given by $$g_{ij}=e_i'\cdot e_j'.$$
If one considers a vector field $X:\mathbb{R}^n\to \mathbb{R}^n$ one can write $X=X^ie_i'$. If we now wish to differentiate $X$, we have to take into account the change of the components $X^i$ and of the basis $e'_j$, which is no longer rigid. That is
$$\frac{\partial X }{\partial y^j}=\frac{\partial (X^i e'_i)}{\partial y^j}=\frac{\partial X^i}{\partial y^j}e_i'+X^i\frac{\partial e'_i}{\partial y^j}.$$
One can write $\frac{\partial e'_i}{\partial y^j}$ as linear combination of the $e_i'$ s, i.e. for some functions $\Gamma_{ij}^k$
$$ \frac{\partial e'_i}{\partial y^j}=\Gamma_{ij}^k e_k'.$$
Upon taking inner product of this equation with $e_l'$, one sees that these coefficients are given by
$$\Gamma_{ij}^k=g^{kl}e_l'\cdot \frac{\partial e'_i}{\partial y^j}.$$
Now, in index notation, the covariant derivative of $X^i$ is given by the
$$\nabla_jX^i=\frac{\partial X^i}{\partial y^j}+\Gamma_{jk}^iX^k.$$
This is of the form $\frac {{\mathcal D}f(x)}{dx}=\frac {df(x)} {dx} +\delta f(x)$, but $f$ must be a vector field (or higher rank tensor), otherwise the covariant and ordinary derivatives concide.