If we pick the subset $ U = \{ f\in L^1(\Omega) | f\geq 0\ \text{in $\Omega$}\} $ of $ L^1(\Omega) $, the square root is well-defined on functions in $ U $. Furthermore, it is clear that for any $ f \in U $, we have $$ \int_\Omega |\sqrt{f}|^2 = \int_\Omega |f| < \infty $$ thus $ \sqrt{f} \in L^2(\Omega) $ for all $ f \in U $.
Making the inequality in $ U $ strict instead, let's call this set $ U' $, intuitively we should be able to consider derivatives of $ \sqrt{\cdot} : U' \to L^2(\Omega) $. For the usual $ \sqrt{\cdot} : \mathbb{R}_+ \to \mathbb{R}_+ $, it is not hard to see that the Fréchet derivative at $ x $ in direction $ h $ is given by $ h(2\sqrt{x})^{-1} $, so I'm guessing that this should carry over.
Considering $$ \left\|\sqrt{x+h} - \sqrt{x} - \frac{h}{2\sqrt{x}}\right\|_{L^2(\Omega)} $$ we can rewrite it by division and multiplication by $ \sqrt{x+h}+\sqrt{x} $ in the first two terms, which leaves us with $$ \left\|\frac{h}{\sqrt{x+h} + \sqrt{x}} - \frac{h}{2\sqrt{x}}\right\|_{L^2(\Omega)}. \tag{$\ast$} $$ Now comes the hard part, where I'm unsure. We would like to show that ($\ast$) is $ \mathcal{O}\left(\|h\|_{L^1(\Omega)}^2\right) $; or at least $ \mathcal{O}\left(\|h\|_{L^1(\Omega)}^p\right) $, where $ p > 1 $.