ALBERT

Description: The relationships between trachytes and peralkaline rhyolites (i.e. pantellerites and comendites), which occur in many continental rift systems, oceanic islands and continental intraplate settings, is unclear. To fill this gap, we have performed phase equilibrium experiments on two representative metaluminous trachytes from Pantelleria to determine both their pre-eruptive equilibration conditions (pressure, temperature, H2O content and redox state) and liquid lines of descent. Experiments were performed in the temperature range 750–950 C, pressure 0 5–1 5 kbar and fluid saturation conditions with XH2O [¼H2O/(H2OþCO2)] ranging between zero and unity. Redox conditions were fixed below the nickel–nickel oxide buffer (NNO). The results show that at 950 C and melt water contents (H2Omelt) close to saturation, trachytes are at liquidus conditions at all pressures. Clinopyroxene is the liquidus phase, being followed by iron-rich olivine and alkali feldspar. Comparison of experimental and natural phases (abundances and compositions) yields the following pre-eruptive conditions: P¼160 5 kbar, T¼925625 C, H2Omelt¼261wt %, and fO2 between NNO– 0 5 and NNO– 2. A decrease in temperature from 950 C to 750 C, as well as of H2Omelt, promotes a massive crystallization of alkali feldspar to over 80 wt %. Iron-bearing minerals show gradual iron enrichment when T and fO2 decrease, trending towards the compositions of the phenocrysts of natural pantellerites. Despite the metaluminous character of the bulk-rock compositions, residual glasses obtained after 80 wt % crystallization evolve toward comenditic compositions, owing to profuse alkali feldspar crystallization, which decreases the Al2O3 of the melt, leading to a consequent increase in the peralkalinity index [PI¼molar (Na2OþK2O)/Al2O3]. This is the first experimental demonstration that peralkaline felsic derivatives can be produced by low-pressure fractional crystallization of metaluminous mafic magmas. Our results show that the pantelleritic magmas of basalt–trachyte–rhyolite igneous suites require at least 95 wt % of parental basalt crystallization, consistent with trace element evidence. Redox conditions, through their effect on Fe–Ti oxide stabilities, control the final iron content of the evolving melt.

Description: Published

Description: 559- 588

Description: 2V. Struttura e sistema di alimentazione dei vulcani

Description: JCR Journal

Description: © The Author(s), 2018. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Frontiers in Marine Science 5 (2018): 90, doi:10.3389/fmars.2018.00090.

Description: Sea turtles inhabiting coastal environments routinely encounter anthropogenic hazards, including fisheries, vessel traffic, pollution, dredging, and drilling. To support mitigation of potential threats, it is important to understand fine-scale sea turtle behaviors in a variety of habitats. Recent advancements in autonomous underwater vehicles (AUVs) now make it possible to directly observe and study the subsurface behaviors and habitats of marine megafauna, including sea turtles. Here, we describe a “smart” AUV capability developed to study free-swimming marine animals, and demonstrate the utility of this technology in a pilot study investigating the behaviors and habitat of leatherback turtles (Dermochelys coriacea). We used a Remote Environmental Monitoring UnitS (REMUS-100) AUV, designated “TurtleCam,” that was modified to locate, follow and film tagged turtles for up to 8 h while simultaneously collecting environmental data. The TurtleCam system consists of a 100-m depth rated vehicle outfitted with a circular Ultra-Short BaseLine receiver array for omni-directional tracking of a tagged animal via a custom transponder tag that we attached to the turtle with two suction cups. The AUV collects video with six high-definition cameras (five mounted in the vehicle nose and one mounted aft) and we added a camera to the animal-borne transponder tag to record behavior from the turtle's perspective. Since behavior is likely a response to habitat factors, we collected concurrent in situ oceanographic data (bathymetry, temperature, salinity, chlorophyll-a, turbidity, currents) along the turtle's track. We tested the TurtleCam system during 2016 and 2017 in a densely populated coastal region off Cape Cod, Massachusetts, USA, where foraging leatherbacks overlap with fixed fishing gear and concentrated commercial and recreational vessel traffic. Here we present example data from one leatherback turtle to demonstrate the utility of TurtleCam. The concurrent video, localization, depth and environmental data allowed us to characterize leatherback diving behavior, foraging ecology, and habitat use, and to assess how turtle behavior mediates risk to impacts from anthropogenic activities. Our study demonstrates that an AUV can successfully track and image leatherback turtles feeding in a coastal environment, resulting in novel observations of three-dimensional subsurface behaviors and habitat use, with implications for sea turtle management and conservation.

Description: This research was funded by National Oceanic and Atmospheric Administration Grant #NA16NMF4720074 to the Massachusetts Division of Marine Fisheries under the Species Recovery Grants to States program. Additional funding was provided by Jean Tempel, Hydroid Inc., and over 100 Project WHOI donors.

Description: © The Author(s), 2018. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Frontiers in Marine Science 5 (2018): 241, doi:10.3389/fmars.2018.00241.

Description: Cryptophyte algae are globally distributed photosynthetic flagellates found in freshwater, estuarine, and neritic ecosystems. While cryptophytes can be highly abundant and are consumed by a wide variety of protistan predators, few studies have sought to quantify in situ grazing rates on their populations. Here we show that autumnal grazing rates on in situ communities of cryptophyte algae in Chesapeake Bay are high throughout the system, while growth rates, particularly in the lower bay, were low. Analysis of the genetic diversity of cryptophyte populations within dilution experiments suggests that microzooplankton may be selectively grazing the fastest-growing members of the population, which were generally Teleaulax spp. We also demonstrate that potential grazing rates of ciliates and dinoflagellates on fluorescently labeled (FL) Rhodomonas salina, Storeatula major, and Teleaulax amphioxeia can be high (up to 149 prey predator−1 d−1), and that a Gyrodinium sp. and Mesodinium rubrum could be selective grazers. Potential grazing was highest for heterotrophic dinoflagellates, but due to its abundance, M. rubrum also had a high overall impact. This study reveals that cryptophyte algae in Chesapeake Bay can experience extremely high grazing pressure from phagotrophic protists, and that this grazing likely shapes their community diversity.

Description: The authors thank the National Science Foundation (OCE 1031718 and 1436169) for providing support for this research.

Description: © The Author(s), 2018. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Frontiers in Marine Science 5 (2018): 168, doi:10.3389/fmars.2018.00168.

Description: North Atlantic right whales (Eubalaena glacialis) are highly endangered and frequently exposed to a myriad of human activities and stressors in their industrialized habitat. Entanglements in fixed fishing gear represent a particularly pervasive and often drawn-out source of anthropogenic morbidity and mortality to the species. To better understand both the physiological response to entanglement, and to determine fundamental parameters such as acquisition, duration, and severity of entanglement, we measured a suite of biogeochemical markers in the baleen of an adult female that died from a well-documented chronic entanglement in 2005 (whale Eg2301). Steroid hormones (cortisol, corticosterone, estradiol, and progesterone), thyroid hormones (triiodothyronine (T3) and thyroxine (T4)), and stable isotopes (δ13C and δ15N) were all measured in a longitudinally sampled baleen plate. This yielded an 8-year profile of foraging and migration behavior, stress response, and reproduction. Stable isotopes cycled in annual patterns that reflect the animal's north-south migration behavior and seasonally abundant zooplankton diet. A progesterone peak, lasting approximately 23 months, was associated with the single known calving event (in 2002) for this female. Estradiol, cortisol, corticosterone, T3, and T4 were also elevated, although variably so, during the progesterone peak. This whale was initially sighted with a fishing gear entanglement in September 2004, but the hormone panel suggests that the animal first interacted with the gear as early as June 2004. Elevated δ15N, T3, and T4 indicate that Eg2301 potentially experienced increased energy expenditure, significant lipid catabolism, and thermal stress approximately 3 months before the initial sighting with fishing gear. All hormones in the panel (except cortisol) were elevated above baseline by September 2004. This novel study illustrates the value of using baleen to reconstruct recent temporal profiles and as a comparative matrix in which key physiological indicators of individual whales can be used to understand the impacts of anthropogenic activity on threatened whale populations.

Description: The Woods Hole Oceanographic Institution's Ocean Life Institute and Marine Mammal Center funded this study and NL was supported by a Postdoctoral Fellowship from Baylor University.

Description: © The Author(s), 2018. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Frontiers in Microbiology 9 (2018): 772, doi:10.3389/fmicb.2018.00772.

Description: © The Author(s), 2018. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Frontiers in Marine Science 5 (2018): 25, doi:10.3389/fmars.2018.00025.

Description: Basking shark (Cetorhinus maximus) populations are considered “vulnerable” globally and “endangered” in the northeast Atlantic by the International Union for the Conservation of Nature (IUCN). Much of our knowledge of this species comes from surface observations in coastal waters, yet recent evidence suggests the majority of their lives may be spent in the deep ocean. Depth preferences of basking sharks have significantly limited movement studies that used pop-up satellite archival transmitting (PSAT) tags as conventional light-based geolocation is impossible for tagged animals that spend significant time below the photic zone. We tagged 57 basking sharks with PSAT tags in the NW Atlantic from 2004 to 2011. Many individuals spent several months at meso- and bathy-pelagic depths where accurate light-level geolocation was impossible during fall, winter and spring. We applied a newly-developed geolocation approach for the PSAT data by comparing three-dimensional depth-temperature profile data recorded by the tags to modeled in situ oceanographic data from the high-resolution HYbrid Coordinate Ocean Model (HYCOM). Observation-based likelihoods were leveraged within a state-space hidden Markov model (HMM). The combined tracks revealed that basking sharks moved from waters around Cape Cod, MA to as far as the SE coast of Brazil (20°S), a total distance of over 17,000 km. Moreover, 59% of tagged individuals with sufficient deployment durations (〉250 days) demonstrated seasonal fidelity to Cape Cod and the Gulf of Maine, with one individual returning to within 60 km of its tagging location 1 year later. Tagged sharks spent most of their time at epipelagic depths during summer months around Cape Cod and in the Gulf of Maine. During winter months, sharks spent extended periods at depths of at least 600 m while moving south to the Sargasso Sea, the Caribbean Sea, or the western tropical Atlantic. Our work demonstrates the utility of applying advances in oceanographic modeling to understanding habitat use of highly migratory, often meso- and bathy-pelagic, ocean megafauna. The large-scale movement patterns of tagged sharks highlight the need for international cooperation when designing and implementing conservation strategies to ensure that the species recovers from the historical effects of over-fishing throughout the North Atlantic Ocean.

Description: We gratefully acknowledge funding from the US National Science Foundation (OCE 0825148), the National Aeronautics and Space Administration (NNS06AA96G), the Massachusetts Environmental Trust, and the Federal Aid in Sport Fish Restoration Program. CB was funded by the Martin Family Society of Fellows for Sustainability Fellowship at the Massachusetts Institute of Technology, the Grassle Fellowship and Ocean Venture Fund at the Woods Hole Oceanographic Institution, and the NASA Earth and Space Science Fellowship. Funding for the development of HYCOM has been provided by the National Ocean Partnership Program and the Office of Naval Research.

Description: © The Author(s), 2018. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Frontiers in Microbiology 9 (2018): 840, doi:10.3389/fmicb.2018.00840.

Description: Earth’s subsurface environment is one of the largest, yet least studied, biomes on Earth, and many questions remain regarding what microorganisms are indigenous to the subsurface. Through the activity of the Census of Deep Life (CoDL) and the Deep Carbon Observatory, an open access 16S ribosomal RNA gene sequence database from diverse subsurface environments has been compiled. However, due to low quantities of biomass in the deep subsurface, the potential for incorporation of contaminants from reagents used during sample collection, processing, and/or sequencing is high. Thus, to understand the ecology of subsurface microorganisms (i.e., the distribution, richness, or survival), it is necessary to minimize, identify, and remove contaminant sequences that will skew the relative abundances of all taxa in the sample. In this meta-analysis, we identify putative contaminants associated with the CoDL dataset, recommend best practices for removing contaminants from samples, and propose a series of best practices for subsurface microbiology sampling. The most abundant putative contaminant genera observed, independent of evenness across samples, were Propionibacterium, Aquabacterium, Ralstonia, and Acinetobacter. While the top five most frequently observed genera were Pseudomonas, Propionibacterium, Acinetobacter, Ralstonia, and Sphingomonas. The majority of the most frequently observed genera (high evenness) were associated with reagent or potential human contamination. Additionally, in DNA extraction blanks, we observed potential archaeal contaminants, including methanogens, which have not been discussed in previous contamination studies. Such contaminants would directly affect the interpretation of subsurface molecular studies, as methanogenesis is an important subsurface biogeochemical process. Utilizing previously identified contaminant genera, we found that ∼27% of the total dataset were identified as contaminant sequences that likely originate from DNA extraction and DNA cleanup methods. Thus, controls must be taken at every step of the collection and processing procedure when working with low biomass environments such as, but not limited to, portions of Earth’s deep subsurface. Taken together, we stress that the CoDL dataset is an incredible resource for the broader research community interested in subsurface life, and steps to remove contamination derived sequences must be taken prior to using this dataset.

Description: We wish to acknowledge the support of the Sloan Foundation and the Deep Carbon Observatory and the Department of Energy, Office of Fossil Energy (Colwell).

Description: © The Author(s), 2018. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Frontiers in Earth Science 6 (2018): 147, doi:10.3389/feart.2018.00147.

Description: Silicic effusive eruptions in deep submarine environments have not yet been directly observed and very few modern submarine silicic lavas and domes have been described. The eruption of Havre caldera volcano in the Kermadec arc in 2012 provided an outstanding database for research on deep submarine silicic effusive eruptions because it produced 15 rhyolite (70–72 wt.% SiO2) lavas and domes with a total volume of ∼0.21 km3 from 14 separate seafloor vents. Moreover, in 2015, the seafloor products were observed, mapped and sampled in exceptional detail (1-m resolution) using AUV Sentry and ROV Jason2 deployed from R/V Roger Revelle. Vent positions are strongly aligned, defining NW-SE and E-W trends along the southwestern and southern Havre caldera margin, respectively. The alignment of the vents suggests magma ascent along dykes which probably occupy faults related to the caldera margin. Four vents part way up the steeply sloping southwestern caldera wall at 1,200–1,300 m below sea level (bsl) and one on the caldera rim (1,060 m bsl) produced elongate lavas. On the steep caldera wall, the lavas consist of narrow tongues that have triangular cross-section shapes. Two of the narrow-tongue segments are connected to wide lobes on the flat caldera floor at ∼1,500 m bsl. The lavas are characterized by arcuate surface ridges oriented perpendicular to the propagation direction. Eight domes were erupted onto relatively flat sea floor from vents at ∼1,000 m bsl along the southern and southwestern caldera rim. They are characterized by steep margins and gently convex-up upper surfaces. With one exception, the domes have narrow spines and deep clefts above the inferred vent positions. One dome has a relatively smooth upper surface. The lavas and domes all consist of combinations of coherent rhyolite and monomictic rhyolite breccia. Despite eruption from deep-water vents (most 〉900 m bsl), the Havre 2012 rhyolite lavas and domes are very similar to subaerial rhyolite lavas and domes in terms of dimensions, volumes, aspect ratio, textures and morphology. They show that lava morphology was strongly controlled by the pre-existing seafloor topography: domes and wide lobes formed where the rhyolite was emplaced onto flat sea floor, whereas narrow tongues formed where the rhyolite was emplaced on the steep slopes of the caldera wall.

Description: This research was funded by an Australian Research Council Postdoctoral fellowship to RJC (DP110102196 and DE150101190), and National Science Foundation grants OCE1357443 and OCE1357216. FI was supported by a Tasmanian Government Postgraduate Award.

Description: © The Author(s), 2018. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Frontiers in Marine Science 5 (2018): 49, doi:10.3389/fmars.2018.00049.

Description: Species inhabiting deep-sea hydrothermal vents are strongly influenced by the geological setting, as it provides the chemical-rich fluids supporting the food web, creates the patchwork of seafloor habitat, and generates catastrophic disturbances that can eradicate entire communities. The patches of vent habitat host a network of communities (a metacommunity) connected by dispersal of planktonic larvae. The dynamics of the metacommunity are influenced not only by birth rates, death rates and interactions of populations at the local site, but also by regional influences on dispersal from different sites. The connections to other communities provide a mechanism for dynamics at a local site to affect features of the regional biota. In this paper, we explore the challenges and potential benefits of applying metacommunity theory to vent communities, with a particular focus on effects of disturbance. We synthesize field observations to inform models and identify data gaps that need to be addressed to answer key questions including: (1) what is the influence of the magnitude and rate of disturbance on ecological attributes, such as time to extinction or resilience in a metacommunity; (2) what interactions between local and regional processes control species diversity, and (3) which communities are “hot spots” of key ecological significance. We conclude by assessing our ability to evaluate resilience of vent metacommunities to human disturbance (e.g., deep-sea mining). Although the resilience of a few highly disturbed vent systems in the eastern Pacific has been quantified, these values cannot be generalized to remote locales in the western Pacific or mid Atlantic where disturbance rates are different and information on local controls is missing.

Description: LM was supported by NSF OCE 1356738 and DEB 1558904. SB was supported by the NSF DEB 1558904 and the Investment in Science Fund at Woods Hole Oceanographic Institution. MB was supported by the Austrian Science Fund grants P20190-B17 and P16774-B03. LL was supported by NSF OCE 1634172 and the JM Kaplan Fund. MN was supported by NSF DEB 1558904. Y-JW was supported by a Korean Institute of Ocean Science and Technology (KIOST) grant PM60210.

Description: © The Author(s), 2018. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Frontiers in Marine Science 5 (2018): 158, doi:10.3389/fmars.2018.00158.

Description: In autumn 2015, several sources reported observations of large amounts of gelatinous material in a large north Norwegian fjord system, either caught when trawling for other organisms or fouling fishing gear. The responsible organism was identified as a physonect siphonophore, Nanomia cara, while a ctenophore, Beroe cucumis, and a hydromedusa, Modeeria rotunda, were also registered in high abundances on a couple of occasions. To document the phenomena, we have compiled a variety of data from concurrent fisheries surveys and local fishermen, including physical samples, trawl catch, and acoustic data, photo and video evidence, and environmental data. Because of the gas-filled pneumatophore, characteristic for these types of siphonophores, acoustics provided detailed and unique insight to the horizontal and vertical distribution and potential abundances (~0.2–20 colonies·m−3) of N. cara with the highest concentrations observed in the near bottom region at ~320 m depth in the study area. This suggests that these animals were retained and accumulated in the deep basins of the fjord system possibly blooming here because of favorable environmental conditions and potentially higher prey availability compared to the shallower shelf areas to the north. Few cues as to the origin and onset of the bloom were found, but it may have originated from locally resident siphonophores. The characteristics of the deep-water masses in the fjord basins were different compared to the deep water outside the fjord system, suggesting no recent deep-water import to the fjords. However, water-masses containing siphonophores (not necessarily very abundant), may have been additionally introduced to the fjords at intermediate depths, with the animals subsequently trapped in the deeper fjord basins. The simultaneous observations of abundant siphonophores, hydromedusae, and ctenophores in the Lyngen-Kvænangen fjord system are intriguing, but difficult to provide a unified explanation for, as the organisms differ in their biology and ecology. Nanomia and Beroe spp. are holopelagic, while M. rotunda has a benthic hydroid stage. The species also have different trophic ecologies and dietary preferences. Only by combining information from acoustics, trawling, genetics, and local fishermen, were the identity, abundance, and the vertical and horizontal distribution of the physonect siphonophore, N. cara, established.

