In Hamiltonian systems, conservation laws are closely related to symmetries of evolutionary equations.
In case of modified Boussinesq hierarchy this relationship is especially tight as its entire infinite set of
conservation laws forms a single involutive orbit of a simple one-parameter symmetry group.
We discuss some geometric properties of this symmetry and show how
its properties ensure involutivity of conservation laws.

Recall that the modified Boussinesq system is formed by the following set of partial differential equations
ut = cvxx + uxv + uvx
vt = − cuxx + uux + kvvx
where u = u(x, t), v = v(x, t) are smooth functions on ℝ^{2}
subjected to zero boundary conditions u(±∞, t) = v(±∞, t) = 0, while c and k are some real constants.
In cases k = − 1 and k = 3 modified Boussinesq system has non-trivial
bi-Hamiltonian structure that drastically simplifies analysis of the system in these sectors.
The first case is described in [2],[5],[6],
while in the present paper we focus on the second sector and show that in case k = 3
bi-Hamiltonian structure of modified Boussinesq system is related to non-Noether symmetry [1]
of equations (1).
Thus in case k = 3 modified Boussinesq equations
ut = cvxx + uxv + uvx
vt = − cuxx + uux + 3vvx
can be rewritten in bi-Hamiltonian form
ut = W(dh ∧ du) = Ŵ(dĥ ∧ du)
vt = W(dh ∧ dv) = Ŵ(dĥ ∧ dv)
where W and Ŵ are compatible Poison bivector fields, i.e.
[W , W] = [W , Ŵ] = [Ŵ , Ŵ] = 0
defined as follows
W = ∫− ∞+ ∞ ½(A ∧ Ax + B ∧ Bx)dx
Ŵ = ∫− ∞+ ∞ (uB ∧ Ax + vB ∧ Bx − cAx ∧ Bx)dx
Note that A, B are vector fields that for every smooth functional R = R(u) are defined
via variational derivatives
A(R) = δRδu, B(R) = δRδv.
Corresponding Hamiltonians in bi-Hamiltonian realization (3) are
h = ½∫− ∞+ ∞ (u^{2}v + v^{3} + 2cuvx)dx
ĥ = ½∫− ∞+ ∞ (u^{2} + v^{2})dx
This bi-Hamiltonian structure is related to symmetry of equations (2), but before we proceed let
us remind that symmetry of evolutionary equations is given by the group of transformations
(u , v) ↦ (g(u) , g(v))
which commutes with time evolution
ddtg(u) = g(ut), ddtg(v) = g(vt)
In case of continuous one-parameter groups of transformation
g(u) = e^{zLE}(u) = u + zLEu + ½z^{2}(LE)^{2}u + ⋯
g(v) = e^{zLE}(v) = v + zLEv + ½z^{2}(LE)^{2}v + ⋯
generated by some vector field E, relation (9) gives rise to the following
conditions for the generator of symmetry E
E(u)t = cE(v)xx + E(u)xv + uE(v)x + uxE(v) + E(u)vx
E(v)t = − cE(u)xx + uE(u)x + 3vE(v)x + E(u)ux + 3E(v)vx
Among solutions of equations (11) there is one important vector field —
the generator of non-Noether symmetry which has the following form
E = ∫− ∞+ ∞
{[xuv + 2t(u^{3} + 3uv^{2} + 6cvvx − 2c^{2}uxx)]Ax − cxvAxx
+ (xuux + xvvx)B +
[xu^{2} + 2xv^{2} + 2t(5v^{3} + 3u^{2}v − 6cvux − 2c^{2}vxx)]Bx
+ cxuBxx}dx
Applying one-parameter group of transformations
g(z) = e^{zLE}
generated by the vector field E to the centre of Poisson algebra
which in our case is formed by functional
J = ∫− ∞+ ∞ (ku + mv)dx
where k, m are arbitrary constants, produces one-parameter family of functions
J(z) = e^{zLE}J = J + zLEJ
+ ½(zLE)^{2}J + ⋯
(actually this is the orbit of non-Noether symmetry group that passes centre of Poisson algebra).
It is interesting that the functionals (LE)^{m}J are in involution.

Theorem.
The orbit (15) of the non-Noether symmetry group
generated by the vector field (12) is involutive
{J(x) , J(y)} = 0 ∀x, y ∈ ℝ
and the functionals
J^{(m)} = (LE)^{m}J
form Lenard scheme with respect to bi-Hamiltonian structure (5)
and produce involutive sequence of conservation laws of the modified Boussinesq hierarchy.

Proof.
The theorem follows from simple geometric properties of the vector field E.
In particular taking the Lie derivative of Poisson bivector field
W along E one gets the second Poisson bivector involved in bi-Hamiltonian system (5)
Ŵ = [E , W]
while the Lie derivative of Ŵ along E vanishes [E , Ŵ] = 0
These properties ensure that the functionals (17) are in involution
(the Poisson bracket of arbitrary two conservation laws from infinite family (17) vanishes)
{J^{(k)} , J^{(m)}} = 0 k, m = 0, 1, 2 ...
Indeed, by applying m-th order Lie derivative (LE)^{m} to the relation
W(dJ^{(0)}) = 0
which reflects the fact that J^{(0)} belongs to the centre of Poisson algebra,
its easy to prove that the functionals (17) form Lenard scheme
W(dJ^{(m + 1)}) = − (1 + m)[E , W](dJ^{(m)})
with respect to bi-Hamiltonian system (5)
From the other hand it is well known [4] that functionals involved in Lenard scheme are
in involution. In the same time calculating the functional
J^{(2)} = (LE)^{2}J^{(0)} = m∫− ∞+ ∞ (u^{2}v + v^{3} + 2cuvx)dx = 2mH
gives rise to Hamiltonian of the modified Boussinesq system and
functionals J^{(m)} being in involution with Hamiltonian must be conservation laws.
∎

By calculating Lie derivatives of J^{(0)} along the vector field E one can
get explicit form of the conservation laws of the modified Boussinesq system:
J^{(0)} = ∫− ∞+ ∞ (ku + mv)dx
J^{(1)} = LEJ^{(0)} = m2∫− ∞+ ∞(u^{2} + v^{2})dx
J^{(2)} = (LE)^{2}J^{(0)} = m∫− ∞+ ∞ (u^{2}v + v^{3} + 2cuvx)dx
J^{(3)} = (LE)^{3}J^{(0)} =
3m4∫− ∞+ ∞ (u^{4} + 5v^{4} + 6u^{2}v^{2}
− 12cv^{2}ux + 4c^{2}ux^{2} + 4c^{2}vx^{2})dx
J^{(m)} = (LE)^{m}J^{(0)} = LEJ^{(m − 1)}

Summary.
The fact that the infinite sequence of conservation laws of modified Boussinesq hierarchy
form single orbit of the one-parameter non-Noether symmetry group indicates that
non-Noether symmetries may play an important role in analysis of certain integrable
models where they drastically simplify calculation of conservation laws and shed more
light on geometric origin of integrable hierarchies. Basic results of the paper can be extended
to the case of periodic boundary conditions u(− ∞) = u(+ ∞) and v(− ∞) = v(+ ∞)
when the modified Boussinesq equations can be considered as bi-Hamiltonian system on a loop space
[4]. Note however that in the periodic case the symmetry (12) does not seem to
preserve boundary conditions.