Investigation of young water inflow in karst aquifers using SF6–CFC–3H/He–85Kr–39Ar and stable isotope components

Karst aquifers are known for their wide distribution of water transfer velocities. From this observation, a multiple geochemical tracer approach seems to be particularly well suited to provide a significant assessment of groundwater flows, but the choice of adapted tracers is essential. In this study, several common tracers in karst aquifers such as physicochemical parameters, major ions, stable isotopes, and d13C to more specific tracers such as dating tracers – 14C, 3H, 3H–3He, CFC-12, SF6 and 85Kr, and 39Ar – were used, in a fractured karstic carbonated aquifer located in Burgundy (France). The information carried by each tracer and the best sampling strategy are compared on the basis of geochemical monitoring done during several recharge events and over longer time periods (months to years). This study’s results demonstrate that at the seasonal and recharge event time scale, the variability of concentrations is low for most tracers due to the broad spectrum of groundwater mixings. The tracers used traditionally for the study of karst aquifers, i.e., physicochemical parameters and major ions, efficiently describe hydrological processes such as the direct and differed recharge, but require being monitored at short time steps during recharge events to be maximized. From stable isotopes, tritium, and Cl� contents, the proportion of the fast direct recharge by the largest porosity was estimated using a binary mixing model. The use of tracers such as CFC-12, SF6, and 85Kr in karst aquifers provides additional information, notably an estimation of apparent age, but they require good preliminary knowledge of the karst system to interpret the results suitably. The CFC-12 and SF6 methods efficiently determine the apparent age of baseflow, but it is preferable to sample the groundwater during the recharge event. Furthermore, these methods are based on different assumptions such as regional enrichment in atmospheric SF6, excess air, and flow models among others. 85Kr and 39Ar concentrations can potentially provide a more direct estimation of groundwater residence time. Conversely, the 3H–3He method is inefficient in the karst aquifer for dating due to 3He degassing.

Active regulation of cerebral venous tone: simultaneous arterial and venous transcranial Doppler sonography during a Valsalva manoeuvre

The aim of this study was to analyse the cerebral venous outflow in relation to the arterial inflow during a Valsalva manoeuvre (VM). In 19 healthy volunteers (mean age 24.1 +/- 2.6 years), the middle cerebral artery (MCA) and the straight sinus (SRS) were insonated by transcranial Doppler sonography. Simultaneously the arterial blood pressure was recorded using a photoplethysmographic method. Two VM of 10 s length were performed per participant. Tracings of the variables were then transformed to equidistantly re-sampled data. Phases of the VM were analysed regarding the increase of the flow velocities and the latency to the peak. The typical four phases of the VM were also found in the SRS signal. The relative flow velocity (FV) increase was significantly higher in the SRS than in the MCA for all phases, particularly that of phase IV (p < 0.01). Comparison of the time latency of the VM phases of the MCA and SRS only showed a significant difference for phase I (p < 0.01). In particular, there was no significant difference for phase IV (15.8 +/- 0.29 vs. 16.0 +/- 0.28 s). Alterations in venous outflow in phase I are best explained by a cross-sectional change of the lumen of the SRS, while phases II and III are compatible with a Starling resistor. However, the significantly lager venous than the arterial overshoot in phase IV may be explained by the active regulation of the venous tone.

Control of ventricular unloading using an electrocardiogram-synchronized Thoratec paracorporeal ventricular assist device

OBJECTIVE: Current pulsatile ventricular assist devices operate asynchronous with the left ventricle in fixed-rate or fill-to-empty modes because electrocardiogram-triggered modes have been abandoned. We hypothesize that varying the ejection delay in the synchronized mode yields more precise control of hemodynamics and left ventricular loading. This allows for a refined management that may be clinically beneficial. METHODS: Eight sheep received a Thoratec paracorporeal ventricular assist device (Thoratec Corp, Pleasanton, Calif) via ventriculo-aortic cannulation. Left ventricular pressure and volume, aortic pressure, pulmonary flow, pump chamber pressure, and pump inflow and outflow were recorded. The pump was driven by a clinical pneumatic drive unit (Medos Medizintechnik AG, Stolberg, Germany) synchronously with the native R-wave. The start of pump ejection was delayed between 0% and 100% of the cardiac period in 10% increments. For each of these delays, hemodynamic variables were compared with baseline data using paired t tests. RESULTS: The location of the minimum of stroke work was observed at a delay of 10% (soon after aortic valve opening), resulting in a median of 43% reduction in stroke work compared with baseline. Maximum stroke work occurred at a median delay of 70% with a median stroke work increase of 11% above baseline. Left ventricular volume unloading expressed by end-diastolic volume was most pronounced for copulsation (delay 0%). CONCLUSIONS: The timing of pump ejection in synchronized mode yields control over left ventricular energetics and can be a method to achieve gradual reloading of a recoverable left ventricle. The traditionally suggested counterpulsation is not optimal in ventriculo-aortic cannulation when maximum unloading is desired.

