Solution:
Rewrite the expression as
$$\sqrt{a} + \sqrt{b} + \sqrt{c} + \sqrt{d} = \frac{1}{2}\sqrt{4a} + \frac{1}{\sqrt{2}}\sqrt{2b} + \frac{\sqrt{3}}{2}\sqrt{\frac{4}{3}c} + \sqrt{d} $$
Then apply the Cauchy-Schwarz inequality:
$$\left(\frac{1}{2}\sqrt{4a} + \frac{1}{\sqrt{2}}\sqrt{2b} + \frac{\sqrt{3}}{2}\sqrt{\frac{4}{3}c} + \sqrt{d}\right)^2 \leq \left(\frac{1}{4} + \frac{1}{2} + \frac{3}{4} + 1\right)\left(4a + 2b + \frac{4}{3}c+d\right) = \frac{10}{4}\left(\frac{12a + 6b + 4c+ 3d}{3}\right) = \frac{10}{12}(6a+ 2(a+b)+(a+b+c)+ 3(a+b+c+d)) \leq \frac{10}{12}(6+2\cdot 5 + 14 + 3\cdot30)= \frac{10}{12}\cdot 120 = 100 $$
$$\implies \sqrt{a} + \sqrt{b} + \sqrt{c} + \sqrt{d} \leq 10 $$
Explaning the magical solution:
This sounds as a magical solution and you may ask "How have I thought in the right coefficients?".Well, here the answer gets a little messy, but my idea was to use Cauchy-Schwarz inequality once and some suitable changes of variables, until be clear which coefficients I should choose. Let be $\alpha,\beta, \gamma>0$, then
$$E =(\sqrt{a} + \sqrt{b} + \sqrt{c} + \sqrt{d})^2 = \left(\frac{1}{\alpha}\sqrt{\alpha^2 a} + \frac{1}{\beta}\sqrt{\beta^2b} + \frac{1}{\gamma}\sqrt{\gamma^2c} + \sqrt{d}\right)^2 $$
$$\leq \left(\frac{1}{\alpha^2}+ \frac{1}{\beta^2} + \frac{1}{\gamma^2}+ 1\right)(\alpha^2a+ \beta^2b + \gamma^2c + d ) \label{1}\tag{1}$$
where in the last inequality we used the Cauchy-Schwarz inequality. The idea now is use the facts that $a\leq 1$, $a + b\leq 5$, $a + b + c\leq 14$ and $a + b + c + d\leq 30$ considering a clever change of variables. So, let's define $x,y,z \geq 0$ such that
$$\gamma^2 = z+1$$
$$\beta^2 = y + \gamma^2 = y + z +1 $$
$$\alpha^2 = x + \beta^2 = x +y + z + 1 $$
Then, considering the second term in \ref{1}:
$$\alpha^2a+ \beta^2b + \gamma^2c + d = (x+ y + z +1)a + (y+ z +1)b + (z+1) c + d $$
$$ = ax + (a+b)y + (a+b+c)z + a+b+c+d \leq x + 5y + 14z + 30 $$
However, it would be better to let all as function of $\alpha, \beta$ and $\gamma$. We have
$$x + 5y + 14z+ 30 = \alpha^2 - \beta^2 + 5\beta^2 - 5\gamma^2 + 14\gamma^2 - 14 + 30 $$
$$= \alpha^2 + 4\beta^2 + 9\gamma^2 + 16 $$
Thus,
$$E \leq \left(\frac{1}{\alpha^2}+ \frac{1}{\beta^2} + \frac{1}{\gamma^2}+ 1\right)(\alpha^2 + 4\beta^2 + 9\gamma^2 + 16) \label{2}\tag{2}$$
Notice that $\alpha,\beta,\gamma$ are positive arbitrary numbers. So we just need to find $\alpha,\beta,\gamma$ such that \ref{2} is less or equal to 100. I tried some values and was unsuccessful. Then, I had the idea to let \ref{2} as symmetric as possible. Firstly, I rewrote it as
$$\left(\frac{1}{\alpha^2}+ \frac{1}{\beta^2} + \frac{1}{\gamma^2}+ 1\right)(\alpha^2 + 4\beta^2 + 9\gamma^2 + 16) = \left(\frac{4}{\alpha^2}+ \frac{4}{\beta^2} + \frac{4}{\gamma^2}+ 4\right)\left( \left(\frac{\alpha}{2}\right)^2 + 4\left(\frac{\beta}{2}\right)^2 + 9\left(\frac{\gamma}{2}\right)^2 + 4\right) $$
Let be $s_0 = \frac{\alpha}{2}$, $w_0 = \frac{\beta}{2}$ and $t_0 = \frac{\gamma}{2}$. So we get,
$$ E \leq \left(\frac{1}{s_0^2}+ \frac{1}{w_0^2} + \frac{1}{t_0^2}+ 4\right)(s_0^2 + 4w_0^2 + 9t_0^2 + 4) $$
But we still can make better. Let be $s_0= s$, $w_0 = w/\sqrt{2}$ and $t_0 = t/\sqrt{3} $. Then,
$$E \leq \left(\frac{1}{s^2}+ \frac{2}{w^2} + \frac{3}{t^2}+ 4\right)(s^2 + 2w^2 + 3t^2 + 4) $$
One can easily see that $s= w = t=1$ provides
$$ E \leq 100 $$
Then, we show that
$$ \sqrt{a} + \sqrt{b} + \sqrt{c} + \sqrt{d} = \sqrt{E} \leq 10 $$
Notice that we need to check that this choice of variables is consistent with the suppositions on $x,y,z$, that is, that they are non-negatives. If we substitute to find the values of $\alpha, \beta$ and $\gamma$, we will find exactly the very same values of the coefficients proposed for the "magical solution".