ALBERT

Description: Summary On 24th August 2016 at 01:36 UTC a ML6.0 earthquake struck several villages in central Italy, among which Accumoli, Amatrice and Arquata del Tronto. The earthquake was recorded by about 350 seismic stations, causing 299 fatalities and damage with macroseismic intensities up to 11. The maximum acceleration was observed at Amatrice station (AMT) reaching 916 cm/s2 on E-W component, with epicentral distance of 15 km and Joyner and Boore distance to the fault surface (RJB) of less than a kilometre. Motivated by the high levels of observed ground motion and damage, we generate broadband seismograms for engineering purposes by adopting a hybrid method. To infer the low frequency seismograms, we considered the kinematic slip model by Tinti et al. (2016). The high frequency seismograms were produced using a stochastic finite-fault model approach based on dynamic corner-frequency. Broadband synthetic time series were therefore obtained by merging the low and high frequency seismograms. Simulated hybrid ground motions were compared both with the observed ground motions and the ground-motion prediction equations (GMPEs), to explore their performance and to retrieve the region-specific parameters endorsed for the simulations. In the near-fault area we observed that hybrid simulations have a higher capability to detect near source effects and to reproduce the source complexity than the use of GMPEs. Indeed, the general good consistency found between synthetic and observed ground motion (both in the time and frequency domain), suggests that the use of regional-specific source scaling and attenuation parameters together with the source complexity in hybrid simulations improves ground motion estimations. To include the site effect in stochastic simulations at selected stations, we tested the use of amplification curves derived from HVRSs (horizontal-to-vertical response spectra) and from HVSRs (horizontal-to-vertical spectral ratios) rather than the use of generic curves according to NTC-18 Italian seismic design code. We generally found a further reduction of residuals between observed and simulated both in terms of time histories and spectra.

Description: Summary We explore here the benefits of using constraints from seismic tomography in gravity data inversion and how inverted density distributions can be improved by doing so. The methodology is applied to a real field case in which we reconstruct the density structure of the Pyrenees along a southwest-northeast transect going from the Ebro basin in Spain to the Arzacq basin in France. We recover the distribution of densities by inverting gravity anomalies under constraints coming from seismic tomography. We initiate the inversion from a prior density model obtained by scaling a pre-existing compressional seismic velocity Vp model using a Nafe-Drake relationship : the Vp model resulting from a full-waveform inversion of teleseismic data. Gravity data inversions enforce structural similarities between Vp and density by minimizing the norm of the cross-gradient between the density and Vp models. We also compare models obtained from 2.5D and 3D inversions. Our results demonstrate that structural constraints allow us to better recover the density contrasts close to the surface and at depth, without degrading the gravity data misfit. The final density model provides valuable information on the geological structures and on the thermal state and composition of the western region of the Pyrenean lithosphere.

Description: Summary The relatively short duration of the early stages of subduction results in a poor geological record, limiting our understanding of this critical stage. Here, we utilize a 2D numerical model of incipient subduction, that is the stage after a plate margin has formed with a slab tip that extends to a shallow depth and address the conditions under which subduction continues or fails. We assess energy budgets during the evolution from incipient subduction to either a failed or successful state, showing how the growth of potential energy, and slab pull, is resisted by the viscous dissipation within the lithosphere and the mantle. The role of rheology is also investigated, as deformation mechanisms operating in the crust and mantle facilitate subduction. In all models, the onset of subduction is characterized by high lithospheric viscous dissipation and low convergence velocities, whilst successful subduction sees the mantle become the main area of viscous dissipation. In contrast, failed subduction is defined by the lithospheric viscous dissipation exceeding the lithospheric potential energy release rate and velocities tend towards zero. We show that development of a subduction zone depends on the convergence rate, required to overcome thermal diffusion and to localise deformation along the margin. The results propose a minimum convergence rate of ∼ 0.5 cm yr−1 is required to reach a successful state, with 100 km of convergence over 20 Myr, emphasizing the critical role of the incipient stage.

Description: Summary Although many studies have revealed that the atmospheric effects of electromagnetic wave propagation (including ionospheric and tropospheric water vapor) have serious impacts on Interferometric Synthetic Aperture Radar (InSAR) measurement results, atmospheric corrections have not been thoroughly and comprehensively investigated in many well-known cases of InSAR focal mechanism solutions, which means there is no consensus on whether atmospheric effects will affect the InSAR focal mechanism solution. Moreover, there is a lack of quantitative assessment on how much the atmospheric effect affects the InSAR focal mechanism solution. In this paper, we emphasized that it was particularly important to assess the impact of InSAR ionospheric and tropospheric corrections on the underground nuclear explosion modeling quantitatively. Therefore, we investigated the 4th North Korea (NKT-4) underground nuclear test using ALOS-2 liters-band SAR images. Because the process of the underground nuclear explosion was similar to the volcanic magma source activity, we modeled the ground displacement using the Mogi model. Both the ionospheric and tropospheric phase delays in the interferograms were investigated. Furthermore, we studied how the ionosphere and troposphere phase delays could bias the estimation of Mogi source parameters. The following conclusions were drawn from our case study: The ionospheric delay correction effectively mitigated the long-scale phase ramp in the full-frame interferogram, the standard deviation decreased from 1.83 cm to 0.85 cm compared to the uncorrected interferogram. The uncorrected estimations of yield and depth were 8.44 kt and 370.33 m, respectively. Compared to the uncorrected estimations, the ionospheric correction increased the estimation of yield and depth to 9.43 kt and 385.48 m while the tropospheric correction slightly raised them to 8.78 kt and 377.24 m. There were no obvious differences in the location estimations among the four interferograms. When both corrections were applied, the overall standard deviation was 1.16 cm, which was even larger than the ionospheric corrected interferogram. We reported the source characteristics of NKT-4 based on the modeling results derived from the ionospheric corrected interferogram. The preferred estimation of NKT-4 was a Mogi source located at 129°04′22.35‘E, 41°17′54.57″N buried at 385.48 m depth. The cavity radius caused by the underground explosion was 22.02 m. We reported the yield estimation to be 9.43 kt. This study showed that for large-scale natural deformation sources such as volcanoes and earthquakes, atmospheric corrections would be more significant, but even if the atmospheric signal did not have much complexity, the corrections should not be ignored.

Description: SUMMARY Self-consistent modelling of magmatic systems is challenging as the melt continuously changes its chemical composition upon crystallization, which may affect the mechanical behaviour of the system. Melt extraction and subsequent crystallization create new rocks while depleting the source region. As the chemistry of the source rocks changes locally due to melt extraction, new calculations of the stable phase assemblages are required to track the rock evolution and the accompanied change in density. As a consequence, a large number of isochemical sections of stable phase assemblages are required to study the evolution of magmatic systems in detail. As the state-of-the-art melting diagrams may depend on nine oxides as well as pressure and temperature, this is a 10-D computational problem. Since computing a single isochemical section (as a function of pressure and temperature) may take several hours, computing new sections of stable phase assemblages during an ongoing geodynamic simulation is currently computationally intractable. One strategy to avoid this problem is to pre-compute these stable phase assemblages and to create a comprehensive database as a hyperdimensional phase diagram, which contains all bulk compositions that may emerge during petro-thermomechanical simulations. Establishing such a database would require repeating geodynamic simulations many times while collecting all requested compositions that may occur during a typical simulation and continuously updating the database until no additional compositions are required. Here, we describe an alternative method that is better suited for implementation on large-scale parallel computers. Our method uses the entries of an existing preliminary database to estimate future required chemical compositions. Bulk compositions are determined within boundaries that are defined manually or through principal component analysis in a parameter space consisting of clustered database entries. We have implemented both methods within a massively parallel computational framework while utilizing the Gibbs free energy minimization program Perple_X. Results show that our autonomous approach increases the resolution of the thermodynamic database in compositional regions that are most likely required for geodynamic models of magmatic systems.

Description: Transformation of refractory cratonic mantle into more fertile lithologies is the key to the fate of cratonic lithosphere. This process has been extensively studied in the eastern North China Craton (NCC) while that of its western part is still poorly constrained. A comprehensive study of newly-found pyroxenite xenoliths from the Langshan area, in the northwestern part of this craton is integrated with a regional synthesis of pyroxenite and peridotite xenoliths to constrain the petrogenesis of the pyroxenites and provide an overview of the processes involved in the modification of the deep lithosphere. The Langshan pyroxenites are of two types, high-Mg# [Mg2+/(Mg2++Fe2+)*100 = ∼ 90, atomic ratios] olivine-bearing websterites with high equilibration temperatures (880 ∼ 970 oC), and low-Mg# (70 ∼ 80) plagioclase-bearing websterites with low equilibration temperatures (550 ∼ 835 oC). The high-Mg# pyroxenites show trade-off abundances of olivine and orthopyroxene, highly depleted bulk Sr-Nd (ƐNd = +11.41, 87Sr/86Sr = ∼0.7034) and low clinopyroxene Sr isotopic ratios (mean 87Sr/86Sr = ∼0.703). They are considered to reflect the reaction of mantle peridotites with silica-rich silicate melts derived from the convective mantle. Their depletion in fusible components (e.g., FeO, TiO2 and Na2O) and progressive exhaustion of incompatible elements suggest melt extraction after their formation. The low-Mg# pyroxenites display layered structures, convex-upward rare earth element patterns, moderately enriched bulk Sr-Nd isotopic ratios (ƐNd = -14.20 ∼ -16.74, 87Sr/86Sr = 0.7070 ∼ 0.7078) and variable clinopyroxene Sr-isotope ratios (87Sr/86Sr = 0.706-0.711). They are interpreted to be crustal cumulates from hypersthene-normative melts generated by interaction between the asthenosphere and heterogeneous lithospheric mantle. Combined with studies on regional peridotite xenoliths, it is shown that the thinning and refertilization of the lithospheric mantle was accompanied by crustal rejuvenation and that such processes occurred ubiquitously in the northwestern part of the NCC. A geodynamic model is proposed for the evolution of the deep lithosphere, which includes long-term mass transfer through a mantle wedge into the deep crust from the Paleozoic to the Cenozoic, triggered by subduction of the Paleo-Asian ocean and the Late Mesozoic lithospheric extension of eastern Asia.

Description: Fishermen are known to try to avoid fishing in stormy weather, as storms pose a physical threat to fishers, their vessels, and their gear. In this article, a dataset and methods are developed to investigate the degree to which fishers avoid storms, estimate storm aversion parameters, and explore how this response varies across vessel characteristics and across regions of the United States. The data consist of vessel-level trip-taking decisions from six federal fisheries across the United States combined with marine storm warning data from the National Weather Service. The estimates of storm aversion can be used to parameterize predictive models. Fishers’ aversion to storms decreases with increasing vessel size and increases with the severity of the storm warning. This information contributes to our understanding of the risk-to-revenue trade-off that fishers evaluate every time they consider going to sea, and of the propensity of fishers to take adaptive actions to avoid facing additional physical risk.

Description: Sardine Sardinops sagax is an ecologically and economically important Clupeid found off the entire South African coast that includes both coastal upwelling and western boundary current systems. Although the management of the sardine fisheries historically assumed a single, panmictic population, the existence of three, semi-discrete subpopulations has recently been hypothesized. We conducted otolith δ18O and microstructure analyses to investigate nursery habitat temperatures and early life growth rates, respectively, of sardine collected from three biogeographic regions around South Africa’s coast to test that hypothesis. Analyses indicated that for both summer- and winter-captured adults and summer-captured juveniles, fishes from the west coast grew significantly slower in water that was several degrees cooler than those from the south and east coasts. This suggests that mixing of sardines between regions, particularly the west and other coasts, is relatively limited and supports the hypothesis of semi-discrete subpopulations. However, the west-south differences disappeared in the results for winter-captured juveniles, suggesting that differences in early life conditions between regions may change seasonally, and/or that all or most winter-captured juveniles originated from the west coast. Further elucidating the interactions between South African sardine subpopulations and the mechanisms thereof is important for sustainable harvesting of this species.

Description: The orogenic development after the continental collision between Laurussia and Gondwana, led to two contrasting associations of mantle-derived magmatic rocks on the territory of the Bohemian Massif: (i) a 340–310 Ma lamprophyre-lamproite orogenic association and (ii) a 300–275 Ma lamprophyre association of anorogenic affinity. Major types of potassic mantle-derived magmatic rocks recognised in the orogenic and anorogenic associations include: (i) calc-alkaline to alkaline lamprophyres, (ii) alkaline “orthopyroxene minettes” (and geochemically related rocks), and (iii) peralkaline lamproites. These three types significantly differ with respect to mineral, whole-rock and Sr–Nd–Pb–Li isotope composition, and spatial distribution. The calc-alkaline lamprophyres occur throughout the entire Saxo-Thuringian and Moldanubian zones, whereas the different types of malte-derived potassic rocks are spatially restricted to particular zones. Rocks of the Carboniferous lamprophyre-lamproite orogenic association are characterised by variable negative εNd(i) and variably radiogenic Sr(i), whereas the rocks of the Permian lamprophyre association of anorogenic affinity are characterised by positive εNd(i) and relatively young depleted-mantle Nd-model ages reflecting increasing input from upwelling asthenospheric mantle. The small variation in the Pb isotopic composition of post-collisional potassic mantle-derived magmatic rocks (of both the orogenic and anorogenic series) implies that the Pb budget of the mantle beneath the Bohemian Massif is dominated by the same crust-derived material, which itself may include material derived from several sources. The source rocks of “orthopyroxene minettes” are characterised by isotopically light (“eclogitic”) Li and strongly radiogenic (crustal) Sr and may have been metasomatised by high-pressure fluids along the edge of a subduction zone. In contrast, the strongly Al2O3 and CaO depleted mantle source of the lamproites is characterised by isotopically heavy Li and high SiO2 and extreme K2O contents. This mantle source may have been metasomatised predominantly by melts. The mantle source of the lamprophyres may have undergone metasomatism by both fluids and melts.

Description: SUMMARY We present a statistical rock physics inversion of the elastic and electrical properties to estimate the petrophysical properties and quantify the associated uncertainty. The inversion method combines statistical rock physics modeling with Bayesian inverse theory. The model variables of interest are porosity and fluid saturations. The rock physics model includes the elastic and electrical components and can be applied to the results of seismic and electromagnetic inversion. To describe the non-Gaussian behaviour of the model properties, we adopt non-parametric probability density functions to sample multimodal and skewed distributions of the model variables. Different from machine learning approach, the proposed method is not completely data-driven but is based on a statistical rock physics model to link the model parameters to the data. The proposed method provides pointwise posterior distributions of the porosity and CO2 saturation along with the most-likely models and the associated uncertainty. The method is validated using synthetic and real data acquired for CO2 sequestration studies in different formations: the Rock Springs Uplift in Southwestern Wyoming and the Johansen formation in the North Sea, offshore Norway. The proposed approach is validated under different noise conditions and compared to traditional parametric approaches based on Gaussian assumptions. The results show that the proposed method provides an accurate inversion framework where instead of fitting the relationship between the model and the data, we account for the uncertainty in the rock physics model.