Description: The work was funded by the Ministry of Fisheries and Coastal Affairs through the Institute of Marine Research (IMR), while the Research Council of Norway (RCN) is thanked for the financial support through the project The Arctic Ocean Ecosystem—(SI_ARCTIC, RCN 228896). AH was supported by the Norwegian Taxonony Initiative (NTI 70184233) and ForBio Research School funding (RCN 248799 and NTI 70184215).

Description: © The Author(s), 2018. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Frontiers in Microbiology 9 (2018): 1201, doi:10.3389/fmicb.2018.01201.

Description: Interactions between microorganisms and algae during bloom events significantly impacts their physiology, alters ambient chemistry, and shapes ecosystem diversity. The potential role these interactions have in bloom development and decline are also of particular interest given the ecosystem impacts of algal blooms. We hypothesized that microbial community structure and succession is linked to specific bloom stages, and reflects complex interactions among taxa comprising the phycosphere environment. This investigation used pyrosequencing and correlation approaches to assess patterns and associations among bacteria, archaea, and microeukaryotes during a spring bloom of the dinoflagellate Alexandrium catenella. Within the bacterial community, Gammaproteobacteria and Bacteroidetes were predominant during the initial bloom stage, while Alphaproteobacteria, Cyanobacteria, and Actinobacteria were the most abundant taxa present during bloom onset and termination. In the archaea biosphere, methanogenic members were present during the early bloom period while the majority of species identified in the late bloom stage were ammonia-oxidizing archaea and Halobacteriales. Dinoflagellates were the major eukaryotic group present during most stages of the bloom, whereas a mixed assemblage comprising diatoms, green-algae, rotifera, and other microzooplankton were present during bloom termination. Temperature and salinity were key environmental factors associated with changes in bacterial and archaeal community structure, respectively, whereas inorganic nitrogen and inorganic phosphate were associated with eukaryotic variation. The relative contribution of environmental parameters measured during the bloom to variability among samples was 35.3%. Interaction analysis showed that Maxillopoda, Spirotrichea, Dinoflagellata, and Halobacteria were keystone taxa within the positive-correlation network, while Halobacteria, Dictyochophyceae, Mamiellophyceae, and Gammaproteobacteria were the main contributors to the negative-correlation network. The positive and negative relationships were the primary drivers of mutualist and competitive interactions that impacted algal bloom fate, respectively. Functional predictions showed that blooms enhance microbial carbohydrate and energy metabolism, and alter the sulfur cycle. Our results suggest that microbial community structure is strongly linked to bloom progression, although specific drivers of community interactions and responses are not well understood. The importance of considering biotic interactions (e.g., competition, symbiosis, and predation) when investigating the link between microbial ecological behavior and an algal bloom’s trajectory is also highlighted.

Description: This work was supported by NSFC (41476092, 41741015), S&T Projects of Shenzhen Science and Technology Innovation Committee (JCYJ20150831192329178, JCYJ20170817160708491, and JCYJ20170412171959157), Key Research and Development Plan of Ministry of Science and Technology of China (2017YFC1403600), as well as by the Woods Hole Center for Oceans and Human Health through the National Science Foundation (Grant OCE-1314642), and National Institute of Environmental Health Sciences (Grant 1-P01-ES021923-01).

Description: © The Author(s), 2018. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Frontiers in Physiology 9 (2018): 838, doi: 10.3389/fphys.2018.00838.

Description: Bottlenose dolphins (Tursiops truncatus) are highly versatile breath-holding predators that have adapted to a wide range of foraging niches from rivers and coastal ecosystems to deep-water oceanic habitats. Considerable research has been done to understand how bottlenose dolphins manage O2 during diving, but little information exists on other gases or how pressure affects gas exchange. Here we used a dynamic multi-compartment gas exchange model to estimate blood and tissue O2, CO2, and N2 from high-resolution dive records of two different common bottlenose dolphin ecotypes inhabiting shallow (Sarasota Bay) and deep (Bermuda) habitats. The objective was to compare potential physiological strategies used by the two populations to manage shallow and deep diving life styles. We informed the model using species-specific parameters for blood hematocrit, resting metabolic rate, and lung compliance. The model suggested that the known O2 stores were sufficient for Sarasota Bay dolphins to remain within the calculated aerobic dive limit (cADL), but insufficient for Bermuda dolphins that regularly exceeded their cADL. By adjusting the model to reflect the body composition of deep diving Bermuda dolphins, with elevated muscle mass, muscle myoglobin concentration and blood volume, the cADL increased beyond the longest dive duration, thus reflecting the necessary physiological and morphological changes to maintain their deep-diving life-style. The results indicate that cardiac output had to remain elevated during surface intervals for both ecotypes, and suggests that cardiac output has to remain elevated during shallow dives in-between deep dives to allow sufficient restoration of O2 stores for Bermuda dolphins. Our integrated modeling approach contradicts predictions from simple models, emphasizing the complex nature of physiological interactions between circulation, lung compression, and gas exchange.

Description: AF (N00014-17-1-2756), PT (N000141512553) and FHJ (N00014-14-1-0410) were supported by the Office of Naval Research, and FHJ by an AIASCOFUND fellowship from Aarhus Institute of Advanced Studies, Aarhus University, under EU's FP7 program (Agreement No. 609033). PT received funding from the MASTS pooling initiative (The Marine Alliance for Science and Technology for Scotland) and their support is gratefully acknowledged.

Description: © The Author(s), 2018. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Conservation Physiology 6 (2018): coy049, doi:10.1093/conphys/coy049.

Description: Male baleen whales have long been suspected to have annual cycles in testosterone, but due to difficulty in collecting endocrine samples, little direct evidence exists to confirm this hypothesis. Potential influences of stress or adrenal stress hormones (cortisol, corticosterone) on male reproduction have also been difficult to study. Baleen has recently been shown to accumulate steroid hormones during growth, such that a single baleen plate contains a continuous, multi-year retrospective record of the whale’s endocrine history. As a preliminary investigation into potential testosterone cyclicity in male whales and influences of stress, we determined patterns in immunoreactive testosterone, two glucocorticoids (cortisol and corticosterone), and stable-isotope (SI) ratios, across the full length of baleen plates from a bowhead whale (Balaena mysticetus), a North Atlantic right whale (Eubalaena glacialis) and a blue whale (Balaenoptera musculus), all adult males. Baleen was subsampled at 2 cm (bowhead, right) or 1 cm (blue) intervals and hormones were extracted from baleen powder with methanol, followed by quantification of all three hormones using enzyme immunoassays validated for baleen extract of these species. Baleen of all three males contained regularly spaced peaks in testosterone content, with number and spacing of testosterone peaks corresponding well to SI data and to species-specific estimates of annual baleen growth rate. Cortisol and corticosterone exhibited some peaks that co-occurred with testosterone peaks, while other glucocorticoid peaks occurred independent of testosterone peaks. The right whale had unusually high glucocorticoids during a period with a known entanglement in fishing gear and a possible disease episode; in the subsequent year, testosterone was unusually low. Further study of baleen testosterone patterns in male whales could help clarify conservation- and management-related questions such as age of sexual maturity, location and season of breeding, and the potential effect of anthropogenic and natural stressors on male testosterone cycles.

Description: This work was supported by (1) the Arizona Board of Regents Technology Research Initiative Fund; (2) the Center for Bioengineering Innovation at Northern Arizona University; (3) the Greenland Institute of Natural Resources; (4) the Woods Hole Oceanographic Institution Ocean Life Institute and (5) Fisheries and Ocean Canada’s (DFO) Priorities and Partnership Strategic Initiatives Fund and Oceans Protection Plan.

Description: © The Author(s), 2018. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Frontiers in Marine Science 5 (2018): 13, doi:10.3389/fmars.2018.00013.

Description: Identifying putative mixotrophic protist species in the environment is important for understanding their behavior, with the recovery of these species in culture essential for determining the triggers of feeding, grazing rates, and overall impact on bacterial standing stocks. In this project, mixotroph abundances determined using tracer ingestion in water and sea ice samples collected in the Ross Sea, Antarctica during the summer of 2011 were compared with data from the spring (Ross Sea) and fall (Arctic) to examine the impacts of bacterivory/mixotrophy. Mixotrophic nanoplankton (MNAN) were usually less abundant than heterotrophs, but consumed more of the bacterial standing stock per day due to relatively higher ingestion rates (1–7 bacteria mixotroph−1 h−1 vs. 0.1–4 bacteria heterotroph−1 h−1). Yet, even with these high rates observed in the Antarctic summer, mixotrophs appeared to have a smaller contribution to bacterivory than in the Antarctic spring. Additionally, putative mixotroph taxa were identified through incubation experiments accomplished with bromodeoxyuridine-labeled bacteria as food, immunoprecipitation (IP) of labeled DNA, and amplification and high throughput sequencing of the eukaryotic ribosomal V9 region. Putative mixotroph OTUs were identified in the IP samples by taxonomic similarity to known phototroph taxa. OTUs that had increased abundance in IP samples compared to the non-IP samples from both surface and chlorophyll maximum (CM) depths were considered to represent active mixotrophy and include ones taxonomically similar to Dictyocha, Gymnodinium, Pentapharsodinium, and Symbiodinium. These OTUs represent target taxa for isolation and laboratory experiments on triggers for mixotrophy, to be combined with qPCR to estimate their abundance, seasonal distribution and potential impact.

Description: This work was supported by National Science Foundation Grants OPP-0838955 (RG) and OPP-0838847 (RS).

Description: © The Author(s), 2018. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Frontiers in Marine Science 5 (2018): 61, doi:10.3389/fmars.2018.00061.

Description: The distribution of dissolved iron (Fe), total organic Fe-binding ligands, and siderophores were measured between the surface and 400 m at Station ALOHA, a long term ecological study site in the North Pacific Subtropical Gyre. Dissolved Fe concentrations were low throughout the water column and strong organic Fe-binding ligands exceeded dissolved Fe at all depths; varying from 0.9 nmol L−1 in the surface to 1.6 nmol L−1 below 150 m. Although Fe does not appear to limit microbial production, we nevertheless found siderophores at nearly all depths, indicating some populations of microbes were responding to Fe stress. Ferrioxamine siderophores were most abundant in the upper water column, with concentrations between 0.1 and 2 pmol L−1, while a suite of amphibactins were found below 200 m with concentrations between 0.8 and 11 pmol L−1. The distinct vertical distribution of ferrioxamines and amphibactins may indicate disparate strategies for acquiring Fe from dust in the upper water column and recycled organic matter in the lower water column. Amphibactins were found to have conditional stability constants (log KcondFeL1,Fe′) ranging from 12.0 to 12.5, while ferrioxamines had much stronger conditional stability constants ranging from 14.0 to 14.4, within the range of observed L1 ligands by voltammetry. We used our data to calculate equilibrium Fe speciation at Station ALOHA to compare the relative concentration of inorganic and siderophore complexed Fe. The results indicate that the concentration of Fe bound to siderophores was up to two orders of magnitude higher than inorganic Fe, suggesting that even if less bioavailable, siderophores were nevertheless a viable pathway for Fe acquisition by microbes at our study site. Finally, we observed rapid production of ferrioxamine E by particle-associated bacteria during incubation of freshly collected sinking organic matter. Fe-limitation may therefore be a factor in regulating carbon metabolism and nutrient regeneration in the mesopelagic.

Description: This work was funded by the Woods Hole Oceanographic Postdoctoral Fellowship for RaB, the Simons Foundation (Award 329108), and the National Science Foundation (OCE-1356747).

Description: © The Author(s), 2018. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Frontiers in Marine Science 5 (2018): 170, doi:10.3389/fmars.2018.00170.

Description: Understanding the mechanisms of coral calcification is critical for accurately projecting coral reef futures under ocean acidification and warming. Recent suggestions that calcification is primarily controlled by organic molecules and the biological activity of the coral polyp imply that ocean acidification may not affect skeletal accretion. The basis for these suggestions relies heavily on correlating the presence of organic matter with the orientation and disorder of aragonite crystals in the skeleton, carrying the assumption that organic matter observed in the skeleton was produced by the polyp to control calcification. Here we use Raman spectroscopy to test whether there are differences in organic matter content between coral skeleton and abiogenic aragonites precipitated from seawater, both before and after thermal annealing (heating). We measured the background fluxorescence and intensity of C-H bonding signals in the Raman spectra, which are commonly attributed to coral polyp-derived skeletal organic matrix (SOM) and have been used to map its distribution. Surprisingly, we found no differences in either fluorescence or C-H bonding between abiogenic aragonite and coral skeleton. Annealing reduced the molecular disorder in coral skeleton, potentially due to removal of organic matter, but the same effect was also observed in the abiogenic aragonites. The presence of organic molecules in the abiogenic aragonites is further supported by measurements of N content and δ15N. Together, our data suggest that some of what has been interpreted in previous studies as polyp-derived SOM may actually be seawater-sourced organic matter or some other signal not unique to biogenic aragonite. Finally, we create a high-resolution Raman map of a Pocillopora skeleton to demonstrate how patterns of fluorescence and elevated calcifying fluid aragonite saturation state (ΩAr) along centers of calcification are consistent with both biological and physico-chemical controls. Our aim is to advance discussion on biological mediation of calcification and the implications for coral resilience in a high-CO2 world.

Description: This study was supported by an ARC Laureate Fellowship (FL120100049) awarded to Professor Malcolm McCulloch and the ARC Centre of Excellence for Coral Reef Studies (CE140100020).

Description: © The Author(s), 2018. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Frontiers in Cellular Neuroscience 12 (2018): 156, doi:10.3389/fncel.2018.00156.

Description: Electrical synapses are ubiquitous in interneuron networks. They form intercellular pathways, allowing electrical currents to leak between coupled interneurons. I explored the impact of electrical coupling on the integration of excitatory signals and on the coincidence detection abilities of electrically-coupled cerebellar basket cells (BCs). In order to do so, I quantified the influence of electrical coupling on the rate, the probability and the latency at which BCs generate action potentials when stimulated. The long-lasting simultaneous suprathreshold depolarization of a coupled cell evoked an increase in firing rate and a shortening of action potential latency in a reference basket cell, compared to its depolarization alone. Likewise, the action potential probability of coupled cells was strongly increased when they were simultaneously stimulated with trains of short-duration near-threshold current pulses (mimicking the activation of presynaptic granule cells) at 10 Hz, and to a lesser extent at 50 Hz, an effect that was absent in non-coupled cells. Moreover, action potential probability was increased and action potential latency was shortened in response to synaptic stimulations in mice lacking the protein that forms gap junctions between BCs, connexin36, relative to wild-type (WT) controls. These results suggest that electrical synapses between BCs decrease the probability and increase the latency of stimulus-triggered action potentials, both effects being reverted upon simultaneous excitation of coupled cells. Interestingly, varying the delay at which coupled cells are stimulated revealed that the probability and the speed of action potential generation are facilitated maximally when a basket cell is stimulated shortly after a coupled cell. These findings suggest that electrically-coupled interneurons behave as coincidence and sequence detectors that dynamically regulate the latency and the strength of inhibition onto postsynaptic targets depending on the degree of input synchrony in the coupled interneuron network.

Description: This work was supported by the laboratory of Brain Physiology at Paris Descartes University (UMR8118), the Centre National de la Recherche Scientifique, the Agence Nationale de la Recherche Grant INterneuron NETwork (INNET), the Laboratory of Cellular and Systemic Neurophysiology, Institute for Physiology I at the University of Freiburg, and the Grass foundation.

Description: © The Author(s), 2018. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Environmental Epigenetics 4 (2018): dvy005, doi:10.1093/eep/dvy005.

Description: There is growing evidence that environmental toxicants can affect various physiological processes by altering DNA methylation patterns. However, very little is known about the impact of toxicant-induced DNA methylation changes on gene expression patterns. The objective of this study was to determine the genome-wide changes in DNA methylation concomitant with altered gene expression patterns in response to 3, 3’, 4, 4’, 5-pentachlorobiphenyl (PCB126) exposure. We used PCB126 as a model environmental chemical because the mechanism of action is well-characterized, involving activation of aryl hydrocarbon receptor, a ligand-activated transcription factor. Adult zebrafish were exposed to 10 nM PCB126 for 24 h (water-borne exposure) and brain and liver tissues were sampled at 7 days post-exposure in order to capture both primary and secondary changes in DNA methylation and gene expression. We used enhanced Reduced Representation Bisulfite Sequencing and RNAseq to quantify DNA methylation and gene expression, respectively. Enhanced reduced representation bisulfite sequencing analysis revealed 573 and 481 differentially methylated regions in the liver and brain, respectively. Most of the differentially methylated regions are located more than 10 kilobases upstream of transcriptional start sites of the nearest neighboring genes. Gene Ontology analysis of these genes showed that they belong to diverse physiological pathways including development, metabolic processes and regeneration. RNAseq results revealed differential expression of genes related to xenobiotic metabolism, oxidative stress and energy metabolism in response to polychlorinated biphenyl exposure. There was very little correlation between differentially methylated regions and differentially expressed genes suggesting that the relationship between methylation and gene expression is dynamic and complex, involving multiple layers of regulation.

Description: This work was supported by the National Institute of Health Outstanding New Environmental Scientist Award to NA (NIH R01ES024915) and Woods Hole Center for Oceans and Human Health [National Institutes of Health (NIH) grant P01ES021923 and National Science Foundation Grant OCE-1314642 to M. Hahn, J. Stegeman, NA and SK].

Description: © The Author(s), 2018. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Frontiers in Marine Science 5 (2018): 362, doi:10.3389/fmars.2018.00362.

Description: Major changes to Arctic marine ecosystems have resulted in longer growing seasons with increased phytoplankton production over larger areas. In the Chukchi Sea, the high productivity fuels intense benthic denitrification creating a nitrogen (N) deficit that is transported through the Arctic to the Atlantic Ocean, where it likely fuels N fixation. Given the rapid pace of environmental change and the potentially globally significant N deficit, we conducted experiments aimed at understanding phytoplankton and microbial N utilization in the Chukchi Sea. Ship-board experiments tested the effect of nitrate (NO3-) additions on both phytoplankton and heterotrophic prokaryote abundance, community composition, photophysiology, carbon fixation and NO3- uptake rates. Results support the critical role of NO3- in limiting summer phytoplankton communities to small cells with low production rates. NO3- additions increased particulate concentrations, abundance of large diatoms, and rates of carbon fixation and NO3- uptake by cells 〉1 μm. Increases in the quantum yield and electron turnover rate of photosystem II in +NO3- treatments suggested that phytoplankton in the ambient dissolved N environment were N starved and unable to build new, or repair damaged, reaction centers. While some increases in heterotrophic prokaryote abundance and production were noted with NO3- amendments, phytoplankton competition or grazers likely dampened these responses. Trends toward a warmer more stratified Chukchi Sea will likely enhance summer oligotrophic conditions and further N starve Chukchi Sea phytoplankton communities.

Description: Fieldwork and analysis for the ICESCAPE program was supported by Ocean Biology and Biogeochemistry Program of the National Aeronautic and Space Administration under Grant No. NNX10AF42G to KA.

