Temporal trends in δ¹⁸O composition of precipitation in Germany: Insights from time series modelling and trend analysis.

Temporal and spatial variations of stable oxygen (18O) and hydrogen (2H) isotope measurements in precipitation act as important proxies for changing hydro-meteorological and regional and global climate patterns. Temporal trends in time series of the stable isotope composition in precipitation were rarely observed, and they are poorly understood. These might be a result of a lack of proper trend detection tools and effort for exploring trend processes. Here, we investigate temporal trends of δ18O in precipitation at 17 observation stations in Germany between 1978 and 2009. We test if significant trends in the isotope time series from different models can be observed. Mann-Kendall trend tests are applied on the isotope series, using general multiplicative seasonal autoregressive integrate moving average (ARIMA) models, which account for first and higher order serial correlations. Effects of temperature, precipitation, and geographic parameters on isotope trends are also investigated in the proposed models. To benchmark our proposed approach, the ARIMA results are compared with a trend-free pre-whitening procedure, the state of the art method for removing the first order autocorrelation in environmental trend studies. Moreover, we further explore whether higher order serial correlations in isotope series affects our trend results. Overall, three out of the 17 stations show significant changes when higher order autocorrelation are adjusted, and four show a significant trend when temperature and precipitation effects are considered. The significant trends in the isotope time series generally occur only at low elevation stations. Higher order autoregressive processes are shown to be important in the isotope time series analysis. Results suggest that the widely used trend analysis with only the first order autocorrelation adjustment may not adequately take account of the high order autocorrelated processes in the stable isotope series. The investigated time series analysis method including higher autocorrelation and external climate variable adjustments is shown to be a better alternative