Description: SUMMARY The horizontal-to-vertical spectral ratio (HVSR) of ambient noise is commonly used to infer a site's resonance frequency (${f_{0,site}}$). HVSR calculations are performed most commonly using the Fourier amplitude spectrum obtained from a single merged horizontal component (e.g. the geometric mean component) from a three-component sensor. However, the use of a single merged horizontal component implicitly relies on the assumptions of azimuthally isotropic seismic noise and 1-D surface and subsurface conditions. These assumptions may not be justified at many sites, leading to azimuthal variability in HVSR measurements that cannot be accounted for using a single merged component. This paper proposes a new statistical method to account for azimuthal variability in the peak frequency of HVSR curves (${f_{0,HVSR}}$). The method uses rotated horizontal components at evenly distributed azimuthal intervals to investigate and quantify azimuthal variability. To ensure unbiased statistics for ${f_{0,HVSR}}$ are obtained, a frequency-domain window-rejection algorithm is applied at each azimuth to automatically remove contaminated time windows in which the ${f_{0,HVSR}}$ values are statistical outliers relative to those obtained from the majority of windows at that azimuth. Then, a weighting scheme is used to account for different numbers of accepted time windows at each azimuth. The new method is applied to a data set of 114 HVSR measurements with significant azimuthal variability in ${f_{0,HVSR}}$, and is shown to reliably account for this variability. The methodology is also extended to the estimation of a complete lognormal-median HVSR curve that accounts for azimuthal variability. To encourage the adoption of this statistical approach to accounting for azimuthal variability in single-station HVSR measurements, the methods presented in this paper have been incorporated into hvsrpy, an open-source Python package for HVSR processing.

Description: SUMMARY This study focuses in the analysis of the internal structure of the upper 3 km of Los Humeros (LH) caldera and the relation of electrical and hydrothermal anomalies. For this purpose, we measured, processed and interpreted 78 broad-band magnetotelluric (MT) soundings. We performed a 3-D inversion of the data set (ModEM) using all MT soundings, although only half of the available frequencies per sounding due to limited computed power. We also carried out the 2-D inversions (NLCG) of the invariant determinant along two orthogonal profiles (EW and NS) crossing the caldera structure; their comparison yields similar resistivity and structural models results. The resistivity modelling is complemented with the results of a joint 3-D inversion of an accurate gravity database of 720 stations, and total field aeromagnetic data (SGM) from the caldera crater. The combined results provide novel details about the structure of the shallow geothermal reservoir of the resurgence caldera complex hosting the active hydrothermal system. Density and resistivity models show the existence of a composed crater basin structure separated by an EW high-density structure; the northern basin is associated to the LH crater, whereas the southern basin associates to the emergent Los Potreros (LP) caldera basin. The magnetization model indicates that there is a common source for the magnetic volcanic products observed at the caldera surface, and that the LP fault is the more magnetized fault of the geothermal system. The propylic zoning under the geothermal field, which according to the MT model results has resistivities above ∼100 Ω-m, was extrapolated using this and additional criteria to obtain the distribution of other hypothetical propylitic zones of hydrothermal potential.

Description: SUMMARY For the time stationary global geomagnetic field, a new modelling concept is presented. A Bayesian non-parametric approach provides realistic location dependent uncertainty estimates. Modelling related variabilities are dealt with systematically by making little subjective a priori assumptions. Rather than parametrizing the model by Gauss coefficients, a functional analytic approach is applied. The geomagnetic potential is assumed a Gaussian process to describe a distribution over functions. A priori correlations are given by an explicit kernel function with non-informative dipole contribution. A refined modelling strategy is proposed that accommodates non-linearities of archeomagnetic observables: First, a rough field estimate is obtained considering only sites that provide full field vector records. Subsequently, this estimate supports the linearization that incorporates the remaining incomplete records. The comparison of results for the archeomagnetic field over the past 1000 yr is in general agreement with previous models while improved model uncertainty estimates are provided.

Description: SUMMARY Rapid development of time-lapse seismic monitoring instrumentations has made it possible to collect dense time-lapse data for tomographically retrieving time-lapse (even continuous) images of subsurface changes. While traditional time-lapse full waveform inversion (TLFWI) algorithms are designed for sparse time-lapse surveys, they lack of effective temporal constraint on time-lapse data, and, more importantly, lack of the uncertainty estimation of the TLFWI results that is critical for further interpretation. Here, we propose a new data assimilation TLFWI method, using hierarchical matrix powered extended Kalman filter (HiEKF) to quantify the image uncertainty. Compared to existing Kalman filter algorithms, HiEKF allows to store and update a data-sparse representation of the cross-covariance matrices and propagate model errors without expensive operations involving covariance matrices. Hence, HiEKF is computationally efficient and applicable to 3-D TLFWI problems. Then, we reformulate TLFWI in the framework of HiEKF (termed hereafter as TLFWI-HiEKF) to predict time-lapse images of subsurface spatiotemporal velocity changes and simultaneously quantify the uncertainty of the inverted velocity changes over time. We demonstrate the validity and applicability of TLFWI–HiEKF with two realistic CO2 monitoring models derived from Frio-II and Cranfield CO2 injection sites, respectively. In both 2-D and 3-D examples, the inverted high-resolution time-lapse velocity results clearly reveal a continuous velocity reduction due to the injection of CO2. Moreover, the accuracy of the model is increasing over time by assimilating more time-lapse data while the standard deviation is decreasing over lapsed time. We expect TLFWI-HiEKF to be equipped with real-time seismic monitoring systems for continuously imaging the distribution of subsurface gas and fluids in the future large-scale CO2 sequestration experiments and reservoir management.

Description: SUMMARY Seven years after the beginning of a massive wastewater injection project in eastern Colombia, local earthquake activity increased significantly. The field operator and the Colombian Geological Survey immediately reinforced the monitoring of the area. Our analysis of the temporal evolution of the seismic and injection data together with our knowledge of the geological parameters of the region indicate that the surge of seismicity is being induced by the re-injection of produced water into the same three producing reservoirs. Earthquake activity began on known faults once disposal rates had reached a threshold of ∼2 × 106 m3 of water per month. The average reservoir pressure had remained constant at 7.6 MPa after several years of production, sustained by a large, active aquifer. Surface injection pressures in the seismically active areas remain below 8.3 MPa, a value large enough to activate some of the faults. Since faults are mapped throughout the region and many do not have seismicity on them, we conclude that the existence of known faults is not the only control on whether earthquakes are generated. Stress conditions of these faults are open to future studies. Earthquakes are primarily found in four clusters, located near faults mapped by the operator. The hypocentres reveal vertical planes with orientations consistent with focal mechanisms of these events. Stress inversion of the focal mechanisms gives a maximum compression in the direction ENE-WSW, which is in agreement with borehole breakout measurements. Since the focal mechanisms of the earthquakes are consistent with the tectonic stress regime, we can conclude that the seismicity is resulting from the activation of critically stressed faults. Slip was progressive and seismic activity reached a peak before declining to few events per month. The decline in seismicity suggests that most of the stress has been relieved on the main faults. The magnitude of a large majority of the recorded earthquakes was lower than 4, as the pore pressure disturbance did not reach the mapped large faults whose activation might have resulted in larger magnitude earthquakes. Our study shows that a good knowledge of the local fault network and conditions of stress is of paramount importance when planning a massive water disposal program. These earthquakes indicate that while faults provide an opportunity to dispose produced water at an economically attractive volume–pressure ratio, the possibility of induced seismicity must also be considered.

Description: To get the best result for seismic imaging using primary reflections, data with densely-spaced sources and receivers are ideally preferred. However, dense acquisition can sometimes be hindered by various obstacles, like platforms or complex topography. Such areas with large data gaps may deter exploration or monitoring, as conventional imaging strategies would either provide poor seismic images or turn out to be very expensive. Surface-related multiples travel along different paths compared to primaries, illuminating a wider subsurface area and hence making them valuable in case of data with large gaps. We propose different strategies of using surface-related multiples to get around the problem of imaging in the case of a large data gap. Conventional least-squares imaging methods that incorporate surface-related multiples do so by re-injecting the measured wavefield in the forward-modelling process, which makes it still sensitive to missing data. We introduce a ‘non-linear’ inversion approach in which the surface multiples are modelled from the original source field. This makes the method less dependent on the receiver geometry, therefore, effectively exploiting the information from surface multiples in cases of limited illumination. However, such an approach is sensitive to the knowledge of the source properties. Therefore, we propose a ‘hybrid’ method that combines the non-linear imaging method with the conventional ‘linear’ multiple imaging method, which further improves our imaging result. We test the methods on numerical as well as field data. The results indicate substantial removal of artefacts in the image derived from linear imaging methods due to incomplete data, by exploiting the surface multiples to a maximum extent.

Description: The Lu–Hf isotope system and Sr–Nd–Hf–Os isotope systematics of mantle rocks are capable of unravelling the early processes in collision belts, especially in a hot subduction context where the Sm–Nd and U–Pb systems in crustal rocks are prone to resetting owing to high temperatures and interaction with melts during exhumation. To improve models of the Devonian–Carboniferous evolution of the Bohemian Massif, we investigated in detail mafic and ultramafic rocks (eclogite, pyroxenite, and peridotite) from the ultrahigh-pressure and ultrahigh-temperature Kutná Hora Crystalline Complex (KHCC: Úhrov, Bečváry, Doubrava, and Spačice localities). Petrography, multiphase solid inclusions, major and trace element compositions of rocks and minerals, and radiogenic isotopic data document contrasting sources and protoliths as well as effects of subduction-related processes for these rocks. The Úhrov peridotite has a depleted composition corresponding to the suboceanic asthenospheric mantle, whereas Bečváry and Doubrava peridotites represent lithospheric mantle that underwent melt refertilization by basaltic and SiO2-undersaturated melts, respectively. Multiphase solid inclusions enclosed in garnet from Úhrov and Bečváry peridotites represent trapped H2O ± CO2-bearing metasomatizing agents and Fe–Ti-rich melts. The KHCC eclogites either formed by high-pressure crystal accumulation from mantle-derived basaltic melts (Úhrov) or represent a fragment of mid-ocean ridge basalt-like gabbroic cumulate (Spačice) and crustal-derived material (Doubrava) both metamorphosed at high P–T conditions. The Lu–Hf age of 395 ± 23 Ma obtained for the Úhrov peridotite reflects garnet growth related to burial of the asthenospheric mantle during subduction of the oceanic slab. By contrast, Spačice and Doubrava eclogites yield younger Lu–Hf ages of ∼350 and 330 Ma, respectively, representing mixed ages as demonstrated by the strong granulite-facies overprint and trace element zoning in garnet grains. We propose a refined model for the Early Variscan evolution of the Bohemian Massif starting with the subduction of the oceanic crust (Saxothuringian ocean) and associated oceanic asthenospheric mantle (Úhrov) beneath the Teplá–Barrandian at ≥380 Ma, which was responsible for melt refertilization of the associated mantle wedge (Bečváry, Doubrava). This was followed by continental subduction (∼370–360 Ma?) accompanied by the oceanic slab break-off and incorporation of the upwelling asthenospheric mantle into the Moldanubian lithospheric mantle and subsequent coeval exhumation of mantle and crustal rocks at ∼350–330 Ma.

Description: Summary Virtual Deep Seismic Sounding (VDSS) uses the arrival time of post-critical SsPmp relative to the direct S wave to infer Moho depth at the Pmp reflection point. Due to the large offset between the virtual source and the receiver, SsPmp is more sensitive to lateral variations of structures than near-vertical phases such as Ps, which is used to construct conventional P receiver functions. However, the way post-critical SsPmp is affected by lateral variations in lithospheric structure is not well understood, and previous studies largely assumed a 1D structure when analyzing SsPmp waveforms. Here we present synthetic tests with various 2D models to show that lateral variations in lithospheric structures, from the lithosphere-asthenosphere boundary to sedimentary basins, profoundly affect travel time, phase, and amplitude of post-critical SsPmp, and that a 1D approximation is usually inappropriate when analyzing 2D data. Despite these strong effects we show, with synthetic examples and the ChinArray data from the Ordos Block in northern China, that a simple ray-theory-based back-projection method can retrieve the geometry of the crust-mantle boundary given array observations in cases with moderate lateral variations in the crust-mantle boundary and/or the lithosphere-asthenosphere boundary. The success of our back-projection method indicates that ray-theory approximations are sufficient in modeling SsPmp travel times in the presence of moderate lateral heterogeneity. In contrast, we show that the ray theory is generally insufficient in modeling SsPmp phase shifts in a strongly heterogeneous lithosphere due to non-planar down-going P waves incident at the crust-mantle boundary. Nonetheless, our results demonstrate the feasibility of direct imaging of the crust-mantle boundary with post-critical SsPmp even in the presence of 2D variations of lithospheric structure.

Description: Progress towards ecosystem-based fisheries management calls for useful tools to prioritize actions. To select suitable methods for local circumstances, evaluating approaches used in other jurisdictions can be a cost-effective first step. We tested Productivity Susceptibility Analysis (PSA) to assess the potential vulnerability of the marine fish community in the Skagerrak–Kattegat (Eastern North Sea) to possible interactions with all Swedish fisheries operating in the area. This analysis combines attributes for a species productivity with attributes related to the susceptibility to capture to quantify a single score for vulnerability: high, medium, or low risk. Results indicate that demersal trawl and gillnet fisheries were associated with the highest risk levels if interaction occurs, i.e. having the highest prevalence of species with potentially high vulnerability to the fisheries. Mixed results were seen when comparing the assessment results with available data. The main benefit of utilizing PSA in the area is the comprehensiveness of the assessment, including data-deficient fisheries and species. Drawbacks include potential overestimation of actual risks. Overall, together with available data, PSA in the studied area provides a comprehensive map of potential risks for further actions and may progress a science-based, precautionary management of the area.

Description: While fisheries provide food and employment for hundreds of millions of people, they also can have significant impact on biodiversity. We explore the potential of area-based fisheries management to simultaneously maintain biodiversity and high levels of sustainable food production. We used two illustrative examples of fisheries that have different gear types, areas, and species to evaluate the trade-off between biodiversity and harvest. We calculate the optimal effort by gear and area that maximizes a weighted objective function of biodiversity and harvest, ranging from 100% of the weight on harvest to 100% on biodiversity. We found for both case studies that the trade-off was highly convex, with win–win solutions allowing for high levels of both fishery harvest and conservation. This is achieved by reducing or eliminating fishing effort that negatively impacts high conservation value species while maintaining fishing effort with gears and in areas where there is low conservation impact. We suggest that, in most fisheries, such situations can be found and that effective area-based management can provide for high levels of biodiversity protection and food production.