Persistent Left Superior Vena Cava in Cardiac Congenital Surgery

Persistent left superior vena cava (LSVC) is a relatively frequent finding in congenital cardiac malformation. The scope of the study was to analyze the timing of diagnosis of persistent LSVC, the timing of diagnosis of associated anomalies of the coronary sinus, and the global impact on morbidity and mortality of persistent LSVC in children with congenital heart disease after cardiac surgery. Retrospective analysis of a cohort of children after cardiac surgery on bypass for congenital heart disease. Three hundred seventy-one patients were included in the study, and their median age was 2.75 years (IQR 0.65-6.63). Forty-seven children had persistent LSVC (12.7 %), and persistent LSVC was identified on echocardiography before surgery in 39 patients (83 %). In three patients (6.4 %) with persistent LSVC, significant inflow obstruction of the left ventricle developed after surgery leading to low output syndrome or secondary pulmonary hypertension. In eight patients (17 %), persistent LSVC was associated with a partially or completely unroofed coronary sinus and in two cases (4 %) with coronary sinus ostial atresia. Duration of mechanical ventilation was significantly shorter in the control group (1.2 vs. 3.0 days, p = 0.04), whereas length of stay in intensive care did not differ. Mortality was also significantly lower in the control group (2.5 vs. 10.6 %, p = 0.004). The results of study show that persistent LSVC in association with congenital cardiac malformation increases the risk of mortality in children with cardiac surgery on cardiopulmonary bypass. Recognition of a persistent LSVC and its associated anomalies is mandatory to avoid complications during or after cardiac surgery.

Drowning--post-mortem imaging findings by computed tomography

The aim of this study was to identify the classic autopsy signs of drowning in post-mortem multislice computed tomography (MSCT). Therefore, the post-mortem pre-autopsy MSCT- findings of ten drowning cases were correlated with autopsy and statistically compared with the post-mortem MSCT of 20 non-drowning cases. Fluid in the airways was present in all drowning cases. Central aspiration in either the trachea or the main bronchi was usually observed. Consecutive bronchospasm caused emphysema aquosum. Sixty percent of drowning cases showed a mosaic pattern of the lung parenchyma due to regions of hypo- and hyperperfused lung areas of aspiration. The resorption of fresh water in the lung resulted in hypodensity of the blood representing haemodilution and possible heart failure. Swallowed water distended the stomach and duodenum; and inflow of water filled the paranasal sinuses (100%). All the typical findings of drowning, except Paltau's spots, were detected using post-mortem MSCT, and a good correlation of MSCT and autopsy was found. The advantage of MSCT was the direct detection of bronchospasm, haemodilution and water in the paranasal sinus, which is rather complicated or impossible at the classical autopsy.

Angioscopy for experimental evaluation of optional IVC filters

PURPOSE: To demonstrate the feasibility of direct angioscopic visualization of an optional inferior vena cava (IVC) filter in situ and during retrieval. MATERIALS AND METHODS: Angioscopy was used for direct visualization of optional IVC filters in six sheep. Cavograms were obtained before the filters were retrieved. After successful filter retrieval, segmental IVC perfusion was performed to evaluate filter retrieval-related damage to the IVC wall. Therefore, all branch vessels were ligated before the IVC segment was flushed with normal saline solution until it was fully distended. Then, the inflow was terminated and the IVC segment observed for deflation. Subsequently, the IVC was harvested en bloc, dissected, and inspected macroscopically. RESULTS: The visibility of IVC filters at angioscopy was excellent. During the retrieval procedure, filter collapse and retraction into the sheath were clearly demonstrated. Angioscopy provided additional information to that obtained with cavography, demonstrating adherent material in three filters. Three filters in place for more than 2 months could not be retrieved because the filter legs were incorporated into the IVC wall. After filter retrieval, there was no perforation at segmental IVC perfusion. At macroscopic inspection of the IVC lumen, a small piece of detached endothelium was found in one animal. CONCLUSION: Angioscopy enabled the direct evaluation of optional IVC filters in situ and during retrieval. Compared with cavography, angioscopy provided additional information about the filter in situ and the retrieval procedure. Future applications of this technique could include studies of filter migration, compression, and clot-trapping efficacy.