Description: Author Posting. © The Authors, 2018. This article is posted here by permission of The Royal Astronomical Society for personal use, not for redistribution. The definitive version was published in Geophysical Journal International 215 (2018): 1072–1087, doi:10.1093/gji/ggy203.

Description: An earthquake rupture process can be kinematically described by rupture velocity, duration and spatial extent. These key kinematic source parameters provide important constraints on earthquake physics and rupture dynamics. In particular, core questions in earthquake science can be addressed once these properties of small earthquakes are well resolved. However, these parameters of small earthquakes are poorly understood, often limited by available data sets and methodologies. The Incorporated Research Institutions for Seismology Community Wavefield Experiment in Oklahoma deployed ∼350 three-component nodal stations within 40 km2 for a month, offering an unprecedented opportunity to test new methodologies for resolving small earthquake finite source properties in high resolution. In this study, we demonstrate the power of the nodal data set to resolve the variations in the seismic wavefield over the focal sphere due to the finite source attributes of an M2 earthquake within the array. The dense coverage allows us to tightly constrain rupture area using the second moment method even for such a small earthquake. The M2 earthquake was a strike-slip event and unilaterally propagated towards the surface at 90 per cent local S-wave speed (2.93 km s−1). The earthquake lasted ∼0.019 s and ruptured Lc ∼70 m and Wc ∼45 m. With the resolved rupture area, the stress-drop of the earthquake is estimated as 7.3 MPa for Mw 2.3. We demonstrate that the maximum and minimum bounds on rupture area are within a factor of two, much lower than typical stress-drop uncertainty, despite a suboptimal station distribution. The rupture properties suggest that there is little difference between the M2 Oklahoma earthquake and typical large earthquakes. The new three-component nodal systems have great potential for improving the resolution of studies of earthquake source properties.

Description: WF is currently supported by the Postdoctoral Scholar Program at the Woods Hole Oceanographic Institution, with funding provided by the Weston Howland Jr. Postdoctoral Scholarship. JM was partially supported by SCEC grant #17177 at Woods Hole Oceanographic Institution. This research was supported by the Southern California Earthquake Center (Contribution No. 8014). SCEC is funded by NSF Cooperative Agreement EAR-1033462 and USGS Cooperative Agreement G12AC20038.

Description: Author Posting. © The Authors, 2018. This article is posted here by permission of The Royal Astronomical Society for personal use, not for redistribution. The definitive version was published in Geophysical Journal International 215 (2018): 942–958, doi:10.1093/gji/ggy316.

Description: Surface waves recorded by global arrays have proven useful for locating tectonic earthquakes and in detecting slip events depleted in high frequency, such as glacial quakes. We develop a novel method using an aggregation of small- to continental-scale arrays to detect and locate seismic sources with Rayleigh waves at 20–50 s period. The proposed method is a hybrid approach including first dividing a large aperture aggregate array into Delaunay triangular subarrays for beamforming, and then using the resolved surface wave propagation directions and arrival times from the subarrays as data to formulate an inverse problem to locate the seismic sources and their origin times. The approach harnesses surface wave coherence and maximizes resolution of detections by combining measurements from stations spanning the whole U.S. continent. We tested the method with earthquakes, glacial quakes and landslides. The results show that the method can effectively resolve earthquakes as small as ∼M3 and exotic slip events in Greenland. We find that the resolution of the locations is non-uniform with respect to azimuth, and decays with increasing distance between the source and the array when no calibration events are available. The approach has a few advantages: the method is insensitive to seismic event type, it does not require a velocity model to locate seismic sources, and it is computationally efficient. The method can be adapted to real-time applications and can help in identifying new classes of seismic sources.

Description: WF is currently supported by the Postdoctoral Scholar Program at the Woods Hole Oceanographic Institution, with funding provided by the Weston Howland Jr. Postdoctoral Scholarship. This work was supported by National Science Foundation grant EAR-1358520 at Scripps Institution of Oceanography, UC San Diego.

Description: © The Author(s), 2018. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Frontiers in Microbiology 9 (2018): 560, doi:10.3389/fmicb.2018.00560.

Description: The observation of significant concentrations of soluble Mn(III) complexes in oxic, suboxic, and some anoxic waters has triggered a re-evaluation of the previous Mn paradigm which focused on the cycling between soluble Mn(II) and insoluble Mn(III,IV) species as operationally defined by filtration. Though Mn(II) oxidation in aquatic environments is primarily bacterially-mediated, little is known about the effect of Mn(III)-binding ligands on Mn(II) oxidation nor on the formation and removal of Mn(III). Pseudomonas putida GB-1 is one of the most extensively investigated of all Mn(II) oxidizing bacteria, encoding genes for three Mn oxidases (McoA, MnxG, and MopA). P. putida GB-1 and associated Mn oxidase mutants were tested alongside environmental isolates Pseudomonas hunanensis GSL-007 and Pseudomonas sp. GSL-010 for their ability to both directly oxidize weakly and strongly bound Mn(III), and to form these complexes through the oxidation of Mn(II). Using Mn(III)-citrate (weak complex) and Mn(III)-DFOB (strong complex), it was observed that P. putida GB-1, P. hunanensis GSL-007 and Pseudomonas sp. GSL-010 and mutants expressing only MnxG and McoA were able to directly oxidize both species at varying levels; however, no oxidation was detected in cultures of a P. putida mutant expressing only MopA. During cultivation in the presence of Mn(II) and citrate or DFOB, P. putida GB-1, P. hunanensis GSL-007 and Pseudomonas sp. GSL-010 formed Mn(III) complexes transiently as an intermediate before forming Mn(III/IV) oxides with the overall rates and extents of Mn(III,IV) oxide formation being greater for Mn(III)-citrate than for Mn(III)-DFOB. These data highlight the role of bacteria in the oxidative portion of the Mn cycle and suggest that the oxidation of strong Mn(III) complexes can occur through enzymatic mechanisms involving multicopper oxidases. The results support the observations from field studies and further emphasize the complexity of the geochemical cycling of manganese.

Description: This work was funded by grants from the Chemical Oceanography program of the National Science Foundation (OCE-1558738 and OCE-1558692).

Description: © The Author(s), 2018. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Frontiers in Physiology 9 (2018): 886, doi:10.3389/fphys.2018.00886.

Description: Diving mammals have evolved a suite of physiological adaptations to manage respiratory gases during extended breath-hold dives. To test the hypothesis that offshore bottlenose dolphins have evolved physiological adaptations to improve their ability for extended deep dives and as protection for lung barotrauma, we investigated the lung function and respiratory physiology of four wild common bottlenose dolphins (Tursiops truncatus) near the island of Bermuda. We measured blood hematocrit (Hct, %), resting metabolic rate (RMR, l O2 ⋅ min-1), tidal volume (VT, l), respiratory frequency (fR, breaths ⋅ min-1), respiratory flow (l ⋅ min-1), and dynamic lung compliance (CL, l ⋅ cmH2O-1) in air and in water, and compared measurements with published results from coastal, shallow-diving dolphins. We found that offshore dolphins had greater Hct (56 ± 2%) compared to shallow-diving bottlenose dolphins (range: 30–49%), thus resulting in a greater O2 storage capacity and longer aerobic diving duration. Contrary to our hypothesis, the specific CL (sCL, 0.30 ± 0.12 cmH2O-1) was not different between populations. Neither the mass-specific RMR (3.0 ± 1.7 ml O2 ⋅ min-1 ⋅ kg-1) nor VT (23.0 ± 3.7 ml ⋅ kg-1) were different from coastal ecotype bottlenose dolphins, both in the wild and under managed care, suggesting that deep-diving dolphins do not have metabolic or respiratory adaptations that differ from the shallow-diving ecotypes. The lack of respiratory adaptations for deep diving further support the recently developed hypothesis that gas management in cetaceans is not entirely passive but governed by alteration in the ventilation-perfusion matching, which allows for selective gas exchange to protect against diving related problems such as decompression sickness.

