Origins and evolution of rhyolitic magmas in the central Snake River Plain: insights from coupled high-precision geochronology,oxygen isotope,and hafnium isotope analyses of zircon

We present new high-precision CA-ID-TIMS and in situ U–Pb ages together with Hf and O isotopic analyses (analyses performed all on the same grains) from four tuffs from the 15?10 Ma Bruneau–Jarbidge center of the Snake River Plain and from three rhyolitic units from the Kimberly borehole in the neighboring 10?6 Ma Twin Falls volcanic center. We find significant intrasample diversity in zircon ages (ranges of up to 3 Myr) and in δ¹⁸O (ranges of up to 6‰) and ε_Hf (ranges of up to 24 ε units) values. Zircon rims are also more homogeneous than the associated cores, and we show that zircon rim growth occurs faster than the resolution of in situ dating techniques. CA-ID-TIMS dating of a subset of zircon grains from the Twin Falls samples reveals complex crystallization histories spanning 10⁴–10⁶ years prior to some eruptions, suggesting that magma genesis was characterized by the cyclic remelting of buried volcanic rocks and intrusions associated with previous magmatic episodes. Age-dependent trends in zircon isotopic compositions show that rhyolite production in the Yellowstone hotspot track is driven by the mixing of mantle-derived melts (normal δ¹⁸O and ε_Hf) and a combination of Precambrian basement rock (normal δ¹⁸O and ε_Hf down to ??60) and shallow Mesozoic and Cenozoic age rocks, some of which are hydrothermally altered (to low δ¹⁸O values) by earlier stages of Snake River Plain magmatism. These crustal melts hybridize with juvenile basalts and rhyolites to produce the erupted rhyolites. We also observe that the Precambrian basement rock is only an important component in the erupted magmas in the first eruption at each caldera center, suggesting that the accumulation of new intrusions quickly builds an upper crustal intrusive body which is isolated from the Precambrian basement and evolves towards more isotopically juvenile and lower-δ¹⁸O compositions over time.     

Severe soil erosion during a 3‐day exceptional rainfall event: combining modelling and field data for a fallow cereal field

Exceptional rainfall events cause significant losses of soil, although few studies have addressed the validation of model predictions at field scale during severe erosive episodes. In this study, we evaluate the predictive ability of the enhanced Soil Erosion and Redistribution Tool (SERT‐2014) model for mapping and quantifying soil erosion during the exceptional rainfall event (~235 mm) that affected the Central Spanish Pyrenees in October 2012. The capacity of the simulation model is evaluated in a fallow cereal field (1.9 ha) at a high spatial scale (1 × 1 m). Validation was performed with field‐quantified rates of soil loss in the rills and ephemeral gullies and also with a detailed map of soil redistribution. The SERT‐2014 model was run for the six rainfall sub‐events that made up the exceptional event, simulating the different hydrological responses of soils with maximum runoff depths ranging between 40 and 1017 mm. Predicted average and maximum soil erosion was 11 and 117 Mg ha^?1 event^?1, respectively. Total soil loss and sediment yield to the La Reina gully amounted to 16.3 and 9.0 Mg event^?1. These rates are in agreement with field estimations of soil loss of 20.0 Mg event^?1. Most soil loss (86%) occurred during the first sub‐event. Although soil accumulation was overestimated in the first sub‐event because of the large amount of detached soil, the enhanced SERT‐2014 model successfully predicted the different spatial patterns and values of soil redistribution for each sub‐event. Further research should focus on stream transport capacity. Copyright © 2014 John Wiley & Sons, Ltd.     

Effect of galvanic distortions on the series and parallel magnetotelluric impedances and comparison with other responses

The series and parallel impedances of the magnetotelluric tensor are appraised in relation to their relative immunity to galvanic electric distortions. The distorted responses are modeled using the Groom-Bailey decomposition of the tensor in terms of twist, shear, statics and strike direction. These four parameters and the undistorted responses are normally considered as unknowns, and are obtained from field data through the solution of an inverse problem. In the present work we use the decomposition as a forward model to simulate distorted sounding curves. Starting with undistorted 2-d te and tm responses, the tensor is distorted by assuming arbitrary values of twist, shear, static and strike direction. By default, both series and parallel responses are immune to the strike direction because they are invariants under rotation. In addition, series responses are immune to twist and shear and parallel responses only to twist. The dependence of the latter on shear is in the form of a real factor that shifts downwards the amplitude curves. On the other hand, the effect of statics on both series and parallel responses is more complicated than that on the impedance tensor because it cannot be accounted for by a simple shift of the curves. On the whole, there is a positive balance on the part of the series and parallel impedances over the te and tm responses because some of the distortions are filtered out by the invariants. It is shown that invariance is not sufficient to be immune to any of the distortions. The example chosen is Eggers’ eigenvalues, which are immune only to the by-the-fault strike direction. Invariance is not necessary either, as evidenced by the phase tensor, whose elements depend on strike but are immune to all distortions. The derivations are illustrated using soundings from the synthetic coprod2s1 and field-recorded coprod2 and bc87 data sets.     

