There are many unavoidable array errors in practical scenarios, which would drastically increase the sidelobe level (SLL) and distort the performance of radar systems accordingly. In this paper, an effective beampattern synthesis approach is proposed to realize a low peak sidelobe level (PSLL) and wide-nulling in the presence of array errors, thus improving the consequent performance of the radar. In particular, the covariance matrix of the sidelobe region (CMSR) is incorporated into the optimization. Considering the randomness of array errors, the statistical mean method is adopted to obtain a more accurate calculation of the CMSR in the presence of array errors based on a Monte Carlo trial. To efficiently and effectively solve the optimization problem, an advanced metaheuristic algorithm, i.e., the water cycle algorithm (WCA), is adopted when seeking the corresponding optimal weight vectors. Numerical results are provided and discussed to demonstrate the effectiveness of the proposed approach including the results based on a wideband linearly spaced magneto-electric (ME) dipole array.