Description: Funding for this project was provided by the Office of Naval Research (ONR YIP Award No. N000141410563, and Dolphin Quest, Inc. FHJ was supported by the Office of Naval Research (Award No. N00014-1410410) and an AIAS-COFUND fellowship from Aarhus Institute of Advanced Studies under the FP7 program of the EU (Agreement No. 609033).

Description: Author Posting. © The Author(s), 2018. This article is posted here by permission of Oxford University Press for personal use, not for redistribution. The definitive version was published in Geophysical Journal International 214 (2018): 2224–2235, doi:10.1093/gji/ggy201.

Description: The key kinematic earthquake source parameters: rupture velocity, duration and area, shed light on earthquake dynamics, provide direct constraints on stress drop, and have implications for seismic hazard. However, for moderate and small earthquakes, these parameters are usually poorly constrained due to limitations of the standard analysis methods. Numerical experiments by Kaneko and Shearer demonstrated that standard spectral fitting techniques can lead to roughly one order of magnitude variation in stress-drop estimates that do not reflect the actual rupture properties even for simple crack models. We utilize these models to explore an alternative approach where we estimate the rupture area directly. For the suite of models, the area averaged static stress drop is nearly constant for models with the same underlying friction law, yet corner-frequency-based stress-drop estimates vary by a factor of 5–10 even for noise-free data. Alternatively, we simulated inversions for the rupture area as parametrized by the second moments of the slip distribution. A natural estimate for the rupture area derived from the second moments is A = πLcWc, where Lc and Wc are the characteristic rupture length and width. This definition yields estimates of stress drop that vary by only 10 per cent between the models but are slightly larger than the true area averaged values. We simulate inversions for the second moments for the various models and find that the area can be estimated well when there are at least 15 available measurements of apparent duration at a variety of take-off angles. The improvement compared to azimuthally averaged corner-frequency-based approaches results from the second moments accounting for directivity and removing the assumption of a circular rupture area, both of which bias the standard approach. We also develop a new method that determines the minimum and maximum values of rupture area that are consistent with a particular data set at the 95 per cent confidence level. For the Kaneko and Shearer models with 20+ randomly distributed observations and ∼10 per cent noise levels, we find that the maximum and minimum bounds on rupture area typically vary by a factor of two and that the minimum stress drop is often more tightly constrained than the maximum.

Description: This work was supported by USGS NEHRP Award G17AP00029. The research was supported by the Southern California Earthquake Center (SCEC; Contribution No. 8013). SCEC is funded by NSF Cooperative Agreement EAR-1033462 and USGS Cooperative Agreement G12AC20038. YK was supported by both public funding from the Government of New Zealand and the Royal Society of New Zealand’s Rutherford Discovery Fellowship.

Description: © The Author(s), 2018. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Frontiers in Marine Science 5 (2018): 273, doi:10.3389/fmars.2018.00273.

Description: Mixotrophic flagellates can comprise significant proportions of plankton biomass in marine ecosystems. Despite the growing recognition of the importance of this ecological strategy, and the identification of major environmental factors controlling phagotrophic behavior (light and nutrients), the physiological and molecular mechanisms underlying mixotrophic behavior are still unclear. In this study, we performed RNA-Seq transcriptomic analysis for two mixotrophic prasinophytes, Micromonas polaris and Pyramimonas tychotreta, under dissolved nutrient regimes that altered their ingestion of bacteria prey. Though the strains examined were polar isolates, both belong to genera with widespread distribution. Our aim was to characterize the transcriptomes of these two non-model phytoflagellates, identify transcripts consistent with phagotrophic activity and assess their differential expression in response to nutrient stress. De novo assembly of the transcriptomes yielded large numbers of novel coding transcripts with no known match within public databases. A summary of the transcripts by Gene Ontology terms showed many expected expression patterns, including genes involved in photosynthetic pathways and enzymes implicated in nutrient uptake pathways. Searches of KEGG databases identified several genes associated with intra-cellular digestive pathways actively transcribed in both prasinophytes. Differential expression analysis showed a larger response in P. tychotreta, where 23,373 genes were up-regulated and 1,752 were down-regulated in the low nutrient treatment when phagotrophy was enhanced. In contrast, in M. polaris, low nutrient treatments resulted in up-regulation of 314 transcripts while down-regulating 371. With respect to phagotrophic-related expression, 37 genes were co-expressed in both P. tychotreta and M. polaris, and although the response was less pronounced in M. polaris, it is consistent with differences in observed ingestion behavior. This study presents the first genomic data for Pyramimonas tychotreta, and also contributes to the limited available data for Micromonas polaris. Furthermore, it provides insight into the presence of genes associated with phagocytosis within the Prasinophyceae and contributes to the understanding of potential target genes required for the construction of a complete model of gene regulation of phagocytic behavior in algae.

Description: The Owlsnest Super-Computing Cluster at Temple University is funded by a National Science Foundation Grant CNS-09-58854. The CUNY HPCC is operated by the College of Staten Island and funded, in part, by grants from the City of New York, State of New York, CUNY Research Foundation, and National Science Foundation Grants CNS-0958379, CNS-0855217, and ACI 1126113. Support for this work was also supplied by National Science Foundation grants PLR-1341362 (RG), PLR-1603538 (RS), and PLR-1603833 (RG).

Description: © The Author(s), 2018. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Frontiers in Environmental Science 6 (2018): 100, doi:10.3389/fenvs.2018.00100.

Description: Determining how microbial communities organize and function at the ecosystem level is essential to understanding and predicting how they will respond to environmental change. Mathematical models can be used to describe these communities, but properly representing all the biological interactions in extremely diverse natural microbial ecosystems in a mathematical model is challenging. We examine a complementary approach based on the maximum entropy production (MEP) principle, which proposes that systems with many degrees of freedom will likely organize to maximize the rate of free energy dissipation. In this study, we develop an MEP model to describe biogeochemistry observed in Siders Pond, a phosphate limited meromictic system located in Falmouth, MA that exhibits steep chemical gradients due to density-driven stratification that supports anaerobic photosynthesis as well as microbial communities that catalyze redox cycles involving O, N, S, Fe, and Mn. The MEP model uses a metabolic network to represent microbial redox reactions, where biomass allocation and reaction rates are determined by solving an optimization problem that maximizes entropy production over time, and a 1D vertical profile constrained by an advection-dispersion-reaction model. We introduce a new approach for modeling phototrophy and explicitly represent oxygenic photoautotrophs, photoheterotrophs and anoxygenic photoautotrophs. The metabolic network also includes reactions for aerobic organoheterotrophic bacteria, sulfate reducing bacteria, sulfide oxidizing bacteria and aerobic and anaerobic grazers. Model results were compared to observations of biogeochemical constituents collected over a 24 h period at 8 depths at a single 15 m deep station in Siders Pond. Maximizing entropy production over long (3 day) intervals produced results more similar to field observations than short (0.25 day) interval optimizations, which support the importance of temporal strategies for maximizing entropy production over time. Furthermore, we found that entropy production must be maximized locally instead of globally where energy potentials are degraded quickly by abiotic processes, such as light absorption by water. This combination of field observations and modeling results indicate that natural microbial systems can be modeled by using the maximum entropy production principle applied over time and space using many fewer parameters than conventional models.

Description: Primary funding for this project was from NSF GG grant EAR-1451356 to JV and JH, with additional support from Gordon and Betty Moore Foundation grant GBMF 3297. JV also received support from NSF Grants OCE-1637630 and OCE-1558710 and Simons Foundation grant 549941. The NSF Center for Dark Energy Biosphere Investigations (C-DEBI; OCE-0939564) also supported the participation of JH.

Description: Author Posting. © The Authors, 2018. This article is posted here by permission of The Royal Astronomical Society for personal use, not for redistribution. The definitive version was published in Geophysical Journal International 215 (2018): 713–735, doi:10.1093/gji/ggy313.

Description: Gas flux in volcanic conduits is often associated with long-period oscillations known as seismic tremor (Lesage et al.; Nadeau et al.). In this study, we revisit and extend the ‘magma wagging’and ‘whirling’models for seismic tremor, in order to explore the effects of gas flux on the motion of a magma column surrounded by a permeable vesicular annulus (Jellinek & Bercovici; Bercovici et al.; Liao et al.). We find that gas flux flowing through the annulus leads to a Bernoulli effect, which causes waves on the magma column to become unstable and grow. Specifically, the Bernoulli effects are associated with torques and forces acting on the magma column, increasing its angular momentum and energy. As the displacement of the magma column becomes large due to the Bernoulli effect, frictional drag on the conduit wall decelerates the motions of the column, restoring them to small amplitude. Together, the Bernoulli effect and the damping effect contribute to a self-sustained wagging-and-whirling mechanism that help explain the longevity of long-period seismic tremor.

Description: This work was supported by National Science Foundation grants EAR-1344538 and EAR-1645057

Description: © The Author(s), 2018. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Frontiers in Earth Science 6 (2018): 88, doi:10.3389/feart.2018.00088.

Description: Shallow seamounts at ocean island hotspots and in other settings may record emergence histories in the form of submarine erosional terraces. Exposure histories are valuable for constraining paleo-elevations and sea levels in the absence of more traditional markers, such as drowned coral reefs. However, similar features can also be produced through primary volcanic processes, which complicate the use of terraced seamounts as an indicator of paleo-shorelines. In the western Galápagos Archipelago, we utilize newly collected bathymetry along with seafloor observations from human-occupied submersibles to document the location and depth of erosional terraces on seamounts near the islands of Santiago, Santa Cruz, Floreana, Isabela, and Fernandina. We directly observed erosional features on 22 seamounts with terraces. We use these observations and bathymetric analysis to develop a framework to identify terrace-like morphologic features and classify them as either erosional or volcanic in origin. From this framework we identify 79 erosional terraces on 30 seamounts that are presently found at depths of 30 to 300 m. Although intermittent subaerial connectivity between the islands has been hypothesized, the depths of these erosional terraces in the Santiago region are the first direct evidence of paleo-connectivity in the modern archipelago. Collectively, the terraces have non-randomly distributed depths. We suggest that peaks in the distribution of terrace depths likely represent long durations of exposure (i.e., sea-level still or lowstands). By comparing these peaks to those of subsidence adjusted sea-level curves, we identify the average subsidence rate that best reproduces the observed terrace distributions. These rates are 0.2–0.4 m/ka for this portion of the central Galápagos, since the formation of the seamounts, consistent with previous independent estimates. Using these subsidence rates and evidence for erosional terraces at depths up to 300 m, we conclude that all islands in the central archipelago have been intermittently connected starting between 435 and 900 ka. Individual island pairs have likely been repeatedly subaerially connected for short intervals since that time.

Description: This project was carried out with financial support from the NSF (OCE-1634685 to SS and OCE-1634952 to VW) and the Dalio Explore Fund.

Description: © The Author(s), 2018. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Frontiers in Molecular Neuroscience 11 (2018): 427, doi:10.3389/fnmol.2018.00427.

Description: Electrical signaling is a cardinal feature of the nervous system and endows it with the capability of quickly reacting to changes in the environment. Although synaptic communication between nerve cells is perceived to be mainly chemically mediated, electrical synaptic interactions also occur. Two different strategies are responsible for electrical communication between neurons. One is the consequence of low resistance intercellular pathways, called “gap junctions”, for the spread of electrical currents between the interior of two cells. The second occurs in the absence of cell-to-cell contacts and is a consequence of the extracellular electrical fields generated by the electrical activity of neurons. Here, we place present notions about electrical transmission in a historical perspective and contrast the contributions of the two different forms of electrical communication to brain function.

Description: This research was supported by National Institutes of Health grants DC03186, DC011099, NS055726, NS085772 and NS0552827 to AP.

Description: Author Posting. © The Authors, 2018. This article is posted here by permission of The Royal Astronomical Society for personal use, not for redistribution. The definitive version was published in Geophysical Journal International 215 (2018): 460–473, doi:10.1093/gji/ggy152.

Description: In this work, we present a new methodology to predict grain-size distributions from geophysical data. Specifically, electric conductivity and magnetic susceptibility of seafloor sediments recovered from electromagnetic profiling data are used to predict grain-size distributions along shelf-wide survey lines. Field data from the NW Iberian shelf are investigated and reveal a strong relation between the electromagnetic properties and grain-size distribution. The here presented workflow combines unsupervised and supervised machine-learning techniques. Non-negative matrix factorization is used to determine grain-size end-members from sediment surface samples. Four end-members were found, which well represent the variety of sediments in the study area. A radial basis function network modified for prediction of compositional data is then used to estimate the abundances of these end-members from the electromagnetic properties. The end-members together with their predicted abundances are finally back transformed to grain-size distributions. A minimum spatial variation constraint is implemented in the training of the network to avoid overfitting and to respect the spatial distribution of sediment patterns. The predicted models are tested via leave-one-out cross-validation revealing high prediction accuracy with coefficients of determination (R2) between 0.76 and 0.89. The predicted grain-size distributions represent the well-known sediment facies and patterns on the NW Iberian shelf and provide new insights into their distribution, transition and dynamics. This study suggests that electromagnetic benthic profiling in combination with machine learning techniques is a powerful tool to estimate grain-size distribution of marine sediments.

Description: This work was funded through DFG Research Center/Cluster of Excellence ‘The Ocean in the Earth System’ and was part of MARUM Research Area SD

Description: © The Author(s), 2018. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Frontiers in Microbiology 9 (2018): 189, doi:10.3389/fmicb.2018.00189.

Description: Only select prokaryotes can biosynthesize vitamin B12 (i.e., cobalamins), but these organic co-enzymes are required by all microbial life and can be vanishingly scarce across extensive ocean biomes. Although global ocean genome data suggest cyanobacteria to be a major euphotic source of cobalamins, recent studies have highlighted that 〉95% of cyanobacteria can only produce a cobalamin analog, pseudo-B12, due to the absence of the BluB protein that synthesizes the α ligand 5,6-dimethylbenzimidizole (DMB) required to biosynthesize cobalamins. Pseudo-B12 is substantially less bioavailable to eukaryotic algae, as only certain taxa can intracellularly remodel it to one of the cobalamins. Here we present phylogenetic, metagenomic, transcriptomic, proteomic, and chemical analyses providing multiple lines of evidence that the nitrogen-fixing cyanobacterium Trichodesmium transcribes and translates the biosynthetic, cobalamin-requiring BluB enzyme. Phylogenetic evidence suggests that the Trichodesmium DMB biosynthesis gene, bluB, is of ancient origin, which could have aided in its ecological differentiation from other nitrogen-fixing cyanobacteria. Additionally, orthologue analyses reveal two genes encoding iron-dependent B12 biosynthetic enzymes (cbiX and isiB), suggesting that iron availability may be linked not only to new nitrogen supplies from nitrogen fixation, but also to B12 inputs by Trichodesmium. These analyses suggest that Trichodesmium contains the genus-wide genomic potential for a previously unrecognized role as a source of cobalamins, which may prove to considerably impact marine biogeochemical cycles.

Description: This work was funded by NSF research grants OCE-1260233, OCE-1260490, OCE-1657757, and OCE-143566.

Description: © The Author(s), 2018. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Frontiers in Microbiology 9 (2018): 358, doi:10.3389/fmicb.2018.00358.

Description: Viruses are key players in ocean ecology and biogeochemistry, not only because of their functional roles but also partially due to their sheer abundance (Fuhrman, 1999; Wilhelm and Suttle, 1999). Because viruses cannot replicate without their hosts' machinery, their abundance is inextricably related to that of their (mostly microbial) hosts. The relationship between viral and microbial abundances is thus of great interest.

Description: Author Posting. © Author(s), 2017. This article is posted here by permission of Oxford University Press for personal use, not for redistribution. The definitive version was published in Geophysical Journal International 212 (2018): 1429–1449, doi:10.1093/gji/ggx488.

Description: We conducted detailed analyses of a global array of trenches, revealing systematic intra- and intertrench variations in plate bending characteristics. The intratrench variations of the Manila and Mariana Trenches were analysed in detail as end-member cases of the relatively young (16–36 Ma) and old (140–160 Ma) subducting plates, respectively. Meanwhile, the intertrench variability was investigated for a global array of additional trenches including the Philippine, Kuril, Japan, Izu-Bonin, Aleutian, Tonga-Kermadec, Middle America, Peru, Chile, Sumatra and Java Trenches. Results of the analysis show that the trench relief (W0) and width (X0) of all systems are controlled primarily by the faulting-reduced elastic thickness near the trench axis (Tme) and affected only slightly by the initial unfaulted thickness (TMe) of the incoming plate. The reduction in Te has caused significant deepening and narrowing of trench valleys. For the cases of relatively young or old plates, the plate age could be a dominant factor in controlling the trench bending shape, regardless the variations in axial loadings. Our calculations also show that the axial loading and stresses of old subducting plates can vary significantly along the trench axis. In contrast, the young subducting plates show much smaller values and variations in axial loading and stresses.

Description: This work was supported by Chinese Academy of Sciences Grants (Y4SL021001, QYZDY-SSW-DQC005, YZ201325 and YZ201534), National Natural Science Foundation of China Grants (91628301, U1606401, 41376063 and 41706056) and HKSAR Research Grant Council Grants (24601515, 14313816).

Description: © The Author(s), 2018. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Frontiers in Microbiology 8 (2018): 2669, doi:10.3389/fmicb.2017.02669.

Description: Molecular methods for microbial community characterization have uncovered environmental and plant-associated factors shaping phyllosphere communities. Variables undetectable using bulk methods can play an important role in shaping plant-microbe interactions. Microscale analysis of bacterial dynamics in the phyllosphere requires imaging techniques specially adapted to the high autoflouresence and 3-D structure of the leaf surface. We present an easily-transferable method (Leaf-FISH) to generate high-resolution tridimensional images of leaf surfaces that allows simultaneous visualization of multiple bacterial taxa in a structurally informed context, using taxon-specific fluorescently labeled oligonucleotide probes. Using a combination of leaf pretreatments coupled with spectral imaging confocal microscopy, we demonstrate the successful imaging bacterial taxa at the genus level on cuticular and subcuticular leaf areas. Our results confirm that different bacterial species, including closely related isolates, colonize distinct microhabitats in the leaf. We demonstrate that highly related Methylobacterium species have distinct colonization patterns that could not be predicted by shared physiological traits, such as carbon source requirements or phytohormone production. High-resolution characterization of microbial colonization patterns is critical for an accurate understanding of microbe-microbe and microbe-plant interactions, and for the development of foliar bacteria as plant-protective agents.

Description: Funding was provided by the J. Unger Vetleson Foundation to SS.

Description: 〈span〉〈div〉Summary〈/div〉The strength of the lithosphere plays a key role in the formation and evolution of tectonic plate boundaries. Localized lithospheric deformation associated with plate tectonics requires a mechanism for weakening across the entire width of the lithosphere, including the strongest cold ductile region. We explore the microphysics of weakening of lithospheric materials, and in particular the coupled evolution of mineral grain size and intragranular defects and their control on lithospheric strength. We propose a model for the interaction between grain-boundaries and dislocation density to reduce the net free energy of grains during dynamic recrystallization (DRX). The driving forces for DRX arise from heterogeneity in dislocation density and grain boundary curvature. Our model shows that grain growth driven by variation in grain boundary curvature can be impeded by variation in dislocation density; this occurs because as the grains grow, to minimize their surface energy, their dislocation density and associated internal energy may increase and offset the driving forces for grain growth. The correlation between grain size and dislocation density can for example arise because the dislocation accumulation in smaller grains is suppressed due to the large stress that is needed to bend and elongate a short dislocation (as dictated by the small grain size), while the larger grains can have long dislocations and reach a steady state dislocation density dictated by the applied stress. In a lithospheric setting, slower grain growth means that it would require less mechanical work to establish weak localized shear zones through grain damage, and retard the healing of previously damaged zones. Furthermore, the competition of two different time-scales - that of grain growth and the dislocation kinetics - can lead to oscillating behavior over 1 to 10 years as the grain size and dislocation density advance towards their steady states. These oscillations are likely to have an effect on the rheology of lithospheric rocks, e.g. their strengthening and weakening through time, and have a potential application to geological processes such as postseismic creep in ductile shear zones.〈/span〉

Description: 〈span〉〈div〉Summary〈/div〉Locating and monitoring passive seismic sources provides us important information for studying subsurface rock deformation, injected fluid migration, regional stress conditions as well as fault rupture mechanism. In this paper, we present a novel passive-source monitoring approach using vector-based elastic time-reversal imaging. By solving the elastic wave equation using observed multicomponent records as boundary conditions, we first compute back-propagated elastic wavefields in the subsurface. Then, we separate the extrapolated wavefields into compressional (P-wave) and shear (S-wave) modes using the vector Helmholtz decomposition. A zero-lag cross-correlation imaging condition is applied to the separated pure-mode vector wavefields to produce passive-source images. We compare imaging results using three implementations, i.e., dot-product, energy and power. Numerical experiments demonstrate that the power imaging condition gives us the highest resolution and is less sensitive to the presence of random noises. To capture the propagation of microseismic fracture and earthquake rupture, we modify the traditional zero-lag cross-correlation imaging condition by summing the multiplication of the separated P- and S-wavefields within local time windows, which enables us to capture the temporal and spatial evolution of earthquake rupture. 2D and 3D numerical examples demonstrate that the proposed method is capable of accurately locating point sources, as well as delineating dynamic propagation of hydraulic fracture and earthquake rupture.〈/span〉

Description: 〈span〉〈div〉Summary〈/div〉Magnitudes of differential stress in the lithosphere, especially in the crust, are still disputed. Earthquake-based stress drop estimates indicate median values 〈 10 MPa whereas the lateral variation of gravitational potential energy per unit area (〈span〉GPE〈/span〉) across significant relief indicates stress magnitudes of ca. 100 MPa in average across a 100 km thick lithosphere between the Indian lowland and the Tibetan plateau. These standard 〈span〉GPE〈/span〉-based stress estimates correspond to membrane stresses, because they are associated with a deformation that is uniform with depth. We show here with new analytical results that lateral variations in 〈span〉GPE〈/span〉 can also cause bending moments and related bending stresses of several hundreds of MPa. Furthermore, we perform two-dimensional thermo-mechanical numerical simulations (1) to evaluate estimates for membrane and bending stresses based on 〈span〉GPE〈/span〉 variations, (2) to quantify minimum crustal stress magnitudes that are required to maintain the topographic relief between Indian lowland and Tibetan plateau for ca. 10 Ma and (3) to quantify the corresponding relative contribution of crustal strength to the total lithospheric strength. The numerical model includes viscoelastoplastic deformation, gravity and heat transfer. The model configuration is based on density fields from the CRUST1.0 data set and from a geophysically and petrologically constrained density model based on 〈span〉in situ〈/span〉 field campaigns. The numerical results indicate that values of differential stress in the upper crust must be 〉 ca. 180 MPa, corresponding to a friction angle of ca. 10°, to maintain the topographic relief between lowland and plateau for 〉 10 Ma. The relative contribution of crustal strength to total lithospheric strength varies considerably laterally. In the region between lowland and plateau and inside the plateau the depth-integrated crustal strength is approximately equal to the depth-integrated strength of the mantle lithosphere. Simple analytical formulae predicting the lateral variation of depth-integrated stresses agree with numerically calculated stress fields, which show both the accuracy of the numerical results and the applicability of simple, rheology-independent, analytical predictions to highly variable, rheology-dependent, stress fields. Our results indicate that (1) crustal strength can be locally equal to mantle lithosphere strength and that (2) crustal stresses must be at least one order of magnitude larger than median stress drops in order to support the plateau relief over a duration of ca. 10 Ma.〈/span〉

Description: 〈span〉〈div〉Summary〈/div〉Scanning magnetometers are increasingly used to characterize the magnetization of mineral grains in rock samples. Up-scaling this measurement technique to large numbers of individual particles is hampered by the intrinsic non-uniqueness of potential-field inversion. Here it is shown that this problem can be circumvented by adding tomographic information that determines the location of the possible field sources. Standard potential theory is used to prove a uniqueness theorem which completely characterizes the mathematical background of the corresponding source-localized inversion. It exactly resolves under which conditions a potential field measurement on a surface can be uniquely decomposed into signals from the different source regions. The intrinsic non-uniqueness of potential field inversion prevents that the source distribution inside the tomographically outlined regions can be recovered, but the potential field of each region is uniquely defined. For scanning magnetometers in rock magnetism, this result implies that magnetic dipole vectors of large numbers of individual magnetic particles can be reliably reconstructed from surface scans of the magnetic field, if the particle positions are independently determined. This provides an incentive to improve scanning methods for future paleomagnetic applications.〈/span〉

Description: 〈span〉〈div〉Summary〈/div〉The most common earthquake forecasting models assume that the magnitude of the next earthquake is independent from the past. This feature severely limits the capability to forecast large earthquakes with high probabilities. Here we investigate empirically on the magnitude-independence assumption, exploring if: i) background and triggered earthquakes have the same frequency-magnitude distribution, ii) variations of seismicity in the space-time-magnitude domain encode some information on the future earthquakes size. For this purpose, and to verify the stability of the findings, we consider seismic catalogues covering different space-time-magnitude windows, such as the Alto Tiberina Near Fault Observatory (TABOO), the California and Japanese seismic catalogues. Our approach is inspired by the nearest-neighbour method proposed by Baiesi & Paczuski (2004) and elaborated by Zaliapin et al. (2008) to distinguish between triggered and background earthquakes. Here we implement the same metric-based correlation to identify the precursory seismicity of any triggered earthquake; this allows us to analyse, for each triggered earthquake, the space-time-magnitude distribution of the seismicity that likely contributed to its occurrence. Our results show that the magnitude-independence assumption holds reasonably well in all catalogues, with a remarkable exception that is consistent with a previous independent study; this departure from the magnitude-independence assumption shows that larger events tend to nucleate at a higher distance from the ongoing sequence. We also notice that the reliability of this assumption may depend on the spatial scale considered; it holds for seismic catalogues of large areas, but we identify possible departures in small areas, reflecting different ways to release locally seismic energy. Finally, we come across an important issue that may lead to misleading results in similar studies, i.e., if a seismic catalogue appears overall complete above a fixed magnitude threshold, it may still yield spurious signals into the analysis. Specifically, we show that some significant departures from the magnitude-independence assumption do not survive when considering spatiotemporal variations of the magnitude of completeness.〈/span〉

Description: 〈span〉〈div〉Summary〈/div〉Vp/Vs models provide important complementary information to Vp and Vs models, relevant to lithology, rock damage, partial melting, water saturation, etc. However, seismic tomography using body-wave traveltime data from local or regional earthquakes does not constrain Vp/Vs well due to the different resolution of Vp and Vs models, with the Vp models usually better constrained than Vs. Since surface-wave data are most sensitive to Vs, which leads to complementary strengths with respect to body-wave data, we jointly invert body- and surface-wave data to better resolve the Vp/Vs models. In order to show the robustness of our joint inversion method, we compare the results to other approaches, including dividing Vp by Vs models and Vp/Vs parameterization with body-wave or both body- and surface-wave data, using synthetic data and real data from the southern California plate boundary region. We confirm that Vp/Vs models from body-wave inversion obtained by division tend to include artifacts, even when both Vp and Vs models seem reasonable. The joint inversion with Vp/Vs parameterization is found to improve the Vp/Vs ratio model significantly and the resultant Vp/Vs model shows more geologically consistent features, such as high Vp/Vs along fault traces at shallow depths likely indicating fault-related damage. The Vp/Vs model also exhibits contrasts at intermediate depths along the Peninsular Range compositional boundary, and high Vp/Vs in the lower crust near the Salton Sea region correlated with high heat flow and may indicate partial melting. The improved Vp/Vs as well as individual Vp and Vs models are useful for earthquake relocation, high-resolution Moho depth imaging, and interpretation of other data and tectonic evolution in the region.〈/span〉

Description: 〈span〉〈div〉Summary〈/div〉Interpretation of surface fault scarps and palaeoseismic trenches is a key component of estimating fault slip rates, earthquake recurrence rates and maximum magnitudes for hazard assessments. Often these analyses rely on the assumption that successive earthquakes all breached the surface and that the ruptures are recorded topographically, or by the deposits exposed in a trench. The M〈sub〉〈span〉w〈/span〉〈/sub〉7.2 1992 Suusamyr earthquake, Kyrgyzstan, is an apparently problematic case for such analyses because its ruptures show significant displacement but are only mapped as having broken the surface along small, disparate portions of the fault. Here we present the results of surveys conducted along the Suusamyr Fault to establish whether that is the case. Two sets of ruptures were identified following the earthquake. They are unusually short for their displacement and are separated by a 25 km gap. Using satellite imagery, high-resolution digital elevation models and palaeoseismic trenching we first reassess the distribution of the 1992 ruptures and then reconstruct the Holocene earthquake record to establish the extent to which the 1992 earthquake is representative of the rupture behaviour of this fault. We find evidence for at least two prehistoric surface rupturing earthquakes in the Holocene: one ∼3 ka and one 〉8 ka that, along with the modern event, gives recurrence intervals of ∼3 kyr and ∼5 kyr. Within spatial gaps between segments of the 1992 ruptures there are clear prehistoric surface ruptures and the ruptures in each prehistoric earthquake were discontinuous. We conclude that there is significant variability in the surface rupture pattern of successive earthquakes on the Suusamyr Fault, with implications for the completeness of palaeoseismic records obtained from thrust scarps.〈/span〉

Description: 〈span〉〈div〉Summary〈/div〉We provide a six-component (6-C) polarization model for 〈span〉P〈/span〉-, 〈span〉SV〈/span〉-, 〈span〉SH〈/span〉-, Rayleigh-, and Love-waves both inside an elastic medium as well as at the free surface. It is shown that single-station 6-C data comprised of three components of rotational motion and three components of translational motion provide the opportunity to unambiguously identify the wave type, propagation direction, and local 〈span〉P〈/span〉- and 〈span〉S〈/span〉-wave velocities at the receiver location by use of polarization analysis. To extract such information by conventional processing of three-component (3-C) translational data would require large and dense receiver arrays. The additional rotational components allow the extension of the rank of the coherency matrix used for polarization analysis. This enables us to accurately determine the wave type and wave parameters (propagation direction and velocity) of seismic phases, even if more than one wave is present in the analysis time window. This is not possible with standard, pure-translational 3-C recordings. In order to identify modes of vibration and to extract the accompanying wave parameters, we adapt the multiple signal classification algorithm (MUSIC). Due to the strong nonlinearity of the MUSIC estimator function, it can be used to detect the presence of specific wave types within the analysis time window at very high resolution. We show how the extracted wavefield properties can be used, in a fully automated way, to separate the wavefield into its different wave modes using only a single 6-C recording station. As an example, we apply the method to remove surface wave energy while preserving the underlying reflection signal and to suppress energy originating from undesired directions, such as side-scattered waves.〈/span〉