Description: To determine the failure characteristics and evolution regularity of the floor strata during fully mechanised top-coal cave mining in typical deep and extra-thick seams in western China, the 61303 working face of an Ordos mine was selected as the engineering research object. A comprehensive monitoring method combining a BOTDR (Brillouin Optical Time-Domain Reflection) distributed fiber strain test and a borehole resistivity CT (Computerised Tomography) test was adopted. The results show that floor stress distribution of the deep-buried and extra-thick coal seam is significantly affected by the different depths of rock lithology. At the interface of the rock strata with a large difference in the elastic modulus, phenomena such as the asynchronism of strata movement and obvious differences in failure easily occur. The failure depth of the floor strata in the 61303 working face is approximately 15.90 m, and the influence depth of the floor disturbance is approximately 32.70 m. Under the influence of the mining pressure, floor stress distribution and crack evolution have obvious spatial and temporal effects. In different inclined boreholes, the data captured by the cable have different values and the fracture locations of the cable also differ. Compared with a single borehole, multiple boreholes with different inclinations, directions and locations can provide more comprehensive and reliable data trends. The knowledge obtained by this monitoring can provide reference information for the study of floor damage under similar conditions and the formulation of technical measures such as those that prevent mine water disasters.

Description: Summary The Mw 6.0 earthquake that hit central Italy on 24 August 2016 caused an abnormally high level of destruction in the town of Amatrice. In order to clarify the role of site response in causing such a disaster, a series of ambient noise recordings acquired in the aftermath of the event are analysed here to identify site resonance properties from the ratios H/V between horizontal and vertical amplitudes of ground motion. Although the noise data acquisition was limited by the emergency management activities, the use of a new analysis technique, which isolates the contribution of Rayleigh waves to the noise wavefield and averages instantaneous estimates of H/V ratios, provided more stable results compared to the standard Nakamura's technique based on mean spectral ratios. The results demonstrated the occurrence of significant resonance phenomena, but without an obvious correlation with the spatial distribution of damage severity. It is apparent that the damage severity was also influenced by some additional local factors related to building vulnerability. Moreover, the time series analysis revealed seasonal variations in the Rayleigh wave ellipticity curves likely related to the water content changes in the surface deposits and their influence on the Poisson coefficient. Finally, the new method proved capable of recognizing time-varying directions of Rayleigh wave propagation. This capability could be exploited to support other passive seismic methods (e.g. ReMi), whose results reliability is limited by the lack of control on wave direction origin.

Description: SUMMARY This paper explores the effects on tsunami simulations of the level of detail of the bathymetric grid in use. For this purpose, we expand available bathymetric data sets of the Pacific Basin in spherical harmonics. For realistic scenarios of tsunamis generated by earthquake dislocations, we conclude that an expansion to a maximum degree lmax = 40, corresponding to wavelengths of 1000 km, is sufficient to reproduce the main features of the tsunami wavefield synthesized in deep water, that is, without considering final shoaling and interaction with coastal features.

Description: SUMMARY Full-waveform inversion (FWI) has demonstrated increasing success in estimating medium properties, but its computational cost still poses challenges in moving towards high-resolution imaging of targets at depth. Here, we propose a target-oriented FWI method that inverts for the medium parameters confined within an arbitrary region of interest. Our method is novel in terms of both local wavefield modelling and data redatuming, in order to build a target-oriented objective function which is sensitive to the target medium only without further assumptions about the medium outside. Based on the convolution-type representation theorem, our local forward modelling operator propagates wavefields within the target medium only while providing full acoustic coupling between the target medium and the surrounding geology. A key requirement of our local FWI method is that the subsurface wavefields surrounding and inside the target be as accurate as possible. As such, the subsurface wavefields are retrieved by the Marchenko method, which can redatum the single-sided surface reflection data to the target zone while preserving both primary and multiple reflections, with minimal a priori knowledge of the full-domain medium. Given a sufficiently accurate initial velocity macromodel, our numerical examples show that our local FWI method resolves the reservoir zone of a 2-D Barrett Unconventional P-wave velocity model much more efficiently than the conventional full-domain FWI without significantly sacrificing accuracy. Our method may further enable FWI approaches to high-resolution imaging of subsurface targets.

Description: The Changning–Menglian orogenic belt (CMOB) in the southeastern Tibetan Plateau is an important link between the Longmu Co–Shuanghu suture (LCSS) in the northern Tibetan Plateau and the Chiang Mai–Inthanon and Bentong–Raub sutures in Thailand and Peninsular Malaysia. These belts and sutures are generally regarded as containing the remnants of the oceanic crust of the Palaeo-Tethys that formed by seafloor spreading as a result of the separation of Gondwana- and Eurasia-derived blocks during the Middle Cambrian. In this paper we report the first discovery of abundant unaltered and retrograde eclogites that occur as irregular lenses and blocks in metasedimentary rocks of the CMOB, and these eclogites form an elongate and almost north–south-trending high-pressure (HP)–ultrahigh-pressure (UHP) metamorphic belt that is ∼200 km long and ∼50 km wide. The newly discovered phengite/talc/epidote–glaucophane eclogites, lawsonite–talc–phengite eclogites, dolomite/magnesite–kyanite eclogites and phengite–kyanite-bearing retrograde eclogites have enriched (E-) and normal mid-ocean ridge basalt (N-MORB)-like affinities and mainly positive as well as some negative whole-rock εNd values (–4·34 to +7·89), which suggest an enriched and depleted oceanic lithosphere source for their protoliths. Magmatic zircons separated from the epidote–glaucophane, magnesite–kyanite and (phengite–kyanite-bearing) retrograde eclogites gave protolith ages of 317–250 Ma, which fit well within the time frame of the opening of the Palaeo-Tethys during the Middle Cambrian and its closure during the Triassic. Abundant metamorphic zircons in the eclogites indicate a Triassic metamorphic event related to the subduction of the Palaeo-Tethys oceanic crust from 235 to 227 Ma. Taking into account previous isotopic age data, we now establish the periods of Early–Middle Triassic (246–227 Ma) and Late Triassic (222–209 Ma) as the ages of subduction and exhumation of the Palaeo-Tethyan oceanic crust, respectively. Thermodynamic modelling revealed that the eclogites record distinct HP–UHP peak metamorphic conditions of 23·0–25·5 kbar and 582–610 °C for the phengite–glaucophane eclogites, 24·0–25·5 kbar and 570–586 °C for the talc–glaucophane eclogites, 29·0–31·0 kbar and 675–712 °C for the dolomite–kyanite eclogites, and 30·0–32·0 kbar and 717–754 °C for the magnesite–kyanite eclogites. These P–T estimates and geochronological data indicate that the Palaeo-Tethys oceanic slab was subducted to different mantle depths from 75 km down to 95 km, forming distinct types of eclogite with a variety of peak eclogite-facies mineral assemblages. The eclogites consistently record clockwise metamorphic P–T–t paths characterized by a heating–compression prograde loop under a low geothermal gradient of 5–10 °C km–1, indicating the rapid subduction of cold oceanic crust at a rate of 4·5–6·0 km Ma–1, followed by isothermal or cooling–decompressive retrogression and exhumation at an average rate of 3·2–4·2 km Ma–1. The newly discovered eclogites of the CMOB with their signatures of ocean-crust subduction are petrologically, geochemically and geochronologically comparable with those of the LCSS, providing powerful support for the idea that a nearly 2000 km long HP–UHP eclogite belt extends from the northern Tibetan Plateau to the southeastern Tibetan Plateau, and that it represents the main boundary suture of the Palaeo-Tethyan domain. These results have far-reaching implications for the tectonic framework and complex metamorphic evolution of the Palaeo-Tethyan domain.

Description: SUMMARY Determining rock microstructure remains challenging, since a proper rock-physics model is needed to establish the relation between pore microstructure and elastic and transport properties. We present a model to estimate pore microstructure based on porosity, ultrasonic velocities and permeability, assuming that the microstructure consists on randomly oriented stiff equant pores and penny-shaped cracks. The stiff pore and crack porosity varying with differential pressure is estimated from the measured total porosity on the basis of a dual porosity model. The aspect ratio of pores and cracks and the crack density as a function of differential pressure are obtained from dry-rock P- and S-wave velocities, by using a differential effective medium model. These results are used to invert the pore radius from the matrix permeability by using a circular pore model. Above a crack density of 0.13, the crack radius can be estimated from permeability, and below that threshold, the radius is estimated from P-wave velocities, taking into account the wave dispersion induced by local fluid flow between pores and cracks. The approach is applied to experimental data for dry and saturated Fontainebleau sandstone and Chelmsford Granite.

Description: Multiple-layered coal seams widely exist in main coal mining areas of China. When these coal seams are exploited, the pillar mining method is always employed. This leads to many coal pillars left in the upper coal seams as a protective barrier. As a result, these residual pillars will not only cause the loss of coal resources but also could trigger environmental issues and a serious of mine disasters. A theoretical model was built to analyse the effect of the residual pillars. From the theoretical model, it was found that four stress concentration areas were formed by the upper residual coal pillars. To address the issues of the residual coal pillars, Datong Coal Mine Group has developed an innovative technology of the roof cutting with a chainsaw. A new protective coal seam mining method using chainsaw roof-cutting technology is introduced. A numerical model is constructed to analyse the mining pressure distribution law in working face within the lower layer coal seam. From the numerical simulation, the new protective layer mining method could reduce about 15.2% of the advancing stress, which contributes a lot to controlling the mining pressure within the lower layer. The field measurement showed that the hydraulic support utilised at the site was at lower pressure levels, which proves the new protective seam mining method can significantly reduce the working face pressure.

Description: Recent application of Fourier transform near infra-red spectroscopy (FT-NIRS) to predict age in fish otoliths has gained attention among fisheries managers as a potential alternative to costly production ageing of managed species. We assessed the age prediction capability of FT-NIRS scans in whole otoliths from red snapper, Lutjanus campechanus, collected from the US Gulf of Mexico and US Atlantic Ocean (South Atlantic). Otoliths were scanned with an FT-NIR spectrometer and resulting spectral signatures were regressed with traditionally estimated ages via partial least squares regression to produce calibration models, which were validated for predictive capability against test sets of otoliths. Calibration models successfully predicted age with R2 ranging 0.94–0.95, mean squared error ≤1.8 years, and bias 31 years were not well predicted, possibly due to light attenuation in the thickest otoliths. Our results suggest that FT-NIRS can improve efficiency in production ageing for fisheries management while maintaining data quality standards.

Description: Replacive symplectites (vermicular intergrowths of two or more minerals) are an important feature of layered igneous intrusions, recording evidence of late-stage reactions between interstitial liquid and crystals. They are common throughout the Layered Series of the 564 Ma Sept Iles layered intrusion in Quebec, Canada, and fall into three types: oxy-symplectites, ‘Type I’ symplectites, and ‘Type II’ symplectites. Oxy-symplectites are comprised of magnetite and orthopyroxene, nucleate on olivine primocrysts, and form via the reaction Olivine + O2 → Orthopyroxene + Magnetite; Type I symplectites (of which there are 3 distinct categories) are comprised of anorthitic plagioclase with pyroxene, amphibole, or olivine vermicules, grow from primocryst oxide grains, and replace primary plagioclase; and Type II symplectites (of which there are 2 distinct categories) are comprised of anorthitic plagioclase with orthopyroxene ± amphibole vermicules, grow from primocryst olivine grains, and replace primocryst plagioclase. Rare symplectites composed of biotite and plagioclase are also present. Symplectite growth occurred at 700-1030 °C with pressure constraints of 1-2 kbar. We propose that Type I symplectites, and some Type II symplectites, formed from interaction of primocrysts with residual Fe-rich liquid as a consequence of differential loss of an immiscible Si-rich liquid conjugate from the crystal mush. However, redistribution and concentration of hydrous fluids in incompletely solidified rock, or an increase in water activity of the interstitial melt, may be more plausible processes responsible for the formation of replacive symplectites comprising abundant hydrous mineral assemblages.

Description: This article investigates the diet of the snow crab (Chionoecetes opilio) and its feeding intensity in the Barents Sea. Data show that snow crab has a diverse diet that includes almost all types of benthic invertebrates living in the Barents Sea. There are differences between the diets of females and males and of juveniles and adults. Juveniles and females typically occupy shallow areas with communities of bivalve molluscs, while males typically live deeper on slopes and depressions where polychaetes and crustaceans are the most abundant groups. Stomach contents were analysed to determine the species composition and frequency of occurrence of various benthic taxa. Consumption of food was estimated and compared with data from the Russian seas of the Pacific region. The total annual consumption of macrozoobenthos by snow crab was calculated in accordance with its current distribution in the Barents Sea. Snow crab consumes at least 30 000 tonnes of benthos annually, which amounts to 0.1–0.2% of the total macrozoobenthic biomass in the investigated area. The population of snow crab causes the largest impact on the benthic communities in the northeastern part of the Barents Sea and near the south side of the Novaya Zemlya archipelago.

Description: Until the late 19th century, extensive beds of flat oyster Ostrea edulis populated the Central North Sea, which have vanished after intensive fisheries. At present, various initiatives are being carried out to investigate the potential to restore this former key species in the area. This historical ecological study contributes by delineating the former oyster bed area and through an assessment of its limits against known gradients in the North Sea. Extensive data from historical maps, texts, and ship-based surveys were used to synthesize our knowledge on the former beds. It was revealed that the area with oyster beds covered ∼6.2% of the total North Sea bottom, with a delineation that could partly be explained by hydrodynamic and temperature gradients. The position and extent of the area are notably different from the area that is used in recent feasibility studies on the restoration of North Sea oyster beds. The offshore oysters lived on muddy sand in relatively cold conditions, and there are several indications that their reproductive rate was low. The apparent disappearance of cold water adapted flat oysters will challenge restoration projects. This study provides indispensable information for the future restoration of flat oyster beds in the North Sea.

