Robust multiferroicity and magnetic modulation of the ferroelectric imprint field in heterostructures comprising epitaxial HfZrO and Co