Description: 〈span〉〈div〉Summary〈/div〉The Sentinel-1 mission comprises two synthetic aperture radar satellites, each with a 12 day orbital repeat, orbiting 6 days apart within a narrow tube. The mission design promises the ability to respond quickly to earthquakes with InSAR, and to facilitate production of interferograms with good interferometric correlation globally. We report on our efforts to study global seismicity using Sentinel-1 Interferometric Wide-Swath data between April 2015 and December 2016. We select 35 potentially detectable terrestrial earthquakes in the range 5.5 ≤ 〈span〉Mw〈/span〉 ≤ 7.8 on the basis of their locations, depths and magnitudes, and process the first post-event interferogram with the shortest possible time-span for each using the ISCE software. We evaluate each interferogram for earthquake deformation signals by visual inspection. We can identify deformation signals attributable to earthquakes in 18 of these interferograms (51%); a further six interferograms (17%) have ambiguous interferometric phase affected by tropospheric noise. 11 events (31%) could not be identified from their interferograms. The majority of these failed detections were due to interferogram decorrelation, particularly apparent for earthquakes that occurred between 15°N and 15°S, where climate conditions promote dense vegetation. The majority of the ambiguous interferograms are affected by tropospheric noise, suggesting that techniques to mitigate such noise could improve detection performance. The largest event we do not detect with Sentinel-1 data is a 〈span〉Mw〈/span〉7.0 earthquake that occurred in Vanuatu in April 2016; we also fail to detect the 2016 〈span〉Mw〈/span〉6.2 Kurayoshi earthquake in one out of two possible 24-day interferograms. We propose these as upper and lower estimates on the magnitude of completeness for earthquakes studied with Sentinel-1 data; to raise the magnitude of completeness we suggest that more frequent (e.g. six day) recurrence may be necessary in low latitude areas.〈/span〉

Description: 〈span〉〈div〉Summary〈/div〉The correct estimation of site-specific attenuation is crucial for the assessment of seismic hazard. Downhole instruments provide in this context valuable information to constrain attenuation directly from data. In this study, we apply an interferometric approach to this problem by deconvolving seismic motions recorded at depth with those recorded at the surface. In doing so, incident and surface-reflected waves can be separated. We apply this technique not only to earthquake data but also to recordings of ambient vibrations. We compute the transfer function between incident and surface-reflected waves in order to infer frequency dependent quality factors for S-waves. The method is applied to a 87 m deep borehole sensor and a co-located surface instrument situated at a hard-rock site in West Bohemia/Vogtland, Germany. We show that the described method provides comparable attenuation estimates using either earthquake data or ambient noise for frequencies between 5-15 Hz. Moreover, a single hour of noise recordings seems to be sufficient to yield stable deconvolution traces and quality factors, thus, offering a fast and easy way to derive attenuation estimates from borehole recordings even in low to mid seismicity regions.〈/span〉

Description: 〈span〉〈div〉Summary〈/div〉Injection of CO〈sub〉2〈/sub〉 into tight reservoirs produces both gravity change and ground deformation, which provide great opportunities for more accurate coupled inverse modelling. In this study, we incorporate signals generated from several synthetic models to estimate the CO〈sub〉2〈/sub〉 distribution in the reservoir. A relationship is found that connects density variations to volumetric changes associated with injected CO〈sub〉2〈/sub〉, taking advantage of a common set of model parameters for both gravitational and geo-mechanical inverse modelling. This is achieved by assuming that the injected CO〈sub〉2〈/sub〉 increases pressure in the reservoir, which in turn generates extra porosity that is then filled in by the CO〈sub〉2〈/sub〉 mass in the generated space. Tikhonov regularization, supported by the Generalized Cross Validation (GCV) technique for finding the optimized model, is used to solve the ill-posed inverse problems. The results indicate that with a combination of gravity and ground deformation monitoring, the uncertainty and ambiguity in gravimetric modelling due to high levels of noise is mitigated by implementing highly accurate ground deformation measurements, which normally have a higher signal to noise ratio.〈/span〉

Description: 〈span〉〈div〉Summary〈/div〉As the number of seismic sensors grows, it is becoming increasingly difficult for analysts to pick seismic phases manually and comprehensively, yet such efforts are fundamental to earthquake monitoring. Despite years of improvements in automatic phase picking, it is difficult to match the performance of experienced analysts. A more subtle issue is that different seismic analysts may pick phases differently, which can introduce bias into earthquake locations. We present a deep-neural-network-based arrival-time picking method called ”PhaseNet” that picks the arrival times of both P and S waves. Deep neural networks have recently made rapid progress in feature learning, and with sufficient training, have achieved super-human performance in many applications. PhaseNet uses three-component seismic waveforms as input and generates probability distributions of P arrivals, S arrivals, and noise as output. We engineer PhaseNet such that peaks in the probability distributions provide accurate arrival times for both P and S waves. PhaseNet is trained on the prodigious available data set provided by analyst-labeled P and S arrival times from the Northern California Earthquake Data Center. The dataset we use contains more than seven hundred thousand waveform samples extracted from over thirty years of earthquake recordings. We demonstrate that PhaseNet achieves much higher picking accuracy and recall rate than existing methods when applied to the waveforms of known earthquakes, which has the potential to increase the number of S-wave observations dramatically over what is currently available. This will enable both improved locations and improved shear wave velocity models.〈/span〉

Description: 〈span〉〈div〉Summary〈/div〉The low frequency earthquakes (LFEs) that constitute tectonic tremor are often inferred to be slow: to have durations of 0.2 to 0.5 s, a factor of 10 to 100 longer than those of typical 〈span〉MW〈/span〉 1-2 earthquakes. Here we examine LFEs near Parkfield, CA in order to assess several proposed explanations for LFEs’ long durations. We determine LFE rupture areas and location distributions using a new approach, similar to directivity analysis, where we examine how signals coming from various locations within LFEs’ finite rupture extents create differences in the apparent source time functions recorded at various stations. We use synthetic ruptures to determine how much the LFE signals recorded at each station would be modified by spatial variations of the source-station travel time within the rupture area given various possible rupture diameters, and then compare those synthetics with the data. Our synthetics show that the methodology can identify inter-station variations created by heterogeneous slip distributions or complex rupture edges, and thus lets us estimate LFE rupture extents for unilateral or bilateral ruptures. To obtain robust estimates of the sources’ similarity across stations, we stack signals from thousands of LFEs, using an empirical Green’s function approach to isolate the LFEs’ apparent source time functions from the path effects. Our analysis of LFEs in Parkfield implies that LFEs’ apparent source time functions are similar across stations at frequencies up to 8 to 16 Hz, depending on the family. The inter-station coherence observed at these relatively high frequencies, or short wavelengths (down to 0.2 to 0.5 km), suggest that LFEs in each of the 7 families examined occur on asperities. They are clustered in patches with sub-1-km diameters. The individual LFEs’ rupture diameters are estimated to be smaller than 1.1 km for all families, and smaller than 0.5 km and 1 km for the two shallowest families, which were previously found to have 0.2-s durations. Coupling the diameters with the durations suggests that it is possible to model these 〈span〉MW〈/span〉 1-2 LFEs with earthquake-like rupture speeds: around 70% of the shear wave speed. However, that rupture speed matches the data only at the edge of our uncertainty estimates for the family with highest coherence. The data for that family are better matched if LFEs have rupture velocities smaller than 40% of the shear wave speed, or if LFEs have different rupture dynamics. They could have long rise times, contain composite sub-ruptures, or have slip distributions that persist from event to event.〈/span〉

Description: 〈span〉〈div〉Summary〈/div〉On June 24, 2015, a 230,000 cubic metre landslide slid into the triangle bayou at the intersection of the Yangtze and Daning Rivers and generated a river tsunami that ran up 6.2 metres on the opposite shoreline at Wushan town. The slope failure and resulting waves killed two people and damaged many shipping facilities. Based on field surveys and eyewitness observations, we apply the ‘Tsunami Squares’ method to model the Hongyanzi landslide and its generated waves. Landslide simulations indicate a maximum impact velocity of ∼16 m/s that matches well with an eyewitness video. The computed post-slide mass stopped on the near riverbed with a shape fitting the observed geological profile. Tsunami simulations reveal a large region of wave impacts that coincide with the observed runup heights. The successful reproduction of the dynamics of this landslide-generated river tsunami emphasizes the capacity and efficiency of Tsunami Squares modeling in emergency reaction and risk assessment.〈/span〉

Description: 〈span〉〈div〉Summary〈/div〉Improvement of global 3D Earth density and velocity models is based in part on measurements of Earth’s normal mode eigenfrequencies and splitting function coefficients. Despite many methods developed inconsistency in measurements still exists and it is difficult to understand which results are more precise, that is which methods introduce less systematic biases in the measurements. Therefore, the main goal of this study is to test the performances of typically used techniques in low-frequency normal mode studies: the optimal sequence estimation stacking method and the autoregressive method in the frequency domain, where validation of the estimates is performed with the phasor walkout method. Motivations for their utilization are their easy and fast implementation and their accurate performances when it comes to eigenfrequency estimates. For this purpose, we first perform the analysis with synthetic seismograms in order to evaluate how the station distributions and noise levels impact the estimates of eigenfrequencies and structure coefficients. Synthetic seismograms are calculated for a 3D realistic Earth model, which includes Earth’s rotation as well as ellipticity and other lateral heterogeneities. They were computed by means of normal mode summation and a perturbation theory for modes up to 1 mHz. The three methods above are also applied to long-period seismometer and superconducting gravimeter data recorded after six earthquakes of magnitude greater than 8.3. Finally, our study shows that the optimal sequence estimation is sensitive to the station distribution under the noise influence, while the autoregressive method for frequency estimation gives us reasonable estimates within the estimated error bars. Moreover, we present new estimates of eigenfrequencies and Q-factors for 〈sub〉0〈/sub〉S〈sub〉2, 0〈/sub〉S〈sub〉3, 2〈/sub〉S〈sub〉1〈/sub〉 and 〈sub〉3〈/sub〉S〈sub〉1〈/sub〉 multiplets. A new value for the 〈span〉c〈/span〉〈sub〉20〈/sub〉 structure coefficient of 〈sub〉0〈/sub〉S〈sub〉2〈/sub〉 multiplet −0.7233 ± 0.0623 μHz is obtained.〈/span〉

Description: 〈span〉〈div〉Summary〈/div〉We propose a novel approach to compute the gravity field due to density anomaly in both outside and inside of the solid Earth with high accuracy and efficiency. The high accuracy comes from the direct employment of the analytic gravitation solution between any point on a two-dimensional (2D) plane in the horizontal direction and individual mass cubes. The high computational efficiency comes from two aspects: 1, the application of the highly efficient 2D discrete convolution algorithm; and 2, a newly developed algorithm for the optimized computation of the weight coefficient matrix. Numerical examples for applying to compare with analytical solutions demonstrated its excellent accuracy. Comparison with other state-of-the-art gravity modeling algorithm has proved that this algorithm has superior performance in both accuracy and efficiency. Application to analyze real topography demonstrated the practicality. This algorithm will be an attractive candidate for carrying out the forward modeling step in geophysical inversion problems with the claimed and proved advantages.〈/span〉

Description: 〈span〉〈div〉Summary〈/div〉We assess the stress field that drives lithospheric deformation, focusing on the active New Zealand plate boundary zone between the Australian and Pacific plates. Here there is a rich database for the horizontal velocity field and crustal structure, used to derive gravitational potential energy (GPE). We solve the stress balance equations, in the context of a thin sheet model of a viscously deforming lithosphere, for characteristic deviatoric stresses and viscosities, defined as the integral of these with depth (divided by the layer thickness), using the stress method of Flesch 〈span〉et al.〈/span〉 (2001), where the input parameters are the fields of strain rate and GPE that apply to the sheet. Synthetic tests show that the stress method is able to resolve the stress field to high (20 per cent) levels of noise in the input strain rates, and the mean stress is a very robust feature of the inversions, regardless of noise levels. We invert for the stress and viscosity fields in New Zealand, calculating the field of GPE from topography/bathymetry and crustal data, and the field of strain rates from either a long term (multi millennial) velocity field inferred from the rate and pattern of Quaternary faulting, or a short term (decadal) velocity field directly observed with decadal Global Positioning System (GPS) measurements. In addition, we consider the effect of shear stresses on the subducted plate interface along the Hikurangi Margin (5–15 MPa), or regionally. We explore the effect of GPE on the inversion results by calculating these for a range of deforming layers (ie thin sheet with 35–150 km thickness), in effect sampling the lithospheric strength in the crust and mantle. The results show that the derived stress magnitudes (square root second invariant of the stress deviator) are in the range 0–35 MPa, with mean values of 13 ± 1 MPa for all models, comparable to typical earthquake stress drops. Gravitationally induced stresses account for approximately half of the full deviatoric stress. Effective characteristic viscosities are 0.5 - 5 × 10〈sup〉21〈/sup〉 Pas in the deforming zone, with an approximate inverse relation between strain rate and viscosity, most likely controlled by thermal structure and/or lithology.〈/span〉

Description: 〈span〉〈div〉Summary〈/div〉I present a source-independent fracture imaging method to use passive seismic data for mapping subwavelength natural fractures. Unlike conventional source-dependent imaging that often adopts reflection-type seismic imaging with known source that is not available in passive seismic surveys, the proposed fracture imaging approach relies on the transmission and diffraction data without the need for source information. I assume that passive seismic data can be decomposed into two types of data: primary transmission wave data and diffraction (coda) wave data. The imaging formula states that primary waves should coincide with coda waves at scatterer points at the time of scattering. Instead of generating source wavefields in the conventional imaging method, the proposed method only need to propagate transmission wave data and diffraction wave data from the receiver arrays and apply an imaging condition to produce an image of fractures. This imaging procedure can be used for processing P-wave or S-wave. In synthetic examples, I evaluate the proposed method in several aspects: inaccurate source location, inaccurate velocity model, sparse receivers and irregular receiver spacing, elastic data, and joint surface and borehole acquisitions. I found that the proposed approach performed well (or even better) comparable to source-dependent fracture imaging when assuming exact source information is known. With perturbed source locations with random shifts (e.g., estimated source location with errors), however, fractures were missing in the source-dependent fracture imaging results but the proposed approach was not influenced. In the presence of velocity errors and sparse and irregular receiver spacing, the proposed method produces better fracture images than the source-dependent imaging results.〈/span〉

Description: 〈span〉〈div〉SUMMARY〈/div〉We introduce a formalism for estimating local spatial averages of the core-mantle boundary (CMB) radial magnetic field and its time derivatives, based on magnetic field observations collected by low-Earth-orbit satellites. This provides a useful alternative to conventional core field modelling based on global spherical harmonic basis functions, where noise in the polar regions maps into all harmonics, and model regularization and spectral truncation are required. A powerful perspective offered by the proposed technique is formal appraisal of the spatial resolution and variance of the resulting field averages. We use the Green’s functions for the Neumann boundary value problem to link the satellite observations to the radial magnetic field on the CMB and estimate field averages using a modified Backus-Gilbert inversion approach. Our approach builds on the Subtractive Optimally Localized Averages (SOLA) method developed in helioseismology, that seeks averaging kernels as close as possible to a chosen target kernel. We are able to account for both internal and external field sources and can easily incorporate data error covariance information, for example describing along-track serial error correlation. As a proof of concept we present a global map collecting local estimates of the radial main field (MF) constructed on a grid at the CMB with one degree spacing in latitude and longitude, derived from one month of three component vector magnetic field data collected by the 〈span〉Swarm〈/span〉 satellite trio, using data from dark and geomagnetically quiet times. Using sums and differences of the field components taken along track and in the east-west direction we obtain estimates with spatial resolution kernel widths varying between 18 and 54 degrees depending on the latitude, and a standard deviation of approximately 10μT (i.e. 5% of the mean CMB field amplitude). The morphology of our CMB radial field map agrees well with results from conventional spherical harmonic field models. In a second application, we determine local estimates of the average rate of change, or secular variation (SV), of the radial field at the CMB, initially considering two year time windows, and performing the analysis on data collected by either the 〈span〉Swarm〈/span〉 or CHAMP satellites. We obtain stable local estimates of the SV at the CMB, and present maps of estimates with averaging kernel widths of approximately 42, 33 and 30 degrees on the equator, with corresponding standard derivations of 0.25, 2.5 and 5 μT/yr. By subtracting SV estimates constructed at different epochs we are able to calculate the local aggregated secular acceleration (SA) and to study its time changes. Differencing SV estimates 2 years apart, and considering an averaging kernel width of 42 degrees on the equator, we obtain SA maps very similar to those found in the CHAOS-6-x7 field model truncated at SH degree 10. Using our approach we are able to directly control the width of the spatial averaging kernel and the length of the time window, enabling us to directly study the robustness of the inferred SA. Pushing to higher resolution in time, considering one year differences of SV estimates constructed using one year windows, we are able to track the evolution of coherent SA structures in time-longitude plots at the equator. At 25° W in mid 2007 we find a distinctive SA ’cross-over’ event, with strong, oppositely signed and adjacent, SA features rapidly changing sign within a year. Our method is well suited for studying such spatio-temporally localized SA events at high resolution; there will be further opportunities for such investigations as the time series of data provided by the 〈span〉Swarm〈/span〉 mission lengthens.〈/span〉

Description: 〈span〉〈div〉Summary〈/div〉The evolution of the ratio between P- and S-wave velocity (〈span〉V〈/span〉〈sub〉P〈/sub〉/〈span〉V〈/span〉〈sub〉S〈/sub〉) with increasing fluid-saturated porosity is computed for isotropic rocks containing spheroidal pores. The ratio 〈span〉V〈/span〉〈sub〉P〈/sub〉/〈span〉V〈/span〉〈sub〉S〈/sub〉 is shown to either decrease or increase with increasing porosity, depending on the aspect ratio α of the pores, fluid to solid bulk modulus ratio ζ, and Poisson’s ratio ν〈sub〉0〈/sub〉 of the solid constituents of the rock. A critical initial Poisson’s ratio ν〈sub〉0, crit〈/sub〉 is computed, separating cases where 〈span〉V〈/span〉〈sub〉P〈/sub〉/〈span〉V〈/span〉〈sub〉S〈/sub〉 increases (if ν〈sub〉0〈/sub〉 〈 ν〈sub〉0, crit〈/sub〉) or 〈span〉decreases〈/span〉 (if ν〈sub〉0〈/sub〉 〉 ν〈sub〉0, crit〈/sub〉) with increasing porosity. For thin cracks and highly compressible fluids, ν〈sub〉0, crit〈/sub〉 is approximated by 0.157 ζ/α, whereas for spherical pores ν〈sub〉0, crit〈/sub〉 is given by 0.2 + 0.8ζ. When ν〈sub〉0〈/sub〉 is close to ν〈sub〉0, crit〈/sub〉, the evolution of 〈span〉V〈/span〉〈sub〉P〈/sub〉/〈span〉V〈/span〉〈sub〉S〈/sub〉 with increasing fluid-saturated porosity is near neutral and depends on subtle changes in pore shape and fluid properties. This regime is found to be relevant to partially dehydrated serpentinites in subduction zones (porosity of aspect ratio near 0.1 and ζ in the range 0.01–0.1), and makes detection of these rocks and possibly elevated fluid pressures difficult from 〈span〉V〈/span〉〈sub〉P〈/sub〉/〈span〉V〈/span〉〈sub〉S〈/sub〉 only.〈/span〉