Advanced chronic heart failure: A position statement from the Study Group on Advanced Heart Failure of the Heart Failure Association of the European Society of Cardiology

Therapy has improved the survival of heart failure (HF) patients. However, many patients progress to advanced chronic HF (ACHF). We propose a practical clinical definition and describe the characteristics of this condition. Patients that are generally recognised as ACHF often exhibit the following characteristics: 1) severe symptoms (NYHA class III to IV); 2) episodes with clinical signs of fluid retention and/or peripheral hypoperfusion; 3) objective evidence of severe cardiac dysfunction, shown by at least one of the following: left ventricular ejection fraction<30%, pseudonormal or restrictive mitral inflow pattern at Doppler-echocardiography; high left and/or right ventricular filling pressures; elevated B-type natriuretic peptides; 4) severe impairment of functional capacity demonstrated by either inability to exercise, a 6-minute walk test distance<300 m or a peak oxygen uptake<12-14 ml/kg/min; 5) history of >1 HF hospitalisation in the past 6 months; 6) presence of all the previous features despite optimal therapy. This definition identifies a group of patients with compromised quality of life, poor prognosis, and a high risk of clinical events. These patients deserve effective therapeutic options and should be potential targets for future clinical research initiatives.

Intrapulmonary 3He gas distribution depending on bolus size and temporal bolus placement

OBJECTIVE: Dynamic ventilation (3)He-MRI is a new method to assess pulmonary gas inflow. As differing airway diameters throughout the ventilatory cycle can influence gas inflow this study intends to investigate the influence of volume and timing of a He gas bolus with respect to the beginning of the tidal volume on inspiratory gas distribution. MATERIALS AND METHODS: An ultrafast 2-dimensional spoiled gradient echo sequence (temporal resolution 100 milliseconds) was used for dynamic ventilation (3)He-MRI of 11 anesthetized and mechanically ventilated pigs. The applied (3)He gas bolus was varied in volume between 100 and 200 mL. A 150-mL bolus was varied in its application time after the beginning of the tidal volume between 0 and 1200 milliseconds. Signal kinetics were evaluated using an in-house developed software after definition of parameters for the quantitative description of (3)He gas inflow. RESULTS: The signal rise time (time interval between signal in the parenchyma reaches 10% and 90% of its maximum) was prolonged with increasing bolus volume. The parameter was shortened with increasing delay of (3)He application after the beginning of the tidal volume. Timing variation as well as volume variation showed no clear interrelation to the signal delay time 10 (time interval between signal in the trachea reaches 50% of its maximum and signal in the parenchyma reaches 10% of its maximum). CONCLUSIONS: Dynamic ventilation (3)He-MRI is able to detect differences in bolus geometry performed by volume variation. Pulmonary gas inflow as investigated by dynamic ventilation (3)He-MRI tends to be accelerated by an increasing application delay of a (3)He gas bolus after the beginning of the tidal volume.

Diastolic dysfunction of the cardiac allograft and maximal exercise capacity

