I have some little problem to give a conception to $\frac{\partial }{\partial \theta}$ and $\frac{\partial }{\partial \varphi}$ on manifold (like $\frac{\partial }{\partial x}$ as well). For example, on the sphere, since $$(x,y,z)=(\sin \varphi\cos\theta,\sin\varphi\sin\theta,\cos\varphi),$$ we get $$\frac{\partial }{\partial \theta}=(-\sin\varphi\sin\theta,\sin\varphi\cos\theta,0)\quad \text{and}\quad \frac{\partial }{\partial \varphi}=(\cos\varphi\cos\theta,\cos\varphi\sin\theta,-\sin\varphi).$$
What is written in my book, it's that is a basis of the tangent space of the sphere (ok, why not... I can imagine that $T_p(\mathbb S^2)$ is a vector space of dimension 2, and locally homeomorphic to $\mathbb R^2$, and above we have a basis). But now, why the notation $\frac{\partial }{\partial \theta}$ and $\frac{\partial }{\partial \varphi}$ ? What would be for example $\frac{\partial f}{\partial \theta}$ of a function $f=f(\theta,\varphi)$ or as well $\frac{\partial f}{\partial \varphi}$ ? How would you compute these derivative on the $\mathbb S^2$ ?
Consider a $2$-sphere with coordinates given by latitude ($\phi$) and longitude ($\theta$). Except at the poles, this gives a local $2$ dimensional coordinate system (each point has many coordinates, but that won't matter).
Through each point, there is a unique curve on which longitude is constant, and through each point except the poles there is a nontrivial curve on which latitude is constant.
In general, partial derivatives of a function $f$ can be viewed as derivatives along curves in which many coordinates are constant, with respect to the change in the coordinate named in the partial derivative. In particular, the partial derivatives $\partial f /\partial \phi$ and $\partial f/ \partial \theta$ can be viewed as derivatives of $f$ along the curves from the previous paragraph. These derivatives can be computed, among other ways, by parameterizing the curves and using the usual formulas for derivatives of a function along parameterized curves.
If you want to get a basis for the tangent space at each point other than the poles, you can use the unit tangent vectors of the curves from the second paragraph, which are basically what you computed in the question.