Description: 〈span〉〈div〉Summary〈/div〉Secondary microseismic sources emit seismic waves over long time spans. Reoccurring signals with similar slowness and frequency therefore arrive at seismic arrays. Blind source separation techniques can be used to identify and isolate such reoccurring signals from other signals and from diffuse seismic noise. Along these lines, we use non-negative matrix factorization as blind source separation technique to decompose continuous seismic array records. We model the recorded energy as a mixture of a few components with static slowness-frequency and time dependent amplitudes. Components and amplitudes are fitted to optimally explain the recorded seismic energy over time. These components represent secondary microseismic signals with quasi-static slowness-frequency vector and fluctuating amplitude. Each fitted component reveals the geographical origin (through the slowness-frequency vector) and time evolution of an active secondary microseism with high precision because it is separated from other signals and diffuse seismic noise. Furthermore, relative travel times can be automatically extracted for the signals that correspond to a specific component that can potentially be used in tomographic studies. We show two examples of seismic signals that were extracted with this technique, one focusing on P-waves from the typhoons Goni and Atsani, and another showing secondary microseism PKP signals from typhoon Glenda.〈/span〉

Description: 〈span〉〈div〉Summary〈/div〉Locating microseismic events is essential for many areas of seismology including volcano and earthquake monitoring and reservoir engineering. Due to the large number of microseismic events in these settings, an automated seismic location method is required to perform real time seismic monitoring. The measurement environment requires a precise and noise-resistant event location method for seismic monitoring. In this paper, we apply Multichannel Coherency Migration (MCM) to automatically locate microseismic events of induced and volcano-tectonic seismicity using sparse and irregular monitoring arrays. Compared to other migration-based methods, in spite of the often sparse and irregular distribution of the monitoring arrays, the MCM can show better location performance and obtain more consistent location results with the catalogue obtained by manual picking. Our MCM method successfully locates many triggered volcano-tectonic events with local magnitude smaller that 0, which demonstrates its applicability on locating very small earthquakes. Our synthetic event location example at a carbon capture and storage site shows that continuous and coherent drilling noise in industrial settings will pose great challenges for source imaging. However, automatic quality control techniques including filtering in the frequency domain and weighting are used to automatically select high quality data, and can thus effectively reduce the effects of continuous drilling noise and improve source imaging quality. The location performance of the MCM method for synthetic and real microseismic datasets demonstrates that the MCM method can perform as a reliable and automatic seismic waveform analysis tool to locate microseismic events.〈/span〉

Description: 〈span〉〈div〉Summary〈/div〉The unification of local height systems has been a classical geodetic problem for a long time, the main challenges of which are the estimation of offsets between different height systems and the correction of tilts along the levelling lines. It has been proposed to address these challenges with clock networks. The latest generation of optical clocks as well as the dedicated frequency links, e.g., optical fibres, are now approaching to deliver the comparison of frequencies at the level of 1.0 × 10〈sup〉−18〈/sup〉. It corresponds to an accuracy of about 1.0 cm in height difference. Clock networks can thus serve as a powerful tool to connect local height systems. To verify the idea, we carried out simulations using the EUVN/2000 (European Unified Vertical Network) as a priori input. Four local height systems were simulated from the EUVN/2000 by introducing individual offsets and tilts, and were re-unified by using measurements in clock networks. The results demonstrate the great potential of clock networks for height system unification. In case that the offsets between different height systems and tilts along national levelling lines in both longitudinal and latitudinal directions are considered, three or four clocks measurements for each local region are sufficient for the unification. These clocks are to be interconnected and should be properly arranged so that they can sense the levelling tilts where necessary. Our results also indicate that even clocks with one magnitude poorer accuracy than the desired ones can still unify the height systems to some extent, but it may cause a shift for the re-unified system.〈/span〉

Description: 〈span〉〈div〉Summary〈/div〉We investigate slip-distribution models of the 2011 Tohoku earthquake, with a particular focus on diffracted tsunamis and uplift-induced waves along the back-arc region of the Japanese Island Arc.The 2011 Tohoku earthquake produced a large amplitude tsunami that diffracted around Kyushu Island before reaching Korea. At the same time, this earthquake co-seismically induced short-period small-amplitude sea waves in the East Sea. We performed tsunami simulations using seven fault models of the Tohoku earthquake to examine whether the models can accurately reproduce the observed waveforms in the open sea of the western Pacific Ocean, the South Sea of Korea, and the coast of the East Sea. For each fault model, we investigate tsunami features due to geomorphological characteristics of the Korean Peninsula in the Korea offshore. To determine which slip distribution model shows a good performance in the tsunami simulations, we set three criteria; the delay time between observations and synthetic waveforms, the normalized mean residual, and the normalized RMS misfit. Depending on the study region, all models show varying degrees of accuracy. The fault dimensions and the amount of slip have a larger effect on the RMS misfit then the slip distribution patterns of the fault models for observations along the Korean coast and the western coast of Japan.〈/span〉

Description: 〈span〉〈div〉SUMMARY〈/div〉The study of radiation patterns in poroelastic media allows us to visually explore the possibility of reconstructing crucial model properties describing reservoir rocks, and to examine the coupling effects between different parameters during full-waveform inversion (FWI). In this paper, we derive analytical formulae for the radiation patterns of single parameter perturbations in fluid-saturated porous media by deriving scattered wavefields based on plane-wave theory and the far-field approximation. We illustrate these scattered wavefields via their radiation patterns expressed as a function of the angle between the incident and scattered waves. To simplify the algebra, we consider poroelastic waves at seismic (low) frequencies, where the fast compressional wave and shear wave are propagating modes but the slow compressional wave is severely dispersive. To verify our derivation of the analytical radiation patterns, we also compute them numerically by perturbing one parameter at a single point, keeping the other parameters fixed at their background values. We find that all the analytical radiation patterns match the wavefronts of the numerically computed scattered wavefields, well indicating that our derivations are correct. Parameters such as the solid density, fluid density, viscosity of the fluid, and intrinsic permeability, have similar radiation patterns and thus show strong coupling effects. Therefore, we anticipate difficulties in recovering these parameters in a multi-parameter FWI procedure. In an attempt to mitigate these trade-offs, we analyze different parameterizations which result in different radiation patterns. As the patterns that we observe are similar to those of the elastic case, we anticipate that parameter separation might be easier when inverting for velocities than for moduli.〈/span〉

Description: 〈span〉〈div〉Abstract〈/div〉In the Sea of Marmara, areas of gas seepage or cold seeps are tightly related to the faults system and understanding the spatial and temporal dynamics in gas-related processes is crucial for geohazard mitigation. Although acoustic surveys proved to be efficient in detecting and locating cold seeps, temporal variability or trends in the gas-related processes are still poorly understood. Two arrays of 10 ocean bottom seismometers were deployed in the western part of the Sea of Marmara in 2011 and 2014, respectively. In addition to the local seismic events, the instruments recorded a large number of short duration events and long-lasting tremors. Short duration events are impulsive signals with duration 〈 1 s, amplitude well above the noise level and a frequency spectrum with one or two narrow peaks. They are not correlated from one site to another, suggesting a very local source. Tremors consist of sequences of clustered impulsive signals lasting for minutes to more than an hour with a multi-peak frequency spectrum. Based on evidence of known seepage and by analogy with volcanic and hydrothermal models, we suggest that short duration events and tremors are associated with gas migration and seepage. There is a relationship between tremors associated with gas emission and the local seismicity, although not systematic. Rather than triggering gas migration out of the seabed, locally strong earthquakes act as catalysts when gas is already present or gas emission is already initiated.〈/span〉

Description: 〈span〉〈div〉SUMMARY〈/div〉Long-period Rayleigh wave Horizontal to Vertical amplitude (H/V) ratios at a station provide information about local earth structure that is complementary to phase velocity. However, a number of studies (Ferreira and Woodhouse, 2007; Lin 〈span〉et al.〈/span〉, 2012) have observed that significant scatter appears in these measurements making it difficult to use H/V ratio measurements to resolve earth structure. Some of the scatter in these measurements has been attributed to local geological structure while some has remained unaccounted for. Most Global Seismographic Network (GSN) stations contain two nearby high-quality broadband seismometers (e.g. in the same vault, but on different piers or in different boreholes). For each broadband sensor in the IRIS/USGS component of the GSN, we estimate H/V ratios of fundamental mode Rayleigh waves using M 〉 6.5 earthquakes from 2001 to 2018 (around 19,000 measurements). We compute these ratios at a number of discrete periods (25, 50, 75, 100, and 150 s) and find that for well-isolated Rayleigh waves (windows where the correlation coefficients between radial and the phase-shifted vertical components are greater than 0.9) significant scatter in H/V ratios occurs between co-located sensors (greater than 25 per cent at 100 s period). This suggests the scatter in H/V ratio measurements can be at least partially attributed to extremely local phenomena such as sensor emplacement in the vault. We also find that H/V ratios can vary as a function of event back-azimuth, indicating that care must be taken when computing average ratios for a station, as a large number of events from a given region could bias H/V ratio measurements at a station.〈/span〉

Description: 〈span〉〈div〉SUMMARY〈/div〉We analyse noise characteristics of spatial and temporal correlation of 260 continuous GPS sites from Crustal Movement Observation Network of China (CMONOC). These data sets were mainly collected between 2010 and 2016, with an average of 6 yr of position time-series. In the functional analysis, a clear regional dependence of seasonal movements has been observed and other significant periodic signals are detected nearby the GPS draconitic period and its harmonics. The distribution of these periodic signals shows a spatial correlation, along with non-negligible local inconsistencies. In the stochastic analysis, impacts of the periodicities on the noise assessment have been investigated and Maximum likelihood estimation is used to study noise properties for deseasonalized residual time-series having the seasonal signals removed and filtered residual time-series having other periodic signals removed further. We demonstrate that for both solutions, the flicker noise is thought to be the dominant time-correlated noise and velocity uncertainties may be underestimated 8–10 times if assuming a pure white noise. Ignoring the periodicities could bias the estimation of noise amplitude, spectral index and velocity uncertainty. After removing the periodic signals, the median flicker noise magnitude shows an average of ∼10 per cent reductions and the noise process shifts closer towards white noise. A correlation between the index variations and RMS variations has been observed, indicating that the index varies more significantly for sites with more periodic signals removed. Besides, the spectral index in the vertical component has a better spatial correlation than that in the horizontal and the spatial distribution of the index of deseasonalized solutions seems to correlate well with the amplitudes of seasonal signals, probably implying common sources of spatial variations of these characteristics. Furthermore, analyses of intersite correlations indicate that the correlation induced from the periodic signals displays a similar pattern to the deseasonalized solutions, confirming that the period pattern is spatially correlated and can induce time-series correlations. Finally, the stochastic processes of the common mode noise are predominantly featured by spatially correlated flicker and white noise over a wide range, consistent with previous results.〈/span〉

Description: 〈span〉〈div〉SUMMARY〈/div〉As temperature increases with depth and the creep resistance of rock decreases exponentially, a high-viscosity sub-lithospheric layer, just beneath the ‘elastic’ lithosphere is expected to exist. Depending on the temperature profile, a low-viscosity asthenosphere may also exist if the temperature deeper down gets high enough. Since the temperature profile is expected to change laterally – especially from below the oceans to cratonic areas underneath continents, rock properties of the lithosphere, high-viscosity sub-lithosphere and low-viscosity asthenosphere are expected to change laterally. Our aim is to constrain sub-lithospheric properties (depth, thickness and viscosity), lateral lithospheric thickness variations and asthenospheric properties using observed GIA data. A Coupled Laplace-Finite Element Method is used to compute gravitationally self-consistent sea level with time-dependent coastline and rotational feedback in addition to changes in deformation, gravity and the state of stress. We start with the VM5a-ICE-6G_C model combination and then modify the lithospheric, sub-lithospheric and asthenospheric properties (including lateral thickness variation) while keeping the mantle viscosities the same as VM5a. Through this study, we confirm that the sub-lithospheric and asthenospheric properties can significantly affect the predicted global relative sea level (RSL), present-day gravity rate-of-change (g-dot) and uplift rate (u-dot) in Laurentia and Fennoscandia. In addition, incorporating the elastic lithosphere with lateral thickness variation, sub-lithosphere and asthenosphere can improve the fit to global RSL, but the predicted peak values of g-dot and u-dot in Laurentia may decrease slightly but not significant enough to affect the fit to the observed data. Our results prefer an elastic lithosphere that has maximum thickness of 140 km under continental cratons but reduces to 60 km underneath the oceans. The results preferred depth of the asthenospheric bottom is around 190–200 km with asthenospheric viscosity around 10〈sup〉20〈/sup〉Pa s. Finally, we show that the best laterally heterogeneous mantle model we found in previous publication when combined with the lithosphere with lateral thickness variaion gives the best fit to global RSL and peak g-dot and u-dot in Laurentia simultaneously.〈/span〉

Description: 〈span〉〈div〉SUMMARY〈/div〉We implement an analytical model based on flexural deflection of a thin elastic disc to investigate the magnitude of lithospheric decompression caused by deglaciation at upper crustal magmatic reservoirs. Considering a published numerical climate model describing the space–time evolution of deglaciation after the Last Glacial Maximum (LGM) along the Southern Volcanic Zone (SVZ) of the Andes, we demonstrate that changes in pressure at upper crustal levels (〈10 km depth) at the scale of several hundred years are of the order of 10–100 MPa. Total decompression and decompression rate (300–150 kPa yr〈sup〉−1〈/sup〉) are 1–2 orders of magnitude larger than values previously estimated by other authors who assume that glacial loads are supported by an elastic half-space, that is, of infinite elastic thickness. The large decompression caused by flexural unbending of an elastic plate of finite thickness as assumed here can easily surpass the tensile strength of rocks (5–20 MPa), creating adequate conditions for failure of the reservoir walls, dike propagation inside and outside the reservoir and the eventual collapse of the reservoir accompanying an explosive eruption. We apply our results to the analysis of post-glacial eruptions of SVZ volcanoes, which erupted large volumes (〉10 km〈sup〉3〈/sup〉) of mafic ignimbrites hundreds to thousands of years after deglaciation onset. We show that this time lag is necessary to achieve a decompression of several tens of megapascals at depths of several kilometres that are consistent with the location of magmatic reservoirs as estimated by independent petrologic, seismic and/or geodetic studies. Moreover the northward increase of this time lag is in agreement with a smaller size of the Andean ice cap in the north than in the south during the LGM. For wet, volatile-rich magmas typical of subduction zones, the effect of large decompression at upper crustal reservoirs caused by flexural unbending of the lithosphere after deglaciation could play a major role in promoting large explosive eruptions through devolatization of the magma, during past deglaciation events as demonstrated here for the LGM along the SVZ and current accelerated ice retreat caused by climate change over large segments of subduction-related arcs at higher latitudes.〈/span〉

Description: 〈span〉〈div〉SUMMARY〈/div〉The unification of local height systems has been a classical geodetic problem for a long time, the main challenges of which are the estimation of offsets between different height systems and the correction of tilts along the levelling lines. It has been proposed to address these challenges with clock networks. The latest generation of optical clocks as well as the dedicated frequency links, for example optical fibres, are now approaching to deliver the comparison of frequencies at the level of 1.0 × 10〈sup〉−18〈/sup〉. It corresponds to an accuracy of about 1.0 cm in height difference. Clock networks can thus serve as a powerful tool to connect local height systems. To verify the idea, we carried out simulations using the EUVN/2000 (European Unified Vertical Network) as 〈span〉apriori〈/span〉 input. Four local height systems were simulated from the EUVN/2000 by introducing individual offsets and tilts, and were reunified by using measurements in clock networks. The results demonstrate the great potential of clock networks for height system unification. In case that the offsets between different height systems and tilts along national levelling lines in both longitudinal and latitudinal directions are considered, three or four clocks measurements for each local region are sufficient for the unification. These clocks are to be interconnected and should be properly arranged so that they can sense the levelling tilts where necessary. Our results also indicate that even clocks with one magnitude poorer accuracy than the desired ones can still unify the height systems to some extent, but it may cause a shift for the reunified system.〈/span〉

Description: 〈span〉〈div〉SUMMARY〈/div〉Waveform backprojection (BP) is a key technique of earthquake-source imaging, which has been widely used for extracting information of earthquake source evolution that cannot be obtained by kinematic source inversion. The technique enjoys considerable popularity, owing to the simplicity of its implementation and the robustness of its processing, but the physical meaning of BP images has remained elusive. In this study, we reviewed the mathematical representation of BP and hybrid BP (HBP) methods, following the pioneering work of Fukahata 〈span〉et al〈/span〉. (〈a href="http://academic.oup.com/gji#bib15"〉2014〈/a〉), to clarify the physical implications of BP images. We found that signal intensity in BP and HBP images is scaled with the amplitude of the Green’s function that corresponds to a unit-step slip, which results in the signal intensity being depth dependent. We propose variants of BP and HBP, which we call kinematic BP and HBP, respectively, to relate the BP signal intensity to slip motion of an earthquake by modifying the normalizing factors used in the original BP and HBP methods. The original BP and HBP images remain useful for assessing the spatiotemporal strength of the wave radiation, which scales with the amplitude of the Green’s function, whereas the kinematic BP and HBP methods are suitable for imaging the slip motion that is responsible for the high-frequency radiation produced during the source-rupture process.〈/span〉

Description: 〈span〉〈div〉SUMMARY〈/div〉 〈strong〉Moment Tensors of hydraulically induced AEs:〈/strong〉 Hydraulic fracturing is an important technique in the development of enhanced geothermal systems and unconventional resources. Although the fracture modes induced by hydraulic fracturing influence the recovery efficiency of the resources, the current understanding of this relationship is insufficient. In this study, we considered the acoustic emissions (AEs) induced during hydraulic fracturing under uniaxial loading conditions in the laboratory, and applied a moment tensor analysis by carefully correcting the coupling condition and directivity of AE transducers. Experiments were conducted for two types of Kurokami–jima granite samples: those with a rift plane perpendicular (Type H) or parallel (Type V) to the expected direction of fracture propagation (i.e. along the loading axis). In the experiments, both sample types experienced a significant number of shear, tensile and compressive events. The dominant fracture mode for Type H samples is found to be tensile events in which the fracture plane is parallel to the loading axis, whereas for Type V samples, shear events are dominant. This difference suggests that the dominant fracture modes induced by hydraulic fracturing are highly dependent on the relationship between the direction of fracture propagation and orientation of pre-existing weak planes.〈/span〉

