Modern diatom, cladocera, chironomid, and chrysophyte cyst assemblages as quantitative indicators for the reconstruction of past environmental conditions in the Alps. I. Climate

Diatom, chrysophyte cyst, benthic cladocera, planktonic cladocera, and chironomid assemblages were studied in the surface sediments of 68 small lakes along an altitudinal gradient from 300 to 2350 m in Switzerland. In addition, 43 environmental variables relating to the physical limnology, geography, catchment characteristics, climate, and water chemistry were recorded or measured for each lake. The explanatory power of each of these predictor variables for the different biological data-sets was estimated by a series of canonical correspondence analyses (CCA) and the statistical significance of each model was assessed by Monte Carlo permutation tests. A minimal set of environmental variables was found for each biological data-set by a forward-selection procedure within CCA. The unique, independent explanatory power of each set of environmental variables was estimated by a series of CCAs and partial CCAs. Inference models or transfer functions for mean summer (June, July, August) air temperature were developed for each biological data-set using weighted-averaging partial least squares or partial least squares. The final transfer functions, after data screening, have root mean squared errors of prediction, as assessed by leave-one-out cross-validation, of 1.37 °C (chironomids), 1.60 °C (benthic cladocera), 1.62 °C (diatoms), 1.77 °C (planktonic cladocera), and 2.23 °C (chrysophyte cysts).     

Rescaling species optima estimated by weighted averaging

The common practice of linear deshrinking in weighted averaging is known to be equivalent to a linear rescaling of the estimated species optima. In published lists of species optima, the use of rescaling is recommended, as it allows values derived from different data sets to be compared and used for new inferences, assuming that taxonomic consistency is assured. Rescaling optima is also shown to influence WA estimates of species tolerances. Non-linear rescaling is also discussed, in the form of cubical rescaling and weighted averaging-partial least squares (WA-PLS). The use of a different deshrinking equation in a small data set did lead to similar prediction errors, probably because of the small size of the data set used.     

Geostatistical integration of near-surface geophysical data

Accurate statics calculation and near‐surface related noise removal require a detailed knowledge of the near‐surface velocity field. Conventional seismic surveys currently are not designed to provide this information, and 3D high‐resolution reflection/refraction acquisition is not feasible for large survey areas. Satellite images and vibrator plate attributes are dense low‐cost data, which can be used in spatially extrapolating velocities from sparse uphole data by geostatistics. We tested this approach in two different areas of Saudi Arabia and found that the optimal recipe depends on the local geology.     

The near-surface information gap for time and depth imaging

The shallowest few hundred metres of the earth cannot be adequately imaged by conventional seismic when tuned for deeper targets. Adding independent measurements (such as uphole or shallow refraction surveys) reduces this information gap, but in some arid areas (such as Saudi Arabia) the near‐surface complexities are not well resolved, even in this way. The joint tomographic inversion of different wave types can contribute to reducing these ambiguities further, by complementing the different penetration ranges and propagation directions of reflected, refracted and diving waves. Here, we demonstrate the weakness of diving waves when used alone, and the value of complementing them by available reflected and refracted arrivals.     

The build up of the isotopic signal in skeletons of the stony coral Porites lutea

The build up of the isotopic signal in corals was followed by sampling the newly formed skeleton at a monthly resolution for a period of two years in order to establish the interrelations between the calcification processes and the skeletal isotopic composition. We deployed two underwater sampling schemes, which provide a monitor of the changes in water temperature and δ¹⁸O and in the corresponding newly accreted skeleton of undisturbed Porites lutea colonies under natural conditions and four transplanted colonies, which maintained the genetic identity throughout the experiment. The results indicate that δ¹⁸O of the newly accreted skeleton does not correlate with ambient temperature although the seasonal temperature variability at the site (winter to summer) is in the order of 6 °C and δ¹⁸O of seawater is constant throughout the year. In contrast to the newly formed surface skeleton, the isotopic compositions of the deep and older parts of the skeleton show the predicted annual isotopic pattern with highly significant correlation between δ¹⁸O_s and SST. The transformation between temperature-independent to temperature-dependent isotopic signal occurs several months after the skeleton was formed at the surface. The position of the skeleton in relation to the open sea may generate the difference between δ¹⁸O_s of the surface skeleton and that of the skeleton previously accreted further down the tissue layer. Our data support the general model of a multi-step skeletogenesis process, where the temperature independent skeleton is entails the first step, the production of skeletal scaffold, and the environmental temperature signature is captured by the next two other steps: the thickening and the periodic abrupt uplift occurring at the depth of the tissue layer. However, re-examination and development of the current isotopic models for coral calcification are required in order to explain the observed different temperature dependency during the growth’s sequence.     

