Background: Nuclear matter fourth-order symmetry energy Esym,4(ρ) may significantly influence the properties of neutron stars such as the core-crust transition density and pressure as well as the proton fraction at high densities. The magnitude of Esym,4(ρ) is, however, largely uncertain. Purpose: Based on systematic analyses of several popular non-relativistic energy density functionals with mean-field approximation, we estimate the value of the Esym,4(ρ) at nuclear normal density ρ0 and its density dependence, and explore the correlation between Esym,4(ρ0) and other macroscopic quantities of nuclear matter properties. Method: We use the empirical values of some nuclear macroscopic quantities to construct model parameter sets by Monte Carlo method for four different energy density functionals with mean-field approximation, namely, the conventional Skyrme-Hartree-Fock (SHF) model, the extended SkyrmeHartree-Fock (eSHF) model, the Gogny-Hartree-Fock (GHF) model, and the momentum-dependent interaction (MDI) model. With the constructed samples of parameter sets, we can estimate the density dependence of Esym,4(ρ) and analyze the correlation of Esym,4(ρ0) with other macroscopic quantities.Results: The value of Esym,4(ρ0) is estimated to be 1.02 ± 0.49 MeV for the SHF model, 1.02 ± 0.50 MeV for the eSHF model, 0.70 ± 0.60 MeV for the GHF model, and 0.74 ± 0.63 MeV for the MDI model. Moreover, our results indicate that the density dependence of Esym,4(ρ) is model dependent, especially at higher densities. Furthermore, we find that the Esym,4(ρ0) has strong positive (negative) correlation with isoscalar (isovector) nucleon effective mass m * s,0 (m * v,0 ) at ρ0. In particular, for the SHF and eSHF models, the Esym,4(ρ) is completely determined by the isoscalar and isovector nucleon effective masses m * s (ρ) and m * v (ρ), and the analytical expression is given. Conclusions: In the mean-field models, the magnitude of Esym,4(ρ0) is generally less than 2 MeV, and its density dependence depends on models, especially at higher densities. Esym,4(ρ0) is strongly correlated with m * s,0 and m * v,0 .