Description: 〈span〉〈div〉SUMMARY〈/div〉We present a sensitivity analysis aimed at testing whether observables related to glacial isostatic adjustment can support or refute the occurrence of a low viscosity melt-rich layer (MRL) above the mantle transition zone, as required by the ‘transition-zone water-filter’ model (Bercovici & Karato 2003). In total, 1600 model runs were performed sampling a range of MRL thicknesses (1, 10 and 20 km) and viscosities (10〈sup〉15〈/sup〉–10〈sup〉19〈/sup〉 Pa·s), plausible viscosity values in the upper and lower mantle regions and four distinct ice histories. To determine decay time constraints, we consider relative sea level (RSL) data from two sites [Ångerman River (ÅR), Sweden and Richmond Gulf (RG), Canada] and use a new method of observational sea level data correction. Comparing model output of postglacial decay times and ${\skew{6}\dot{J}_2}$ to observational constraints, we find numerous possible solutions, largely as a result of parameter trade-off. The investigated observables are sensitive to the existence of an MRL and reasonable variations in its thickness and viscosity. The magnitude and nature of this sensitivity varies between the two data types as well as the adopted background viscosity structure. Decay time results from either considered location do not strictly support or exclude MRL existence. However, both locations offer MRL viscosity requirements for given thicknesses, with ÅR being more restrictive. RG constraints allow MRL viscosities as low as 10〈sup〉16〈/sup〉 Pa·s (10 km) and 10〈sup〉17〈/sup〉 Pa·s (20 km). ÅR results narrow these permitted viscosity ranges to 10〈sup〉18〈/sup〉 Pa·s or greater for both 10 and 20 km MRL thicknesses. In the case of a 1 km thick MRL, ÅR constraints permit the viscosity to be as low as 10〈sup〉17〈/sup〉 Pa·s, whereas those of RG permit any MRL viscosity. The decay time observations are satisfied by only a small subset of ‘background’ mantle viscosities (regardless of the MRL properties), none of which support a spherically symmetric solution of Earth viscosity. Finally, comparing model output to the observed ${\skew{6}\dot{J}_2}$ value did not provide any constraints on MRL properties. However, our results show that this observable has a strong preference for viscosity values in the lower mantle that are equal to or greater than 10〈sup〉22〈/sup〉 Pa·s.〈/span〉

Description: 〈span〉〈div〉SUMMARY〈/div〉We investigated an 〈span〉M〈/span〉〈sub〉w〈/sub〉 ∼ 6.2 earthquake doublet on the border of the USA and Canada using ALOS2 Light-of-Sight displacements and GPS measurements. We selected three L-band ALOS-2 interfergorams with temporal baselines of one yr to extract coseismic deformation maps, in which master and slave images were both acquired in July. A subpixel-based alignment and another range spectral splitting techniques under the GAMMA InSAR software framework were applied to improve the interferometric coherence and reduce the effects of phase anomalies in two of the three interferometric pairs due to either ionospheric delay or a potential focusing issues in the generation of the ALOS2 SLC data. The updated interferograms convincingly reveal deformation fringe patterns produced by the two earthquakes. We conducted a nonlinear geophysical inversion to estimate the geometric parameters of the earthquakes with the InSAR and GPS measurements. The best-fitting model shows that a thrust faulting on a reverse fault and left-lateral strike-slip faulting on a nearly vertical fault with the centroid depths of 9.3±0.6 and 8.4±0.7 km, respectively, are most likely responsible for the earthquake doublet. The eastern Denali fault (EDF) and Duke River fault are major active faults in the region and the earthquake doublet could be due to reactivation of the part of the two faults system.〈/span〉

Description: 〈span〉〈div〉SUMMARY〈/div〉In this paper, we present a series of mathematical abstractions for seismologically relevant wave equations discretized using finite-element methods, and demonstrate how these abstractions can be implemented efficiently in computer code. Our motivation is to mitigate the combinatorial complexity present when considering geophysical waveform modelling and inversion, where a variety of spatial discretizations, material models, and boundary conditions must be considered simultaneously. We accomplish this goal by first considering three distinct classes of abstract mathematical models: (1) those representing the physics of an underlying wave equation, (2) those describing the discretization of the chosen equation onto a finite-dimensional basis and (3) those describing any spatial transforms. A full representation of the discrete wave equation can then be constructed using a hierarchical nesting of models from each class. Additionally, each class is functionally orthogonal to the others, and with certain restrictions models within one class can be interchanged independently from changes in another. We then show how this recasting of the relevant equations can be implemented concisely in computer software using an abstract object-oriented design, and discuss how recent developments in the numerical and computational sciences can be naturally incorporated. This builds to a set of results where we demonstrate how the developments presented can lead to an implementation capable of multiphysics waveform simulations in completely unstructured domains, on both hypercubical and simplical spectral-element meshes, in both two and three dimensions, while remaining concise, efficient and maintainable.〈/span〉

Description: 〈span〉〈div〉SUMMARY〈/div〉Traveltime approximation plays an important role in seismic data processing, for example, anisotropic parameter estimation and seismic imaging. By exploiting seismic traveltimes, it is possible to improve the accuracy of anisotropic parameter estimation and the resolution of seismic imaging. Conventionally, the traveltime approximations in anisotropic media are obtained by expanding the anisotropic eikonal equation in terms of the anisotropic parameters and the elliptically anisotropic eikonal equation based on perturbation theory. Such an expansion assumes a small perturbation and weak anisotropy. In a realistic medium, however, the assumption of small perturbation likely breaks down. We present a retrieved zero-order deformation equation that creates a map from the anisotropic eikonal equation to a linearized partial differential equation system based on the homotopy analysis method. By choosing the linear and nonlinear operators in the retrieved zero-order deformation equation, we develop new traveltime approximations that allow us to compute the traveltimes for a medium of arbitrarily strength anisotropy. A comparison of the traveltimes and their errors from the homotopy analysis method and from the perturbation method suggests that the traveltime approximations provide a more reliable result in strongly anisotropic media.〈/span〉

Description: 〈span〉〈div〉SUMMARY〈/div〉Guided waves in a water layer overlaying an elastic half-space are known as normal modes. They are often present in seismic recordings at long offsets in shallow-water environment and generally considered coherent noise. The normal modes, however, carry important information about the near-surface and, as demonstrated by a number of authors, can be used to obtain the shallow velocity model. There is a growing evidence that the latter needs not to be isotropic due to various geological reasons. Motivated by that, we consider the normal-mode propagation in case the elastic half-space exhibits orthorhombic anisotropy. We derive the period equation that describes the normal-mode phase velocity dispersion. To simplify the complicated expression, we present acoustic and ellipsoidal orthorhombic approximations. We also outline the approach towards the group velocity and group azimuth calculation and apply it to the ellipsoidal case to obtain concise and intuitive expressions. Using numerical test, we study the relation between phase and group domains in elastic orthorhombic case. The deviation between velocities and azimuths in these domains is the strongest for low frequencies and it rapidly decreases with increasing frequency. For higher frequencies, the anisotropy effects of the underlaying half-space are barely detectable since the observed signal is composed mainly of the direct acoustic wave, resulting in the two domains being nearly indistinguishable.〈/span〉

Description: 〈span〉〈div〉SUMMARY〈/div〉Understanding why Earth’s lithosphere is divided into several plates while other terrestrial bodies have unbroken lids is a long-standing challenge, often addressed with the help of numerical modelling. A key mechanism defining the transition between these two convective regimes is the formation of shear zones that cut through the entire lithosphere in regions with high stresses. Here we present a modelling study in which lithospheric stresses resulting from small-scale convection in the upper mantle are analysed. We perform model simulations that include elasticity and a free surface and evaluate how these physical complexities affect stress distribution inside the lithosphere, which in turn controls the depths of yielding and the possible initiation of subduction. We show that the spatial distribution of stress is significantly altered by the presence of elastic deformation only when the model lithosphere acts as a thick plate capable of bending. Whether or not this is the case depends on the viscosity model. For an Arrhenius viscosity limited by a cut-off value that produces an essentially rigid lid, flexure dominates the observed lithospheric stress pattern in simulations with a free surface. The amplitudes of the stress are, when a free surface is assumed but elasticity is neglected, largely overestimated. Including both a free surface and elasticity results in stresses with maximum amplitudes close to those observed in the traditional models with a viscous rheology and a free-slip upper boundary, suggesting that having no additional complexity is, in a way, better than employing just a free surface. We also demonstrate how the use of impermeable free-slip side boundaries can result in the formation of unnatural, laterally locked convection cells, and bias the results of a parametric study. For each point in the parameter space, we perform several simulations with slightly different initial temperature fields in order to statistically eliminate the occurrence of locked states.〈/span〉

Description: 〈span〉〈div〉SUMMARY〈/div〉The design of an array configuration is an important task in array seismology during experiment planning. Often the array response function (ARF), which depends on the relative position of array stations and frequency content of the incoming signals, is used as the array design criterion. In practice, additional constraints and parameters have to be taken into account, for example, land ownership, site-specific noise levels or characteristics of the seismic sources under investigation. In this study, a flexible array design framework is introduced that implements a customizable scenario modelling and optimization scheme by making use of synthetic seismograms. Using synthetic seismograms to evaluate array performance makes it possible to consider additional constraints. We suggest to use synthetic array beamforming as an array design criterion instead of the ARF. The objective function of the optimization scheme is defined according to the monitoring goals, and may consist of a number of subfunctions. The array design framework is exemplified by designing a seven-station small-scale array to monitor earthquake swarm activity in Northwest Bohemia/Vogtland in central Europe. Two subfunctions are introduced to verify the accuracy of horizontal slowness estimation; one to suppress aliasing effects due to possible secondary lobes of synthetic array beamforming calculated in horizontal slowness space and the other to reduce the event’s mislocation caused by miscalculation of the horizontal slowness vector. Subsequently, a weighting technique is applied to combine the subfunctions into one single scalar objective function to use in the optimization process.〈/span〉

Description: 〈span〉〈div〉SUMMARY〈/div〉We develop a new method for measuring ellipticity of Rayleigh waves from ambient noise records by degree-of-polarization (DOP) analysis. The new method, named DOP-E, shows a good capability to retrieve accurate ellipticity curves separated from incoherent noise. In order to validate the method we perform synthetic tests simulating noise in a 1-D earth model. We also perform measurements on real data from Antarctica and Northern Italy. Observed curves show a good fit with measurements from earthquake records and with theoretical ellipticity curves. The inversion of real data measurements for 〈span〉vS〈/span〉 structure shows a good agreement with previous models. In particular, the shear-wave structure beneath Concordia station shows no evidence of a significant layer of liquid water at the base of the ice. The new method can be used to measure ellipticity at high frequency and therefore it will allow the imaging of near-surface structure, and possibly of temporal changes in subsurface properties. It promises to be useful to study near-surface processes in a wide range of geological settings, such as volcanoes, fault zones and glaciers.〈/span〉

Description: 〈span〉〈div〉SUMMARY〈/div〉The interaction of subducted oceanic lithosphere with the discontinuities of the mantle transition zone (MTZ) provides insight into the composition and temperature of the subducted slab as well as potential melting of the slab or the surrounding mantle and loss of volatiles from the slab. Detailed mapping of the structure of the MTZ will help to better understand how slabs transport material and volatiles into the mantle and how phase transitions affect the slab dynamics. Here we use a dense network of seismic stations in northern Anatolia to image the structure of the MTZ discontinuities in detail using 〈span〉P〈/span〉-wave receiver functions. With a station spacing of about 7 km and a surface footprint of ∼35 km × ∼70 km, analysing receiver functions calculated from teleseismic earthquakes that occurred during an ∼18-month deployment produced clear images of where the MTZ interacts with the Tethys/Cyprus slabs that either lie flat on the 660-km discontinuity or pass into the lower mantle. We observe an undulating 660-km discontinuity depressed by up to 30 km and a slightly depressed (1–2 km) 410-km discontinuity, apparently undisturbed by the slab. The MTZ is thickened to ∼270 km as result of the cool slab in the MTZ influencing the 660-km discontinuity and includes an arrival at ∼520-km depth likely from the top of a flat lying slab or a discontinuity related to a solid–solid phase transition in the olivine component of the mantle. We find evidence for low-velocity zones both above and below the 410-km discontinuity and above the 660-km discontinuity. The low-velocity zones around the 410-km discontinuity might be the result of hydration of the MTZ from the slab and upward convection of MTZ material into the upper mantle. The origin of the low-velocity zone around the 660-km discontinuity is less clear and could be related to sedimentation of subducted mid-ocean ridge basalts. The small footprint of the seismic array provides accurate information on the structure of the MTZ in an area influenced by subduction and shows small-scale changes in MTZ structure that might be lost in studies covering larger areas with sparser sampling.〈/span〉

Description: 〈span〉〈div〉SUMMARY〈/div〉We perform non-linear time-series analysis on a harmonic tremor seismogram recorded at 830 m away from the centre of the crater during the 2011 eruption at Shinmoedake, Japan. We found features suggesting the existence of period doubling bifurcation in the harmonic tremor signal, implying that the harmonic tremor might be generated by a non-linear process. In order to quantify the non-linearity in the harmonic tremor signal, we measure the correlation dimension 〈span〉D〈/span〉 and the maximal Lyapunov exponent λ. For one short but stable segment of the harmonic tremor seismogram, we obtained 〈span〉D〈/span〉 = 1.12 and λ = 0.03 s〈sup〉−1〈/sup〉. This result implies that the stable oscillation of the harmonic tremor is predominantly a limit cycle with small amounts of chaos present. We then use surrogate data analysis to check that our measurements of 〈span〉D〈/span〉 and λ do not include any false positive detection of non-linearity. Limit cycles imply that the harmonic tremor is generated by self-sustained oscillations. We show that the autonomous Julian tremor model is able to exhibit period doubling bifurcation. We also show that the non-autonomous Julian tremor model with a step like increase and decrease in input pressure is able to exhibit oscillations of varying amplitude while keeping a constant frequency spectrum. Both phenomena were observed at Shinmoedake. We also demonstrate that the non-autonomous Julian tremor model with a transient input pressure is able to exhibit long period-like events in addition to harmonic tremor-like events, implying that the same non-linear mechanism could be responsible for the generation of both type of events.〈/span〉

Description: 〈span〉〈div〉SUMMARY〈/div〉Reservoir models are numerical representations of the subsurface petrophysical properties such as porosity, volume of minerals and fluid saturations. These are often derived from elastic models inferred from seismic inversion in a two-step approach: first, seismic reflection data are inverted for the elastic properties of interest (such as density, 〈span〉P〈/span〉-wave and 〈span〉S〈/span〉-wave velocities); these are then used as constraining properties to model the subsurface petrophysical variables. The sequential approach does not ensure a proper propagation of uncertainty throughout the entire geo-modelling workflow as it does not describe a direct link between the observed seismic data and the resulting petrophysical models. Rock physics models link the two domains. We propose to integrate seismic and rock physics modelling into an iterative geostatistical seismic inversion methodology. The proposed method allows the direct inference of the porosity, volume of shale and fluid saturations by simultaneously integrating well-logs, seismic reflection data and rock physics model predictions. Stochastic sequential simulation is used as the perturbation technique of the model parameter space, a calibrated facies-dependent rock physics model links the elastic and the petrophysical domains and a global optimizer based on cross-over genetic algorithms ensures the convergence of the methodology from iteration to iteration. The method is applied to a 3-D volume extracted from a real reservoir data set of a North Sea reservoir and compared to a geostatistical seismic AVA.〈/span〉

Description: 〈span〉〈div〉SUMMARY〈/div〉We investigate slip-distribution models of the 2011 Tohoku earthquake, with a particular focus on diffracted tsunamis and uplift-induced waves along the backarc region of the Japanese Island Arc. The 2011 Tohoku earthquake produced a large amplitude tsunami that diffracted around Kyushu Island before reaching Korea. At the same time, this earthquake coseismically induced short-period small-amplitude sea waves in the East Sea. We performed tsunami simulations using seven fault models of the Tohoku earthquake to examine whether the models can accurately reproduce the observed waveforms in the open sea of the western Pacific Ocean, the South Sea of Korea, and the coast of the East Sea. For each fault model, we investigate tsunami features due to geomorphological characteristics of the Korean Peninsula in the Korea offshore. To determine which slip distribution model shows a good performance in the tsunami simulations, we set three criteria; the delay time between observations and synthetic waveforms, the normalized mean residual, and the normalized RMS misfit. Depending on the study region, all models show varying degrees of accuracy. The fault dimensions and the amount of slip have a larger effect on the RMS misfit then the slip distribution patterns of the fault models for observations along the Korean coast and the western coast of Japan.〈/span〉

Description: 〈span〉〈div〉Abstract〈/div〉Economically important deposits of nickel-copper and platinum group elements tend to be associated with small mafic-ultramafic intrusions characterized by compositional layering. Our rock magnetic study of an exploration drill core recovered from a typical small intrusion from the recently discovered Kun-Man'e ore field have found that concentration-independent hysteresis parameters of the studied rocks generally co-vary with petrological characteristics, such as major element abundances and compositions of olivine and clinopyroxene. Although monoclinic pyrrhotite and defect-poor magnetite both contribute to magnetic properties of the samples, we argue that hysteresis parameters are not affected significantly by pyrrhotite content. Rather, the variations in the hysteresis parameters reflect variations in median grain-size of magnetite crystals. The correlation between the rock-magnetic grain-size indicators and geochemical crystallization temperature indicators suggests that the size of magnetite crystals is controlled by cooling regime. Compared to conventional laborious crystal size distribution (CSD) studies, rock magnetic measurements can provide a quick and sensitive way to evaluate CSD of magmatic magnetite in mafic and ultramafic rocks.〈/span〉

Description: 〈span〉〈div〉SUMMARY〈/div〉The Reykjanes geothermal system is a high-temperature seawater system situated in SW-Iceland. Interferometric analysis of the Sentinel-1 satellite synthetic aperture radar (InSAR) data has been used to determine a time series of ground deformation induced by geothermal utilization between April 2015 and October 2017. Surface displacements have been estimated at coherent pixels, indicating a steady and linear subsidence within a sub-circular bowl centered on the well field at a maximum near-vertical rate of about 25 mm/yr, together with horizontal contraction. The average line-of-sight (LOS) displacement rates from ascending and descending tracks are inverted to determine the characteristics of the deformation source at depth, modeling the geothermal reservoir as a body of simple geometry within an elastic half space. The results indicate a deformation source at about 1 km depth contracting at a rate of (0.7–0.9) × 10〈sup〉5〈/sup〉 m〈sup〉3〈/sup〉/yr during the 2015–2017 period. Using pressure and temperature monitoring data at 900 m depth as well as an analysis of the reservoir structure and rock properties, we infer that the recent estimated volume change can be attributed to processes in the steam cap situated in the topmost part of the geothermal reservoir, in the 800–1200 m depth range. Processes involve a combination of compaction under pressure decrease and/or thermal contraction due to cooling of the rocks within or near the steam cap. The steam cap expanded in response to a sudden pressure drop resulting from the increase in extraction of geothermal fluids for a new power plant in 2006.〈/span〉

