Stacking of discontinuous regional palaeoclimate records: Chironomid-based summer temperatures from the Alpine region

Since multi-site reconstructions are less affected by site-specific climatic effects and artefacts, regional palaeotemperature reconstructions based on a number of sites can provide more robust estimates of centennial- to millennial-scale temperature trends than individual, site-specific records. Furthermore, reconstructions based on multiple records are necessary for developing continuous climate records over time scales longer than covered by individual sequences. Here, we present a procedure for developing such reconstructions based on relatively short (centuries to millennia), discontinuously sampled records as are typically developed when using biotic proxies in lake sediments for temperature reconstruction. The approach includes an altitudinal correction of temperatures, an interpolation of individual records to equal time intervals, a stacking procedure for sections of the interval of interest that have the same records available, as well as a splicing procedure to link the individual stacked records into a continuous reconstruction. Variations in the final, stacked and spliced reconstruction are driven by variations in the individual records, whereas the absolute temperature values are determined by the stacked segment based on the largest number of records. With numerical simulations based on the NGRIP δ18O record, we demonstrate that the interpolation and stacking procedure provides an approximation of a smoothed palaeoclimate record if based on a sufficient number of discontinuously sampled records. Finally, we provide an example of a stacked and spliced palaeotemperature reconstruction 15000–90 calibrated 14C yr BP based on six chironomid records from the northern and central Swiss Alps and eastern France to discuss the potential and limitations of this approach.

Alpine lacustrine varved record reveals summer temperature as main control of glacier fluctuations over the past 2250 years

Glacier fluctuations are a key indicator of changing climate. Their reconstruction beyond historical times unravels glacier variability and its forcing factors on long time scales, which can considerably improve our understanding of the climate–glacier relationship. Here, we present a 2250-year-long reconstruction of particle-mass accumulation rates recorded in the lacustrine sediments of Lake Trüebsee (Central Swiss Alps) that are directly related to glacier extent, thus reflecting a continuous record of fluctuations of the upstream-located Titlis Glacier. Mass accumulation rate values show strong centennial to multi-centennial fluctuations and reveal 12 well-pronounced periods of enhanced values corresponding to times of maximum extent of the neighboring Lower Grindelwald Glacier. This result supports previous studies of proglacial lake sediments that documented high mass accumulation rate values during glacier advances. The strong variability in the Lake Trüebsee mass accumulation rate record thus represents a highly sensitive paleoclimatic archive, which mirrors rapid and pronounced feedbacks of Titlis Glacier to climatic changes over the past 2250years. The comparison of our data with independent paleo-temperature reconstructions from tree rings suggests that variations in mean summer temperature were the primary driving factor of fluctuations of Titlis Glacier. Also, advances of Titlis Glacier occurred during the grand solar minima (Dalton, Maunder, Spörer, Wolf) of the last millennium. This relation of glacier extent with summer temperature reveals strong evidence that the mass balance of this Alpine glacier is primarily controlled by the intensity of glacier melting during summer.

The role of human-induced fire and sweet chestnut (Castanea sativa Mill.) cultivation on the long-term landscape dynamics of the southern Swiss Alps

Changes in fire occurrence during the last decades in the southern Swiss Alps make knowledge on fire history essential to understand future evolution of the ecosystem composition and functioning. In this context, palaeoecology provides useful insights into processes operating at decadal-to-millennial time scales, such as the response of plant communities to intensified fire disturbances during periods of cultural change. We provide a high-resolution macroscopic charcoal and pollen series from Guèr, a well-dated peat sequence at mid-elevation (832 m.a.s.l.) in southern Switzerland, where the presence of local settlements is documented since the late Bronze Age and the Iron Age. Quantitative fire reconstruction shows that fire activity sharply increased from the Neolithic period (1–3 episodes/1000 year) to the late Bronze and Iron Age (7–9 episodes/1000 year), leading to extensive clearance of the former mixed deciduous forest (Alnus glutinosa, Betula, deciduous Quercus). The increase in anthropogenic pollen indicators (e.g. Cerealia-type, Plantago lanceolata) together with macroscopic charcoal suggests anthropogenic rather than climatic forcing as the main cause of the observed vegetation shift. Fire and controlled burning were extensively used during the late Roman Times and early Middle Ages to promote the introduction and establishment of chestnut (Castanea sativa) stands, which provided an important wood and food supply. Fire occurrence declined markedly (from 9 to 5–6 episodes/1000 year) during late Middle Ages because of fire suppression, biomass removal by human population, and landscape fragmentation. Land-abandonment during the last decades allowed forest to partly re-expand (mainly Alnus glutinosa, Betula) and fire frequency to increase.