BACKGROUND: Peak oxygen uptake (peak Vo(2)) is an established integrative measurement of maximal exercise capacity in cardiovascular disease. After heart transplantation (HTx) peak Vo(2) remains reduced despite normal systolic left ventricular function, which highlights the relevance of diastolic function. In this study we aim to characterize the predictive significance of cardiac allograft diastolic function for peak Vo(2). METHODS: Peak Vo(2) was measured using a ramp protocol on a bicycle ergometer. Left ventricular (LV) diastolic function was assessed with tissue Doppler imaging sizing the velocity of the early (Ea) and late (Aa) apical movement of the mitral annulus, and conventional Doppler measuring early (E) and late (A) diastolic transmitral flow propagation. Correlation coefficients were calculated and linear regression models fitted. RESULTS: The post-transplant time interval of the 39 HTxs ranged from 0.4 to 20.1 years. The mean age of the recipients was 55 +/- 14 years and body mass index (BMI) was 25.4 +/- 3.9 kg/m(2). Mean LV ejection fraction was 62 +/- 4%, mean LV mass index 108 +/- 22 g/m(2) and mean peak Vo(2) 20.1 +/- 6.3 ml/kg/min. Peak Vo(2) was reduced in patients with more severe diastolic dysfunction (pseudonormal or restrictive transmitral inflow pattern), or when E/Ea was > or =10. Peak Vo(2) correlated with recipient age (r = -0.643, p < 0.001), peak heart rate (r = 0.616, p < 0.001) and BMI (r = -0.417, p = 0.008). Of all echocardiographic measurements, Ea (r = 0.561, p < 0.001) and Ea/Aa (r = 0.495, p = 0.002) correlated best. Multivariate analysis identified age, heart rate, BMI and Ea/Aa as independent predictors of peak Vo(2). CONCLUSIONS: Diastolic dysfunction is relevant for the limitation of maximal exercise capacity after HTx.

Holocene climatic change and the development of the lake-effect snowbelt in Michigan, USA

Lake-effect snow is an important constraint on ecological and socio-economic systems near the North American Great Lakes. Little is known about the Holocene history of lake-effect snowbelts, and it is difficult to decipher how lake-effect snowfall abundance affected ecosystem development. We conducted oxygen-isotope analysis of calcite in lake-sediment cores from northern Lower Michigan to infer Holocene climatic variation and assess snowbelt development. The two lakes experience the same synoptic-scale climatic systems, but only one of them (Huffman Lake) receives a significant amount of lake-effect snow. A 177-cm difference in annual snowfall causes groundwater inflow at Huffman Lake to be 18O-depleted by 2.3‰ relative to O'Brien Lake. To assess when the lake-effect snowbelt became established, we compared calcite-δ18O profiles of the last 11,500 years from these two sites. The chronologies are based on accelerator-mass-spectrometry 14C ages of 11 and 17 terrestrial-plant samples from Huffman and O'Brien lakes, respectively. The values of δ18O are low at both sites from 11,500 to 9500 cal yr BP when the Laurentide Ice Sheet (LIS) exerted a dominant control over the regional climate and provided periodic pulses of meltwater to the Great Lakes basin. Carbonate δ18O increases by 2.6‰ at O'Brien Lake and by 1.4‰ at Huffman Lake between 9500 and 7000 cal yr BP, suggesting a regional decline in the proportion of runoff derived from winter precipitation. The Great Lakes snowbelt probably developed between 9500 and 5500 cal yr BP as inferred from the progressive 18O-depletion at Huffman Lake relative to O'Brien Lake, with the largest increase of lake-effect snow around 7000 cal yr BP. Lake-effect snow became possible at this time because of increasing contact between the Great Lakes and frigid arctic air. These changes resulted from enhanced westerly flow over the Great Lakes as the LIS collapsed, and from rapidly rising Great Lakes levels during the Nipissing Transgression. The δ18O difference between Huffman and O'Brien lakes declines after 5500 cal yr BP, probably because of a northward shift of the polar vortex that brought increasing winter precipitation to the entire region. However, δ18O remains depleted at Huffman Lake relative to O'Brien Lake because of the continued production of lake-effect snow.

Influence of left ventricular relaxation on the pressure half time of aortic regurgitation

