This paper investigates a scalable recursive state estimation method for nonlinear two-dimensional (2-D) stochastic complex networks. In these networks, the switching topology is modelled by random variables generated by the Bernoulli distribution. Considering the Extended Kalman Filter structure, we employ the Taylor series expansion to handle the nonlinear functions in the model. Hence, high-order terms emerging in the linearization process are modelled by norm-bounded parameter uncertainties. Due to the uncertainty of such a system, we obtain an upper bound for the estimation error covariance. Furthermore, we obtain an upper bound for the cross covariance to reduce the proposed filter's computational complexity and improve scalability. Within the established setting, next, the proposed filter is designed in such a way that the presented upper bounds are minimized. Finally, a numerical example is provided to illustrate the effectiveness of the proposed state estimation algorithm.