LU Vel (GJ 375): A M3.5Ve Flare and Double-Lined Spectroscopic Binary

High resolution echelle spectroscopic observations taken with the FEROS spectrograph at the 2.2 m telescope ESO confirm the binary nature of the flare M3.5V star LU Vel (GJ 375, RE J0958-462) previously reported by Christian and Mathioudakis (2002). Emission of similar intensity from both components is detected in the Balmer, Na i D₁&D₂, He i D₃, Ca ii H&K, and Ca ii IRT lines. We have determined precise radial velocities by cross correlation with radial velocity standard stars, which have allowed us to obtain for the first time the orbital solution of the system. The binary consists of two near-equal M3.5V components with an orbital period shorter than 2 days. We have analyzed the behaviour of the chromospheric activity indicators (variability and possible flares). In addition, we have determined its rotational velocity and kinematics.     

A Practical Procedure for the Measurement of Fragmentation by Blasting by Image Analysis

Summary. The operation of a digital image analysis system in a limestone quarry is described. The calibration of the system, required in order to obtain moderately reliable fragmentation values, is done from muckpile sieving data by tuning the image analysis software settings so that the fragmentation curve obtained matches as close as possible the sieving. The sieving data have also been used to extend the fragment size distribution curves measured to sizes below the system’s optical resolution and to process the results in terms of fragmented rock, discounting the material coming from a loose overburden (natural fines) that is cast together with the fragmented rock. Automatic and manual operation modes of the image analysis are compared. The total fragmentation of a blast is obtained from the analysis of twenty photographs; a criterion for the elimination of outlier photographs has been adopted using a robust statistic. The limitations of the measurement system due to sampling, image processing and fines corrections are discussed and the errors estimated whenever possible. An analysis of consistency of the results based on the known amount of natural fines is made. Blasts with large differences in the amount of fines require a differentiated treatment, as the fine sizes tend to be the more underestimated in the image analysis as they become more abundant; this has been accomplished by means of a variable fines adjustment factor. Despite of the unavoidable errors and the large dispersion always associated with large-scale rock blasting data, the system is sensitive to relative changes in fragmentation.     

A Parametric Approach for Dealing with Compositional Rounded Zeros

In this work, a parametric approach for replacing data below the detection limit, also known as rounded zeros, in compositional data sets is proposed. Compositional rounded zeros correspond to small proportions of some whole that cannot be reliably detected by the analytical instruments under given operating conditions. This kind of zeros appear frequently in the data collection process in geosciences. They must be treated in an adequate way before some multivariate analysis can be applied. Our procedure results from a modification of the Expectation-Maximization (EM) algorithm and is based on the additive log-ratio transformation. Its coherence with the nature of compositional data and with basic operations in the simplex sample space is checked. Using real data sets, we find that this approach improves other parametric and non-parametric techniques for compositional rounded zeros.     

Exact cosmological solutions in Brans and Dicke's scalar-tensor theory,I

The exact solution is sought for the cosmological equations of Brans and Dicke's scalar-tensor theory when a power law exists between the gravitational constant and the radius of curvature of the universe. For the space of negative curvature no solution is possible. On the contrary for a closed space the gravitational constant and the radius of curvature increase linearly with respect to the age of the universe. The parameter of the scalar-tensor theory is necessarily negative and can be determined by the present values of the mass-density of the universe, the Hubble-constant and the gravitational constant. The solution has no analogy in Einstein's theory with vanishing cosmological constant, even when the deviations from Einstein's values of the solar relativistic effects are small.     

Analyzing ice dynamics using Sentinel-1 data at the Solheimajoküll Glacier,Iceland

ABSTRACT

The climate in southern Iceland has warmed over the last 70 years, resulting in accelerated glacier dynamics at the Solheimajoküll glacier. In this study, we compare glacier terminus locations from 1973 to 2018, to changes in climate across the study area, and we derive ice-surface velocities (2015–2018) from satellite remote-sensing imagery (Sentinel-1) using the offset-tracking method. There have been two regional temperature trends in the study period: cooling (1973–1979) and warming (1980–2018). Our results indicate a time lag of about 20 years between the onset of glacier retreat (?53 m/year since 2000) and the inception of the warming period. Seasonally, the velocity time series suggest acceleration during the summer melt season since 2016, whereas glacier velocities during accumulation months were constant. The highest velocities were observed at high elevations where the ice-surface slope is the steepest. We tested several scenarios to assess the hydrological time response to glacier accelerations, with the highest correlations being found between one and 30 days after the velocity estimates. Monthly correlation analyses indicated inter-annual and intra-annual variability in the glacier dynamics. Additionally, we investigate the linkage between glacier velocities and meltwater outflow parameters as they provide useful information about internal processes in the glacier. Velocity estimates positively correlate with water level and negatively correlate with water conductivity between April and August. There is also a disruption in the correlation trend between water conductivity and ice velocity in June, potentially due to a seasonal release of geothermal water.