Description: Summary Relative location of microearthquakes that occurred at Mt. Pollino (Italy) from 2011 to 2013 have been analyzed with the aim of a detailed imaging of the geometry of active faults. We identified 27 clusters composed of a number of earthquakes from 9 to 33, with local magnitude in the range 0.6–2.7. The relative location shows that the distribution of hypocenters in each cluster is characterized by extension from few tens of meters to at most 350 m. For each cluster the hypocenter distribution was fitted by a plane to infer the fault orientation, and results were compared with the fault plane solutions corresponding to the focal mechanism of earthquakes of the same cluster. The comparison shows a good agreement in most of the cases. The relative location analysis, generally applied to earthquakes with similar waveform, has been improved to permit also the relative location of earthquakes characterized by not similar signals. To achieve this purpose a modified procedure that overcome the condition of very similar waveforms has been applied to estimate the time delay between first pulses of the master events. The relative location of master events of all clusters shows a precise imaging of the relative position of all analyzed sources and allows also to follow with high accuracy the evolution in time of the seismic swarm within the selected periods. The hypocenter position of master events and the nearly parallel fitting planes of any clusters suggest that most of the analyzed earthquakes were produced by different patches of the same fault. The final results depict a main fault plane characterized by NW-SE strike, dip of about 35–45 degrees, and depth between 4.5 and 6.5 km b.s.l.. Focal mechanisms, used also to evaluate the local stress field, are mostly of normal type with few strike slip solutions for the shallowest events. This result is in good agreement with the local tectonic stress regime that is characterized by predominant NE-SW trans-tension, as inferred from structural, seismological and geophysical data.

Description: SUMMARY The source characteristics of offshore seismic events, especially regular (or fast) and slow earthquakes, can provide key information on their source physics and frictional conditions at the plate boundary. Due to strong 3-D heterogeneities in offshore regions, such as those relating to sea water, accretionary prism and small-scale velocity heterogeneity, conventional methods using a 1-D earth model may mis-estimate source parameters such as the duration and radiation energy. Estimations could become severe inaccuracies for small offshore seismic events because high-frequency (〉1 Hz) seismograms, which are strongly affected by 3-D heterogeneities, are only available for analysis because of their signal-to-noise ratio. To investigate the effects of offshore heterogeneities on source parameter estimation for small seismic events, we analysed both observed and simulated high-frequency seismograms southeast off the Kii Peninsula, Japan, in the Nankai subduction zone. Numerical simulations of seismic wave propagation using a 3-D velocity structure model clarified the effects of each heterogeneity. Comparisons between observations and model simulations demonstrated that the thick low-velocity accretionary prism has significant effects on high-frequency seismic wave propagation. Especially for shallow low-frequency tremors occurring at depths just below the accretionary prism toe, seismogram durations are significantly broader than an assumed source duration, even for stations with epicentral distances of approximately 10 km. Spindle-shape seismogram envelopes were observed even at such close stations. Our results suggest that incorporating 3-D heterogeneities is necessary for practical estimation of source parameters for small offshore events.

Description: SUMMARY We present an updated joint tomographic method to simultaneously invert receiver function waveforms and surface wave dispersions for a 3-D S-wave velocity (Vs) model. By applying this method to observations from ∼900 seismic stations and with a priori Moho constraints from previous studies, we construct a 3-D lithospheric S-wave velocity model and crustal-thickness map for the central–east Tibetan plateau. Data misfit/fitting shows that the inverted model can fit the receiver functions and surface wave dispersions reasonably well, and checkerboard tests show the model can retrieve major structural information. The results highlight several features. Within the plateau crustal thickness is 〉60 km and outwith the plateau it is ∼40 km. Obvious Moho offsets and lateral variations of crustal velocities exist beneath the eastern (Longmen Shan Fault), northern (central–east Kunlun Fault) and northeastern (east Kunlun Fault) boundaries of the plateau, but with decreasing intensity. Segmented high upper-mantle velocities have varied occurrences and depth extents from south/southwest to north/northeast in the plateau. A Z-shaped upper-mantle low-velocity channel, which was taken as Tibetan lithospheric mantle, reflecting deformable material lies along the northern and eastern periphery of the Tibetan plateau, seemingly separating two large high-velocity mantle areas that, respectively, correspond to the Indian and Asian lithospheres. Other small high-velocity mantle segments overlain by the Z-shaped channel are possibly remnants of cold microplates/slabs associated with subductions/collisions prior to the Indian–Eurasian collision during the accretion of the Tibetan region. By integrating the Vs structures with known tectonic information, we derive that the Indian slab generally underlies the plateau south of the Bangong–Nujiang suture in central Tibet and the Jinsha River suture in eastern Tibet and west of the Lanchangjiang suture in southeastern Tibet. The eastern, northern, northeastern and southeastern boundaries of the Tibetan plateau have undergone deformation with decreasing intensity. The weakly resisting northeast and southeast margins, bounded by a wider softer channel of uppermost mantle material, are two potential regions for plateau expansion in the future.

Description: SUMMARY We present unambiguous evidence that the Mahogany salt body, located in the Northern part of the Gulf of Mexico, is seismically anisotropic. Evidence of anisotropy comes from shear wave splitting data obtained from a vertical seismic profile VSP. The data set consists of 48 vertically aligned receivers in a borehole drilled through the salt body. Splitting analysis is performed on shear wave phases that are converted from compressional waves at the top and bottom of the salt body. The phase converted at the top of the salt layer shows a clear signature of seismic anisotropy, while the phase at the base of the salt layer shows negligible splitting. We investigate the possibility of rock salt halite LPO as a cause of the observed anisotropy. A finite element geomechanical salt deformation model of the Mahogany salt body is developed, where deformation history is used as an input to the texture plasticity simulation program VPSC. Assuming a halite salt body, a full elasticity model is then calculated and used to create a synthetic VSP splitting data set. The comparison between the synthetic and real VSP data set shows that LPO of rock salt can explain the observed anisotropy remarkably well. This is the strongest evidence to date of seismic anisotropy in a deforming salt structure. Furthermore, for the first time, we are able to demonstrate clear evidence that deforming halite is the most likely cause of this anisotropy, combining data set analysis and synthetic full wave form modelling based on calculated rock salt elasticities. Neglecting anisotropy in seismic processing in salt settings could lead to potential imaging errors, for example the deformation models show an averaged delta parameter of δ = –0.06, which would lead in a zero offset reflection setting to a depth mismatch of 6.2 per cent. Our work also show how observations of salt anisotropy can be used to probe characteristics of salt deformation.

Description: Summary About a decade ago, noise-based monitoring became a key tool in seismology. One of the tools is Passive Image Interferometry (PII), which uses Noise Correlation Functions (NCF) to retrieve seismic velocity variations. Most studies apply PII to vertical components recording oceanic low-frequent ambient noise ( 1 Hz) on three three-component sensors. With environmental sensors inside the subsurface and in the air, we are able to connect observed velocity variations with environmental parameters. Temperatures below 0○C correlate well with strong shear wave velocity increases. The temperature sensors inside the ground suggest that a frozen layer of less than 5cm thickness causes apparent velocity increases above 2%, depending on the channel pair. The observations indicate that the different velocity variation retrieved from the different channel pairs are due to different surface wave responses inherent in the channel pairs. With dispersion curve modelling in a 1D medium we can verify that surfaces waves of several tens of meters wavelength experience a velocity increase of several percent due to a centimeters thick frozen layer. Moreover, the model verifies that Love waves show larger velocity increases than Rayleigh waves. The findings of this study provide new insights for monitoring with PII. A few days with temperature below 0○C can already mask other potential targets (e.g. faults or storage sites). Here we suggest to use vertical components, which is less sensitive to the frozen layer at the surface. If the target is the seasonal freezing, like in permafrost studies, we suggest to use three-component sensors in order to retrieve the Love wave response. This opens the possibility to study other small scale processes at the shallow subsurface with surface wave responses.

Description: SUMMARY The behaviour of tsunami waves at any location depends on the local morphology of the coasts, the encountered bathymetric features, and the characteristics of the source. However, the importance of accurately modelling the geometric properties of the causative fault for simulations of seismically induced tsunamis is rarely addressed. In this work, we analyse the effects of using two different geometric models of the subduction interface of the Calabrian Arc (southern Italy, Ionian Sea) onto the simulated tsunamis: a detailed 3-D subduction interface obtained from the interpretation of a dense network of seismic reflection profiles, and a planar interface that roughly approximates the 3-D one. These models can be thought of as representing two end-members of the level of knowledge of fault geometry. We define three hypothetical earthquake ruptures of different magnitudes (Mw 7.5, 8.0, 8.5) on each geometry. The resulting tsunami impact is evaluated at the 50-m isobath in front of coastlines of the central and eastern Mediterranean. Our results show that the source geometry imprint is evident on the tsunami waveforms, as recorded at various distances and positions relative to the source. The absolute differences in maximum and minimum wave amplitudes locally exceed one metre, and the relative differences remain systematically above 20 per cent with peaks over 40 per cent. We also observe that tsunami energy directivity and focusing due to bathymetric waveguides take different paths depending on which fault is used. Although the differences increase with increasing earthquake magnitude, there is no simple rule to anticipate the different effects produced by these end-member models of the earthquake source. Our findings suggest that oversimplified source models may hinder our fundamental understanding of the tsunami impact and great care should be adopted when making simplistic assumptions regarding the appropriateness of the planar fault approximation in tsunami studies. We also remark that the geological and geophysical 3-D fault characterization remains a crucial and unavoidable step in tsunami hazard analyses.

