Attributing trends in extremely hot days to changes in atmospheric dynamics

This paper presents a method for attributing regional trends in the frequency of extremely hot days (EHDs) to changes in the frequency of the atmospheric patterns that characterize such extraordinary events. The study is applied to mainland Spain and the Balearic Islands for the extended summers of the period 1958–2008, where significant and positive trends in maximum temperature (Tx) have been reported during the second half of the past century. First, the study area was split into eight regions attending to their different temporal variability of the daily Tx series obtained from the Spain02 gridded data set using a clustering procedure. Second, the large-scale atmospheric situations causing EHDs are defined by circulation types (CTs). The obtainment of the CTs differs from the majority of CT classifications proposed in the literature. It is based on regional series and on a previous characterization of the main atmospheric situations obtained using only some days classified as extremes in the different regions. Three different atmospheric fields (SLP, T850, and Z500) from ECMWF reanalysis and analysis data and combinations of them (SLP–T850, SLP–Z500, and T850–Z500) are used to produce six different CT classifications. Subsequently, links between EHD occurrence in the different regions and CT for all days have been established. Finally, a simple model to relate the trends in EHDsfor each region to the changes in the CT frequency appearance has been formulated. Most regions present positive and significant trends in the occurrence of EHDs. The CT classifications using two variables perform better. In particular, SLP–T850 is the best for characterizing the atmospheric situations leading to EHD occurrences for most of the regions. Only a small number of CTs have significant trends in their frequency and are associated with high efficiency causing EHD occurrences in most regions simultaneously, especially in the northern and central regions. Attribution results show that changes in circulation can only explain some part of the regional EHD trends. The percentage of the trend attributable to changes in atmospheric dynamics varies from 15 to 50 %, depends on the region and is sensitive to the selected large-scale variables.

Sensitivity of the WRF model to PBL parametrisations and nesting techniques: evaluation of wind storms over complex terrain

Simulating surface wind over complex terrain is a challenge in regional climate modelling. Therefore, this study aims at identifying a set-up of the Weather Research and Forecasting Model (WRF) model that minimises system- atic errors of surface winds in hindcast simulations. Major factors of the model configuration are tested to find a suitable set-up: the horizontal resolution, the planetary boundary layer (PBL) parameterisation scheme and the way the WRF is nested to the driving data set. Hence, a number of sensitivity simulations at a spatial resolution of 2 km are carried out and compared to observations. Given the importance of wind storms, the analysis is based on case studies of 24 historical wind storms that caused great economic damage in Switzerland. Each of these events is downscaled using eight different model set-ups, but sharing the same driving data set. The results show that the lack of representation of the unresolved topography leads to a general overestimation of wind speed in WRF. However, this bias can be substantially reduced by using a PBL scheme that explicitly considers the effects of non-resolved topography, which also improves the spatial structure of wind speed over Switzerland. The wind direction, although generally well reproduced, is not very sensitive to the PBL scheme. Further sensitivity tests include four types of nesting methods: nesting only at the boundaries of the outermost domain, analysis nudging, spectral nudging, and the so-called re-forecast method, where the simulation is frequently restarted. These simulations show that restricting the freedom of the model to develop large-scale disturbances slightly increases the temporal agreement with the observations, at the same time that it further reduces the overestimation of wind speed, especially for maximum wind peaks. The model performance is also evaluated in the outermost domains, where the resolution is coarser. The results demonstrate the important role of horizontal resolution, where the step from 6 to 2 km significantly improves model performance. In summary, the combination of a grid size of 2 km, the non-local PBL scheme modified to explicitly account for non-resolved orography, as well as analysis or spectral nudging, is a superior combination when dynamical downscaling is aimed at reproducing real wind fields.

Glacial–interglacial temperature change in the tropical West Pacific: A comparison of stalagmite-based paleo-thermometers

In the tropics, geochemical records from stalagmites have so far mainly been used to qualitatively reconstruct changes in precipitation, but several new methods to reconstruct past temperatures from stalagmite material have emerged recently: i) liquide vapor homogenization of fluid inclusion water ii) noble gas concentrations in fluid inclusion water, iii) the partitioning of oxygen isotopes between fluid inclusion water and calcite, and iv) the abundance of the 13C18O16O(‘clumped’) isotopologue in calcite. We present, for the first time, a direct comparison of these four paleo-thermometers by applying them to a fossil stalagmite covering nearly two glaciale interglacial cycles (Marine Isotope Stages (MIS) 12 e 9) and to two modern stalagmites, all from northern Borneo. The temperature estimates from the different methods agree in most cases within errors for both the old and recent samples; reconstructed formation temperatures of the recent samples match within 2-sigma errors with measured cave temperatures. However, slight but systematic deviations are observed between noble gas and liquide vapor homogenization temperatures. Whereas the temperature sensitivity of fluid inclusion d18O and clumped isotopes is currently debated, we find that the calibration of Tremaine et al. (2011) for fluid inclusion d18O and a synthetic calcite-based clumped isotope calibration (Ziegler et al., in prep.) yield temperature estimates consistent with the other methods. All methods (with the potential exception of clumped isotopes) show excellent agreement on the amplitude of glaciale interglacial temperature change, indicating temperature shifts of 4-5 C°. This amplitude is similar to the amplitude of Mg/Ca-based regional sea surface temperature records, when correcting for sea level driven changes in cave elevation. Our reconstruction of tropical temperature evolution over the time period from 440 to 320 thousand years ago (ka) adds support to the view that climate sensitivity to varying greenhouse forcing is substantial also in the deep tropics.

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

Deuterium (δD) and oxygen (δ18O) 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 δ17O excess derived from precise measurement of δ17O and δ18O gives new and additional insights in tracing the hydrological cycle whereas uncertainties surround this proxy. However, 17O excess 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 (δD,δ17O and δ18O) 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 δD,δ18O and δ17O 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 17O excess 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 (Δ of ~ + 10‰ for δD). The isotopic composition of cave drip water exhibits a weak seasonal signal which is shifted by around 8 - 10 months (groundwater residence time) when compared to the precipitation. Moreover, we carried out the first δ17O measurement in speleothem fluid inclusions, as well as the first comparison of the δ17 O 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.

Quantification of water content across a cement-clay interface using high resolution neutron radiography

In many designs for radioactive waste repositories, cement and clay will come into direct contact. The geochemical contrast between cement and clay will lead to mass fluxes across the interface, which consequently results in alteration of structural and transport properties of both materials that may affect the performance of the multi-barrier system. We present an experimental approach to study cement-clay interactions with a cell to accommodate small samples of cement and clay. The cell design allows both in situ measurement of water content across the sample using neutron radiography and measurement of transport parameters using through-diffusion tracer experiments. The aim of the high- resolution neutron radiography experiments was to monitor changes in water content (porosity) and their spatial extent. Neutron radiographs of several evolving cement-clay interfaces delivered quantitative data which allow resolving local water contents within the sample domain. In the present work we explored the uncertainties of the derived water contents with regard to various input parameters and with regard to the applied image correction procedures. Temporal variation of measurement conditions created absolute uncertainty of the water content in the order of ±0.1 (m3/m3), which could not be fully accounted for by correction procedures. Smaller relative changes in water content between two images can be derived by specific calibrations to two sample regions with different, invariant water contents.