  We report the first observation of two resonance-like structures in the pi+
chi_c1 invariant mass distribution near 4.1 GeV/c^2 in exclusive B0-bar --> K-
pi+ chi_c1 decays. From a Dalitz plot analysis in which the pi+ chi_c1 mass
structures are represented by Breit-Wigner resonance amplitudes, we determine
masses and widths of: M_1 = (4051 +-14 +20-41) MeV/c^2, Gamma_1 = (82 +21-17
+47-22) MeV, M_2 = (4248 +44-29 +180-35}) MeV/c^2, and Gamma_2 =(177 +54-39
+316-61) MeV; and product branching fractions of BF(B0-bar --> K- Z+_1,2) x
BF(Z+_1,2 --> pi+ chi_c1) = (3.0 +1.5-0.8 +3.7-1.6) x 10^-5 and (4.0 +2.3-0.9
+19.7-0.5) x 10^-5 respectively. Here the first uncertainty is statistical, the
second is systematic. The significance of each of the pi+ chi_c1 structures
exceeds 5 sigma, including the systematic uncertainty from various fit models.
This analysis is based on 657 x 10^6 BB-bar events collected at the Upsilon(4S)
resonance with the Belle detector at the KEKB asymmetric-energy e+e- collider.
