Reconstruction of paleoenvironments of the Great Salt Lake Basin during the late Cenozoic

This study summarizes the results of micropaleontological, sedimentological, and isotope geochemical analyses of cuttings from five deep wells drilled in the Great Salt Lake (Utah, USA). Spanning the last 5.0 million yrs, our environmental history of the Great Salt Lake distinguishes four intervals based on paleobiological and sedimentological characteristics, using a previously developed tephrochronology for age control. For most of its history, the Great Salt Lake Basin has been occupied by a mixture of marsh, shallow lacustrine and sand flat conditions. In contrast, open lake conditions, typical of the Bonneville cycles and the modern Great Salt Lake apparently have only dominated the basin for the past 0.6-0.8 Ma. The two main structural basins in the study area (the North and South Basins) experienced different lacustrine histories. Large but frequently saline lakes occupied the North Basin after about 0.6 Ma. In the South Basin, ephemeral, saline lacustrine conditions started at 2.1 Ma and developed to full lacustrine conditions at 0.3 Ma. Our paleoenvironmental interpretations are broadly consistent with the aquatic palynological records from the same wells, as well as with the prior core- and outcrop-based lines of evidence. However, the differences in lake history between the North and South Basin have not been previously recognized.     

Configurational thermodynamics of aluminous pyroxenes: A generalized pair approximation

A pair approximation is used to estimate the effects of short-range order on the thermodynamic properties of aluminous clinopyroxenes on the joins diopside (CaMg-Si₂O₆)-jadeite (NaAlSi₂O₆) and diopside-CaTs (CaAl₂SiO₆). The generalized pair approximation is the simplest model for concentrated solutions which includes short-range order. Short-range order is expected to be especially significant in coupled solid solutions, such as aluminous pyroxenes, since atoms of different valence substitute for each other. The calculations show that the random model, in which the configurational entropy is calculated as if atoms on each crystallographic site mix randomly, is appropriate as a first approximation. The excess entropy relative to the random model behaves regularly, is always negative, and becomes more negative as temperature decreases or the ordering energies increase. The excess entropy relative to the random model can be modeled reasonably well with a simple power series, or Margules-type, formulation. In contrast, the excess entropy relative to a molecular model, in which the ideal activity is assumed to be equal to some mole fraction, is irregular, can be positive or negative, and even changes in sign with variations in temperature and composition. The configurational enthalpy is positive at high temperatures, and becomes negative with decreasing temperature or increasing ordering energy. The mixing enthalpy can have non-configurational contributions, in addition to the effective short-range configurational contributions considered explicitly. The pair approximation predicts an ordering transition from C2/c to P2₁/n for CaTs and diopside-CaTs solutions at moderate to low temperatures, respectively. A field where C2/c orders to C2 is also found. A higher order approximation, different relative ordering energies, or quantitative consideration of strain contributions is required to account for the C2/c to P2/n transition in omphacites. There is no justification for molecular models, in which the configurational entropy is calculated as if endmember “molecules” were mixing in the crystal, in either concentrated or dilute solutions. Molecular models do not represent limiting ordered states for coupled solid solutions.     

Crustal structure of the Levant Basin, eastern Mediterranean

A seismic refraction/wide-angle reflection experiment was undertaken in the Levant Basin, eastern Mediterranean. Two roughly east–west profiles extend from the continental shelf of Israel toward the Levant Basin. The northern profile crosses the Eratosthenes Seamount and the southern profile crosses several distinct magnetic anomalies. The marine operation used 16 ocean bottom seismometers deployed along the profiles with an air gun array and explosive charges as energy sources. The results of this study strongly suggest the existence of oceanic crust under portions of the Levant Basin and continental crust under the Eratosthenes Seamount. The seismic refraction data also indicate a large sedimentary sequence, 10–14 km thick, in the Levant Basin and below the Levant continental margin. Assuming the crust is of Cretaceous age, this gives a fairly high sedimentation rate. The sequence can be divided into several units. A prominent unit is the 4.2 km/s layer, which is probably composed of the Messinian evaporites. Overlying the evaporitic layer are layers composed of Plio–Pleistocene sediments, whose velocity is 2.0 km/s. The refraction profiles and gravity and magnetic models indicate that a transition from a two layer continental to a single-layer oceanic crust takes place along the Levant margin. The transition in the structure along the southern profile is located beyond the continental margin and it is quite gradual. The northern profile, north of the Carmel structure, presents a different structure. The continental crust is much thinner there and the transition in the crustal structure is more rapid. The crustal thinning begins under western Galilee and terminates at the continental slope. The results of the present study indicate that the Levant Basin is composed of distinct crustal units and that the Levant continental margin is divided into at least two provinces of different crustal structure.     

Comment on G. Wansard and F. Mezquita, The response of ostracode shell chemistry to seasonal change in a Mediterranean freshwater spring environment

In a study of the minor element chemistry of ostracode shells, Wansard and Mezquita (2001) concluded that the Mg/Ca and Sr/Ca of the ambient water controlled the respective ratios in ostracode calcite. Contrary to their conclusion that minor element chemistry is not influenced by temperature, we find a very strong statistical correlation (r = 0.97) between temperature and shell chemistry in their data. This discrepancy apparently arises due to their use of a partition coefficient model of shell chemistry that masks a number of strong correlations in their data. We argue here that the partition coefficient model is not appropriate for this biologically-mediated carbonate, and that a broader range of possible controls on shell chemistry needs to be tested.     

Slowness domain offset and traveltime approximations in layered vertical transversely isotropic media