Establishing the skill of climate field reconstruction techniques for precipitation with pseudoproxy experiments

This study aims at assessing the skill of several climate field reconstruction techniques (CFR) to reconstruct past precipitation over continental Europe and the Mediterranean at seasonal time scales over the last two millennia from proxy records. A number of pseudoproxy experiments are performed within the virtual reality ofa regional paleoclimate simulation at 45 km resolution to analyse different aspects of reconstruction skill. Canonical Correlation Analysis (CCA), two versions of an Analog Method (AM) and Bayesian hierarchical modeling (BHM) are applied to reconstruct precipitation from a synthetic network of pseudoproxies that are contaminated with various types of noise. The skill of the derived reconstructions is assessed through comparison with precipitation simulated by the regional climate model. Unlike BHM, CCA systematically underestimates the variance. The AM can be adjusted to overcome this shortcoming, presenting an intermediate behaviour between the two aforementioned techniques. However, a trade-off between reconstruction-target correlations and reconstructed variance is the drawback of all CFR techniques. CCA (BHM) presents the largest (lowest) skill in preserving the temporal evolution, whereas the AM can be tuned to reproduce better correlation at the expense of losing variance. While BHM has been shown to perform well for temperatures, it relies heavily on prescribed spatial correlation lengths. While this assumption is valid for temperature, it is hardly warranted for precipitation. In general, none of the methods outperforms the other. All experiments agree that a dense and regularly distributed proxy network is required to reconstruct precipitation accurately, reflecting its high spatial and temporal variability. This is especially true in summer, when a specifically short de-correlation distance from the proxy location is caused by localised summertime convective precipitation events.

Can sympatric speciation by disruptive sexual selection explain rapid evolution of cichlid diversity in Lake Victoria?

Rapid speciation can occur on ecological time scales and interfere with ecological processes, resulting in species distribution patterns that are difficult to reconcile with ecological theory. The haplochromine cichlids in East African lakes are an extreme example of rapid speciation. We analyse the causes of their high speciation rates. Various studies have identified disruptive sexual selection acting on colour polymorphisms that might cause sympatric speciation. Using data on geographical distribution, colouration and relatedness from 41 species endemic to Lake Victoria, we test predictions from this hypothesis. Plotting numbers of pairs of closely related species against the amount of distributional overlap between the species reveals a bimodal distribution with modes on allopatric and sympatric. The proportion of sister species pairs that are heteromorphic for the traits under disruptive selection is higher in sympatry than in allopatry. These data support the hypothesis that disruptive sexual selection on colour polymorphisms has caused sympatric speciation and help to explain the rapid radiation of haplochromine species flocks.

Giant Electric Field Enhancement in Split Ring Resonators Featuring Nanometer-Sized Gaps

Today's pulsed THz sources enable us to excite, probe, and coherently control the vibrational or rotational dynamics of organic and inorganic materials on ultrafast time scales. Driven by standard laser sources THz electric field strengths of up to several MVm−1 have been reported and in order to reach even higher electric field strengths the use of dedicated electric field enhancement structures has been proposed. Here, we demonstrate resonant electric field enhancement structures, which concentrate the incident electric field in sub-diffraction size volumes and show an electric field enhancement as high as ~14,000 at 50 GHz. These values have been confirmed through a combination of near-field imaging experiments and electromagnetic simulations.

Solvation-Driven Charge Transfer and Localization in Metal Complexes

