A measurement of WZ electroweak (EW) vector boson scattering is presented. The measurement is performed in the leptonic decay modes WZ→ℓνℓ′ℓ′, where ℓ,ℓ′=e,μ. The analysis is based on a data sample of proton-proton collisions at s=13 TeV at the LHC collected with the CMS detector and corresponding to an integrated luminosity of 35.9 fb−1. The WZ plus two jet production cross section is measured in fiducial regions with enhanced contributions from EW production and found to be consistent with standard model predictions. The EW WZ production in association with two jets is measured with an observed (expected) significance of 2.2 (2.5) standard deviations. Constraints on charged Higgs boson production and on anomalous quartic gauge couplings in terms of dimension-eight effective field theory operators are also presented.