Ok, I tried to write something below but it gets technical very quickly and its not very well written (too lazy). So if you cannot handle wrong-prints please don't read the text below as it will only be an offense.
With that said let me see if I can shed some light on this issue. Normally when we are dealing with measurement uncertainty the error reported on a manufacturer sheet is a 95% confidence interval of the error. We are here assuming that the error is normally distributed and the error is 2 standard deviations from the mean.
Let us consider measurement X with uncertainty U i.e. we are seeing X+U in our instrument, where U is N(0,s), where s is the reported error of the instrument halved. Let f denote the function transforming our measurement to an angle, i.e. we are interested in f(X+U)-f(X). But notice that f is not linear and f(X+U)-f(X) is not normally distributed anymore, the author mentions "the linear part of the error" which I guess is referring to the fact that f is linearized as a function F and we can proceed in assuming that F(X+U)-F(X) is normally distributed (linearity). The author transforms the 2*Std(U) to 2*Std(F(X+U)-F(X)), i.e. two times the standard deviation.
Another important point here is that although the presence of one error affects the effect of another error so their transformed error (angle) are not independent anymore (they where independent to begin with). This can essentially be disregarded since the influence on angle is so small that considering them as independent is not a great error.
Now we come to the conclusion of this. Consider the four angular errors which according to the above can be regarded as independent, we wish to compute the standard deviation of this total error, so that we can compute a 95% confidence interval. Due to independence the sum of the variance is the variance of the sum, i.e. we can just add the squares together to get the variance of the sum.
To make this clearer we are considering measurements X1,X2,X3,X4 with uncertainties U1,U2,U3,U4, and denote Std(U1)=s1,…,Std(U4)=s4. We are also considering four functions f1,f2,f3,f4 that transforms a measurement to the angle (given the other values), to be clear f1(X1+U1,X2+U2,X3+U3,X4+U4), but we approximated this with F1(X1+U1,X2,X3,X4), (where F1 is a linear approximation of f1 around the point (X1,X2,X3,X4), in the following we will suppress the X2,X3,X4 depencence). Now what the author calculated was F1(X1+2*s1) which due to linearity is 2*Std(F1(X1+U1)), the same goes for F2,F3,F4. Now what we want to calculate is
2*Std(F1(X1+U1)+F2(X2+U2)+F3(X3+U3)+F4(X4+U4)) = 2*sqrt(Var(F1(X1+U1)+F2(X2+U2)+F3(X3+U3)+F4(X4+U4))
where Var indicates Variance (square of standard deviation). Now using independence we get
2*sqrt(Var(F1(X1+U1)+F2(X2+U2)+F3(X3+U3)+F4(X4+U4)) = 2*sqrt(Var(F1(X1+U1))+Var(F2(X2+U2))+Var(F3(X3+U3))+Var(F4(X4+U4)))
Which is then equivalent to what the author performed. So this means that a 95% confidence interval for the angle error is roughly 0.006 degrees when rounded off.
Quite sure somewhere is an error but the overall idea is sound, we are considering the calculated angular errors as standard deviations of independent random variables and thus squaring, summing and then taking the root gives the standard deviation of the sum.