In any physicochemical process in liquids, the dynamical response of the solvent to the solutes out of equilibrium plays a crucial role in the rates and products: the solvent molecules react to the changes in volume and electron density of the solutes to minimize the free energy of the solution, thus modulating the activation barriers and stabilizing (or destabilizing) intermediate states. In charge transfer (CT) processes in polar solvents, the response of the solvent always assists the formation of charge separation states by stabilizing the energy of the localized charges. A deep understanding of the solvation mechanisms and time scales is therefore essential for a correct description of any photochemical process in dense phase and for designing molecular devices based on photosensitizers with CT excited states. In the last two decades, with the advent of ultrafast time-resolved spectroscopies, microscopic models describing the relevant case of polar solvation (where both the solvent and the solute molecules have a permanent electric dipole and the mutual interaction is mainly dipole−dipole) have dramatically progressed. Regardless of the details of each model, they all assume that the effect of the electrostatic fields of the solvent molecules on the internal electronic dynamics of the solute are perturbative and that the solvent−solute coupling is mainly an electrostatic interaction between the constant permanent dipoles of the solute and the solvent molecules. This well-established picture has proven to quantitatively rationalize spectroscopic effects of environmental and electric dynamics (time-resolved Stokes shifts, inhomogeneous broadening, etc.). However, recent computational and experimental studies, including ours, have shown that further improvement is required. Indeed, in the last years we investigated several molecular complexes exhibiting photoexcited CT states, and we found that the current description of the formation and stabilization of CT states in an important group of molecules such as transition metal complexes is inaccurate. In particular, we proved that the solvent molecules are not just spectators of intramolecular electron density redistribution but significantly modulate it. Our results solicit further development of quantum mechanics computational methods to treat the solute and (at least) the closest solvent molecules including the nonperturbative treatment of the effects of local electrostatics and direct solvent−solute interactions to describe the dynamical changes of the solute excited states during the solvent response.

Comparison of co-located independent ground-based middle atmospheric wind and temperature measurements with numerical weather prediction models

High-resolution, ground-based and independent observations including co-located wind radiometer, lidar stations, and infrasound instruments are used to evaluate the accuracy of general circulation models and data-constrained assimilation systems in the middle atmosphere at northern hemisphere midlatitudes. Systematic comparisons between observations, the European Centre for Medium-Range Weather Forecasts (ECMWF) operational analyses including the recent Integrated Forecast System cycles 38r1 and 38r2, the NASA’s Modern-Era Retrospective Analysis for Research and Applications (MERRA) reanalyses, and the free-running climate Max Planck Institute–Earth System Model–Low Resolution (MPI-ESM-LR) are carried out in both temporal and spectral dom ains. We find that ECMWF and MERRA are broadly consistent with lidar and wind radiometer measurements up to ~40 km. For both temperature and horizontal wind components, deviations increase with altitude as the assimilated observations become sparser. Between 40 and 60 km altitude, the standard deviation of the mean difference exceeds 5 K for the temperature and 20 m/s for the zonal wind. The largest deviations are observed in winter when the variability from large-scale planetary waves dominates. Between lidar data and MPI-ESM-LR, there is an overall agreement in spectral amplitude down to 15–20 days. At shorter time scales, the variability is lacking in the model by ~10 dB. Infrasound observations indicate a general good agreement with ECWMF wind and temperature products. As such, this study demonstrates the potential of the infrastructure of the Atmospheric Dynamics Research Infrastructure in Europe project that integrates various measurements and provides a quantitative understanding of stratosphere-troposphere dynamical coupling for numerical weather prediction applications.

European summer temperatures since Roman times

