Without any assumption on the regime (stationary, transient, periodic,...), KCL comes from electric charge continuity applied to circuit nodes, i.e. (approximately) zero-dimensional volumes.
In general,
$$\dfrac{d}{dt} \int_{V} \rho + \oint_{\partial V} \mathbf{j} \cdot \hat{\mathbf{n}} = 0 \ ,$$
as you shrink the volume to a "small volume", the volume decreases as $\sim L^3$, while its boundary decreases as $\sim L^2$. Taking the limit for $L \rightarrow 0$, the leading term is the boundary term, and the leading term of the equation reads,
$$0 = \oint_{\partial V} \mathbf{j} \cdot \hat{\mathbf{n}} = \sum_k I_k \ ,$$
i.e. KCL.
Of course, if you take a finite-dimenisonal volume $V$, without any other assumptions, you can store net electric charge inside the volume, and the integral form of the charge balance gives you
$$\dfrac{d}{dt}{Q}_V = \sum_k I_k \ .$$