An evolutionary strategy-based error parameterization method that searches for the most ideal error adjustment factors was developed to obtain better assimilation results. Numerical experiments were designed using some classical nonlinear models (i.e., the Lorenz-63 model and the Lorenz-96 model). Crossover and mutation error adjustment factors of evolutionary strategy were investigated in four aspects: the initial conditions of the Lorenz model, ensemble sizes, observation covariance, and the observation intervals. The search for error adjustment factors is usually performed using trial-and-error methods. To solve this difficult problem, a new data assimilation system coupled with genetic algorithms was developed. The method was tested in some simplified model frameworks, and the results are encouraging. The evolutionary strategy-based error handling methods performed robustly under both perfect and imperfect model scenarios in the Lorenz-96 model. However, the application of the methodology to more complex atmospheric or land surface models remains to be tested.