The spatial context is critical when assessing present-day climate anomalies, attributing them to potential forcings and making statements regarding their frequency and severity in a long-term perspective. Recent international initiatives have expanded the number of high-quality proxy-records and developed new statistical reconstruction methods. These advances allow more rigorous regional past temperature reconstructions and, in turn, the possibility of evaluating climate models on policy-relevant, spatio-temporal scales. Here we provide a new proxy-based, annually-resolved, spatial reconstruction of the European summer (June–August) temperature fields back to 755 CE based on Bayesian hierarchical modelling (BHM), together with estimates of the European mean temperature variation since 138 BCE based on BHM and composite-plus-scaling (CPS). Our reconstructions compare well with independent instrumental and proxy-based temperature estimates, but suggest a larger amplitude in summer temperature variability than previously reported. Both CPS and BHM reconstructions indicate that the mean 20th century European summer temperature was not significantly different from some earlier centuries, including the 1st, 2nd, 8th and 10th centuries CE. The 1st century (in BHM also the 10th century) may even have been slightly warmer than the 20th century, but the difference is not statistically significant. Comparing each 50 yr period with the 1951–2000 period reveals a similar pattern. Recent summers, however, have been unusually warm in the context of the last two millennia and there are no 30 yr periods in either reconstruction that exceed the mean average European summer temperature of the last 3 decades (1986–2015 CE). A comparison with an ensemble of climate model simulations suggests that the reconstructed European summer temperature variability over the period 850–2000 CE reflects changes in both internal variability and external forcing on multi-decadal time-scales. For pan-European temperatures we find slightly better agreement between the reconstruction and the model simulations with high-end estimates for total solar irradiance. Temperature differences between the medieval period, the recent period and the Little Ice Age are larger in the reconstructions than the simulations. This may indicate inflated variability of the reconstructions, a lack of sensitivity and processes to changes in external forcing on the simulated European climate and/or an underestimation of internal variability on centennial and longer time scales.

Computational fluid dynamics modelling in cardiovascular medicine.

This paper reviews the methods, benefits and challenges associated with the adoption and translation of computational fluid dynamics (CFD) modelling within cardiovascular medicine. CFD, a specialist area of mathematics and a branch of fluid mechanics, is used routinely in a diverse range of safety-critical engineering systems, which increasingly is being applied to the cardiovascular system. By facilitating rapid, economical, low-risk prototyping, CFD modelling has already revolutionised research and development of devices such as stents, valve prostheses, and ventricular assist devices. Combined with cardiovascular imaging, CFD simulation enables detailed characterisation of complex physiological pressure and flow fields and the computation of metrics which cannot be directly measured, for example, wall shear stress. CFD models are now being translated into clinical tools for physicians to use across the spectrum of coronary, valvular, congenital, myocardial and peripheral vascular diseases. CFD modelling is apposite for minimally-invasive patient assessment. Patient-specific (incorporating data unique to the individual) and multi-scale (combining models of different length- and time-scales) modelling enables individualised risk prediction and virtual treatment planning. This represents a significant departure from traditional dependence upon registry-based, population-averaged data. Model integration is progressively moving towards 'digital patient' or 'virtual physiological human' representations. When combined with population-scale numerical models, these models have the potential to reduce the cost, time and risk associated with clinical trials. The adoption of CFD modelling signals a new era in cardiovascular medicine. While potentially highly beneficial, a number of academic and commercial groups are addressing the associated methodological, regulatory, education- and service-related challenges.

Dynamics of dissipative Bose-Einstein condensation

We resolve the real-time dynamics of a purely dissipative s=1/2 quantum spin or, equivalently, hard-core boson model on a hypercubic d-dimensional lattice. The considered quantum dissipative process drives the system to a totally symmetric macroscopic superposition in each of the S3 sectors. Different characteristic time scales are identified for the dynamics and we determine their finite-size scaling. We introduce the concept of cumulative entanglement distribution to quantify multiparticle entanglement and show that the considered protocol serves as an efficient method to prepare a macroscopically entangled Bose-Einstein condensate.

Biotic responses to rapid climatic changes