Description: 〈span〉〈div〉Abstract〈/div〉The seismic speed and anisotropy are two methods frequently used to image the assembly inside the Earth. We study the crust assembly beneath the Biga Peninsula and the surrounding area in Northwest Turkey using the accelerometer and broadband recordings where short-to-medium period (5-20 s) Love-Rayleigh surface waves are utilized to extract the group-phase speed data (fundamental mode). Single-station and two-station techniques are engaged to understand the detected surface waves for the speed and anisotropy assemblies. The single-station group speeds are inverted in a tomographic approach to attain the two-dimensional group speed diagrams. The least-squares inversion procedure is utilized to find the speed-depth profiles (one-dimensional) under each grid location. The one-dimensional results are cooperatively inferred to attain the three-dimensional appearance of the S-wave speeds below the measured region. This process is reiterated for Love and Rayleigh waves. Isotropic configuration is not sufficient to concurrently describe the present detected Love-Rayleigh surface waves. Vertical transverse isotropic crust assembly is found to better elucidate the detected data showing the Rayleigh-Love discrepancy. Complex arrangement of sills and dykes due to the widespread plutonic and volcanic activity in the region linked to the interaction between the Turkish plate and the African plate (northward subducting) is thought to depict the crust assembly deformations causing the detected long-wavelength vertical transverse isotropy. The mineral orientation within horizontal sills and vertical dykes following the magma flow, which is independent of seismic wavelength, adds to the detected anisotropy. The upper crust vertical transverse isotropy is mostly negative; i.e. SV-wave is faster than SH-wave, which is assumed to be due to the existence of dykes. The middle-to-lower crust vertical transverse isotropy is commonly positive; i.e. SH-wave is faster than SV-wave, which is assumed to be due to the existence of sills. The two-station analyses operating on cross-correlograms give analogous vertical transverse isotropic results to those of the single-station estimates.〈/span〉

Description: 〈span〉〈div〉Summary〈/div〉The representative elementary volume (REV) is a critically important concept in fractured rock investigations as it tells us at what scale the fractured domain can be represented by an anisotropic tensor as opposed to requiring the details of each individual fracture for modeling purposes. Whereas the REV size and corresponding tensor characteristics for the hydraulic conductivity (〈span〉K〈/span〉) in fractured rock have been the subject of numerous previous investigations, no studies to date have focused on the electrical conductivity (σ). This is despite the fact that geoelectrical measurements are arguably the most popular means of geophysically investigating fractured rock, typically via azimuthal resistivity surveying where the observed electrical anisotropy is commonly used to infer hydraulic characteristics. In this paper, we attempt to fill this void and present a systematic numerical study of the impacts of changes in fracture-network properties on the REV size and equivalent tensor characteristics for both the electrical and hydraulic conductivities. We employ a combined statistical and numerical approach where the size of the REV is estimated from the conductivity variability observed across multiple stochastic fracture-network realizations for various domain sizes. Two important differences between fluid and electric current flow in fractured media are found to lead to significant differences in the REV size and tensor characteristics for σ and 〈span〉K〈/span〉; these are the greater importance of the matrix in the electrical case and the single-power instead of cubic dependence of electric current flow upon aperture. Specifically, the REV for the electrical conductivity will always be smaller than that for the hydraulic conductivity, and the corresponding equivalent tensor will exhibit less anisotropy, often with notably different principal orientations. These findings are of key importance for the eventual interpretation of geoelectrical measurements in fractured rock, where we conclude that extreme caution must be taken when attempting to make the link to hydraulic properties.〈/span〉

Description: 〈span〉〈div〉Summary〈/div〉Stress is a crucial factor that affects the permeability, and the influence of stress on permeability has been studied extensively by many investigators. In this article the relationship between confining pressure and permeability is investigated. For this purpose, a void ratio- pressure model based on the theory of incompressible elastic porous media is first proposed, then, based on the model and Kozeny-Carman equation, a new pressure–permeability relationship is proposed. Finally the relationship is utilized to predict the permeability of coal and mudstone under different confining pressure, and these predictions are compared with other models and found to be in better agreement with experimental data reported recently.〈/span〉

Description: 〈span〉〈div〉Abstract〈/div〉We perform nonlinear time series analysis on a harmonic tremor seismogram recorded at 830 m away from the centre of the crater during the 2011 eruption at Shinmoedake, Japan. We found features suggesting the existence of period doubling bifurcation in the harmonic tremor signal, implying that the harmonic tremor might be generated by a nonlinear process. In order to quantify the nonlinearity in the harmonic tremor signal, we measure the correlation dimension 〈span〉D〈/span〉 and the maximal Lyapunov exponent λ. For one short but stable segment of the harmonic tremor seismogram, we obtained 〈span〉D〈/span〉 = 1.12 and λ = 0.03 s〈sup〉−1〈/sup〉. This result implies that the stable oscillation of the harmonic tremor is predominantly a limit cycle with small amounts of chaos present. We then use surrogate data analysis to check that our measurements of 〈span〉D〈/span〉 and λ do not include any false positive detection of nonlinearity. Limit cycles imply that the harmonic tremor is generated by self-sustained oscillations. We show that the autonomous Julian tremor model [〈span〉Julian〈/span〉, 1994] is able to exhibit period doubling bifurcation. We also show that the non-autonomous Julian tremor model with a step like increase and decrease in input pressure is able to exhibit oscillations of varying amplitude while keeping a constant frequency spectrum. Both phenomena were observed at Shinmoedake. We also demonstrate that the non-autonomous Julian tremor model with a transient input pressure is able to exhibit long period-like events in addition to harmonic tremor-like events, implying that the same nonlinear mechanism could be responsible for the generation of both type of events.〈/span〉

Description: 〈span〉〈div〉Summary〈/div〉We present a joint analysis of newly acquired gravity and teleseismic data in the North Tanzanian Divergence, where the lithospheric break-up is at its earliest stage. The impact of a mantle upwelling in more mature branches of the East African Rift has been extensively studied at a lithospheric scale. However, few studies have been completed that relate the deep-seated mantle anomaly detected in broad regional seismic tomography with the surface deformation observed in the thick Archaean Pan-African suture zone located in North Tanzania. Our joint inversion closes the gap between local and regional geophysical studies, providing velocity and density structures from the surface down to ca. 250 km depth with new details. Our results support the idea of a broad mantle upwelling rising up to the lithosphere and creating a thermal modification along its path. However, our study clearly presents an increasing amplitude of the associated anomaly both in velocity and density above 200 km depth, which cannot be solely explained by a temperature rise. We infer from our images the combined impact of melt (2-3%), composition and hydration that accompany the modification of a thick heterogenous cratonic lithosphere are a response to the hot mantle rising. The detailed images we obtained in density and velocity assert that Archaean and Proterozoic units interact with the mantle upwelling to restrict the lithosphere modifications within the Magadi-Natron-Manyara rift arm. The composition and hydration variations associated with those units equilibrate the thermal erosion of the craton root and allow for its stability between 100 and 200 km depth. Above 80 km depth, the crustal part is strongly affected by intruding bodies (melt and gas) which produces large negative anomalies in both velocity and density beneath the main magmatic centers. In addition to the global impact of a superplume, the velocity and density anomaly pattern suggests a 3D distribution of the crust and mantle lithospheric stretching, which is likely to be controlled by inherited fabrics and enhanced by lateral compositional and hydration variations at the Tanzanian craton-orogenic belt boundary.〈/span〉

Description: 〈span〉〈div〉Summary〈/div〉We report palaeomagnetic and K–Ar geochronologic results of two volcanic sequences from Ethiopia. The Belessa section, dated around 29–30 Ma and spanning ∼1 km in thickness, is related to the Oligocene Afro-Arabian traps, whereas the ∼700-m-thick Debre Sina section was emplaced during the Miocene in two periods around 10–11 and 14–15 Ma. We sampled 67 flows of predominantly basaltic rocks near Belessa and 59 rhyolitic to trachybasaltic flows near Debre Sina. From a geodynamic viewpoint, the magnetostratigraphy of the Belessa sequence confirms that the Ethiopian traps were emplaced at a minimum rate of ∼1 m/kyr, with a possible acceleration of the volume of volcanism over time. To provide insight into the evolution of the geomagnetic field in the Afro-Arabian region over the past 30 Myr, we combined our results with previous studies in the same area. Recentred directional distributions were elongated in the meridian plane, in coherence with field models for a dipole-dominated field. The dispersion 〈span〉S〈/span〉 of the virtual geomagnetic poles, representative of the vigour of palaeosecular variation, was approximately 50% higher during the 10–30 Ma interval than during the past 5 Myr. As the reversal frequency 〈span〉f〈/span〉 was more than two times lower during the Early Oligocene than during the Plio-Pleistocene, it appears that 〈span〉S〈/span〉 and 〈span〉f〈/span〉 are uncorrelated in this near-equatorial region. It remains an open question whether this apparent decoupling is ascribable to a local anomaly, is only sporadic in time, or represents a general feature of the geodynamo.〈/span〉

Description: 〈span〉〈div〉Summary〈/div〉The non-linear interaction of ocean surface waves produces coherent infrasound noise—microbaroms—between 0.1 and 0.5 Hz. Microbaroms propagate through the atmosphere over thousands of kilometres due to low absorption and efficient ducting between the ground and the stratopause. These signals are globally and permanently detected by the International Monitoring System (IMS) infrasound network, which has been established to monitor compliance with the Comprehensive Nuclear-Test-Ban Treaty. At the International Data Centre (IDC) in Vienna, where IMS data are routinely processed, microbarom detections appear in overlapping frequency bands, and are treated as false alarms. Therefore, understanding the variability in microbarom detections is essential to support the IDC in the reduction of the false alarm rate. In this study, microbarom amplitudes and the direction of arrivals at the German infrasound station IS26 were modelled. For the simulations, the source was described by an operational ocean wave interaction model, and the signal amplitude was modelled using a semi-empirical attenuation relation. This relation strongly depends on middle atmosphere (MA; i.e. 15–90 km altitude) dynamics; however, vertical temperature and wind profiles, provided by numerical weather prediction (NWP) models, exhibit significant biases and differences when compared with high-resolution lidar soundings in altitudes where infrasound signals propagate. To estimate uncertainties in the modelled amplitude, a fully autonomous lidar for MA temperature measurements was installed at IS26. Temperature and wind perturbations, considering observed biases and deviations, were added to the operational high-resolution atmospheric model analysis produced by the European Centre for Medium-Range Weather Forecasts. Such uncertainties in horizontal winds and temperature strongly impact propagation conditions, explaining almost 97 per cent of the actual detections, compared to 77 per cent when using the direct output of the NWP model only. Incorporating realistic wind and temperature uncertainties in NWP models can thus significantly improve the understanding of microbarom detections as well as the detection capability of a single station throughout the year.〈/span〉

Description: 〈span〉〈div〉Summary〈/div〉The interaction of subducted oceanic lithosphere with the discontinuities of the mantle transition zone (MTZ) provides insight into the composition and temperature of the subducted slab as well as potential melting of the slab or the surrounding mantle and loss of volatiles from the slab. Detailed mapping of the structure of the MTZ will help to better understand how slabs transport material and volatiles into the mantle and how phase transitions affect the slab dynamics. Here we use a dense network of seismic stations in northern Anatolia to image the structure of the MTZ discontinuities in detail using P-wave receiver functions. With a station spacing of about 7 km and a surface footprint of ∼35 km by ∼70 km, analysing receiver functions calculated from teleseismic earthquakes that occurred during an ∼18 month deployment produced clear images of where the mantle transition zone interacts with the Tethys/Cyprus slabs that either lie flat on the 660-km discontinuity or pass into the lower mantle. We observe an undulating 660-km discontinuity depressed by up to 30 km and a slightly depressed (1–2 km) 410-km discontinuity, apparently undisturbed by the slab. The MTZ is thickened to ∼270 km as result of the cool slab in the MTZ influencing the 660-km discontinuity and includes an arrival at ∼520-km depth likely from the top of a flat lying slab or a discontinuity related to a solid-solid phase transition in the olivine component of the mantle. We find evidence for low-velocity zones both above and below the 410-km discontinuity and above the 660-km discontinuity. The low velocity zones around the 410-km discontinuity might be the result of hydration of the MTZ from the slab and upward convection of MTZ material into the upper mantle. The origin of the low velocity zone around the 660-km discontinuity is less clear and could be related to sedimentation of subducted mid-ocean ridge basalts. The small footprint of the seismic array provides accurate information on the structure of the MTZ in an area influenced by subduction and shows small-scale changes in MTZ structure that might be lost in studies covering larger areas with sparser sampling.〈/span〉

Description: 〈span〉〈div〉SUMMARY〈/div〉Long-period Rayleigh wave horizontal to vertical amplitude (H/V) ratios at a station provide information about local earth structure that is complementary to phase velocity. However, a number of studies have observed that significant scatter appears in these measurements making it difficult to use H/V ratio measurements to resolve earth structure. Some of the scatter in these measurements has been attributed to local geological structure while some has remained unaccounted for. Most Global Seismographic Network (GSN) stations contain two nearby high-quality broad-band seismometers (e.g. in the same vault, but on different piers or in different boreholes). For each broad-band sensor in the IRIS/USGS component of the GSN, we estimate H/V ratios of fundamental mode Rayleigh waves using 〈span〉M〈/span〉 〉 6.5 earthquakes from 2001 to 2018 (around 19 000 measurements). We compute these ratios at a number of discrete periods (25, 50, 75, 100 and 150 s) and find that for well-isolated Rayleigh waves (windows where the correlation coefficients between radial and the phase-shifted vertical components are greater than 0.9) significant scatter in H/V ratios occurs between colocated sensors (greater than 25 per cent at 100 s period). This suggests the scatter in H/V ratio measurements can be at least partially attributed to extremely local phenomena such as sensor emplacement in the vault. We also find that H/V ratios can vary as a function of event backazimuth, indicating that care must be taken when computing average ratios for a station, as a large number of events from a given region could bias H/V ratio measurements at a station.〈/span〉

Description: In their recent review, Bailey et al. (2017) make the compelling case that indirect genetic effects (IGEs) should be more widely incorporated into behavioral ecology research. IGEs occur whenever the trait of one individual is affected by the genotype of another “interacting” individual, usually via one or more “social cues.” IGEs create an additional source of heritable variation that can influence an individual’s trait expression above and beyond any direct contributions of its genotype (direct genetic effect [DGE]) and the environment ( Lynch and Walsh 1998 ). When conspecifics interact, as is the case for behavioral interactions, they make up part of an individual’s “social environment” and because these interacting partners are themselves genetically variable, additional heritable variation can be generated. So why does this matter? A key outcome of IGEs is that traits affected by them can evolve not only due to direct genetic variation but also due to evolution of the social environment itself ( Wolf et al. 1998 ).

Description: Behavior is rapidly flexible and highly context-dependent, which poses obvious challenges to researchers attempting to dissect its causes. However, over a century of unresolved debate has also focused on whether the very flexibility and context-dependence of behavior lends it a unique role in the evolutionary origins and patterns of diversity in the Animal Kingdom. Here, we propose that both challenges can benefit from studying how indirect genetic effects (IGEs: the effects of genes expressed in one individual on traits in another individual) shape behavioral phenotypes. We provide a sketch of the theoretical framework that grounds IGEs in behavioral ecology research and focus on recent advances made from studies of IGEs in areas of behavioral ecology such as sexual selection, sexual conflict, social dominance, and parent–offspring interactions. There is mounting evidence that IGEs have important influences on behavioral phenotypes associated with these processes, such as sexual signals and preferences and behaviors which function to manipulate interacting partners. IGEs can also influence both responses to selection and selection itself, and considering IGEs refines evolutionary predictions and provides new perspectives on the origins of seemingly perplexing behavioral traits. A key unresolved question, but one that has dominated the behavioral sciences for over a century, is whether behavior is more likely than other types of traits to contribute to evolutionary change and diversification. We advocate taking advantage of an IGE approach to outline falsifiable hypotheses and a general methodology to rigorously test this frequently proposed, yet still contentious, special role of behavior in evolution.

Description: The notion that behavioral traits “must often take the lead in evolution” because they are “especially plastic” ( West-Eberhard 2003 ) sits somewhat uncomfortably with standard quantitative genetic theory. Under simple models, plasticity actually slows evolution: more plastic traits have greater variance, reduced heritability and so lowered expected selection response. Here, Bailey et al. (2017) highlight the tantalizing possibility that indirect genetic effects (IGEs) could reconcile these contrasting positions. They provide a neat structured list of questions to guide empirical efforts, effectively asking:

Description: The key proposal in Bailey et al. (2018) is that evolutionary processes may be misinterpreted if indirect genetic effects (IGEs) are not taken into consideration. Further, that such processes are a fundamental part of behavior and thus the genetic analysis of behavior and its evolution must take into account IGEs. The authors present a very convincing case for this proposal both from theoretical argument and empirical evidence.

Description: There is growing evidence for consistent among-individual variation in individual sociability (e.g., tendency to be sociable) in a number of species. However, sexes often differ in their social behaviors, as well as the selection pressures which they experience. This may translate into differences in repeatability of sociability, although this has not yet been tested. Here, we investigated whether eastern water dragons ( Intellegama leseurii ) exhibited evidence of consistent among-individual variation (i.e., repeatability) in 4 different measurements of sociability. Specifically, we measured sociability in 4 ways (degree, centrality, proportion of time spent being social, and number of preferences), and tested whether there was evidence for sex differences in the repeatability of these sociability measurements, or whether observed levels of repeatability could be explained by a stable social environment. Our findings provide new evidence for sex differences in social personality: we found that males were significantly repeatable in 3 of 4 sociability measurements (degree, centrality, and proportion of time spent being social), whereas females were not. Further, we found that these differences were not a result of differences in the dynamics of the social environment. We discuss our findings in the context of sexual selection, as well as sex differences in the evolutionary drivers of social behavior.

Description: When foraging in their natural environment, many animals readily complement their personal knowledge with additional social information. To balance the costs and benefits of copying others, animals have to discern situations in which it is more advantageous to use social rather than personal information. Here, we used foraging bumblebees ( Bombus terrestris ) in a controlled laboratory setting and showed that the difficulty of a foraging task affects how the bees weight the 2 types of information. We used artificial flowers to devise easy and difficult discriminatory tasks, and tested the influence of floral and social cues on decision making. When facing an easy discrimination task, foraging bees were likely to rely on personal information and were only marginally affected by social information. By contrast, they prioritized social over personal information when flower discrimination was difficult and therefore the probability of making errors was higher. In summary, bees are able to use social and personal information to optimize foraging decisions in a flexible way.

Description: We propose that indirect genetic effects (IGEs) represent an appealing way to dissect the genetics and evolutionary dynamics of complex phenotypes studied in behavioral ecology. IGEs are complementary to inclusive fitness approaches, but as Kruuk and Wilson (2018) observe, are distinct because they do not require assumptions about relatedness (see also McGlothlin et al. 2014 and McDonald et al. 2017 ). Thus, IGEs open up the opportunity to study any social interaction and any trait expressed in a social interaction (See Roff 2018 ). These features of IGEs have enabled behavioral ecologists to study behavior from a slightly different perspective, as an evolutionary feedback process that links individuals and their environments through genes.

Description: By definition, teaching is both cooperative and costly and so we appreciate Cordony’s (2017) commentary in this spirit. We also agree that the cost of teaching in the superb fairy-wren ( Malurus cyaneus ) system can and should be more fully explored. Kleindorfer, Hoi, et al. (2014) quantified a cost of incubation calling ( Colombelli-Négrel et al. 2012 ) to fairy-wren mothers through the detection of higher egg predation at nests with more incubation calls. The commentary raises 2 main concerns: that 1) Kleindorfer, Hoi, et al. (2014) misinterpreted the operational definition of cost in animal teaching by Caro and Hauser (1992) and its operational interpretation by Thornton and Raihani (2010) , because 2) nest predation is neither an immediate nor a personal cost to the mother.