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 or
der (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.
Predictions of analytical models for hadronisation, namely the dispersive model, the shape function and the single dressed gluon approximation, are compared with moments of hadronic event shape distributions measured in epem annihilation at centre-of
-mass energies between 14 and 209 GeV. In contrast to Monte Carlo models for hadronisation, analytical models require to adjust only two universal parameters, the strong coupling and a second quantity parametrising nonperturbative corrections. The extracted values of as are consistent with the world average and competitive with previous measurements. The variance of event shape distributions is compared with predictions given by some of these models. Limitations of the models, probably due to unknown higher order corrections, are demonstrated and discussed.
