Geostatistical radar-raingauge combination with nonparametric correlograms: methodological considerations and application in Switzerland

Modelling spatial covariance is an essential part of all geostatistical methods. Traditionally, parametric semivariogram models are fit from available data. More recently, it has been suggested to use nonparametric correlograms obtained from spatially complete data fields. Here, both estimation techniques are compared. Nonparametric correlograms are shown to have a substantial negative bias. Nonetheless, when combined with the sample variance of the spatial field under consideration, they yield an estimate of the semivariogram that is unbiased for small lag distances. This justifies the use of this estimation technique in geostatistical applications. Various formulations of geostatistical combination (Kriging) methods are used here for the construction of hourly precipitation grids for Switzerland based on data from a sparse realtime network of raingauges and from a spatially complete radar composite. Two variants of Ordinary Kriging (OK) are used to interpolate the sparse gauge observations. In both OK variants, the radar data are only used to determine the semivariogram model. One variant relies on a traditional parametric semivariogram estimate, whereas the other variant uses the nonparametric correlogram. The variants are tested for three cases and the impact of the semivariogram model on the Kriging prediction is illustrated. For the three test cases, the method using nonparametric correlograms performs equally well or better than the traditional method, and at the same time offers great practical advantages. Furthermore, two variants of Kriging with external drift (KED) are tested, both of which use the radar data to estimate nonparametric correlograms, and as the external drift variable. The first KED variant has been used previously for geostatistical radar-raingauge merging in Catalonia (Spain). The second variant is newly proposed here and is an extension of the first. Both variants are evaluated for the three test cases as well as an extended evaluation period. It is found that both methods yield merged fields of better quality than the original radar field or fields obtained by OK of gauge data. The newly suggested KED formulation is shown to be beneficial, in particular in mountainous regions where the quality of the Swiss radar composite is comparatively low. An analysis of the Kriging variances shows that none of the methods tested here provides a satisfactory uncertainty estimate. A suitable variable transformation is expected to improve this.

Reduced space optimal interpolation of daily rain gauge precipitation in Switzerland

The coarse spacing of automatic rain gauges complicates near-real- time spatial analyses of precipitation. We test the possibility of improving such analyses by considering, in addition to the in situ measurements, the spatial covariance structure inferred from past observations with a denser network. To this end, a statistical reconstruction technique, reduced space optimal interpolation (RSOI), is applied over Switzerland, a region of complex topography. RSOI consists of two main parts. First, principal component analysis (PCA) is applied to obtain a reduced space representation of gridded high- resolution precipitation fields available for a multiyear calibration period in the past. Second, sparse real-time rain gauge observations are used to estimate the principal component scores and to reconstruct the precipitation field. In this way, climatological information at higher resolution than the near-real-time measurements is incorporated into the spatial analysis. PCA is found to efficiently reduce the dimensionality of the calibration fields, and RSOI is successful despite the difficulties associated with the statistical distribution of daily precipitation (skewness, dry days). Examples and a systematic evaluation show substantial added value over a simple interpolation technique that uses near-real-time observations only. The benefit is particularly strong for larger- scale precipitation and prominent topographic effects. Small-scale precipitation features are reconstructed at a skill comparable to that of the simple technique. Stratifying the reconstruction method by the types of weather type classifications yields little added skill. Apart from application in near real time, RSOI may also be valuable for enhancing instrumental precipitation analyses for the historic past when direct observations were sparse.

Ground-based aerosol optical depth trends at three high-altitude sites in Switzerland and Southern Germany from 1995–2010

Ground-based aerosol optical depth (AOD) climatologies at three high-altitude sites in Switzerland (Jungfraujoch and Davos) and Southern Germany (Hohenpeissenberg) are updated and re-calibrated for the period 1995 – 2010. In addition, AOD time-series are augmented with previously unreported data, and are homogenized for the first time. Trend analysis revealed weak AOD trends (λ = 500 nm) at Jungfraujoch (JFJ; +0.007 decade-1), Davos (DAV; +0.002 decade-1) and Hohenpeissenberg (HPB; -0.011 decade-1) where the JFJ and HPB trends were statistically significant at the 95% and 90% confidence levels. However, a linear trend for the JFJ 1995 – 2005 period was found to be more appropriate than for 1995 – 2010 due to the influence of stratospheric AOD which gave a trend -0.003 decade-1 (significant at 95% level). When correcting for a recently available stratospheric AOD time-series, accounting for Pinatubo (1991) and more recent volcanic eruptions, the 1995 – 2010 AOD trends decreased slightly at DAV and HPB but remained weak at +0.000 decade-1 and -0.013 decade-1 (significant at 95% level). The JFJ 1995 – 2005 AOD time-series similarly decreased to -0.003 decade-1 (significant at 95% level). We conclude that despite a more detailed re40 analysis of these three time-series, which have been extended by five years to the end of 2010, a significant decrease in AOD at these three high-altitude sites has still not been observed.

