•± J BC model is examined in the lowest approximation of cluster variation method.•The ground state phase diagrams are calculated.•The sublattice order-parameters are studied in detail to obtain the phase diagrams.•The model gives tetracritical points, reentrant and double reentrant behaviors.•The FM, AFM, FI, PM and chaotic phase regions are observed. The Blume-Emery-Griffiths (BEG) model is considered on the Bethe lattice (BL) in terms of exact recursion relations (ERR) under the effect of a crystal field (D) which was either turned on or off randomly for a given probability. The repulsive case of biquadratic exchange interaction ((K < 0)) between the nearest-neighbor (NN) spins is assumed and the effects of changing the coordination number are also investigated. The thermal variation of the order-parameters, i.e. dipole and quadrupole moments, is examined to obtain the phase diagrams. Very rich phase diagrams with the second- and first-order phase transition lines, tricritical, bicritical and critical end points and the occurrence of reentrant behavior are observed. Different phase regions were explored which include the ferromagnetic (F), paramagnetic (P), staggered quadrupolar (SQ) and ferrimagnetic (FI) phases.