Description: SUMMARY We present results from complementary geological, topographic, seismic and electrical resistivity surveys at the Sagebrush Flat (SGB) site along the Clark fault (CF) strand of the San Jacinto fault zone trifurcation area southeast of Anza, California. Joint interpretation of these data sets, each with unique spatiotemporal sensitivities, allow us to better characterize the shallow (

Description: The 2013 Lushan Mw 6.7 earthquake is the largest blind thrust event ever occurred on the southern segment of the Longmen Shan fault system. It has attracted extensive attention since it occurred 5 yr later following the 2008 Mw 7.8 Wenchuan destructive earthquake in this region. However, its slip distribution is still on debate due to the complex tectonic settings and limited near-field observations. In this study, we added some near-field GPS data, together with previously published GPS data and levelling data, and take consideration of possible coseismic and post-seismic effects caused by the 2008 Wenchuan earthquake, to construct a more accurate horizontal and vertical coseismic surface displacement field associated with the 2013 Lushan earthquake with a better spatial coverage. Then we invert for a refined slip distribution based on a flat-ramp-flat fault suggested by the relocated aftershock sequence and seismic imaging. Our preferred fault plane is striking southwest with 211° and dipping at varying angles of 4°, 35° and 12° separately for such a flat-ramp-flat geometry. The main rupture is roughly characterized by two asperities, including a round disk on the ramp with larger slips and an adjoining oval asperity on the shallow flat with smaller slips. The maximum slip is 1.2 m at 14.3 km focal depth, located at ∼20 km to the northwest of the GCMT epicentre. The released geodetic moment is 1.50 $ imes $ 1019 Nm, equivalent to a Mw 6.7 earthquake. The slips on the fault plane clearly illustrate that this event is dominated by the thrusting and minor striking, which is consistent with its tectonic settings. Furthermore, if we assume the 2013 Mw 6.7 Lushan event to be the characteristic earthquake on the southern section of the Longmen Shan thrust zone, the accumulated strain should not be fully released by this strong event, and a potential seismic risk still exists in this region.

Description: Larval settlement is a key process in the lifecycle of benthic marine organisms; however, little is known on how it could change in reduced seawater pH and carbonate saturation states under future ocean acidification (OA). This is important, as settlement ensures species occur in optimal environments and, for commercially important species such as abalone, reduced settlement could decrease future population success. We investigated how OA could affect settlement success in the New Zealand abalone Haliotis iris by examining: (1) direct effects of seawater at ambient (pHT 8.05) and reduced pHT (7.65) at the time of settlement, (2) indirect effects of settlement substrates (crustose coralline algae, CCA) preconditioned at ambient and reduced pHT for 171 days, and (3) carry-over effects, by examining settlement in larvae reared to competency at ambient and reduced pHT (7.80). We found no effects of seawater pH or CCA incubation on larval settlement success. OA-induced carry-over effects were evident, with lower settlement in larvae reared at reduced pH. Understanding the mechanisms behind these responses is key to fully comprehend the extent to which OA will affect marine organisms and the industries that rely on them.

Description: 2-D full-waveform inversion (FWI) of shallow-seismic wavefields has recently become a novel way to reconstruct S-wave velocity models of the shallow subsurface with high vertical and lateral resolution. In most applications, seismic wave attenuation is ignored or considered as a passive modelling parameter only. In this study, we explore the feasibility and performance of multiparameter viscoelastic 2-D FWI in which seismic velocities and attenuation of P and S waves, respectively, and mass density are inverted simultaneously. Synthetic reconstruction experiments reveal that multiple crosstalks between all viscoelastic material parameters may occur. The reconstruction of S-wave velocity is always robust and of high quality. The parameters P-wave velocity and density exhibit weaker sensitivity and can be reconstructed more reliably by multiparameter viscoelastic FWI. Anomalies in S-wave attenuation can be recovered but with limited resolution. In a field-data application, a small-scale refilled trench is nicely delineated as a low P- and S-wave velocity anomaly. The reconstruction of P-wave velocity is improved by the simultaneous inversion of attenuation. The reconstructed S-wave attenuation reveals higher attenuation in the shallow weathering zone and weaker attenuation below. The variations in the reconstructed P- and S-wave velocity models are consistent with the reflectivity observed in a ground penetrating radar (GPR) profile.

Description: The presence of sets of open fractures is common in most reservoirs, and they exert important controls on the reservoir permeability as fractures act as preferential pathways for fluid flow. Therefore, the correct characterization of fracture sets in fluid-saturated rocks is of great practical importance. In this context, the inversion of fracture characteristics from seismic data is promising since their signatures are sensitive to a wide range of pertinent fracture parameters, such as density, orientation and fluid infill. The most commonly used inversion schemes are based on the classical linear slip theory (LST), in which the effects of the fractures are represented by a real-valued diagonal excess compliance matrix. To account for the effects of wave-induced fluid pressure diffusion (FPD) between fractures and their embedding background, several authors have shown that this matrix should be complex-valued and frequency-dependent. However, these approaches neglect the effects of FPD on the coupling between orthogonal deformations of the rock. With this motivation, we considered a fracture model based on a sequence of alternating poroelastic layers of finite thickness representing the background and the fractures, and derived analytical expressions for the corresponding excess compliance matrix. We evaluated this matrix for a wide range of background parameters to quantify the magnitude of its coefficients not accounted for by the classical LST and to determine how they are affected by FPD. We estimated the relative errors in the computation of anisotropic seismic velocity and attenuation associated with the LST approach. Our analysis showed that, in some cases, considering the simplified excess compliance matrix may lead to an incorrect representation of the anisotropic response of the probed fractured rock.

Description: SUMMARY Seismic full waveform inversion (FWI) is a state-of-the-art technique for estimating subsurface physical models from recorded seismic waveform, but its application requires care because of high non-linearity and non-uniqueness. The final outcome of global convergence from conventional FWI using local gradient information relies on an informative starting model. Bayesian inference using Markov chain Monte Carlo (MCMC) sampling is able to remove such dependence, by a direct extensive search of the model space. We use a Bayesian trans-dimensional MCMC seismic FWI method with a parsimonious dipping layer parametrization, to invert for subsurface velocity models from pre-stack seismic shot gathers that contain mainly reflections. For the synthetic study, we use a simple four-layer model and a modified Marmousi model. A recently collected multichannel off-shore seismic reflection data set, from the Lord Howe Rise (LHR) in the east of Australia, is used for the field data test. The trans-dimensional FWI method is able to provide model ensembles for describing posterior distribution, when the dipping-layer model assumption satisfies the observed data. The model assumption requires narrow models, thus only near-offset data to be used. We use model stitching with lateral and depth constraints to create larger 2-D models from many adjacent overlapping submodel inversions. The inverted 2-D velocity model from the Bayesian inference can then be used as a starting model for the gradient-based FWI, from which we are able to obtain high-resolution subsurface velocity models, as demonstrated using the synthetic data. However, lacking far-offset data limits the constraints for the low-wavenumber part of the velocity model, making the inversion highly non-unique. We found it challenging to apply the dipping-layer based Bayesian FWI to the field data. The approximations in the source wavelet and forward modelling physics increase the multimodality of the posterior distribution; the sampled velocity models clearly show the trade-off between interface depth and velocity. Numerical examples using the synthetic and field data indicate that trans-dimensional FWI has the potential for inverting earth models from reflection waveform. However, a sparse model parametrization and far offset constraints are required, especially for field application.

Description: SUMMARY Parker–Oldenburg's method is perhaps the most commonly used technique to estimate the depth of density interface from gravity data. To account for large density variations reported, for instance, at the Moho interface, between the ocean seawater density and marine sediments, or between sediments and the underlying bedrock, some authors extended this method for variable density models. Parker–Oldenburg's method is suitable for local studies, given that a functional relationship between gravity data and interface geometry is derived for Earth's planar approximation. The application of this method in (large-scale) regional, continental or global studies is, however, practically restricted by errors due to disregarding Earth's sphericity. Parker–Oldenburg's method was, therefore, reformulated also for Earth's spherical approximation, but assuming only a uniform density. The importance of taking into consideration density heterogeneities at the interface becomes even more relevant in the context of (large-scale) regional or global studies. To address this issue, we generalize Parker–Oldenburg's method (defined for a spherical coordinate system) for the depth of heterogeneous density interface. Furthermore, we extend our definitions for gravity gradient data of which use in geoscience applications increased considerably, especially after launching the Gravity field and steady-state Ocean Circulation Explorer (GOCE) gravity-gradiometry satellite mission. For completeness, we also provide expressions for potential. The study provides the most complete review of Parker–Oldenburg's method in planar and spherical cases defined for potential, gravity and gravity gradient, while incorporating either uniform or heterogeneous density model at the interface. To improve a numerical efficiency of gravimetric forward modelling and inversion, described in terms of spherical harmonics of Earth's gravity field and interface geometry, we use the fast Fourier transform technique for spherical harmonic analysis and synthesis. The (newly derived) functional models are tested numerically. Our results over a (large-scale) regional study area confirm that the consideration of a global integration and Earth's sphericty improves results of a gravimetric forward modelling and inversion.

Description: The growing seismic networks and the increasing number of permanent seismic stations can help in improving the physical basis of seismic hazard assessment. For this purpose, the definition of reference site conditions is of great significance. If a reliable estimate of the reference ground motion is known, its modification at any given site can be modelled with respect to that reference site. Since the choice of a well-characterized reference site is not straightforward, mainly due to the high variability in the shallow layers, such choices prove to be affected by large uncertainties. While proxy parameters like the average S-wave velocity over the uppermost 30 m (vS30) might help in characterizing reference site conditions, such parameters are neither available at all sites nor do they allow concluding that the site is not affected by amplification and attenuation effects. In this study, we identify prospective reference sites across Europe in a harmonized and fully data-driven way. All analysis is based on freely available geological and geophysical data and no on-site measurements or site-specific proxies are required. The study accounts for both the influence of amplification and attenuation in a large frequency range. To address the key conceptual issues, we verify our classification based on machine learning techniques in which the influence of the individual site characterization parameters is investigated. Our study indicates that around 250 sites in Europe over more than 2000 investigated are not affected by local site effects and can de facto be considered as reference sites based on the criteria applied.

Description: The S-receiver function (SRF) technique is an effective tool to study seismic discontinuities in the upper mantle such as the mid-lithospheric discontinuity (MLD) and the lithosphere-asthenosphere boundary (LAB). This technique uses deconvolution and aligns traces along the maximum of the deconvolved SV signal. Both of these steps lead to acausal signals, which may cause interference with real signals from below the Moho. Here we go back to the origin of the S-receiver function method and process S-to-P converted waves using S-onset times as the reference time and waveform summation without any filter like deconvolution or bandpass. We apply this ‘causal’ SRF (C-SRF) method to data of the USArray and obtain partially different results in comparison with previous studies using the traditional acausal SRF method. The new method does not confirm the existence of an MLD beneath large regions of the cratonic US. The shallow LAB in the western US is, however, confirmed with the new method. The elimination of the MLD signal below much of the cratonic US reveals lower amplitude but highly significant phases that previously had been overwhelmed by the apparent MLD signals. Along the northern part of the area with data coverage we see relics of Archean or younger north-west directed low-angle subduction below the entire Superior Craton. In the cratonic part of the US we see indications of the cratonic LAB near 200 km depth. In the Gulf Coast of the southern US we image relics of southeast directed shallow subduction, likely of mid-Paleozoic age.

Description: Ambient noise correlation has been used extensively to retrieve traveltimes of surface waves. However, studies of retrieving amplitude information and attenuation from ambient noise are limited. In this study, we develop methods and strategies to extract Rayleigh wave amplitude and attenuation from ambient noise correlation, based on theoretical derivation, numerical simulation, and practical considerations of real seismic data. The synthetic data included a numerical simulation of a highly anisotropic noise source and Earth-like temporally varying strength. Results from synthetic data validate that amplitudes and attenuations can indeed be extracted from noise correlations for a linear array. A temporal flattening procedure is effective in speeding up convergence while preserving relative amplitudes. The traditional one-bit normalization and other types of temporal normalization that are applied to each individual station separately are problematic in recovering attenuation and should be avoided. In this study, we propose an ‘asynchronous’ temporal flattening procedure for real data that does not require all stations to have data at the same time. Furthermore, we present the detailed procedure for amplitude retrieval from ambient noise. Tests on real data suggest attenuations extracted from our noise-based methods are comparable with those from earthquakes. Our study shows an exciting promise of retrieving amplitude and attenuation information from ambient noise correlations and suggests practical considerations for applications to real data.

Description: The themed article set (TS) “Science in support of a nonlinear non-equilibrium world” reflects the challenge presented by the increasing potential for complex non-linear behaviour in marine ecosystems, many of which are undergoing dramatic changes due to anthropogenic perturbations. Marine ecosystems are complex adaptive systems, yet management strategies are often guided by a linear, stable perspective that excludes non-linearities and the possibility for evolution and adaptation. Rapidly increasing amounts of observational data, the interdisciplinary development of powerful mathematical approaches from complexity theory, and the evolving complex adaptive systems approach that includes human behaviour promise to substantially advance the development of management strategies. For these reasons, the ICES Journal of Marine Science solicited contributions to a TS that would take up these issues. In this introduction, I focus on three important areas—causality, prediction, and emergence—where a “non-linear” perspective can advance our understanding and better support sustainable management of ecosystems. I briefly present the nine contributions that are included in the themed set and suggest some ways forward. We hope that these articles serve to convince managers and marine scientists of the benefits of incorporating approaches and results from research on non-linear dynamics.

Description: The island of Pantelleria, located in the Sicily Channel Rift Zone (Italy), has been the site of violent peralkaline silicic magmatism alternating with minor effusive to low-intensity Strombolian eruptions of basaltic composition. The basaltic rock suites exposed on the island were sampled to investigate the plumbing system dynamics through the study of chemical stratigraphy and temporal records of olivine crystals. Our petrographic and geochemical observations, together with the compositional variability of olivine, suggest different evolutionary histories for basaltic magmas erupted over two major periods divided by the ∼45 ka Green Tuff (GT) eruption. Core-to-rim compositional traverses across olivine crystals document different types of zoning. We recognized olivine zones affected by Fo oscillations at very fine scales in the inner cores, rims and/or in intermediate portions of crystals and used them to reconstruct the residence and passage of crystals through different magmatic environments, with P–T–ƒO2 and compositional characteristics constrained by thermodynamic modeling. The sequence of magmatic environments evidenced by olivine zoning indicate that the pre-GT volcanic period was dominated by injection at shallow crustal levels (∼300–200 MPa) of primitive melts, initially moving from a deep storage zone at the crust-mantle boundary. Supply of this magma significantly decreased after the GT eruption, while the dynamics of magma transfer within the upper portion of the plumbing system were greatly enhanced. The diffusive relaxation of olivine zoning provided the timing of storage and migration of a crystal through different environments. For magmas feeding the ancient (〉45 ka) basaltic activity we retrieved transfer histories that are much longer (up to ∼3 years) if compared with those calculated for the post-GT basalts (1–9 months). The compositional and temporal dataset presented in this study supports the idea that the GT eruption and the subsequent collapse of the volcanic edifice could have caused major changes to the internal structural setting of Pantelleria, creating more favorable conditions for the migration of magmas in the upper portions of the plumbing system.

Description: Volcanism following the initiation of subduction is vital to our understanding of this specific magma-generation environment. This setting is represented by the first development of the Izu–Bonin–Mariana arc system as subduction commenced along the Western Pacific margin in the Eocene. A new collection of volcanic rocks recovered from the islands and exposed crustal sections of the Bonin Ridge spans the first 10 Myr of arc evolution. An elemental and radiogenic isotope dataset from this material is presented in conjuction with new 40Ar/39Ar ages and a stratigraphic framework developed by a detailed mapping campaign through the volcanic sections of the Bonin Islands. The dating results reveal that both the locus and type of magmatism systematically changed with time in response to the progressive sinking of the slab until the establishment of steady-state subduction at around 7–8 Ma. Following initial mid-ocean ridge basalt (MORB)-like spreading-related basalt magmatism, volcanic centres migrated away from the trench and changed from high-Si boninite to low-Si boninite or high-Mg andesite, then finally tholeiitic or calcalkaline arc magma. Subducting pelagic sediment combined with Pacific-type igneous ocean crust dominates the slab input to the shallow source of high-Si boninites at 49 Ma, but high-precision Pb isotope data show that this sediment varies in composition along the subducting plate. At around 45 Ma, volcanism switched to low-Si boninite and the pelagic sediment signature was almost entirely replaced by volcanic or volcaniclastic material originating from a HIMU ocean island source. These low-Si boninites are isotopically consistent with a slab component comprising variable proportions of HIMU volcaniclastic rocks and Pacific MORB. In turn, this signature was replaced by a Pacific MORB-dominated flux in the post 45 Ma tholeiite and calcalkaline volcanic rocks. Notably, each change in slab-derived flux coincided with a change in the magma type. Fluctuations in the slab-derived geochemical signature were superimposed on a change in the mantle wedge source from highly depleted harzburgite to a depleted MORB-type mantle-type source. In turn, this may correspond to the increasing depth of the leading edge of the slab through this 5 Myr period.

Description: SUMMARY Distributed acoustic sensing (DAS) acquisition is becoming more and more popular for its dense sampling at a lower cost than seismometers. However, data processing for DAS data is challenging, especially for surface-deployedfibers, in which only the horizontal component of strain variation is effectively recorded. Also, the coupling between the fiber and the Earth is usually poor and the recorded single-component data are noisy. Thus, we introduce data processing strategies dedicated to enhancing the ambient-noise and active-source seismic data recorded by a horizontally deployed tactical fiber-optics cable buried in a sand dune area in Saudi Arabia. We propose a similarity-weighted stacking of randomly selected short-time duration windows to generate virtual common shot gathers (CSG) from the recorded ambient noise. The similarity-weighted stacking only counts the primary contributions of coherent events, while a short-time correlation can suppress the crosstalk usually present in late arrivals. The stacking fold is preserved or even can be increased by generating plenty of random time segments compared to stacking the full recording time. For the recorded active-source data, we skip the interferometric step, but use the envelope of the CSG. The envelope is needed to mitigate the complexity of waveforms, while preserving the slopes of arrivals. Then, we use the wave-equation-based Rayleigh-wave dispersion spectrum inversion, which utilizes all the dispersion modes available and does not require picking the dispersion curve, in estimating the shallow S-wave velocities. The local cross-correlation objective function allows for additional freedom in matching the modelled and observed data, and thus, helps us avoid falling into a local minimum when starting with kinematically poor velocity models.

Description: SUMMARY In this study, new methods are developed to estimate the dissipation factors, inhomogeneity parameters and phase velocities of the reflected waves at the free surface of a poro-viscoelastic solid in which the seismic wave propagation is described by effective Biot theory. The Christoffel equations of an effective Biot medium are solved for a general harmonic plane wave and the three complex velocities obtained corresponding to the shear wave (SV), fast-P wave and slow-P wave, together with their polarizations. Based on the complex form of the energy balance equation in an effective Biot material, expressions are derived for the energy ratios at the free surface. Moreover, the equations for the inhomogeneity parameters are derived as functions of the complex slowness or the unit polarization vectors. Based on the implicit and the explicit dissipation factor expressions, two methods are developed to obtain the dissipation factors, the inhomogeneity parameters and the phase velocities of mode-converted waves. These methods are illustrated by numerical examples which show that the dissipation factors, inhomogeneity parameters and phase velocities of reflected waves can strongly depend on the incidence angle (also reflected angle), the incident wave inhomogeneity parameter and the wave frequency. Ignoring these dependencies and using dissipation factors only valid for homogeneous waves can cause discrepancies in computed phase velocities and dissipation factors for interface generated (reflected/transmitted) inhomogeneous waves.

Description: Oxygen decline poses increasing risks to global shelf communities. This study was conducted to measure species oxygen exposures in situ and to assess risks of low oxygen based on the hypothesis that species risk varies spatially and is dependent on the interaction of shelf oxygen dynamics with habitat pattern, species mobility, and tolerance to low oxygen. Here, we report concomitant observations of oxygen, depth, and habitat for the most common benthic and epibenthic megafauna on the southern California shelf (30–125 m). The study was conducted when oxygen concentrations were the lowest observed within the previous 30 years off southern California (2012–2014), yet hypoxic conditions (

Description: In North America, studies regarding effects of CO2-induced low pH in bivalve aquaculture are largely restricted to the US Pacific coast. Studies on species from the northwest Atlantic are lacking. Furthermore, information on the roles of intergenerational exposure and biological sex in bivalve responses to low pH, particularly in an aquaculture-specific context, is scant. We tested if short-term (1 month) exposure to CO2-induced reductions in pHNBS affected the reproductive development of male and female eastern oysters (Crassostrea virginica) during hatchery-specific reproductive conditioning and whether maternal and/or paternal exposure influenced larval responses. Reduced pH (7.5–7.7) increased the rate of reproductive development in both males and females. There was no indication of intergenerational effects; adult pH conditions did not affect early larval development. In contrast, low pH conditions experienced by gametes during spawning, fertilization, and embryo incubation (48 h) resulted in higher larval survival (+6–8% from control), reduced shell height (−2 to 3 µm), and increased deformities (abnormal shell shape; +3–5%). We suggest that local adaptation to acidic land runoff may account for the positive effects of low pH observed in this study. Bioeconomic assessments are now needed to understand the implications of reduced pH on aquaculture operations in these regions of Atlantic Canada.

Description: Body waves can be extracted from correlation functions computed from seismic records even at teleseismic distances. Here we propose to use P and PcP waves from the secondary microseism frequency band that are propagating between Europe and the Eastern US to image the Core-Mantle Boundary (CMB) and D’ structure beneath the North Atlantic. This study presents the first 3D image of the lower mantle obtained from ocean-generated microseism data. Robustness of our results is evaluated by comparing images produced by propagation in both directions. Our observations reveal complex patterns of lateral and vertical variations of P-wave reflectivity with a particularly strong anomaly extending upward in the lower mantle up to 2600 km deep. We compare these results with synthetic data and associate this anomaly to a Vp velocity increase above the CMB. Our image aims at promoting the study of the lower mantle with microseism noise excitations.

Description: Uniaxial compression tests with combined acousto-optical monitoring techniques are conducted on flawed granite specimens, with the aim of investigating the fracture-related acoustic emission (AE) event rate characteristics at the unstable cracking phase in flawed rocks. The interevent time (IET) function F(τ) is adopted to interpret the AE time-series from damage stress (σcd) to ultimate failure, and photographic data are used to evaluate unstable cracking behaviours in flawed granite. The results show that a high AE event rate is always registered but intermittently interrupted by macrofracturing at the unstable cracking phase. The reversed U-shaped curve relation between the AE event rate and the loading time is documented in unstable flawed granite for the first time. The acoustic quiescence has a mechanismic linkage and quantitative correlation with stress drop, and this synchronous acousto-mechanical behaviour is a typical result of the initiation, growth and coalescence of macrocracks initiated from the flaw tips. Moreover, the reactivation and intensification of fracture process zones (FPZs) by increasing loads are the dominant mechanism triggering unstable crack growth in flawed granite.

Description: SUMMARY Seismoelectric (SE) signals, accompanying seismic wave fields radiated from earthquakes, can be observed on records of magnetotelluric stations. Assuming that these SE signals are generated by electrokinetic coupling we investigate whether they can be used as a ‘pore-space monitoring’-tool. Regarding future field experiments we analyse synthetic SE waveforms calculated for a fully saturated base model consisting of five layers overlying a half-space, resembling the conditions of the Armutlu Peninsula (Turkey). This example site stands for a location with near-surface thermal aquifers exposed to tectonic stress and significant microseismicity. As expected, coseismic SE waves arrive simultaneously with the seismic onsets whereas interface response (IRs) SE waves arrive (shortly) before the generating seismic onsets. Therein, so-called evanescent IRs show a similar moveout as seismic phases and so-called radiation IRs travel with zero slowness. We found that the influence of IRs on the overall SE signal can be identified by envelope analysis of SE time series and by seismoelectric spectral ratios (SESRs) in the frequency domain. For a sensitivity analysis we added an extra layer to the base model with differing porosity, porefluid salinity and permeability values. At near-epicentral distances both trace-envelopes and SESRs are sensitive to the porosity and porefluid salinity changes in the simulated near-surface aquifer. The SESRs’ and SE envelopes’ amplitudes vary in the order of up to some 10 per cent in response to porosity and salinity increases of factor 2 and 100, respectively. In contrast, a decrease of the permeability value by the factor 100 leads to an SESR amplitude variation of less than 1–10 per cent. In the Armutlu model the largest relative changes of SE signals occur near the epicentre where the ratio between coseismic and IR amplitudes is close to 1. For 1–6 km deep source depth the SE detection swell at the earth surface is in the order of magnitude 2–3, depending on the ambient electromagnetic noise and hypocentral distance. This estimate assumes that SE signals are recorded with standard magnetotelluric stations. It can improve if array methods are applied.

Description: SUMMARY Rock thermal conductivity is an essential input parameter for enhanced oil recovery methods design and optimization and for basin and petroleum system modelling. Absence of any effective technique for direct in situ measurements of rock thermal conductivity makes the development of well-log based methods for rock thermal conductivity determination highly desirable. A major part of the existing problem solutions is regression model-based approaches. Literature review revealed that there are only several studies performed to assess the applicability of neural network-based algorithms to predict rock thermal conductivity from well-logging data. In this research, we aim to define the most effective machine-learning algorithms for well-log based determination of rock thermal conductivity. Well-logging data acquired at a heavy oil reservoir together with results of thermal logging on cores extracted from two wells were the basis for our research. Eight different regression models were developed and tested to predict vertical variations of rock conductivity from well-logging data. Additionally, rock thermal conductivity was determined based on Lichtenecker–Asaad model. Comparison study of regression-based and theoretical-based approaches was performed. Among considered machine learning techniques Random Forest algorithm was found to be the most accurate at well-log based determination of rock thermal conductivity. From a comparison of the thermal conductivity—depth profile predicted from well-logging data with the experimental data, and it can be concluded that thermal conductivity can be determined with a total relative error of 12.54 per cent. The obtained results prove that rock thermal conductivity can be inferred from well-logging data for wells that are drilled in a similar geological setting based on the Random Forest algorithm with an accuracy sufficient for industrial needs.

Description: SUMMARY Atmospheric shock waves are a common phenomenon in explosive volcanic eruptions. We consider the motion of a spherical shock wave generated by a point source in the strong shock approximation. The shock front corresponds to discontinuities in the gas velocity, density, pressure and temperature, which are calculated as functions of the energy of the explosion. The problem is solved analytically for the distributions of velocity, density, pressure and temperature in the atmosphere as functions of the distance from the source. The motion of the shock wave being supersonic, the solution is valid for a few seconds after the explosion, corresponding to a distance of few kilometres. The acoustic effect of the shock wave, expressed by the peak sound pressure level, is calculated and may reach hundreds of decibels. The pressure waveform that could be recorded in the vicinity of the volcano is calculated and compared with typical waveforms in weak shock conditions. The change in the refractive index of air due to density inhomogeneity is calculated and the conditions under which a condensation cloud is formed behind the shock front are investigated.

Description: I describe my unlikely path into marine science from a childhood in the Arizona desert and Oregon woods. Without realizing it, I developed a sense of place in nature and the value of open interdisciplinary communication among diverse scientists. My undergraduate education emphasized physiological adaptations to the environment or what might now be considered the “fundamental niche”, and my graduate thinking was inspired by a population/community based evolutionary understanding of how strong interactions define a “realized niche”. I have attempted to define strong interactions in three different ecosystems. This difficult problem is confounded by the loss of natural systems resulting from human impacts. I discuss my frustrations with eroding conservation efforts in a society that is rapidly devaluing nature and consider how we might recover our most fundamental values. I conclude that there is an urgent need to improve field-based teaching of undergraduate non-majors about nature and to be much more effective in our interactions with the general public. If we hope to have our legacy include a liveable world with natural places, we urgently need to act unilaterally to shift some of our values and reward systems towards the challenge of educating the undergraduates and especially the general public.

Description: I present a prescription for computing free-air coseismic and post-seismic gravity changes induced by seismic sources in a viscoelastic earth model. I assume a spherical earth geometry and a 2.5-D calculation, that is, 3-D motions that satisfy the equations of quasi-static equilibrium on a 2-D viscoelastic structure. The prescription permits application to regional gravity computations where a 2-D structure adequately represents the structural heterogeneity. I use a hybrid approach where deformation is computed on a discretized domain and the resulting density perturbations are expanded with spherical harmonics to produce the free-air gravity field. Starting with a solution to the equations of quasi-static displacements in the Laplace transform domain for a given dislocation source, I solve Poisson’s equation using Lagrangian interpolation on spectral element nodes to compute the required deformation quantities that contribute to free-air gravity. A numerical inverse Laplace transform then yields time domain results. This methodology is tested with analytic solutions on a spherically stratified viscoelastic structure, then applied to evaluate the effect of a descending slab of relatively high viscosity on post-seismic gravity in a megathrust faulting setting.

Description: Accounting for movement and mixing in stock assessment is important for managing sustainable fisheries, particularly for highly migratory species. However, many fisheries management approaches continue to use single-stock, single-area models to assess mixed-population stocks that are known to have complex movement dynamics. We evaluated a single-stock, single-area stock assessment model’s performance on fishery pseudodata generated using a spatially complex operating model that incorporates movement and mixing of simulated Atlantic bluefin tuna-like populations. Structural model misspecification produced positively biased perceptions of size and productivity of the smaller western population, based on supplement by the larger eastern population, and negatively biased perceptions of the size and productivity of the eastern population due to net movement of fish out of the eastern stock area. This bias could lead to unintended overexploitation of the smaller western population and potential for foregone yield of the larger eastern population. Our findings provide a greater understanding of the effects of movement and mixing on single-stock, single-area model-based management approaches and emphasize the importance of explicitly considering these dynamics in ensuring the sustainability of highly migratory species like Atlantic bluefin tuna.

Description: An algorithm developed in Tanimoto & Wang for deriving shallow elasticity structure from colocated pressure and seismic instruments is applied to data at nine stations in the Piñon Flat Observatory (PFO) in Southern California. Depth resolution kernels indicate that this approach can recover near-surface shear-modulus structure in the upper 50–100 m of the Earth. Our estimate for Vs30 at the closest station (BPH01) agrees with an independent result by Yong et al. which was based on an active source survey with a geophone array. Vs30 derived from the borehole S-wave speed model by Fletcher et al. also agrees with estimates at nearby locations, considering the size of uncertainties. Derived structures are much smoother than their models, however. This is because our method analyses slow quasi-static deformation of the solid Earth in the frequency range 0.005–0.05 Hz that inherently limits depth resolution. Lack of detailed variations in structures does not seem to be the problem for estimating Vs30 as Vs30 is an averaged quantity from the upper 30 m of the Earth. Also our estimates are not affected by different choices of frequency range in the inversion. However, our method may have some difficulty at stiff, hard-rock sites because ground deformation caused by surface pressure can become small and the key observables in our method, the ratios between seismic amplitude and surface pressure change, become difficult to measure accurately.

Description: Machine learning, a subfield of artificial intelligence, offers various methods that can be applied in marine science. It supports data-driven learning, which can result in automated decision making of de novo data. It has significant advantages compared with manual analyses that are labour intensive and require considerable time. Machine learning approaches have great potential to improve the quality and extent of marine research by identifying latent patterns and hidden trends, particularly in large datasets that are intractable using other approaches. New sensor technology supports collection of large amounts of data from the marine environment. The rapidly developing machine learning subfield known as deep learning—which applies algorithms (artificial neural networks) inspired by the structure and function of the brain—is able to solve very complex problems by processing big datasets in a short time, sometimes achieving better performance than human experts. Given the opportunities that machine learning can provide, its integration into marine science and marine resource management is inevitable. The purpose of this themed set of articles is to provide as wide a selection as possible of case studies that demonstrate the applications, utility, and promise of machine learning in marine science. We also provide a forward-look by envisioning a marine science of the future into which machine learning has been fully incorporated.

Description: In probabilistic Bayesian inversions, data uncertainty is a crucial parameter for quantifying the uncertainties and correlations of the resulting model parameters or, in transdimensional approaches, even the complexity of the model. However, in many geophysical inference problems it is poorly known. Therefore, it is common practice to allow the data uncertainty itself to be a parameter to be determined. Although in principle any arbitrary uncertainty distribution can be assumed, Gaussian distributions whose standard deviation is then the unknown parameter to be estimated are the usual choice. In this special case, the paper demonstrates that a simple analytical integration is sufficient to marginalise out this uncertainty parameter, reducing the complexity of the model space without compromising the accuracy of the posterior model probability distribution. However, it is well known that the distribution of geophysical measurement errors, although superficially similar to a Gaussian distribution, typically contains more frequent samples along the tail of the distribution, so-called outliers. In linearized inversions these are often removed in subsequent iterations based on some threshold criterion, but in Markov chain Monte Carlo (McMC) inversions this approach is not possible as they rely on the likelihood ratios, which cannot be formed if the number of data points varies between the steps of the Markov chain. The flexibility to define the data error probability distribution in McMC can be exploited in order to account for this pattern of uncertainties in a natural way, without having to make arbitrary choices regarding residual thresholds. In particular, we can regard the data uncertainty distribution as a mixture between a Gaussian distribution, which represent valid measurements with some measurement error, and a uniform distribution, which represents invalid measurements. The relative balance between them is an unknown parameter to be estimated alongside the standard deviation of the Gauss distribution. For each data point, the algorithm can then assign a probability to be an outlier, and the influence of each data point will be effectively downgraded according to its probability to be an outlier. Furthermore, this assignment can change as the McMC search is exploring different parts of the model space. The approach is demonstrated with both synthetic and real tomography examples. In a synthetic test, the proposed mixed measurement error distribution allows recovery of the underlying model even in the presence of 6 per cent outliers, which completely destroy the ability of a regular McMC or linear search to provide a meaningful image. Applied to an actual ambient noise tomography study based on automatically picked dispersion curves, the resulting model is shown to be much more consistent for different data sets, which differ in the applied quality criteria, while retaining the ability to recover strong anomalies in selected parts of the model.

Description: SUMMARY In a dislocation problem, a paradoxical discordance is known to occur between an original smooth curve and an infinitesimally discretized curve. To solve this paradox, we have investigated a non-hypersingular expression for the integral kernel (called the stress Green’s function) which describes the stress field caused by the displacement discontinuity. We first develop a compact alternative expression of the non-hypersingular stress Green’s function for general 2-D and 3-D infinite homogeneous elastic media. We next compute the stress Green’s functions on a curved fault and revisit the paradox. We find that previously obtained non-hypersingular stress Green’s functions are incorrect for curved faults, and that smooth and infinitesimally segmented faults are equivalent. Their compatibility bridges the gap between analytical methods featuring curved faults and numerical methods using subdivided flat patches.

Description: The goal of this contribution is to investigate the expected improvement of temporal gravity field determination via a couple of high-low satellite-to-satellite tracking (HLSST) missions. The simulation system is firstly validated by determining monthly gravity field models within situ GRACE GPS tracking data. The general consistency between the retrieved solutions and those developed by other official agencies indicates the good performance of our software. A 5-yr full-scale simulation is then performed using the full error sources including all error components. Analysis of each error component indicates that orbit error is the main contributor to the overall HLSST-derived gravity field model error. The noise level of monthly solution is therefore expected to reduce 90 per cent in terms of RMSE over ocean when the orbit accuracy improves for a magnitude of one order. As for the current HLSST mission consisting of a current GNSS receiver and an accelerometer (10−10 and 10−9 m s–2 noise for sensitive and non-sensitive axes), it is expected to observe monthly (or weekly) gravity solution at the spatial resolution of about 1300 km (or 2000 km). As for satellite constellations, a significant improvement is expected by adding the second satellite with the inclination of 70° and the third satellite with the inclination of 50°. The noise reduction in terms of cumulative geoid height error is approximately 51 per cent (or 62 per cent) when the observations of two (or three) HLSST missions are used. Moreover, the accuracy of weekly solution is expected to improve 40–70 per cent (or 27–59 per cent) for three (or two) HLSST missions when compared to one HLSST mission. Due to the low financial costs, it is worthy to build a satellite constellation of HLSST missions to fill the possible gaps between the dedicated temporal gravity field detecting missions.

Description: Seismoelectric measurements are conducted with a synthetic porous rock sample to model an ocean exploration. Two kinds of seismoelectric coupling signals, that is, the interfacial EM wave signal and the coseismic electric signal, have been recorded by the electrodes buried inside a rock sample instead of those located in the fluid or in the solid region near the interface as performed in previous works. These seismoelectric signals are clearly observed and identified with a high signal-to-noise ratio. The characteristics of the measured interfacial EM wave and coseismic electric signals are analysed with the experimental data. We also simulate the seismoelectric conversion fields and make a comparison between the measured and simulated seismoelectric signals. The result shows that the simulated and measured signals match well for both the interfacial EM wave and the coseismic electric fields accompanying the fast P wave. Our results also show that the amplitudes of seismoelectric signals are in the order of tens to hundreds of microvolts with our experimental system. This confirms that the seismoelectric signals are measurable in the interior of the rocks with current measurement techniques, suggesting the seismoelectric measurement to be a potential method for studying characteristics of the material beneath the seafloor.

Description: High-pressure (〉15 kbar) melts of intermediate–felsic materials have been well studied by experiments, whereas their existence in nature, especially in orogenic belts, is rarely examined. With the aim of identifying and characterizing high-pressure partial melts of intermediate–felsic continental crusts, this study presents comprehensive geochemical and geochronological data for 47 Jurassic granites (166∼157 Ma) from the Sulu orogen. These Sulu Jurassic granites (SJG) consist of quartz, K-feldspar and plagioclase with minor mineral assemblages of biotite ± muscovite ± garnet ± epidote ± allanite. Their low mafic mineral abundance, high SiO2 and Al2O3, and low FeOt + MgO contents show leucogranite-like affinities. They have low Mg#, low Rb/Sr, and mildly peraluminous features, collectively suggesting an intermediate–felsic orthogneissic source. Whole-rock Zr saturation thermometry and Ti-in-zircon thermometry together suggest initial magma temperatures between 695 ± 32 °C and 751 ± 27 °C (1 standard deviation), indicating derivation from water-present melting. The SJG notably feature high Sr contents (average 792 ppm), high Sr/CaO ratios (average 476) as well as inter-correlated low REE concentrations (average ΣREE 87 ppm), low Th concentrations (average 5·1 ppm) and positive Eu anomalies (Eu/Eu* up to 2·94). These characteristics are best explained by partial melting of intermediate–felsic sources under high pressure (〉15 kbar), leaving residuum where feldspar is sparse or absent and allanite is present. Inherited zircon age spectra and Sr–Nd–Pb isotopic compositions suggest that their source components could be mainly the Triassic orthogneisses whose protoliths are from the northern margin of the South China Block, probably in a wedge structure where the exhumed felsic slabs were wedged into the crust of the North China Block in the middle–late Jurassic and formed a stacked thickened crust. The wedge structure was most probably driven by synchronous large-scale strike-slip of the Tanlu fault, as a far-field effect of the oblique subduction of the paleo-Pacific plate. The characteristic chemical features observed in this study may be applied to identifying partial melts with similar petrogenesis elsewhere.

Description: SUMMARY The volcanism of the Eifel volcanic field (EVF), in west-central Germany, is often considered an example of hotspot volcanism given its geochemical signature and the putative mantle plume imaged underneath. EVF's setting in a stable continental area provides a rare natural laboratory to image surface deformation and test the hypothesis of there being a thermally buoyant plume. Here we use Global Positioning System (GPS) data to robustly image vertical land motion (VLM) and horizontal strain rates over most of intraplate Europe. We find a spatially coherent positive VLM anomaly over an area much larger than the EVF and with a maximum uplift of ∼1 mm yr−1 at the EVF (when corrected for glacial isostatic adjustment). This rate is considerably higher than averaged over the Late-Quaternary. Over the same area that uplifts, we find significant horizontal extension surrounded by a radial pattern of shortening, a superposition that strongly suggests a common dynamic cause. Besides the Eifel, no other area in NW Europe shows significant positive VLM coupled with extensional strain rates, except for the much broader region of glacial isostatic adjustment. We refer to this 3-D deformation anomaly as the Eifel Anomaly. We also find an extensional strain rate anomaly near the Massif Central volcanic field surrounded by radial shortening, but we do not detect a significant positive VLM signal there. The fact that the Eifel Anomaly is located above the Eifel plume suggests that the plume causes the anomaly. Indeed, we show that buoyancy forces induced by the plume at the bottom of the lithosphere can explain this remarkable surface deformation. Plume-induced deformation can also explain the relatively high rate of regional seismicity, particularly along the Lower Rhine Embayment.

Description: SUMMARY Seismic body wave traveltime tomography and surface wave dispersion tomography have been used widely to characterize earthquakes and to study the subsurface structure of the Earth. Since these types of problem are often significantly non-linear and have non-unique solutions, Markov chain Monte Carlo methods have been used to find probabilistic solutions. Body and surface wave data are usually inverted separately to produce independent velocity models. However, body wave tomography is generally sensitive to structure around the subvolume in which earthquakes occur and produces limited resolution in the shallower Earth, whereas surface wave tomography is often sensitive to shallower structure. To better estimate subsurface properties, we therefore jointly invert for the seismic velocity structure and earthquake locations using body and surface wave data simultaneously. We apply the new joint inversion method to a mining site in the United Kingdom at which induced seismicity occurred and was recorded on a small local network of stations, and where ambient noise recordings are available from the same stations. The ambient noise is processed to obtain inter-receiver surface wave dispersion measurements which are inverted jointly with body wave arrival times from local earthquakes. The results show that by using both types of data, the earthquake source parameters and the velocity structure can be better constrained than in independent inversions. To further understand and interpret the results, we conduct synthetic tests to compare the results from body wave inversion and joint inversion. The results show that trade-offs between source parameters and velocities appear to bias results if only body wave data are used, but this issue is largely resolved by using the joint inversion method. Thus the use of ambient seismic noise and our fully non-linear inversion provides a valuable, improved method to image the subsurface velocity and seismicity.

Description: Understanding recruitment variability in marine fisheries has benefits for the stock management and dependent fishers’ ability to plan their income. Here, we overview past and new research on the complex recruitment dynamics of redleg banana prawns Penaeus indicus in Australia’s Joseph Bonaparte Gulf to assess themes dating back to the time of Hjort and identify new challenges. During 2015 and 2016, redleg prawn catches and catch per unit effort decreased to anomalously low levels, suggesting a substantial decrease in prawn biomass. We hypothesized that low catches could be explained by temporary drops in sea level and rainfall potentially reducing the ability of postlarvae to reach their nursery ground. We contend that very bad prawn catch years may be predicted using two variables that are possible drivers of recruitment—the January Southern Oscillation Index and the combined January to February cumulative rainfall. However, due to challenges in verifying and defining such environmental relationships for inclusion in a stock assessment, we propose development of a harvest strategy framework to support management recommendations. Our study highlights the increasing role of anthropogenic climate change in exaggerating the impacts of environmental drivers on recruitment processes and the need to also focus on multidisciplinary research.

Description: A career does not follow a straight path. Determination, decision-making, and focus are essential ingredients, as well as a fair amount of flexibility, especially when one is struggling with contradictory signals. Career planning and the necessary decision-making must be learned however, and this may be particularly difficult when negative outcomes are likely and encouragement is rare. Under such circumstances, finding a job that makes one happy could be considered a noteworthy measure of success. However, even after attaining such a position, many tend to compare their own performance and career development with those of the celebrities in the field. This can only result in frustration and insecurity. Furthermore, success in marine science is generally characterized by metrics, together with the manner in which one’s career has advanced through a series of positions occupied in the zig-zag from student life to retirement. For me, a more personal kind of success has been to overcome the fear of failure that arises through constant comparison of my own performance and achievements with those who are perceived as the best in the field. This might be viewed more as social anxiety than fear as I will explain in this article.

Description: SUMMARY If a crystal lattice is subjected to a stress, it becomes distorted and no longer represents the ideal crystal symmetry, and if the stress introduces defects such as dislocations, some of this distortion is preserved after the applied stress is removed. In this study, we investigate lattice distortion in quartz at the micron scale with synchrotron X-ray Laue diffraction. From Laue images the local deviatoric strain tensor is derived and corresponding stresses are calculated based on elastic properties. The method is applied to metasedimentary quartzites from the Bergell Alps that were deformed at conditions of greenschist facies metamorphism. The residual palaeostrain is represented in maps of the deviatoric strain tensor components and with deviatoric strain axis pole figures. Data suggest overall shortening perpendicular to the schistosity plane but with considerable asymmetry relative to foliation and lineation, probably attributed to simple shear. Crystallographic pole figures from Laue diffraction agree with neutron diffraction and EBSD measurements and display quartz c-axes girdle distributions with maxima also perpendicular to schistosity. The method shows promise to be used as a palaeo-piezometer to unravel the stress field during tectonic deformation.

Description: Both equilibrium and dynamic crystallization experiments have been performed on a hydrous haplogranitic melt at 200 MPa to model nucleation and growth mechanisms and simulate pegmatite textures. The equilibrium results provide a reference frame (phase assemblages and compositions, liquidus and solidus temperatures and dependence with the melt H2O concentration) to parametrize the kinetic experiments. The seven H2O-saturated dynamic crystallization experiments followed a specific time-temperature path. After a pre-conditioning step at 800°C, charges were cooled between 3·5 and 7°C/min to 700, 660 and 600°C corresponding to ΔT of 20, 60 and 120°C. Dwell times ranged from 42 up to 1440h. Variable mineral assemblages and textures, and two types of polymineralic assemblages were obtained depending on ΔT and t. For ΔT = 120°C, crystallization is sequential and includes graphic quartz–alkali-feldspar intergrowths characteristic of pegmatite textures. The crystallization sequence reflects nucleation and growth of kinetically-favoured metastable phases and solid solution compositions from the supercooled melt. Early alkali-feldspars are more K-rich than expected at equilibrium and late albites more Na-rich. The K-rich graphic texture progressively evolves to a Na-rich intergrowth texture. Melts also follow a progressive though limited sodic evolution with time. At the interface of growing alkali-feldspars, melts are enriched in SiO2 and depleted in Al2O3, Na2O and, to a lesser extent, K2O. H2O accumulates at the interface reaching concentration levels higher (by 1–2 wt %) than the saturation. Rejection of SiO2 and H2O at the interface controls the effective undercooling in the local melt and promotes rapid textural changes toward larger grain sizes at the front of graphic zones. Textural ripening takes place contemporaneously to sequential crystallization. Growth rates for quartz and alkali-feldspar are tightly grouped, between 7·3 x 10-11 and 1·6 x 10-12 m s-1. Textures from the dynamic crystallization experiments closely resemble natural pegmatites, but layered aplite units have not been reproduced. Our results confirm and strengthen the importance of liquidus undercooling to generate pegmatite textures.

Description: SUMMARY This study establishes the proper framework in which to compare seismic observations with mineral physics constraints for studies of the inner core by determining how the elements of the elasticity tensor are sampled by the normal modes of the Earth. The obtained mapping between the elements of the elasticity tensor and the seismic wave speeds shows that the choice of averaging scheme used to calculate isotropic properties is crucial to understand the composition of the inner core, especially for comparison with the shear wave speed such as that provided in PREM. For example, the appropriate shear wave speed calculated for an Fe-Ni-Si hcp alloy at inner-core conditions differs from the shear wave speed obtained by taking a Reuss average by as much as $27, { m per, cent}$. It is also shown for the first time that by combining the isotropic observations based upon normal-mode characteristic frequencies and anisotropic parameters from their splitting, the five independent elastic parameters (A, C, F, L and N) that fully describe a transversely isotropic inner core can be uniquely constrained. The elastic values based upon a variety of mode-splitting studies are reported, and the differences between models from various research groups are shown to be relatively small ($lt 10, { m per, cent}$). Additionally, an analogous body-wave methodology is developed to approximately estimate the five independent elastic constants from observations of compressional wave traveltime anomalies. The body-wave observations are utilized to consider the depth dependence of inner-core anisotropy, in particular, the structure of the innermost inner core. Finally, we demonstrate that substantial errors may result when attempting to relate seismically observed P and S wave speeds from Debye velocities obtained through nuclear resonant inelastic X-ray scattering. The results of these experiments should be compared directly with the Debye velocity calculated from seismically constrained elastic constants. This manuscript provides a new set of formulae and values of seismic observations of the inner core that can be easily compared against mineral physics constraints for better understanding of the inner-core composition.

Description: SUMMARY The loess from the northern piedmont of the Dabie Mountains is in a transition area between loess from the Chinese Loess Plateau, the Quaternary red soils of southern China and the Xiashu loess. Despite its significant location, the study has been inadequate. In this study, the Guangshan section in the northern piedmont of the Dabie Mountains was selected for investigation. Environmental magnetism, geochemistry, colour reflectance and optical diffuse reflection spectroscopy analyses were applied to detect the magnetic variations in the loess. The results showed that (1) the magnetic minerals consisted mainly of magnetite, maghemite, hematite and goethite, which are the same as those in the Quaternary loess from the Chinese Loess Plateau, the Xiashu loess and the Quaternary red soils of southern China. The average magnetic particles were in the pseudo-single domain, like those of the Chinese Loess Plateau loess. (2) Unit III of the Guangshan section (2.4–4 m), with high chemical index of alteration and low Ba-index, was demonstrated as the most strongly developed palaeosol in the whole section, in agreement with field observations (more Fe-Mn films and weakly vermiculated development). However, it exhibited minimal susceptibility values and the lowest concentration of fine ferrimagnetic minerals. Simultaneously, the unit had low hematite to goethite ratio (Hm/Gt), suggesting that the pedogenic environment was humid; and it also had high values of b* and Gt%, implying that there was more goethite. Therefore, we can conclude that excessive soil moisture and intensive pedogenesis dissolved the fine ferrimagnetic minerals originally produced by pedogenesis and transformed them into goethite. These results could help to trace the palaeoclimatic evolution of the study area and clarify the magnetic variations of loess in different climates throughout China.

Description: SUMMARY In this paper, we present 3-D numerical simulations in which a compositional mantle plume rises through a shallow mantle layer at the rear of a retreating slab. The slab–plume buoyancy flux ratio Bs/Bp is varied from 3 to 30 over nine slab–plume simulations. The plume causes an overall decrease of the slab retreat rate from 15 to 7 % in the corresponding range $B_s/B_p=[3,30]$. The retreat rate decrease occurs in two stages: the first decrease occurs remotely when plume and slab are hundreds of kilometres apart; the second decrease is linked to the slab–plume-head impact. Continuous tracking of key positions along the plume head and conduit together with velocity profiling further shows a very close interplay between the conduit and plume head dynamics. In particular, the combination of the slab– and plume(head)–induced flows at the rear of the slab increases the advection and tilt of the conduit and causes its flaring with height in the direction parallel to the trench. As a result, the conduit source slowly drifts away from the slab by hundreds of kilometres and flares by one and a half times its original radius over 30 Myr of plume head spreading. The conduit tilt and flaring result in an increasingly unbalanced azimuthal distribution of the incoming plume flux from the feeding conduit into the head. These changes in the feeding conditions coupled with the mantle flow at the level of the plume head lead to the asymmetric spreading of plume material beneath the plate in the preferential direction that is parallel to the trench. Upon its arrival against the slab, it can be a front of buoyant material set to subduct along the slab width that has widened up to two and a half times its initial dimension. In nature, it is expected to extend from hundreds to thousands of kilometres depending on the slab–plume buoyancy flux ratio. To our knowledge, this study is the first to highlight the inter-relationships between a plume head and its feeding conditions in the plume source region.

Description: SUMMARY We describe array methods to search for low signal-to-noise ratio (SNR) signals in long-period seismic data using Fourier analysis. This is motivated by published results that find evidence of solar free oscillations in the Earth's seismic hum. Previous work used data from only one station. In this paper, we describe methods for computing spectra from array data. Arrays reduce noise level through averaging and provide redundancy that we use to distinguish coherent signal from a random background. We describe two algorithms for calculating a robust spectrum from seismic arrays, an algorithm that automatically removes impulsive transient signals from data, a jackknife method for estimating the variance of the spectrum, and a method for assessing the significance of an entire spectral band. We show examples of their application to data recorded by the Homestake Mine 3-D array in Lead, SD and the Piñon Flats PY array. These are two of the quietest small aperture arrays ever deployed in North America. The underground Homestake data has exceptionally low noise, and the borehole sensors of the PY array also have very low noise, making these arrays well suited to finding very weak signals. We find that our methods remove transient signals effectively from the data so that even low-SNR signals in the seismic background can be found and tested. Additionally, we find that the jackknife variance estimate is comparable to the noise floor, and we present initial evidence for solar g-modes in our data through the T2 test, a multivariate generalization of Student's t-test.

Description: The rock type most commonly associated with komatiite throughout Earth’s history is tholeiitic basalt. Despite this well-known association, the link between komatiite and basalt formation is still debated. Two models have been suggested: that tholeiitic basalts represent the products of extensive fractional crystallization of komatiite, or that basalts are formed by lower degrees of mantle melting than komatiites in the cooler portions of a zoned plume. We present major and trace element data for tholeiitic basalts (∼7·5 wt% MgO) and dunites (46–48 wt% MgO) from the Palaeoproterozoic Winnipegosis Komatiite Belt (WKB), which, along with previous data for WKB komatiites (17–26 wt% MgO), are utilized to explore the potential links between komatiite and basalt via crystallization processes. The dunites are interpreted as olivine + chromite cumulates that were pervasively serpentinized during metamorphism. Their major and immobile trace element relationships indicate that the accumulating olivine was highly magnesian (Mg# = 0·91–0·92), and that small amounts (4–7 wt% on average) of intercumulus melt were trapped during their formation. These high inferred olivine Mg# values suggest that the dunites are derived from crystallization of komatiite. The tholeiitic basalts have undergone greenschist-facies metamorphism and have typical geochemical characteristics for Palaeoproterozoic basalts, with the exception of high FeO contents. Their REE patterns are similar to Winnipegosis komatiites, although absolute concentrations are higher by a factor of ∼2·5. The ability of thermodynamic modelling with MELTS software to reproduce komatiite liquid lines of descent (LLD) is evaluated by comparison with the crystallization sequence and mineral compositions observed for Winnipegosis komatiites. With minor caveats, MELTS is able to successfully reproduce the LLD. This modelling is extended to higher pressures to simulate crystallization of komatiitic melt in an upper crustal magma chamber. All major and rare earth element characteristics of the tholeiitic basalts can be reproduced by ∼60 % crystallization of a Winnipegosis komatiite-like parental melt at pressures of ∼1·5–2·5 kbar at oxygen fugacities between QFM − 1 and QFM + 1, where QFM is the quartz–fayalite–magnetite buffer. Winnipegosis basalts have low Mg# values that are not in equilibrium with mantle peridotite. They therefore cannot represent primary mantle melts derived from cooler mantle than the komatiites, and require fractional crystallization processes in their formation. Furthermore, their trace element characteristics indicate a depth of melting indistinguishable from that of the Winnipegosis komatiites, and derivation from an identical depleted mantle source. All geochemical and geological evidence is therefore consistent with their derivation from a komatiitic melt, and the presence of a large komatiite-derived dunite body in the WKB provides evidence of extensive fractionation of komatiite in the upper crust. The observed uniform basalt compositions are interpreted as the result of a density minimum in the evolving komatiitic melt at temperatures between clinopyroxene and plagioclase saturation, with efficient extraction of melt from a mixture containing ∼60 % crystals. We conclude that the WKB basalts formed by fractional crystallization of a komatiitic parental melt, and suggest that this model may be more broadly applicable to other localities where komatiites and associated basalts show similar geochemical or isotopic characteristics.

Description: SUMMARY We consider the viscoelastic rheology of the solid Earth under the Antarctic Peninsula due to ice mass loss that commenced after the breakup of the Larsen-B ice shelf. We extend the previous analysis of nearby continuous GPS time-series to include five additional years and the additional consideration of the horizontal components of deformation. They show strong uplift from ∼2002 to 2011 followed by reduced uplift rates to 2018. Modelling the GPS-derived uplift as a viscoelastic response to ongoing regional ice unloading from a new ice model confirms earlier estimates of low upper-mantle viscosities of ∼0.3–3 × 1018 Pa s in this region but allows a wide range of elastic lithosphere thickness. The observed and modelled north coordinate component shows little nonlinear variation due to the location of ice mass change to the east of the GPS sites. However, comparison of the observed and modelled east coordinate component constrains the upper-mantle viscosity to be less than ∼9 × 1018 Pa s, consistent with the viscosity range suggested by the uplift rates alone and providing important, largely independent, confirmation of that result. Our horizontal analysis showed only marginal sensitivity to modelled lithospheric thickness. The results for the horizontal components are sensitive to the adopted plate rotation model, with the estimate based on ITRF2014 suggesting that the sum of residual plate motion and pre-2002 glacial isostatic adjustment is likely less than ∼±0.5 mm yr−1 in the east component.

Description: SUMMARY The Jiaodong Peninsula consists of the Jiaobei massif and the Northern Sulu UHP massif. These are separated by the Wulian suture zone (WSZ), a key region for understanding the collision between the North China Craton (NCC) and South China Block (SCB). To interpret this collisional zone, a broad-band seismic profile of 20 stations was installed across the WSZ. Shear wave splitting analysis of teleseismic data revealed a contrast in the splitting patterns beneath different structural zones of the Jiaodong Peninsula. The anisotropic structures of the Jiaobei massif and Northern Sulu UHP massif can be explained by a single anisotropic layer model with WNW-ESE or E-W oriented fast directions. In the WSZ, splitting parameters exhibit pronounced variation in backazimuths indicating a two-layer anisotropy pattern. The lower layer exhibits a WNW-ESE fast direction consistent with that observed in the other two regions. Because the fast direction is generally parallel to the present-day direction of Pacific plate subduction, the anisotropy most likely represents asthenospheric return flow in the big mantle wedge caused by Pacific plate subduction. The upper layer exhibits an NE fast direction, that is, parallel to faulting associated with the WSZ. The lithosphere may preserve fossilized anisotropy induced by the Late Triassic collision of the NCC and SCB even after subsequent destruction and thinning from the Late Mesozoic to Cenozoic. We infer that the WSZ represents a lithospheric-scale structural boundary between the NCC and SCB.

Description: John A Armstrong, Christopher M J Osborne and Lyndsay Fletcher examine how neural networks can be used to explore the nature and location of solar activity.

Description: Zachary Dickeson and Joel Davis review current morphological observations and future questions about the existence, origin and fate of ancient oceans on Mars.

Description: SUMMARY Observations of changes in the geomagnetic field provide unique information about processes in the outer core where the field is generated. Recent geomagnetic field reconstructions based on palaeomagnetic data show persistent westward drift at high northern latitudes at the core–mantle boundary (CMB) over the past 4000 yr, as well as intermittent occurrence of high-latitude weak or reverse flux patches. To further investigate these features, we analysed time-longitude plots of a processed version of the geomagnetic field model pfm9k.1a, filtered to remove quasi-stationary features of the field. Our results suggest that westward drift at both high northern and southern latitudes of the CMB have been a persistent feature of the field over the past 9000 yr. In the Northern Hemisphere we detect two distinct signals with drift rates of 0.09° and 0.25° yr−1 and dominant zonal wavenumbers of m = 2 and 1, respectively. Comparisons with other geomagnetic field models support these observations but also highlight the importance of sedimentary data that provide crucial information on high-latitude geomagnetic field variations. The two distinct drift signals detected in the Northern Hemisphere can largely be decomposed into two westward propagating waveforms. We show that constructive interference between these two waveforms accurately predicts both the location and timing of previously observed high-latitude weak/reverse flux patches over the past 3–4 millennia. In addition, we also show that the 1125-yr periodicity signal inferred from the waveform interference correlates positively with variations in the dipole tilt over the same time period. The two identified drift signals may partially be explained by the westward motion of high-latitude convection rolls. However, the dispersion relation might also imply that part of the drift signal could be caused by magnetic Rossby waves riding on the mean background flow.

Description: On his journey from Cape Town to Oxford and now Cambridge, the exoplanet hunter has also been responsible for disproving a spurious detection. Interview by Paul Johnson.

Description: Leigh N Fletcher shares the research success and future prospects discussed at the UK Solar System Planetary Atmospheres Community Meeting 2020.

Description: The RAS 200 project in partnership with Girlguiding UK and the UK Space Agency is proving successful for all ages, as Sue Bowler reports.

Description: Sue Bowler chats with the organizers of the 13 March RAS Specialist Discussion Meeting “The near-Sun solar wind at solar minimum” about how they made it work for attendees in person and online – and the challenges and opportunities for our communities in the coming months.

Description: SUMMARY Mechanisms of stress transfer and probabilistic models have been widely investigated to explain earthquake clustering features. However, these approaches are still far from being able to link individual events and to determine the number of earthquakes caused by a single event. An alternative approach based on proximity functions allows us to generate hierarchical clustering trees and to identify pairs of nearest-neighbours between consecutive levels of hierarchy. Then, the productivity of an earthquake is the number of events of the next level to which it is linked. Using a relative magnitude threshold ΔM to account for scale invariance in the triggering process, we show that the ΔM-productivity attached to each event is a random variable that follows an exponential distribution. The exponential rate of this distribution does not depend on the magnitude of triggering events and systematically decreases with depth. These results could now be used to characterize active fault systems and improve epidemic models of seismicity.

Description: SUMMARY The existence of a stably stratified layer underneath the core–mantle boundary (CMB) has been recently revived by corroborating evidences coming from seismic studies, mineral physics and thermal evolution models. Such a layer could find its physical origination either in compositional stratification due to the accumulation of light elements at the top or the core or in thermal stratification due to the heat flux becoming locally subadiabatic. The exact properties of this stably stratified layer, namely its size $mathcal {H}_s$ and the degree of its stratification characterized by the Brunt–Väisälä frequency N, are however uncertain and highly debated. A stable layer underneath the CMB can have crucial dynamical impacts on the geodynamo. Because of the inhibition of the convective motions, a stable layer is expected to primarily act as a low-pass filter on the magnetic field, smoothing out the rapidly varying and small-scale features by skin effect. To investigate this effect more systematically, we compute 70 global geodynamo models varying the size of the stably stratified layer from 0 to 300 km and its amplitude from N/Ω = 0 to N/Ω ≃ 50, Ω being the rotation rate. We show that the penetration of the convective flow in the stably stratified layer is controlled by the typical size of the convective eddies and by the local variations of the ratio N/Ω. Using quantitative measures of the degree of morphological semblance between the magnetic field obtained in numerical models and the geomagnetic field at the CMB, we establish an upper bound for the stable layer thickness $mathcal {H}_slt (N/Omega )^{-1} mathcal {L}_s$, $mathcal {L}_s$ being the horizontal size of the convective flow at the base of the stable layer. This defines a strong geomagnetic constraint on the properties of a stably stratified layer beneath the CMB. Unless unaccounted double-diffusive effects could drastically modify the dynamics of the stable layer, our numerical geodynamo models hence favour no stable stratification atop the core.

Description: Human fisheries provide scavengers with abundant and predictable feeding opportunities that may schedule their behavioural patterns. Using miniaturized global positioning system (GPS) tracking technology, we evaluated how Audouin’s gull (Ichthyaetus audouinii), a Mediterranean endemic seabird that makes extensive use of feeding opportunities provided by fisheries, co-occurred (i.e. presumably interacted) with the most important fishing fleets operating off the NE Iberian Peninsula (i.e. diurnal trawlers and nocturnal purse seiners), both in space and time. Results showed that individuals were able to adapt their distribution and activity patterns to the scheduled routines of these fisheries. Waveform analyses based on co-occurring positions revealed that most interactions with trawlers occurred during the afternoon (16:00 h GMT + 1) when discarding occurs as vessels return to port. In contrast, gull-purse seiner interactions largely occurred at night (between 02:00 and 04:00 h) coinciding with the setting and hauling of the nets. Moreover, we found an individual component in seabird–fishery interactions, showing that there may be differential use of fisheries by individuals within the population. In addition to implications for our understanding of the behavioural ecology of this species, these results may have important management implications, particularly under the current European Union Common Fisheries Policy scenario of largely restricting discards.