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: Submarine landslides are mass movements that transport sediment across the continental shelf to the deep ocean. This phenomenon happens when the shear stress exceeds the frictional resistance of the slope. We analyze a variety of seismic attributes to interpret large submarine slide blocks on the North Slope, Alaska. Results show that the slide blocks appear as mounds with scarps associated with them on the seismic section. The slide blocks vary in size, depending on their distance away from the shelf. The pattern of the slide blocks affects the overlying sedimentation. Geological feature: Submarine slide blocks Seismic appearance: Mound-like steep ramp and scarp characteristics on seismic sections; blocky and irregular features with sharp boundaries on the horizon slices and seismic attributes Features with similar appearance: Mass-transport deposits; Remnant blocks; Reef deposits; Submarine channels; Gullies Formation: Torok Formation Age: Cretaceous Location: North Slope, Alaska Seismic data: Obtained from the Alaska Department of Natural Resources, Division of Oil and Gas, through the tax-credit program ( State of Alaska, 2017 , http://dggs.alaska.gov/gmc/seismic-well-data.php ) Analysis tools: Seismic attributes (such as coherent energy, Sobel-filter similarity, dip magnitude, and dip azimuth) and geobody extraction

Description: Linear features at an acute angle with the flight direction are imaged as a series of aligned circular anomalies in the images of the Area 15 aeromagnetic survey, which covered part of the Brazilian southeastern region. These features are interpolation artifacts, a recurring problem found in airborne magnetic images that cause problems for qualitative and quantitative geophysical-geological interpretation. This imaging problem is attributed to spatial aliasing. By running simulations of magnetic data on a synthetic model, we physically proved that the interpolation artifacts from Area 15 are due to inappropriate survey design. Besides the most common expression of artifacts, we described a geologically non-coherent linear pattern as a new type of artifact. Supported by spectral analyses we shown that Area 15 aliased spectrum is similar to geological high-frequency magnetic features, which constitutes a motive for unearthing the correct geophysical signal. Thus, we made use of four techniques for removing the artifacts. The trend enforcement method was able to partially improve the images, while the inverse interpolation method was ineffective, apparently because Area 15 data are severely aliased. Both constrained coherence diffusion and multi-trend gridding methods were able to significantly reduce the presence of artifacts. Despite the high-frequency attenuation, these tools adequately enhanced the magnetic trends and minimized the artifacts. Therefore, the improved images are better suited for a reliable geological interpretation.

Description: Seismic survey data collected for coal gas exploration show that there are many collapse columns distributed in the subsurface of Qinshui Basin, China. The interesting features of the collapse columns are observed by the seismic attributes, including the circular discontinuous patches on the horizon of the Shanxi Formation and multiple parallel discontinuities in vertical profiles of amplitudes. We speculate that the wide presence of these collapse columns are point constraints for the migration and accumulation of coal gas on a large scale. Geological feature: Collapse columns within coal reservoirs Seismic appearance: The coherence illuminates circular/oval discontinuities on the horizon of the Shanxi Formation; the vertical amplitude profiles show cylindrical/funnel-shaped discontinuities. Alternative interpretations: Fault damage zones; velocity pulldown from the overburden Features with similar appearance: Fault-karst in carbonate reservoir; reef pinnacles Formation: Permian Shanxi Formation and Carboniferous Taiyuan Formation Age: Late Permian Location: Qinshui Basin in Shanxi, north-central China Seismic data: Provided by PetroChina Huabei Oilfield Company Contributors: Zonghu Liao, Lin Zhang, and Lianbo Zeng Analysis tools: The seismic amplitude and attribute of coherence from the seismic survey (prestack time migrated)

Description: To quantitatively image fractures with high resolution, we develop an elastic least-squares migration (LSM) algorithm coupled with linear-slip theory, which accurately addresses seismic wave interaction with thin structures. We derive a linearized waveform inversion using the Born approximation to the boundary integral equation for scattered waves, including linear-slip interfaces for P-SV and SH wavefields. Numerical modeling tests assuming a laboratory-scale fracture where a 20 cm long fracture is illuminated by waves with 50 KHz center frequency, show that the proposed LSM successfully estimates the fracture compliances. Furthermore, due to the presence of coupling compliances at the fracture, the results using the proposed LSM show better images than those using the conventional LSM estimating Lamé constants. We also numerically illustrate that the proposed LSM can be successfully applied to dipole acoustic borehole logging data with 3 KHz center frequency for single-well reflection imaging of a 10 m long, dipping fracture embedded in a random background. Finally, we apply the LSM to laboratory experimental data, measuring PP reflections from a fluid-filled fracture. We confirm that the estimated fracture compliances correspond well to those estimated by earlier AVO inversion. Furthermore, the LSM resolves the spatially varying fracture compliances due to local filling of water in the fracture. Because the linear-slip theory can be applied to thin structures in a wide range of scales, high-resolution imaging results and estimated fracture compliance distributions will be crucial to further address small-scale properties at fractures, joints and geological faults.

Description: The Chang 7 oil layer from the upper Triassic Yanchang Formation is an important layer for hydrocarbon exploration. Most studies on the Chang 7 oil layer have focused on the source rocks, while research on the sandstone is still inadequate, especially on the petrography and geochemical characteristics. Using seven sandstone samples of the Chang 7 oil layer in the Yanhe profile, the grain-size analysis, major elements, trace elements, and rare earth elements were tested. The results find that the sandstone of fine-grained sediments of the Chang 7 oil layer is dominated by arkose with a minor number of lithic arkose. The range of grain size (Mz) is from 2.72 to 3.92 Φ, and the C value and M value of the sandstone samples suggest characteristics of turbidity deposition. The Al/Si ratios of all of the samples imply high clay mineral content. The results of trace and rare earth elements demonstrate the reducing condition, freshwater, and cold and dry weather. The provenance of the sandstone samples is mainly from island arc acidic volcanic rock, and the type of provenance is mixed with sedimentary rock, granite, and alkaline basalt. The tectonic background is continental island arc. This study provides a systematic geologic foundation for the formation of sandstone of Chang 7 oil layer in Ordos Basin.

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: We have synthesized macroscopic data (cores) and microscopic data (thin sections and computed tomography images) to characterize volcanic oil and gas reservoir spaces in the Chagan Sag, Yin-E Basin, China. We also have assessed the controlling factors affecting the formation of these reservoirs. Primary pores significantly improve the reservoir quality because they enable the formation of secondary pores (dissolution pores) by allowing migration, alteration, and filling by inorganic and organic fluids. Consequently, the development of secondary pores modifies the permeability of the volcanic rocks. The formation of volcanic reservoirs is influenced by lithology, diagenic minerals, faults, and fractures. The size and abundance of pores, especially gas pores, varies with the lithology. Gas pores are consistently larger and more abundant in basalts than in andesites and tuffs. Moreover, the carbonates and chlorites readily dissolve in the vesicles and fractures, thus resulting in the formation of secondary pores. However, quartzes are not susceptible to dissolution in this study area. Further, faults and fractures are necessary for the formation of reservoir spaces and for oil and gas migration. The tomographic images provide 3D distributions of fractures, vesicles, and matrix dissolution pores. The results suggest that fractures connecting isolated gas pores provide favorable spaces for oil accumulation.

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: Interpretation of recent, high-quality seismic data in the Gulf of Mexico (GOM) has led to competing hypotheses regarding the basin’s rift to drift transition. Some studies suggest a fault-controlled mechanism that ultimately results in mantle exhumation prior to seafloor spreading. Others suggest voluminous magmatic intrusion accommodates the terminal extension phase and results in the extrusion of volcanic seaward dipping reflectors (SDRs). Whereas it has been generally accepted that the plate motions between the rift and drift phases of the GOM are nearly perpendicular to each other, it has not been greatly discussed if the breakup mechanism plays a role in accommodating the transition in plate motion. We have developed a plate kinematic and crustal architecture hypothesis to address the transition from rift to drift in the GOM. We support the proposition of a fault-controlled breakup mechanism, in which slip on a detachment between the crust and mantle may have exhumed the mantle. However, we stress that this mechanism is not exclusive of synrift magmatism, though it does imply that SDRs observed in the GOM are not in this case indicative of a volcanic massif separating attenuated continental and normal oceanic crust. We support our hypothesis through a geometrically realistic 2D potential field model, which includes a magnetic seafloor spreading model constrained by recent published seismic data and analog rock properties. The 2D model suggests that magnetic anomalies near the continent-ocean transition may be related to removal of the lower continental crust during a phase of hyperextension prior to breakup, ending in mantle exhumation. The kinematics of breakup, derived from recent satellite gravity data and constrained by our spreading model and the global plate circuit, suggests that this phase of hyperextension accommodated the change in plate motion direction and a diachronous breakup across the GOM.

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 Frade field, located within the Campos Basin in the southeastern Brazilian margin, is a key oil field that produces from Oligo-Miocene turbidite reservoirs that derived their structural positioning due to the presence of an underlying salt diapir. The evolution of the Frade salt structure was examined using well data, selected 2D lines, and a 3D volume that were interpreted in detail focusing on the Aptian evaporite interval and its influence on the overburden. Analysis of the salt-sediment interaction indicated a complex deformation history that included five main stages of deformation, some assisted by tectonic reactivation episodes. (1) Post-Albian reactivation of a nearby north–northwest-south–southeast basement fault caused the Albian carbonate interval to fault, forming a west–northwest-east–southeast shear zone with a dextral strike-slip component. This movement initiated thin-skinned tectonics that offset the Albian carbonates and formed a pull-apart basin that accommodated a thick Late Cretaceous interval, which weakened the overburden and allowed for the initial formation of the Frade salt diapir. (2) Renewed diapir growth thickened and redistributed the Cenomanian-Maastrichtian sedimentary package proximal to the Frade salt anticline. (3) An initial and localized collapse of the Frade salt anticline occurred during early Paleogene extension. (4) Paleogene shortening caused the salt to flow, resulting in salt withdrawal in the southeast and diapir rejuvenation near its present-day apex, forming several inversion structures. In addition, the Paleogene shortening resulted in a low-relief anticlinal structure that rotated the turbidites into geometries favoring hydrocarbon accumulation. (5) A return to an extensional regime occurred during the late Oligocene/early Miocene. The results of this study provide a new insight into the development of strike-slip salt tectonic structures and show for the first time within the Campos Basin an Albian-level pull-apart basin that formed in association with salt tectonics.

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: Faulting processes have created large damage zones with complex structures in the field; however, estimating the width and geometry of such fault structures in the subsurface is challenging due to a lack of data. Seismic attributes (e.g., coherence and variance) from seismic surveys have been used for the characterization of faults, but most cases do not detail the effectiveness of this approach. By using forward modeling and the associated seismic attributes of variance, four fault models of idealized damage zones are characterized and the frequency effect is evaluated on the width estimation of fault damage zones in the subsurface. The main results indicate that (1) the general geometric pattern of damage zones could be identified by using simulated amplitude and seismic variance with main frequencies of 10, 25, and 40 Hz; (2) the estimated widths of damage zones at a low frequency of 10 Hz are larger (up to twofold) than those at frequencies of 25 and 40 Hz; for large damage zones (〉400 m), the width is best estimated by a frequency of 25 Hz; and (3) scattering noise and diffraction around the fault are found in data at a high frequency of 40 Hz, which results in width overestimation of the damage zones by approximately 17%. The internal structures are difficult to distinguish as scattering noise and chaotic reflections dominate seismic signals. More factors that may influence the accuracy of damage zone width estimation via seismic attributes, include the bedding thickness, fracture density, and velocity. An in-depth understanding of this approach is useful in the application of seismic variance to characterize fault damage zones that may significantly control the fluid migration in the subsurface.

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: The petrophysical properties of rocks have an important influence on shale quality. To characterize the difference of petrophysical properties between organic-rich and organic-poor shale reservoirs, we used the low-field nuclear magnetic resonance (NMR) technique and field emission scanning electron microscopy analysis after argon-ion polishing or natural section to measure porosity and permeability from six core samples from well SY6 in the Sangzhi block, Northwest Hunan province. Some information about pore types, pore structure, residual porosity, movable porosity, and permeability based on the T2 spectrums’ difference of organic-rich and organic-poor shale samples were discussed. The shale sample test results show that the main pores size is mesopore, which provide most of shale gas reservoir space. The continuous peaks demonstrated the pores’ connectivity better than the isolated peaks, and shale gas can migrate freely between these connected pores and fractures. The permeability of all samples calculated by the classic Coates model is extremely low, which is not conducive to the migration of shale gas. We evaluated the dominating factors of NMR porosity and permeability and found that the relationships between NMR porosity and permeability and total organic carbon content, quartz minerals, and clay minerals are not clear, which may be a comprehensive influence. The research results have important guiding significance for shale reservoir quality evaluation in this area.

Description: The synchrosqueezing generalized S-transform (SSGST) is commonly used to generate an isofrequency component of a signal by squeezing the decomposed frequency components of the signal. However, for seismic signals, the single-trace process can have a lack of lateral information in the squeezed results and lead to some discontinuous geologic information that will mislead the interpreter. Thus, to improve the stability of SSGST, we have developed a multichannel seismic trace squeezing method. Multichannel SSGST (MSSGST) considers the decomposed frequency components of neighboring traces of the analysis seismic trace and then reconstructs the center trace. Therefore, compared with SSGST, the multichannel processing improves the stability of the squeezing and produces more laterally continuous results that properly follow the geologic phenomenon. The effectiveness of MSSGST is validated using various field data. We use the application to demonstrate the potential of multichannel squeezing to perform well at (1) improving the energy focusing and continuity of the decomposed frequency components, (2) depicting the boundaries of geologic structures, and (3) identifying the thin layers.

Description: We have carried out petrophysically and geologically guided inversions (PGIs) to jointly invert airborne and ground-based gravity data and airborne magnetic data to recover a quasi-geology model of the DO-27 kimberlite pipe in the Tli Kwi Cho (also referred to as TKC) cluster. DO-27 is composed of three main kimberlite rock types in contact with each other and embedded in a granitic host rock covered by a thin layer of glacial till. The pyroclastic kimberlite (PK), which is diamondiferous, and the volcanoclastic kimberlite (VK) have anomalously low density, due to their high porosity, and weak magnetic susceptibility. They are indistinguishable from each other based upon their potential-field responses. The hypabyssal kimberlite (HK), which is not diamondiferous, has been identified as highly magnetic and remanent. Quantitative petrophysical signatures for each rock unit are obtained from sample measurements, such as the increasing density of the PK/VK unit with depth and the remanent magnetization of the HK unit, and are represented as a Gaussian mixture model (GMM). This GMM guides the PGI toward generating a 3D quasi-geology model with physical properties that satisfies the geophysical data sets and the petrophysical signatures. Density and magnetization models recovered individually yield volumes that have physical property combinations that do not conform to any known petrophysical characteristics of the rocks in the area. A multiphysics PGI addresses this problem by using the GMM as a coupling term, but it puts a volume of the PK/VK unit at a location that is incompatible with geologic information from drillholes. To conform to that geologic knowledge, a fourth unit is introduced, PK-minor, which is petrophysically and geographically distinct from the main PK/VK unit. This inversion produces a quasi-geology model that presents good structural locations of the diamondiferous PK unit and can be used to provide a resource estimate or decide the locations of future drillholes.

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: The inadequate sampling of seismic data in the spatial dimension results in migration artifacts. Conventional least-squares reverse time migration (LSRTM) could improve the image quality. However, even LSRTM will not work in some inadequately sampling situations. To mitigate the impact of migration artifacts, we propose a new LSRTM method with a sparse regularization, which takes advantage of the effective sparse representation of the subsurface reflectivity model in the 2D undecimated wavelet transform (UWT) domain. Different from other sparse regularizations, a sparseness constraint in the 2D UWT domain is applied on the angle slices of the image. To efficiently solve the least-squares inversion problem, we employ an inversion scheme using the conjugate gradient method that uses a soft threshold method to achieve sparse constraint in updating the conjugate gradient direction. Compared with the sparse constraint based on discrete wavelet transform (DWT), the threshold in this method is angle-dependent and is determined according to the energy distribution of each angle slice. To avoid over regularization that can lead to instability and increase the number of iterations, we also apply an exponential threshold strategy. Numerical tests on synthetic datasets demonstrate that our method is capable of improving the image quality by enhancing the resolution and suppressing migration artifacts caused by inadequately sampled seismic data. The method can converge more rapidly than the conventional LSRTM. Since this method performs sparse regularization on several slopes, it achieves better performance on enhancing complex structures with discontinuities like the steeply dipping faults compared with the DWT-based regularization.

Description: Fractured reservoirs, as one kind of unconventional reservoirs, have great potential for oil and gas development, and their accurate characterization requires the development of rock-physics models that better simulate real fractured rocks. However, current models focus mainly on the elastic properties of rocks with aligned cracks, while the effects of randomly orienting cracks in transversely isotropic (TI) rocks are poorly studied even though such conditions are frequently encountered in the earth. To address this problem, we have derived models for the elastic properties of rocks with a TI background permeated by 3D inclined cracks and randomly orienting cracks. Then, based on the developed models, we comprehensively study the effects of the two inclination angles (i.e., the dip angle between the cracks and the isotropic plane and the rotation angle between the cracks and the plane normal to the isotropic plane, respectively) of 3D inclined cracks on the elastic properties of TI rocks. We determine that the two angles have significant influences on the elastic coefficients and hence the elastic velocities, and that their influences on the elastic properties are varying in different directions. We further investigate the effects of crack density and aspect ratio of randomly orienting cracks on the elastic properties of the fractured rocks with a TI background. The results show that the increasing crack density and crack aspect ratio reduce the elastic coefficients and velocities for rocks with randomly orienting cracks, in which the relations between compressional-wave velocities and the crack properties (i.e., crack density and crack aspect ratio) are obtained to aid the interpretation of the acquired acoustic exploration data. The proposed new models can greatly improve the modeling capability for the elastic properties of rocks with a TI background permeated by inclined and randomly orienting cracks.

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: Recently, there have been excellent articles from SEG's senior management team in this space. When I received notification that it was my turn to contribute a column, I immediately began thinking up avoidance tactics. I thought about pretending I didn't get the e-mail with my assignment, or saying, “I'm still the new guy. What do I have to write about?” Alternatively, “I'm an accountant. I build spreadsheets. I don't write columns.” As I began to realize that none of that weaseling was likely to work, I told myself, “I need to find a ghostwriter.”

Description: The following Active Members have been selected by the SEG Committee on Nominations and agreed to be candidates for the 2020–2021 board of directors. Voting will end on 31 July with results announced early-August. The new board will begin its term at the end of the 2020 Annual Meeting. Elected District Representatives will serve a two-year term beginning1 August 2020.

Description: Last month, Annual Meeting Steering Committee Vice Chair Cengiz Esmersoy and SEG Executive Director John Koehr described two SEG 2020 Annual Meeting panels devoted to preparing geophysicists for the future needs of our industry. This month, we continue with a similar theme of preparing future generations for earth science careers and informing and enlightening oil and gas (O&G) geophysicists about important work outside of their native profession.

Description: My term on the SEG Board of Directors started at the SEG Annual Meeting in Houston in 2017, immediately after the historic floods caused by Hurricane Harvey. This unprecedented natural event was often described as “apocalyptic” — but it clearly doesn't compare to the events of 2020.

Description: Welcome to the latest installment of Geophysics Bright Spots. For most of us, the pandemic has upended our daily routines. Personally, a minor silver lining in all of this chaos has been the gain of a few extra hours every week. I have used this time to catch up on the fascinating work being done in the field of artificial intelligence and its applications to numerous disciplines. An emerging trend is physics-informed machine learning, which will help us bridge the gap between traditional theoretical approaches and more recent data-driven methodologies, leading to physically plausible and meaningful results. To follow is a list of research that the editors found interesting in the latest issue of Geophysics. I sincerely hope these articles enlighten you and take your mind off some of the chaos surrounding us.

Description: Deva Prasad Ghosh 1935–2019

Description: How can couples best navigate dual careers? How do you balance work and life throughout a career? How can dual-career couples benefit companies? Eve Sprunt, author of A Guide for Dual-Career Couples, answers these questions and shows how management and individuals alike can activate the power of dual-career couples. The following is an excerpt from episode 81 of SEG's podcast, Seismic Soundoff. Listen to the full episode at https://seg.org/podcast/post/9041 .

Description: Numerical Methods of Exploration Seismology: With Algorithms in MATLAB, by Gary F. Margrave and Michael P. Lamoureux, ISBN 978-1-107-17014-8, 2019, Cambridge University Press, 450 p., US$79.99 (print), US$64 (eBook). Plate Tectonics and Great Earthquakes: 50 Years of Earth-shaking Events, by Lynn R. Sykes, ISBN 978-0-231-18688-9, 2019, Colombia University Press, 272 p., US$35 (print), US$35 (eBook). Geologic Structures of the Arctic Basin, by Alexey Piskarev, Victor Poselov, and Valery Kaminsky, ISBN 978-3-319-77741-2, 2019, Springer, 376 p., US$139.99 (print), US$109 (eBook).

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: Lectures are generally used to teach new knowledge and skills, promote reflection, and stimulate further work and learning. Some advantages of lectures are that they are straightforward, efficient, and can be recorded and recycled. As evidenced by more than 5300 attendees in more than 29 countries during 2019, the SEG Distinguished Lecture, Honorary Lecture, and Virtual Near Surface Global Lecture programs are considered to be some of the most effective outreach programs for the Society.

Description: Seismic 4D monitoring technology has not been as widely employed for gas fields as it has for oil. Many gas fields rely on depletion drive, which has a 4D seismic response that can be uncertain and difficult to predict. On the other hand, aquifer-supported gas fields with measurable water ingress have a reasonable chance of success in terms of generating an interpretable 4D amplitude signal. Pluto gas field in the North West Shelf of Australia falls into this category. Following discovery in 2005, Pluto was appraised by five wells, which found a consistent gas gradient and gas-water contact across the entire field and its various reservoirs. Gas production began in 2012. Time-lapse seismic feasibility studies concluded that gas-saturation changes could be observed with a monitor seismic survey acquired three to four years after first gas. The Pluto 4D Monitor 1 survey was acquired at the start of 2016 and revealed both hardening and softening anomalies. Hardening is interpreted as water ingress (expected) and softening as gas expansion (unexpected). The Pluto 4D results provided important insights into reservoir connectivity and discontinuities. Large hardening anomalies at the TR27 (lower) level can be clearly seen in the data, showing avenues for water ingress. More importantly, a large softening anomaly below the original gas-water contact in the TR29 (upper) reservoir is interpreted to be gas expansion into the aquifer created by a U-tubing effect around a possible barrier in the gas leg. This suggests that the entire TR29 reservoir may not be accessed by the producing PLA04 well. Based on this 4D interpretation, the PLA07 well was drilled and completed in 2019 to produce the TR29 gas updip from the gas expansion anomaly and to increase Pluto field recovery.

Description: Muon tomography is applied to realistic density models of a steam-assisted gravity drainage (SAGD) reservoir at 1.25 and 5 years after initial reservoir production. Forward models of muon count and opacity based on the density models are computed, as well as inverse models of the synthetic muon observations for various simulated detector arrays. The results demonstrate that both phases of reservoir development, namely the rising phase and the spreading phase, can be resolved by muon detectors placed 30 m below the bitumen reservoir at 230 m total vertical depth. The total mass change in the reservoir was recovered from the inversion model and differs from the true mass change by 20%–29%. The spatial distribution of density change shows very good agreement in the horizontal direction, while the vertical is less well constrained in this modeled sensor array configuration. The inverse models provide improved insights into reservoir depletion patterns and indicate muon tomography to be an applicable tool for continuous reservoir monitoring. The numerical modeling approach developed herein is able to model a wide range of SAGD reservoir geometries and detector arrays toward planning of optimized monitoring solutions.

Description: Well-based technologies for seismic subsurface monitoring increasingly utilize fiber-optic cables installed in boreholes as distributed acoustic sensing (DAS) systems. A DAS cable allows measuring linear strain of the fiber and can serve as an array of densely spaced seismic receivers. The strain amplitudes recorded by the DAS cable depend on the near-well formation properties (the softer the medium, the larger the strain). Thus, these properties can be estimated by measuring relative variations of the amplitudes of seismic waves propagating along the well. An advantage of such an approach to subsurface characterization and monitoring is that no active seismic source is required. Passive sources such as earthquakes can be utilized. A synthetic data example demonstrates viability of the approach for monitoring of small-scale CO2 injection into an aquifer. Two field DAS data examples based on signal recordings from several distant earthquakes show that the relevant properties of the near-well formation can be estimated with an accuracy of approximately 5%.

Description: Permanent reservoir monitoring is important for cases of induced seismicity in which there may be a risk to people or to the environment. In such cases, accurately locating microearthquakes and assessing their hazard level can help keep production at safe levels. The process can benefit greatly from the use of automation. With the shift toward full-waveform microearthquake location algorithms and workflows, greater accuracy and information can be retrieved compared to that offered by traditional traveltime estimation techniques, but the complexity of these workflows and run-time costs can be higher. Results are presented from an automatic elastic event location and moment tensor inversion workflow that has been highly parallelized on clustered computer hardware. Run times that previously took up to several days to complete using a manually intensive execution of the workflow now can be achieved in approximately 1 hour. Some 180 events recorded at the Groningen gas field and ranging in magnitude from 0.1 to 3.4 MW (ML) have been located and analyzed with the automatic workflow. The results indicate equivalent location accuracy when compared to the manually intensive workflow execution. However, larger errors are noted in the depth positions of some events and in the range and nature of the focal mechanism, as derived from moment tensor inversion. High grading of the manual and automatic results has been performed and used to study the geomechanical behavior of the microearthquakes in the Groningen region, which exhibit mainly dip-slip, double-couple motion, in areas of previous production activity.

Description: The P-Cable high-resolution 3D marine acquisition system tows many short, closely separated streamers behind a small source. It can provide 3D seismic data of very high temporal and spatial resolution. Since the system is containerized and has small dimensions, it can be deployed at short notice and relatively low cost, making it attractive for time-lapse seismic reservoir monitoring. During acquisition of a 3D high-resolution survey in the Gulf of Mexico in 2014, a pair of sail lines were repeated to form a time-lapse seismic test. We processed these in 2019 to evaluate their geometric and seismic repeatability. Geometric repetition accuracy was excellent, with source repositioning errors below 10 m and bin-based receiver positioning errors below 6.25 m. Seismic data comparisons showed normalized root-mean-square difference values below 10% between 40 and 150 Hz. Refinements to the acquisition system since 2014 are expected to further improve repeatability of the low-frequency components. Residual energy on 4D difference seismic data was low, and timing stability was good. We conclude that the acquisition system is well suited to time-lapse seismic surveying in areas where the reservoir and time-lapse seismic signal can be adequately imaged by small-source, short-offset, low-fold data.

Description: A legacy seismic data set from 2001 and 2011 was used for time-lapse interpretation over a Brazilian carbonate reservoir in Field-X of Campos Basin. The acquired 4D seismic data set was noisy and initially deemed to be uninterpretable. Pressure data were limited, and the initial simulation model was poorly calibrated. All of these challenges warranted the need to establish an interdisciplinary interpretation workflow. In this paper, we introduce a tiered integrated approach to optimize data value from multiple sources. The results of interwell connectivity from production data analysis were used as a basis for assignment of 4D signal confidence flags, which were later combined with simulation-to-seismic modeling of history-matched realizations, enabling interpretation efforts of the noisy 4D seismic data set. This integrated approach led to identification of the injected water front and a potential sweet spot for an infill well.

Description: Time-lapse monitoring using 3D distributed acoustic sensing vertical seismic profiles (DAS VSPs) is rapidly maturing as a nonintrusive low-cost solution for target-oriented monitoring in deep water. In a Gulf of Mexico field, DAS fibers deployed in active wells enable detailed tracking of the water flood in two deep reservoirs. Multiple tests in adverse well conditions let us understand the impact of source size and other factors on the spatially dependent quality of time-lapse DAS data and prove that excellent image repeatability is achievable under typical field conditions. Frequent repeat surveys allowed us to predict the timing of water arrival in a producer and to observe new water injection patterns that are important for understanding water-flood performance. Going forward, DAS 4D monitoring is envisioned as a tool that can assist with proactive wells and reservoir management, new well planning, and reservoir model updates.

Description: The papers submitted to this special section demonstrate that the topic of reservoir monitoring is extremely diverse. This diversity is reflected in the wide range of geologic settings covered by these papers — deepwater unconsolidated clastics, more cemented sandstones in onshore fields, and carbonates. Diversity is seen in the range of production scenarios described by these papers — water sweep of oil and gas, thermal recovery using steam-assisted gravity drainage (SAGD), and enhanced recovery using CO2 injection. Moreover, the papers in this section cover much more than time-lapse 3D seismic. Although about half of the submitted papers use 4D surface seismic data to monitor reservoirs, the remainder cover a diversity of methods that include time-lapse vertical seismic profiles (VSPs), repeat well logging using distributed acoustic sensing (DAS), and muon tomography. Even the concept of the “reservoir” is expanded to include monitoring microseismicity that might result from production activity.

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.