BACKGROUND The severity of aortic regurgitation can be estimated using pressure half time (PHT) of the aortic regurgitation flow velocity, but the correlation between regurgitant fraction and PHT is weak. AIM To test the hypothesis that the association between PHT and regurgitant fraction is substantially influenced by left ventricular relaxation. METHODS In 63 patients with aortic regurgitation, subdivided into a group without (n = 22) and a group with (n = 41) left ventricular hypertrophy, regurgitant fraction was calculated using the difference between right and left ventricular cardiac outputs. Left ventricular relaxation was assessed using the early to late diastolic Doppler tissue velocity ratio of the mitral annulus (E/ADTI), the E/A ratio of mitral inflow (E/AM), and the E deceleration time (E-DT). Left ventricular hypertrophy was assessed using the M mode derived left ventricular mass index. RESULTS The overall correlation between regurgitant fraction and PHT was weak (r = 0.36, p < 0.005). In patients without left ventricular hypertrophy, there was a significant correlation between regurgitant fraction and PHT (r = 0.62, p < 0.005), but not in patients with left ventricular hypertrophy. In patients with a left ventricular relaxation abnormality (defined as E/ADTI< 1, E/AM< age corrected lower limit, E-DT >/= 220 ms), no associations between regurgitant fraction and PHT were found, whereas in patients without left ventricular relaxation abnormalities, the regurgitant fraction to PHT relations were significant (normal E/AM: r = 0.57, p = 0.02; E-DT< 220 ms: r = 0.50, p < 0.001; E/ADTI < 1: r = 0.57, p = 0.02). CONCLUSIONS Only normal left ventricular relaxation allows a significant decay of PHT with increasing aortic regurgitation severity. In abnormal relaxation, which is usually present in left ventricular hypertrophy, wide variation in prolonged backward left ventricular filling may cause dissociation between the regurgitant fraction and PHT. Thus the PHT method should only be used in the absence of left ventricular relaxation abnormalities.

In vitro precision of fit of computer-aided designed and computer-aided manufactured titanium screw-retained fixed dental prostheses before and after ceramic veneering

OBJECTIVE To compare the precision of fit of full-arch implant-supported screw-retained computer-aided designed and computer-aided manufactured (CAD/CAM) titanium-fixed dental prostheses (FDP) before and after veneering. The null-hypothesis was that there is no difference in vertical microgap values between pure titanium frameworks and FDPs after porcelain firing. MATERIALS AND METHODS Five CAD/CAM titanium grade IV frameworks for a screw-retained 10-unit implant-supported reconstruction on six implants (FDI tooth positions 15, 13, 11, 21, 23, 25) were fabricated after digitizing the implant platforms and the cuspid-supporting framework resin pattern with a laser scanner (CARES(®) Scan CS2; Institut Straumann AG, Basel, Switzerland). A bonder, an opaquer, three layers of porcelain, and one layer of glaze were applied (Vita Titankeramik) and fired according to the manufacturer's preheating and fire cycle instructions at 400-800°C. The one-screw test (implant 25 screw-retained) was applied before and after veneering of the FDPs to assess the vertical microgap between implant and framework platform with a scanning electron microscope. The mean microgap was calculated from interproximal and buccal values. Statistical comparison was performed with non-parametric tests. RESULTS All vertical microgaps were clinically acceptable with values <90 μm. No statistically significant pairwise difference (P = 0.98) was observed between the relative effects of vertical microgap of unveneered (median 19 μm; 95% CI 13-35 μm) and veneered FDPs (20 μm; 13-31 μm), providing support for the null-hypothesis. Analysis within the groups showed significantly different values between the five implants of the FDPs before (P = 0.044) and after veneering (P = 0.020), while a monotonous trend of increasing values from implant 23 (closest position to screw-retained implant 25) to 15 (most distant implant) could not be observed (P = 0.169, P = 0.270). CONCLUSIONS Full-arch CAD/CAM titanium screw-retained frameworks have a high accuracy. Porcelain firing procedure had no impact on the precision of fit of the final FDPs. All implant microgap measurements of each FDP showed clinically acceptable vertical misfit values before and after veneering. Thus, the results do not only show accurate performance of the milling and firing but show also a reproducible scanning and designing process.

Amplified inception of European Little Ice Age by sea ice-ocean-atmosphere feedbacks

The inception of the Little Ice Age (~1400–1700 AD) is believed to have been driven by an interplay of external forcing and climate system internal variability. While the hemispheric signal seems to have been dominated by solar irradiance and volcanic eruptions, the understanding of mechanisms shaping the climate on a continental scale is less robust. In an ensemble of transient model simulations and a new type of sensitivity experiments with artificial sea ice growth, the authors identify a sea ice–ocean–atmosphere feedback mechanism that amplifies the Little Ice Age cooling in the North Atlantic–European region and produces the temperature pattern suggested by paleoclimatic reconstructions. Initiated by increasing negative forcing, the Arctic sea ice substantially expands at the beginning of the Little Ice Age. The excess of sea ice is exported to the subpolar North Atlantic, where it melts, thereby weakening convection of the ocean. Consequently, northward ocean heat transport is reduced, reinforcing the expansion of the sea ice and the cooling of the Northern Hemisphere. In the Nordic Seas, sea surface height anomalies cause the oceanic recirculation to strengthen at the expense of the warm Barents Sea inflow, thereby further reinforcing sea ice growth. The absent ocean–atmosphere heat flux in the Barents Sea results in an amplified cooling over Northern Europe. The positive nature of this feedback mechanism enables sea ice to remain in an expanded state for decades up to a century, favoring sustained cold periods over Europe such as the Little Ice Age. Support for the feedback mechanism comes from recent proxy reconstructions around the Nordic Seas.