Considering horizontally layered transversely isotropic media with vertical symmetry axis and all types of pure‐mode and converted waves we present a new wide‐angle series approximation for the kinematical characteristics of reflected waves: horizontal offset, intercept time, and total reflection traveltime as functions of horizontal slowness. The method is based on combining (gluing) both zero‐offset and (large) finite‐offset series coefficients. The horizontal slowness is bounded by the critical value, characterised by nearly horizontal propagation within the layer with the highest horizontal velocity. The suggested approximation uses five parameters to approximate the offset, six parameters to approximate the intercept time or the traveltime, and seven parameters to approximate any two or all three kinematical characteristics. Overall, the method is very accurate for pure‐mode compressional waves and shear waves polarised in the horizontal plane and for converted waves. The application of the method to pure‐mode shear waves polarised in the vertical plane is limited due to cusps and triplications. To demonstrate the high accuracy of the method, we consider a synthetic, multi‐layer model, and we plot the normalised errors with respect to numerical ray tracing.     

Slowness-domain kinematical characteristics for horizontally layered orthorhombic media. Part I: Critical slowness match

Kinematical characteristics of reflected waves in anisotropic elastic media play an important role in the seismic imaging workflow. Considering compressional and converted waves, we derive new, azimuthally dependent, slowness-domain approximations for the kinematical characteristics of reflected waves (radial and transverse offsets, intercept time and traveltime) for layered orthorhombic media with varying azimuth of the vertical symmetry planes. The proposed method can be considered an extension of the well-known ‘generalized moveout approximation’ in the slowness domain, from azimuthally isotropic to azimuthally anisotropic models. For each slowness azimuth, the approximations hold for a wide angle range, combining power series coefficients in the vicinity of both the normal-incidence ray and an additional wide-angle ray. We consider two cases for the wide-angle ray: a ‘critical slowness match’ and a ‘pre-critical slowness match’ studied in Parts I and II of this work, respectively. For the critical slowness match, the approximations are valid within the entire slowness range, up to the critical slowness. For the ‘pre-critical slowness match’, the approximations are valid only within the bounded slowness range; however, the accuracy within the defined range is higher. The critical slowness match is particularly effective when the subsurface model includes a dominant high-velocity layer where, for nearly critical slowness values, the propagation in this layer is almost horizontal. Comparing the approximated kinematical characteristics with those computed by numerical ray tracing, we demonstrate high accuracy.     

Traveltime approximation in vertical transversely isotropic layered media

The well‐known asymptotic fractional four‐parameter traveltime approximation and the five‐parameter generalised traveltime approximation in stratified multi‐layer transversely isotropic elastic media with a vertical axis of symmetry have been widely used for pure‐mode and converted waves. The first three parameters of these traveltime expansions are zero‐offset traveltime, normal moveout velocity, and quartic coefficient, ensuring high accuracy of traveltimes at short offsets. The additional parameter within the four‐parameter approximation is an effective horizontal velocity accounting for large offsets, which is important to avoid traveltime divergence at large offsets. The two additional parameters in the above‐mentioned five‐parameter approximation ensure higher accuracy up to a given large finite offset with an exact match at this offset. In this paper, we propose two alternative five‐parameter traveltime approximations, which can be considered extensions of the four‐parameter approximation and an alternative to the five‐parameter approximation previously mentioned. The first three short‐offset parameters are the same as before, but the two additional long‐offset parameters are different and have specific physical meaning. One of them describes the propagation in the high‐velocity layer of the overburden (nearly horizontal propagation in the case of very large offsets), and the other characterises the intercept time corresponding to the critical slowness that includes contributions of the lower velocity layers only. Unlike the above‐mentioned approximations, both of the proposed traveltime approximations converge to the theoretical (asymptotic) linear traveltime at the limit case of very large (“infinite”) offsets. Their accuracy for moderate to very large offsets, for quasi‐compressional waves, converted waves, and shear waves polarised in the horizontal plane, is extremely high in cases where the overburden model contains at least one layer with a dominant higher velocity compared with the other layers. We consider the implementation of the proposed traveltime approximations in all classes of problems in which the above‐mentioned approximations are used, such as reflection and diffraction analysis and imaging.     

Large roots dominate the contribution of trees to slope stability

Tree roots provide surface erosion protection and improve slope stability through highly complex interactions with the soil due to the nature of root systems. Root reinforcement estimation is usually performed by in situ pullout tests, in which roots are pulled out of the soil to reliably estimate the root strength of compact soils. However, this test is not suitable for the scenario where a soil progressively fails in a series of slump blocks – for example, in unsupported soils near streambanks and road cuts where the soil has no compressive resistance at the base of the hillslope. The scenario where a soil is unsupported on its downslope extent and progressively deforms at a slow strain rate has received little attention, and we are unaware of any study on root reinforcement that estimates the additional strength provided by roots in this situation. We therefore designed two complementary laboratory experiments to compare the force required to pull the root out. The results indicate that the force required to pull out roots is reduced by up to 50% when the soil fails as slump blocks compared to pullout tests. We also found that, for slump block failure, roots had a higher tendency to slip than to break, showing the importance of active earth pressure on root reinforcement behaviour, which contributes to reduced friction between soil and roots. These results were then scaled up to a full tree and tree stand using the root bundle and field-measured spatial distributions of root density. Although effects on the force mobilized in small roots can be relevant, small roots have virtually no effect on root reinforcement at the tree or stand scale on hillslopes. When root distribution has a wide range of diameters, the root reinforcement results are controlled by large roots, which hold much more force than small roots. © 2019 John Wiley & Sons, Ltd.