Subduction factory processes beneath the Guguan cross-chain, Mariana Arc: no role for sediments, are serpentinites important?

We need to understand chemical recycling at convergent margins and how chemical interactions between subducted slab and the overlying mantle wedge affect mantle evolution and magmagenesis. This requires distinguishing contributions from recycled individual subducted components as well as those contributed by the mantle. We do this by examining magmatic products generated at different depths above a subduction zone, in an intra-oceanic arc setting. The Guguan cross-chain in the intra-oceanic Mariana arc overlies subducted Jurassic Pacific plate lithosphere at depths of ~125--230 km and erupts mostly basalt. Basalts from rear-arc volcanoes are more primitive than those from the magmatic front, in spite of being derived by lower degrees of melting of less-depleted mantle. Rear-arc magmas also show higher temperatures and pressures of equilibration. Coexisting mineral compositions become more MORB- or OIB-like with increasing height above the subduction zone. Trace element and isotopic variations indicate that the subduction component in cross-chain lavas diminishes with increasing depth to the subduction zone, except for water contents. There is little support for the idea that melting beneath the Mariana Trough back-arc basin depleted the source region of arc magmas, but melting to form rear-arc volcanoes may have depleted the source of magmatic front volcanoes. Enrichments in rear-arc lavas were not caused by sediment melting; the data instead favor an OIB-type mantle that has been modestly affected by subduction zone fluids. Our most important conclusion is that sediment fluids or melts are not responsible for the K--h relationship and other cross-chain chemical and isotopic variations. We speculate that an increasing role for supercritical fluids released from serpentinites interacting with modestly enriched mantle might be responsible for cross-chain geochemical and isotopic variations. Electronic Supplementary Material Supplementary material is available for this article at and is accessible for authorized users.     

Enriched back-arc basin basalts from the northern Mariana Trough: implications for the magmatic evolution of back-arc basins

The composition of basalts erupted at the earliest stages in the evolution of a back-arc basin permit unique insights into the composition and structure of the sub-arc mantle. We report major and trace element chemical data and O-, Sr-, Nd-, and Pb- isotopic analyses for basalts recovered from four dredge hauls and one ALVIN dive in the northern Mariana Trough near 22°N. The petrography and major element chemistry of these basalts (MTB-22) are similar to tholeiites from the widest part of the Trough, near 18°N (MTB-18), except that MTB-22 have slightly more K₂O and slightly less TiO₂. The trace element data exhibit a very strong arc signature in MTB-22, including elevated K, Rb, Sr, Ba, and LREE contents; relatively lowK/Ba and highBa/La andSr/Nd. The Sr- and Nd- isotopic data plot in a field displaced from that of MTB-18 towards Mariana arc lavas, and the Pb-isotopic composition of MTB-22 is indistinguishable from Mariana arc lavas and much more homogeneous than MTB-18. Mixing of 50–90% Mariana arc component with a MORB component is hypothesized. We cannot determine whether this resulted from physical mixing of arc mantle and MORB mantle, or whether the arc component is introduced by metasomatism of MORB-like mantle by fluids released from the subducted lithosphere. The strong arc signature in back-arc melts from the Mariana Trough at 22°N, where the back-arc basin is narrow, supports general models for back-arc basin evolution whereby early back-arc basin basalts have a strong arc component which diminishes in importance relative to MORB as the back-arc basin widens.     

Sidewall effects and sand trap efficiency in a large wind tunnel

Aeolian sand transport is a widespread physical phenomenon on the surface of Earth, as well as on Mars and Titan. Accurate measurements of the components of the transport system are necessary if we are to understand the nature of the physical processes. Sand traps are typically used to measure sediment transport rates, and issues associated with the sampling efficiency of traps and the development of reliable traps have received considerable attention in recent decades. In this study, we measured aeolian transport rate at five distances from a wind tunnel sidewall using a vertically‐segmented sand trap. Total transport rates were determined by weighing the bed sediment before and after each experiment, and with and without a trap installed. The following results were obtained: (1) sand transport increased linearly with the distance away from the sidewall, and the appropriate location to measure maximum transport is within the central 20% of the wind tunnel; (2) current methods overestimate the sampling efficiency of sand traps when comparing trap data to transport rate data obtained by weighing sand moved through the entire tunnel because the effects of the sidewalls in decreasing total transport are neglected; (3) the efficiency of the vertically‐segmented trap that we tested ranged from 11.57% to 31.68% using our revised methods, whereas standard methods caused efficiency to be overestimated by 32–72% of the efficiency; (4) using either method, the efficiency of the trap increased exponentially with shear velocity for the range we used. Copyright © 2017 John Wiley & Sons, Ltd.     

Waves, currents, sediment flux and morphological response in a barred nearshore system

A shore-normal array of seven, bi-directional electromagnetic flowmeters and nine surface piercing, continuous resistance wave staffs were deployed across a multiple barred nearshore at Wendake Beach, Georgian Bay, Canada, and monitored for a complete storm cycle. Time-integrated estimates of total (ITVF) and net (INVF) sediment volume flux together with bed elevation changes were determined using depth-of-activity rods.

The three bars, ranging in height from 0.10 to 0.40 m accreted during the storm (0.03 m), and the troughs were scoured (0.05 m). Sediment reactivation depths reached 0.14 m and 12% of the nearshore control volume was mobilized. However, the INVF value for the storm was less than 1% of the control volume revealing a near balance in sediment volume in the bar system. Landward migration of the inner, crescentic and second, sinuous bars occurred in association with an alongshore migration of the bar form itself; the outermost, straight, shore-parallel bar remained fixed in location.

