The Godbillon-Vey invariant occurs in homology theory, and algebraic topology, when conditions for a co-dimension 1, foliation of a 3D manifold are satisfied. The magnetic Godbillon-Vey helicity invariant in magnetohydrodynamics (MHD) is a higher order helicity invariant that occurs for flows, in which the magnetic helicity density h m = A·B = A·(∇ × A) = 0, where A is the magnetic vector potential and B is the magnetic induction. This paper obtains evolution equations for the magnetic Godbillon-Vey field η = A × B/|A| 2 and the Godbillon-Vey helicity density h gv = η·(∇ × η) in general MHD flows in which either h m = 0 or h m = 0. A conservation law for h gv occurs in flows for which h m = 0. For h m = 0 the evolution equation for h gv contains a source term in which h m is coupled to h gv via the shear tensor of the background flow. The transport equation for h gv also depends on the electric field potential ψ, which is related to the gauge for A, which takes its simplest form for the advected A gauge in which ψ = A · u where u is the fluid velocity. An application of the Godbillon-Vey magnetic helicity to nonlinear force-free magnetic fields used in solar physics is investigated. The possible uses of the Godbillon-Vey helicity in zero helicity flows in ideal fluid mechanics, and in zero helicity Lagrangian kinematics of three-dimensional advection are discussed.