  Event Shape Data from $e^+e^-$ annihilation into hadrons collected by the
JADE experiment at centre-of-mass energies between 14 GeV and 44 GeV are used
to determine the strong coupling $\alpha_S$. QCD predictions complete to
next-to-next-to-leading order (NNLO), alternatively combined with resummed
next-to-leading-log-approximation (NNLO+NLLA) calculations, are used. The
combined value from six different event shape observables at the six JADE
centre-of-mass energies using the NNLO calculations is
  $\alpha_S(M_Z)$= 0.1210 +/- 0.0007(stat.) +/- 0.0021(expt.) +/- 0.0044(had.)
+/- 0.0036(theo.) and with the NNLO+NLLA calculations the combined value is
$\alpha_S$= 0.1172 +/- 0.0006(stat.) +/- 0.0020(expt.) +/- 0.0035(had.) +/-
0.0030(theo.) . The stability of the NNLO and NNLO+NLLA results with respect to
missing higher order contributions, studied by variations of the
renormalisation scale, is improved compared to previous results obtained with
NLO+NLLA or with NLO predictions only. The observed energy dependence of
$\alpha_S$ agrees with the QCD prediction of asymptotic freedom and excludes
absence of running with 99% confidence level.