Qualitative and quantitative changes in fossil flora and fauna have been used in many studies to infer climatic change. Here we ask a different question: how do flora and fauna respond to climatic changes such as rapid warming or cooling? As an independent proxy for paleotemperature we take the ratio of oxygen isotopes in biogenically precipitated lake marl and in ostracod shells. This introductory paper describes the project design and the five sites on an altitudinal transect from 600 m to about 2300 m asl in the western Swiss Alps. As cases of climatic cooling and warming we use the beginning and end of the Younger Dryas as major changes, and the Gerzensee and Preboreal oscillations as minor changes. At the two sites of Gerzensee and Leysin these changes are recorded in stable-isotope ratios, and there the time scales can be derived by correlations to the GRIP ice core (Schwander et al., 2000 and von Grafenstein et al., 2000). Biotic responses to climate changes are treated in individual papers using pollen (Wick, 2000), plant macrofossils (Tobolski and Ammann, 2000), and remains of chironomids (Brooks, 2000), beetles and other insects (Lemdahl, 2000), and chydorid Cladocera (Hofmann, 2000). They are followed by a synthesis focusing on quantification of biotic responses (Ammann et al., 2000). In addition, a reconstruction of summer temperatures for the Allerød and the Younger Dryas at Gerzensee is provided by Lotter et al. (2000).

Stable isotope response to lake eutrophication: Calibration of a high-resolution lacustrine sequence from Baldeggersee, Switzerland

Stable isotope analyses of discrete seasonal layers from a 108-yr annually laminated freeze-core from Baldeg-gersee, a small, eutrophic lake in central Switzerland, provide information on the climatological and environmental factors, including lake eutrophication, that control oxygen and carbon isotopic composition of epilimnic biologically induced calcite precipitate. During the last 100 yr, Baldeggersee has undergone major increases in productivity and eutrophication in response to nutrient loading from agriculture and industrialization in the lake's watershed. Calibration of the isotopic signal in Baldeggersee to historical limnological data quantitatively links evidence of isotopic depletion in the sedimented calcite to trophic state of the lake. δ18O values from the spring/summer “light” sediment layers steadily diverged to more depleted values in response to historical eutrophication: measured δ18O values were up to 21.5‰ more negative than calculated equilibrium δ18O values. Evidence for 13C depletion in the calcite, relative to equilibrium values, is more difficult to ascertain because of an overall dominance of isotopic enrichment in the dissolved inorganic pool as productivity in Baldeggersee increases. A positive association exists between the degree of oxygen-18 depletion and the calcite crystal size. Thus, large amorphous calcite grains can be used as a proxy for recognizing apparent isotopic nonequilibrium in sediment sequences from highly productive lacustrine environments from all geologic time scales. In contrast to the light layers, the oxygen isotopic composition of the calcite in the late summer/fall “dark” sediment layers is unaffected by the apparent isotope nonequilibrium. Oxygen and carbon isotope values from the dark laminae in the Baldeggersee sediment therefore provide environmental and climatological proxies that can be calibrated with known environmental and regional climate data for the last century.

Chronology for the EDML ice core from Dronning Maud Land, Antarctica, over the last 150 000 years

A chronology called EDML1 has been developed for the EPICA ice core from Dronning Maud Land (EDML). EDML1 is closely interlinked with EDC3, the new chronology for the EPICA ice core from Dome-C (EDC) through a stratigraphic match between EDML and EDC that consists of 322 volcanic match points over the last 128 ka. The EDC3 chronology comprises a glaciological model at EDC, which is constrained and later selectively tuned using primary dating information from EDC as well as from EDML, the latter being transferred using the tight stratigraphic link between the two cores. Finally, EDML1 was built by exporting EDC3 to EDML. For ages younger than 41 ka BP the new synchronized time scale EDML1/EDC3 is based on dated volcanic events and on a match to the Greenlandic ice core chronology GICC05 via 10Be and methane. The internal consistency between EDML1 and EDC3 is estimated to be typically ~6 years and always less than 450 years over the last 128 ka (always less than 130 years over the last 60 ka), which reflects an unprecedented synchrony of time scales. EDML1 ends at 150 ka BP (2417 m depth) because the match between EDML and EDC becomes ambiguous further down. This hints at a complex ice flow history for the deepest 350 m of the EDML ice core.