Characterization of microcirculation in multiple sclerosis lesions by dynamic texture parameter analysis (DTPA)

OBJECTIVE Texture analysis is an alternative method to quantitatively assess MR-images. In this study, we introduce dynamic texture parameter analysis (DTPA), a novel technique to investigate the temporal evolution of texture parameters using dynamic susceptibility contrast enhanced (DSCE) imaging. Here, we aim to introduce the method and its application on enhancing lesions (EL), non-enhancing lesions (NEL) and normal appearing white matter (NAWM) in multiple sclerosis (MS). METHODS We investigated 18 patients with MS and clinical isolated syndrome (CIS), according to the 2010 McDonald's criteria using DSCE imaging at different field strengths (1.5 and 3 Tesla). Tissues of interest (TOIs) were defined within 27 EL, 29 NEL and 37 NAWM areas after normalization and eight histogram-based texture parameter maps (TPMs) were computed. TPMs quantify the heterogeneity of the TOI. For every TOI, the average, variance, skewness, kurtosis and variance-of-the-variance statistical parameters were calculated. These TOI parameters were further analyzed using one-way ANOVA followed by multiple Wilcoxon sum rank testing corrected for multiple comparisons. RESULTS Tissue- and time-dependent differences were observed in the dynamics of computed texture parameters. Sixteen parameters discriminated between EL, NEL and NAWM (pAVG = 0.0005). Significant differences in the DTPA texture maps were found during inflow (52 parameters), outflow (40 parameters) and reperfusion (62 parameters). The strongest discriminators among the TPMs were observed in the variance-related parameters, while skewness and kurtosis TPMs were in general less sensitive to detect differences between the tissues. CONCLUSION DTPA of DSCE image time series revealed characteristic time responses for ELs, NELs and NAWM. This may be further used for a refined quantitative grading of MS lesions during their evolution from acute to chronic state. DTPA discriminates lesions beyond features of enhancement or T2-hypersignal, on a numeric scale allowing for a more subtle grading of MS-lesions.

Eemian and Holocene sea-surface conditions in the southern Black Sea: organic-walled dinoflagellate cyst record from core 22-GC3

In order to compare the sea-surface conditions in the Black Sea during the Holocene and Eemian, sapropelic parts of marine core 22-GC3 (42°13.53′N/36°29.55′E, 838 m water depth) were studied for organic-walled dinoflagellate cyst content. The record shows a change from freshwater/brackish assemblages (Pyxidinopsis psilata, Spiniferites cruciformis, and Caspidinium rugosum) to more marine assemblages (Lingulodinium machaerophorum and Spiniferites ramosus complex) during each interglacial, due to the inflow of saline Mediterranean water. The lacustrine–marine transitions in 22-GC3 occurred at ~ 8.3 cal kyr BP during the early Holocene and ~ 128 kyr BP during the early Eemian, slightly later compared to the onset of interglacial conditions on the adjacent land. Dinoflagellate cyst assemblages reveal higher sea-surface salinity (~ 28–30) (e.g. Spiniferites pachydermus, Bitectatodinium tepikiense, and Spiniferites mirabilis) around ~ 126.5–121 kyr BP in comparison to the Holocene (~ 15–20) as well as relatively high sea-surface temperature (e.g. Tuberculodinium vancampoae, S. pachydermus, and S. mirabilis) especially at ~ 127.6–125.3 kyr BP. Establishment of high sea-surface salinity during the Eemian correlates very well with reconstructed relatively high global sea-level and is explained as a combined effect of increased Mediterranean supply and high temperatures at the beginning of the last interglacial. The observed changes in the dinocyst record highlight the importance of nutrients for the composition of the Eemian and Holocene dinocyst assemblages.