The surf zone was highly dissipative throughout the storm (ε = 3.8 × 10²–192 × 10²) and the wave spectrum was dominated by energy at the incident frequency. Spectral peaks at frequencies of the first harmonic and at one quarter that of the incident wave were associated with secondary wave generation just prior to breaking and a standing edge wave, respectively. The former spectral peak was within the 95% confidence band for the spectrum while the latter contributed not more than 10% to the total energy in the surface elevation spectrum even near the shoreline.

During the storm wave height exceeded 2 m (H_s) and periods reached 5 s (T_{p k}): orbital velocities exceeded 0.5 m s⁻¹ (u_{rm s}) and were above the threshold of motion for the medium-to-fine sands throughout the storm. Shore-parallel flows in excess of 0.4 m s⁻¹ were recorded with maxima in the troughs and minima just landward of the bar crest.

The rate and direction of sediment flux is best explained by the interaction of antecedent bed slopes with spatial gradients in the mean and asymmetry of the shore-normal velocity field. These hydrodynamic parameters represent “steady” flows superimposed on the dominantly oscillatory motion and assumed a characteristic spatial pattern from the storm peak through the decay period. Increases spatially in the magnitudes of both the mean flows and flow asymmetries cause an increasing net transport potential (erosion); decreases in these values spatially cause a decreasing net transport potential and thus deposition. These transport potentials are increased or decreased through the gravity potential induced by the local bed slope. Shore-parallel flow was important in explaining sediment flux and morphological change where orbital velocities, mean flows and flow asymmetries were at a minimum.     

Ion association in concentrated NaCl brines from ambient to supercritical conditions: results from classical molecular dynamics simulations

Highly concentrated NaCl brines are important geothermal fluids; chloride complexation of metals in such brines increases the solubility of minerals and plays a fundamental role in the genesis of hydrothermal ore deposits. There is experimental evidence that the molecular nature of the NaCl–water system changes over the pressure–temperature range of the Earth's crust. A transition of concentrated NaCl–H₂O brines to a "hydrous molten salt" at high P and T has been argued to stabilize an aqueous fluid phase in the deep crust.     

The genesis of seismic activity on faults in Central Asia in real time and its variations

We used the Digital Faults geoinformation system that we developed to propose an algorithm for quantitative estimation of seismic activity on faults. The resulting technique was used to study the spatiotemporal patterns in the present-day activity of faults in Central Asia. Fault activity was found to vary at frequencies of a few years and cannot be explained by changes in the regional stress fields. We studied the tendency of seismic events to be localized to areas of dynamic influence due to faults. The active faults were grouped by the criteria of seismicity organization in the influence areas of these faults. It was shown that fault activity and its comparatively high frequency on real time scales are caused by strain waves, which may be generated by interplate and interblock movements in the brittle lithosphere. Judging by the speed of strain waves, the active faults are classified into groups that differ in their geological and geophysical parameters. They can be used to estimate the directions of strain wave fronts and to identify areas of dominant fault activation over intervals of real (geologically speaking) time. We give a map showing active faults in Central Asia, plots of a quantitative index of their seismic activity, and the directivity vectors of strain waves that excite fault activity. The methods we developed for classifying active faults by the quantitative index of seismic activity and for determining the vectors of strain waves that excite fault activity are all tools that significantly expand our possibilities when developing tectonophysical models of the seismic process in earthquake-generating zones of the lithosphere and open new methods for attacking problems in intermediate-term earthquake prediction.     

Influence of averaging interval on shear velocity estimates for aeolian transport modeling

Recent investigations of aeolian transport have focused on increasingly short time scales because of growing recognition that wind unsteadiness is a major factor in the dynamics of sediment transport. However, the statistical reliability of shear velocity (u_*) estimates becomes increasingly uncertain as averaging interval is decreased. This study provides an empirical assessment of the influence of averaging interval on the reliability of u_* estimates. The data consist of 15-min wind-speed profiles (1 Hz sampling) collected at four coastal sites. Each profile was subdivided into progressively shorter fixed-length time intervals, and estimates of u_* and the 95% confidence interval for u_* were determined for each time-block using standard statistical techniques.The logarithmic model accurately represents the measured wind-speed profiles, even with relatively brief averaging intervals. Mean r² values remain robust down to block lengths as short as 10–20 s, typically retaining better than 98% of the r² value found for the full-length data sets. Fewer than 2% of the individual 10-s blocks had r² values less than 0.9. However, mean confidence intervals typically expanded by 70–80% of the full-record value as block length decreased from 900 to 10 s. For highly log-linear profiles, this amounted to an absolute increase from about ±8% to only ±14% of u_*, so that the additional information gained through the use of shorter averaging intervals may outweigh the increase in statistical uncertainty. Nevertheless, given that rates of aeolian transport are generally modeled as a function of u_*³, this increase in uncertainty may be significant for transport modeling. Thus, very short averaging intervals should be used with caution when predicting aeolian sediment flux. It is proposed that transport modeling should incorporate the shear velocity confidence interval as an indicator of the potential error associated with this source of uncertainty.     

Aquila X-1: a low-inclination soft X-ray transient

We have obtained I -band photometry of the neutron star X-ray transient Aql X-1 during quiescence. We find a periodicity at 2.487 cycles d⁻¹, which we interpret as twice the orbital frequency (19.30±0.05 h). Folding the data on the orbital period, we model the light-curve variations as the ellipsoidal modulation of the secondary star. We determine the binary inclination to be 20°–30° (90 per cent confidence) and also determine the 95 per cent upper limits to the radial velocity semi-amplitude and rotational broadening of the secondary star to be 117 and 50 km s⁻¹, respectively.