Impact of US Brown Swiss genetics on milk quality from low-input herds in Switzerland: interactions with grazing intake and pasture type

This study investigated the effect of, and interactions between, contrasting crossbreed genetics (US Brown Swiss [BS] × Improved Braunvieh [BV] × Original Braunvieh [OB]) and feeding regimes (especially grazing intake and pasture type) on milk fatty acid (FA) profiles. Concentrations of total polyunsaturated FAs, total omega-3 FAs and trans palmitoleic, vaccenic, α-linolenic, eicosapentaenoic and docosapentaenoic acids were higher in cows with a low proportion of BS genetics. Highest concentrations of the nutritionally desirable FAs, trans palmitoleic, vaccenic and eicosapentaenoic acids were found for cows with a low proportion of BS genetics (0-24% and/or 25-49%) on high grazing intake (75-100% of dry matter intake) diets. Multivariate analysis indicated that the proportion of OB genetics is a positive driver for nutritionally desirable monounsaturated and polyunsaturated FAs while BS genetics proportion was positive driver for total and undesirable individual saturated FAs. Significant genetics × feeding regime interactions were also detected for a range of FAs.

Triple isotope (δD, δ17O, δ18O) study on precipitation, drip water and speleothem fluid inclusions for a central European cave (NW Switzerland)

Deuterium (dD) and oxygen (d18O) isotopes are powerful tracers of the hydrological cycle and have been extensively used for paleoclimate reconstructions as they can provide information on past precipitation, temperature and atmospheric circulation. More recently, the use of 17Oexcess derived from precise measurement of d17O and d18O gives new and additional insights in tracing the hydrological cycle whereas uncertainties surround this proxy. However, 17Oexcess could provide additional information on the atmospheric conditions at the moisture source as well as about fractionations associated with transport and site processes. In this paper we trace water stable isotopes (dD, d17O and d18O) along their path from precipitation to cave drip water and finally to speleothem fluid inclusions for Milandre cave in northwestern Switzerland. A two year-long daily resolved precipitation isotope record close to the cave site is compared to collected cave drip water (3 months average resolution) and fluid inclusions of modern and Holocene stalagmites. Amount weighted mean dD, d18O and d17O are �71.0‰, �9.9‰, �5.2‰ for precipitation, �60.3‰, �8.7‰, �4.6‰ for cave drip water and �61.3‰, �8.3‰, �4.7‰ for recent fluid inclusions respectively. Second order parameters have also been derived in precipitation and drip water and present similar values with 18 per meg for 17Oexcess whereas d-excess is 1.5‰ more negative in drip water. Furthermore, the atmospheric signal is shifted towards enriched values in the drip water and fluid inclusions (D of ~ þ 10‰ for dD). The isotopic composition of cave drip water exhibits a weak seasonal signal which is shifted by around 8e10 months (groundwater residence time) when compared to the precipitation. Moreover, we carried out the first d17O measurement in speleothem fluid inclusions, as well as the first comparison of the d17O behaviour from the meteoric water to the fluid inclusions entrapment in speleothems. This study on precipitation, drip water and fluid inclusions will be used as a speleothem proxy calibration for Milandre cave in order to reconstruct paleotemperatures and moisture source variations for Western Central Europe.

Timing and nature of the penultimate deglaciation in a high alpine stalagmite from Switzerland

The timing and nature of the penultimate deglaciation, also known as Termination II (T-II), is subject of controversial discussions due to the scarcity of precisely-dated palaeoclimate records. Here we present a new precisely-dated and highly-resolved multi-proxy stalagmite record covering T-II from the high alpine Schafsloch Cave in Switzerland, an area where climate is governed by the North Atlantic. The inception of stalagmite growth at 137.4 ± 1.4 kyr before present (BP) indicates the presence of drip water and cave air temperatures of above 0°C, and is related to a climate-induced change in the thermal state (from cold- to warm-based) of the glacier above the cave. The cessation of stalagmite growth between 133.1 ± 0.7 and 131.9 ± 0.6 kyr BP is most likely related to distinct drop in temperature associated with Heinrich stadial 11. The resumption of stalagmite growth at 131.9 ± 0.6 kyr BP is accompanied by an abrupt increase in temperature and precipitation as indicated by distinct shifts in the oxygen and carbon isotopic composition as well as in trace element concentrations. The mid-point of T-II is around 131.8 ± 0.6 kyr BP in the Schafsloch Cave record is significantly earlier compared to the age of 129.1 ± 0.1 kyr BP in the Sanbao Cave record from China. The different ages between both records can be best explained by the competing effects of insolation and glacial boundary forcing on seasonality and snow cover extent in Eurasia.