In this paper, a novel hybrid semi-implicit meshless weighted least-squares (WLS) scheme H-SIFPM is developed for the first time by coupling the semi-implicit finite point-set method and finite difference method (FDM), and then applied to predict the time-memory nonlinear behavior dominated by a 2D variable-order time-fractional nonlinear Schrödinger equation (TF-NLSE). The proposed H-SIFPM for TF-NLSE is mainly derived from that: firstly, the variable-order time-fractional with Caputo derivative is approximated by a FDM scheme and then imposed to the approximation of WLS as a boundary condition; secondly, the hybrid process is implicitly discretized in temporal level, and an implicit meshless FPM scheme is developed; thirdly, an iterative technique is adopted to treat the implicit format. Subsequently, the estimated error and stability of the proposed scheme is tentatively proved. Moreover, the accuracy and capacity of the proposed H-SIFPM are also tested by solving several examples with constant/variable-order cases, in which a second-order numerical convergent rate is illustrated. We also demonstrate the advantages of the meshless method by solving the problem on complex irregular domain with local refinement. Finally, the nonlinear dispersion behavior dominated by 2D TF-NLSE/TF-GPE with different boundary conditions are predicted, and the FDM results are also presented for comparison. All numerical results show that the proposed H-SIFPM scheme for TF-NLSEs is stable, accurate and flexible.