Let $(M, g)$ be a Riemannian manifold of dimension $n \geq 4$, and let $W$ be the $(0, 4)$ Weyl tensor on $M$. Let $\{e_i : 1 \leq i \leq n\}$ be an orthonormal frame on $M$, that is, a family of smooth vector fields $e_i: U \to T M$, with $U \subseteq M$ being open and nonempty, that form an orthonormal basis of the tangent space at each point in $U$. Furthermore, let $A(t)$ be an arbitrary family of orthogonal transformations on the tangent space $T M$ smoothly dependent on $t \in \mathbb R$ with the additional property that for $t = 0$, they become the identity function on $T M$.
Let's assume we have
$$W(A(t) e_1, A(t) e_2, A(t) e_3, A(t) e_4) = 0 \tag{1}$$
for all $t \in \mathbb R$. I now want to calculate the derivative of this equation with respect to $t$ at $t = 0$. More precisely, I want to derive the equation
$$W(B e_1, e_2, e_3, e_4) + W(e_1, B e_2, e_3, e_4) + B(e_1, e_2, B e_3, e_4) + W(e_1, e_2, e_3, B e_4) = 0 \tag{2}$$
where $B$ is a skew-symmetric matrix. Furthermore, I want to show that for each skew-symmetric B, there exists an orthogonal matrix A with such a property.
I'm honestly a bit lost here on how to approach this. I'm guessing that at some point, the chain rule for differentiation comes to play, and arriving at a skew-symmetric matrix in the tensor components also seems plausible as the derivative of the transformation $A(t)$ (represented as a matrix) always is $B(t) A(t)$ for some skew-symmetric matrix function $B(t)$ (source), and since we're differentiating in $t = 0$ where $A(t)$ is the identity, we have $A'(0) = B(0) = B$ for some skew-symmetric matrix $B$.
(Although I'm not really sure how we can show that each such skew-symmetric matrix $B$ is already the derivative of an orthogonal $A(t)$ with $A(0) = I$?)
But my main problem is that I don't know what I can do with the Weyl-tensor, and how exactly I can differentiate it here. I've searched the web a bit for tensor derivatives, and the closest I could find was the wiki article about covariant derivates for tensor fields where the formula remotely looks like the sum I'm supposed to get, but my main problem here is that I'm supposed to differentiate with respect to a real variable that's somehow within this tensor expression, so I have no idea if it's the covariant derivative or something else that comes into play here, and if should be the covariant derivative, then how do I get there from differentiating a function with a real-valued parameter $t$ as input.