8 Gs seems like a lot. And position error will indeed maginfy second derivatives. I'd first fit a smoothed flight path through 10 points or so, then take second derivatives of the smoothed function. Just run a regression with a quadratic function for the x y and z coordinates separately +/- say 5 points, and use the coefficient on the quadratic term to get local acceleration.
John Cochrane