Numerical integration for the real time production of fundamental ephemerides over a wide time span

A simplified model of the solar system has been developed along with an integration method, enabling to compute planetary and lunar ephemerides to an accuracy better than 1 and 2 milliarcsecs, respectively. On current personal computers, the integration procedure (SOLEX) is fast enough that by using a relatively small ( 20 Kbytes/Cy) database of starting conditions, any epoch in the time interval (up to ±100 Cy) covered by the database can be reached by the integrator in a few seconds. This makes the algorithm convenient for the direct computation of high precision ephemerides over a time span of several millennia.     

Dynamic active earth pressure on cantilever retaining walls

In normal practice, the active earth pressure on cantilever retaining wall is evaluated with different procedures relating to an ideal vertical plane passing through the heel of the wall. If the wall presents a long heel, failure planes do not interfere with the vertical stem, so that the limit Rankine conditions can develop freely in the backfill. The inclination of lateral actions along the ideal plane is assumed to be constant and depends on the geometry of the ground level and on the friction angle φ. The Authors recently proposed a new method to evaluate the active earth pressure coefficient due to seismic loading with a pseudo-static stress plasticity solution. The present paper describes the application of this method to a retaining wall supporting a φ soil backfill with an irregular surface. For two different configurations of wall-soil system, the behaviour is also studied by continuum FDM dynamic analyses, utilising four Italian accelerometric time-histories scaled at the same peak ground acceleration. The comparison between different procedures is also analysed.     

Fractionation of oxygen and hydrogen isotopes in evaporating water

Variations in oxygen and hydrogen isotope ratios of water and ice are powerful tools in hydrology and ice core studies. These variations are controlled by both equilibrium and kinetic isotope effects during evaporation and precipitation, and for quantitative interpretation it is necessary to understand how these processes affect the isotopic composition of water and ice. Whereas the equilibrium isotope effects are reasonably well understood, there is controversy on the magnitude of the kinetic isotope effects of both oxygen and hydrogen and the ratio between them. In order to resolve this disagreement, we performed evaporation experiments into air, argon and helium over the temperature range from 10 to 70 °C. From these measurements we derived the isotope effects for vapor diffusion in gas phase (ε_{diff(HD¹⁶O)} for D/H and ε_{diff(H2¹⁸O)} for ¹⁸O/¹⁶O). For air, the ratio ε_{diff(HD¹⁶O)}/ε_{diff(H2¹⁸O)} at 20 °C is 0.84, in very good agreement with Merlivat (1978) (0.88), but in considerable inconsistency with Cappa et al. (2003) (0.52). Our results support Merlivat’s conclusion that measured ε_{diff(HD¹⁶O)}/ε_{diff(H2¹⁸O)} ratios are significantly different than ratios calculated from simplified kinetic theory of gas diffusion. On the other hand, our experiments with helium and argon suggest that this discrepancy is not due to isotope effects of molecular collision diameters. We also found, for the first time, that the ε_{diff(HD¹⁶O)}/ε_{diff(H2¹⁸O)} ratio tends to increase with cooling. This new finding may have important implications to interpretations of deuterium excess (d-excess = δD − 8δ¹⁸O) in ice core records, because as we show, the effect of temperature on d-excess is of similar magnitude to glacial interglacial variations in the cores.     

Thermal state of the Campi Flegrei caldera inferred from seismic attenuation tomography

A three-dimensional Qp image of the Campi Flegrei caldera between 0 and 3 km of depth has been inferred by the inversion of P rise time and pulse width data of 87 local earthquakes recorded during the last bradiseismic crisis by a local array deployed in the area by the University of Wisconsin. The availability of both thermal measurements in 5 deep boreholes and of a heat flow surface map of the area allowed us to calibrate the local temperature F vs. Qp relationship. The comparison of Qp, Vp and Vp/Vs images, combined with hydrogeological and geochemical data from deep boreholes, allowed us to distinguish some low-Qp anomalies related to the presence of fluids in the rocks from a deep low-Qp anomaly related to the conductive cooling of a magma reservoir. The deep anomaly is located in the same zone where several authors believe that the volcanic and magmatic activity migrated after the Neapolitan Yellow Tuff eruption. Moreover this anomaly includes the area where the existence of a magma chamber at depth between 4 and 5 km was inferred by an active seismic experiment.