ALBERT

Description: Secondary microseismic noise is generated by non-linear interactions between ocean waves at the ocean surface. We present here the theory for computing the site effect of the ocean layer upon body waves generated by noise sources distributed along the ocean surface. By defining the wavefield as the superposition of plane waves, we show that the ocean site effect can be described as the constructive interference of multiply reflected P waves in the ocean that are then converted to either P or SV waves at the ocean–crust interface. We observe that the site effect varies strongly with period and ocean depth, although in a different way for body waves than for Rayleigh waves. We also show that the ocean site effect is stronger for P waves than for S waves. We validate our computation by comparing the theoretical noise body wave sources with the sources inferred from beamforming analysis of the three seismogram components recorded by the Southern California Seismic Network. We use rotated traces for the beamforming analysis, and we show that we clearly detect P waves generated by ocean gravity wave interactions along the track of typhoon Ioke (2006 September). We do not detect the corresponding SV waves, and we demonstrate that this is because their amplitude is too weak.

Description: Published

Description: 1096-1106

Description: 1T. Geodinamica e interno della Terra

Description: JCR Journal

Description: restricted

Description: "With so much media and political criticism of their shortcomings and failures, it is easy to overlook the fact that many governments work pretty well much of the time. Great Policy Successes turns the spotlight on instances of public policy that are remarkably successful. It develops a framework for identifying and assessing policy successes, paying attention not just to their programmatic outcomes but also to the quality of the processes by which policies are designed and delivered, the level of support and legitimacy they attain, and the extent to which successful performance endures over time. The bulk of the book is then devoted to 15 detailed case studies of striking policy successes from around the world, including Singapore's public health system, Copenhagen and Melbourne's rise from stilted backwaters to the highly liveable and dynamic urban centres they are today, Brazil's Bolsa Familia poverty relief scheme, the US's GI Bill, and Germany's breakthrough labour market reforms of the 2000s. Each case is set in context, its main actors are introduced, key events and decisions are described, the assessment framework is applied to gauge the nature and level of its success, key contributing factors to success are identified, and potential lessons and future challenges are identified. Purposefully avoiding the kind of heavy theorizing that characterizes many accounts of public policy processes, each case is written in an accessible and narrative style ideally suited for classroom use in conjunction with mainstream textbooks on public policy design, implementation, and evaluation.

Description: A major uncertainty in determining the mass balance of the Antarctic ice sheet from measurements of satellite gravimetry, and to a lesser extent satellite altimetry, is the poorly known correction for the ongoing deformation of the solid Earth caused by glacial isostatic adjustment (GIA). Although much progress has been made in consistently modeling the ice-sheet evolution throughout the last glacial cycle, as well as the induced bedrock deformation caused by these load changes, forward models of GIA remain ambiguous due to the lack of observational constraints on the ice sheet's past extent and thickness and mantle rheology beneath the continent. As an alternative to forward-modeling GIA, we estimate GIA from multiple space-geodetic observations: Gravity Recovery and Climate Experiment (GRACE), Envisat/ICESat and Global Positioning System (GPS). Making use of the different sensitivities of the respective satellite observations to current and past surface-mass (ice mass) change and solid Earth processes, we estimate GIA based on viscoelastic response functions to disc load forcing. We calculate and distribute the viscoelastic response functions according to estimates of the variability of lithosphere thickness and mantle viscosity in Antarctica. We compare our GIA estimate with published GIA corrections and evaluate its impact in determining the ice-mass balance in Antarctica from GRACE and satellite altimetry. Particular focus is applied to the Amundsen Sea Sector in West Antarctica, where uplift rates of several centimetres per year have been measured by GPS. We show that most of this uplift is caused by the rapid viscoelastic response to recent ice-load changes, enabled by the presence of a low-viscosity upper mantle in West Antarctica. This paper presents the second and final contributions summarizing the work carried out within a European Space Agency funded study, REGINA (www.regina-science.eu).

Description: © The Author(s), 2015. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Bioinformatics 31 (2015): 1872-1874, doi:10.1093/bioinformatics/btv045.

Description: The association of organisms to their environments is a key issue in exploring biodiversity patterns. This knowledge has traditionally been scattered, but textual descriptions of taxa and their habitats are now being consolidated in centralized resources. However, structured annotations are needed to facilitate large-scale analyses. Therefore, we developed ENVIRONMENTS, a fast dictionary-based tagger capable of identifying Environment Ontology (ENVO) terms in text. We evaluate the accuracy of the tagger on a new manually curated corpus of 600 Encyclopedia of Life (EOL) species pages. We use the tagger to associate taxa with environments by tagging EOL text content monthly, and integrate the results into the EOL to disseminate them to a broad audience of users.

Description: The Encyclopedia Of Life Rubenstein Fellows Program [CRDF EOL-33066-13/E33066], the LifeWatchGreece Research Infrastructure [384676-94/GSRT/ NSRF(C&E)] and the Novo Nordisk Foundation Center for Protein Research [NNF14CC0001].

Description: Author Posting. © The Authors, 2019. 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 218(3), (2019): 1822-1837, doi: 10.1093/gji/ggz253.

Description: Joint inversion of multiple electromagnetic data sets, such as controlled source electromagnetic and magnetotelluric data, has the potential to significantly reduce uncertainty in the inverted electrical resistivity when the two data sets contain complementary information about the subsurface. However, evaluating quantitatively the model uncertainty reduction is made difficult by the fact that conventional inversion methods—using gradients and model regularization—typically produce just one model, with no associated estimate of model parameter uncertainty. Bayesian inverse methods can provide quantitative estimates of inverted model parameter uncertainty by generating an ensemble of models, sampled proportional to data fit. The resulting posterior distribution represents a combination of a priori assumptions about the model parameters and information contained in field data. Bayesian inversion is therefore able to quantify the impact of jointly inverting multiple data sets by using the statistical information contained in the posterior distribution. We illustrate, for synthetic data generated from a simple 1-D model, the shape of parameter space compatible with controlled source electromagnetic and magnetotelluric data, separately and jointly. We also demonstrate that when data sets contain complementary information about the model, the region of parameter space compatible with the joint data set is less than or equal to the intersection of the regions compatible with the individual data sets. We adapt a trans-dimensional Markov chain Monte Carlo algorithm for jointly inverting multiple electromagnetic data sets for 1-D earth models and apply it to surface-towed controlled source electromagnetic and magnetotelluric data collected offshore New Jersey, USA, to evaluate the extent of a low salinity aquifer within the continental shelf. Our inversion results identify a region of high resistivity of varying depth and thickness in the upper 500 m of the continental shelf, corroborating results from a previous study that used regularized, gradient-based inversion methods. We evaluate the joint model parameter uncertainty in comparison to the uncertainty obtained from the individual data sets and demonstrate quantitatively that joint inversion offers reduced uncertainty. In addition, we show how the Bayesian model ensemble can subsequently be used to derive uncertainty estimates of pore water salinity within the low salinity aquifer.

Description: We gratefully acknowledge funding support from National Science Foundation grants 1458392 and 1459035. We thank the captain and crew of the R.V. Marcus G. Langseth for a successful cruise and the Marine EM Lab at Scripps Institution of Oceanography for providing the instrumentation. We also thank Chris Armerding, Marah Dahn, John Desanto, Jimmy Elsenbeck, Matt Folsom, Keiichi Ishizu, Jeff Pepin, Charlotte Wiman and Georgie Zelenak for participating in the cruise. We gratefully acknowledge Alberto Malinverno for the idea to use a Monte Carlo scheme to estimate the distribution of pore fluid salinity, and William Menke for many constructive conversations and suggestions.

Description: Author Posting. © The Authors, 2019. 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 218(3), (2019): 2122-2135, doi: 10.1093/gji/ggz272.

Description: We have conducted the first passive Ocean Bottom Seismograph (OBS) experiment near the Challenger Deep at the southernmost Mariana subduction zone by deploying and recovering an array of 6 broad-band OBSs during December 2016–June 2017. The obtained passive-source seismic records provide the first-ever near-field seismic observations in the southernmost Mariana subduction zone. We first correct clock errors of the OBS recordings based on both teleseismic waveforms and ambient noise cross-correlation. We then perform matched filter earthquake detection using 53 template events in the catalogue of the US Geological Survey and find 〉7000 local earthquakes during the 6-month OBS deployment period. Results of the two independent approaches show that the maximum clock drifting was ∼2 s on one instrument (OBS PA01), while the rest of OBS waveforms had negligible time drifting. After timing correction, we locate the detected earthquakes using a newly refined local velocity model that was derived from a companion active source experiment in the same region. In total, 2004 earthquakes are located with relatively high resolution. Furthermore, we calibrate the magnitudes of the detected earthquakes by measuring the relative amplitudes to their nearest relocated templates on all OBSs and acquire a high-resolution local earthquake catalogue. The magnitudes of earthquakes in our new catalogue range from 1.1 to 5.6. The earthquakes span over the Southwest Mariana rift, the megathrust interface, forearc and outer-rise regions. While most earthquakes are shallow, depths of the slab earthquakes increase from ∼100 to ∼240 km from west to east towards Guam. We also delineate the subducting interface from seismicity distribution and find an increasing trend in dip angles from west to east. The observed along-strike variation in slab dip angles and its downdip extents provide new constraints on geodynamic processes of the southernmost Mariana subduction zone.

Description: We express our appreciation to the science parties and crew members of the R/V Shiyan 3 for deployment and collection of the OBS instruments during the Mariana expeditions. This study is supported by the Hong Kong Research Grant Council Grants (No. 14313816), Faculty of Science at CUHK, Chinese Academy of Sciences (No. Y4SL021001, QYZDY-SSW-DQC005, 133244KYSB20180029), the National Natural Science Foundation of China (No. 41890813, 91628301, 41676042, U1701641, 41576041, 91858207 and U1606401), the National Key R&D Program of China (2018YFC0309800 and 2018YFC0310100). Generic Mapping Tools (Wessel & Smith 1991) and PSSAC (developed by Prof Lupei Zhu) are used for data analysis and figure preparation in this study. Constructive comments from Dr Lidong Bie and two anonymous reviewers are helpful in improving the manuscript.

Description: Author Posting. © Oxford University Press, 2016. 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 205 (2016): 728-743, doi:10.1093/gji/ggw044.

Description: While elasticity is a defining characteristic of the Earth's lithosphere, it is often ignored in numerical models of long-term tectonic processes in favour of a simpler viscoplastic description. Here we assess the consequences of this assumption on a well-studied geodynamic problem: the growth of normal faults at an extensional plate boundary. We conduct 2-D numerical simulations of extension in elastoplastic and viscoplastic layers using a finite difference, particle-in-cell numerical approach. Our models simulate a range of faulted layer thicknesses and extension rates, allowing us to quantify the role of elasticity on three key observables: fault-induced topography, fault rotation, and fault life span. In agreement with earlier studies, simulations carried out in elastoplastic layers produce rate-independent lithospheric flexure accompanied by rapid fault rotation and an inverse relationship between fault life span and faulted layer thickness. By contrast, models carried out with a viscoplastic lithosphere produce results that may qualitatively resemble the elastoplastic case, but depend strongly on the product of extension rate and layer viscosity U × ηL. When this product is high, fault growth initially generates little deformation of the footwall and hanging wall blocks, resulting in unrealistic, rigid block-offset in topography across the fault. This configuration progressively transitions into a regime where topographic decay associated with flexure is fully accommodated within the numerical domain. In addition, high U × ηL favours the sequential growth of multiple short-offset faults as opposed to a large-offset detachment. We interpret these results by comparing them to an analytical model for the fault-induced flexure of a thin viscous plate. The key to understanding the viscoplastic model results lies in the rate-dependence of the flexural wavelength of a viscous plate, and the strain rate dependence of the force increase associated with footwall and hanging wall bending. This behaviour produces unrealistic deformation patterns that can hinder the geological relevance of long-term rifting models that assume a viscoplastic rheology.

Description: This work was supported by NSF grants OCE-11-54238 (JAO, MDB), EAR-10-10432 (MDB) and OCE-11-55098 (GI), as well as a WHOI Deep Exploration Institute grant and start-up support from the University of Idaho (EM).

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

Description: Reproduction of mysticete whales is difficult to monitor, and basic parameters, such as pregnancy rate and inter-calving interval, remain unknown for many populations. We hypothesized that baleen plates (keratinous strips that grow downward from the palate of mysticete whales) might record previous pregnancies, in the form of high-progesterone regions in the sections of baleen that grew while the whale was pregnant. To test this hypothesis, longitudinal baleen progesterone profiles from two adult female North Atlantic right whales (Eubalaena glacialis) that died as a result of ship strike were compared with dates of known pregnancies inferred from calf sightings and post-mortem data. We sampled a full-length baleen plate from each female at 4 cm intervals from base (newest baleen) to tip (oldest baleen), each interval representing ∼60 days of baleen growth, with high-progesterone areas then sampled at 2 or 1 cm intervals. Pulverized baleen powder was assayed for progesterone using enzyme immunoassay. The date of growth of each sampling location on the baleen plate was estimated based on the distance from the base of the plate and baleen growth rates derived from annual cycles of stable isotope ratios. Baleen progesterone profiles from both whales showed dramatic elevations (two orders of magnitude higher than baseline) in areas corresponding to known pregnancies. Baleen hormone analysis shows great potential for estimation of recent reproductive history, inter-calving interval and general reproductive biology in this species and, possibly, in other mysticete whales.

Description: This work was supported by the Eppley Foundation for Research, the National Oceanographic and Atmospheric Administration Marine Mammal Health and Stranding Program and the Woods Hole Oceanographic Institution Ocean Life Institute.

Description: © The Author(s), 2014. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Genome Biology and Evolution 7 (2015): 457-464, doi:10.1093/gbe/evu289.

Description: The uptake and transport of vitamin B12 (cobalamin; Cbl) in mammals involves a refined system with three evolutionarily related transporters: transcobalamin 1 (Tcn1), transcobalamin 2 (Tcn2), and the gastric intrinsic factor (Gif). Teleosts have a single documented binder with intermediate features to the human counterparts. Consequently, it has been proposed that the expansion of Cbl binders occurred after the separation of Actinopterygians. Here, we demonstrate that the diversification of this gene family took place earlier in gnathostome ancestry. Our data indicates the presence of single copy orthologs of the Sarcopterygii/Tetrapoda duplicates Tcn1 and Gif, and Tcn2, in Chondrichthyes. In addition, a highly divergent Cbl binder was found in the Elasmobranchii. We unveil a complex scenario forged by genome, tandem duplications and lineage-specific gene loss. Our findings suggest that from an ancestral transporter, exhibiting large spectrum and high affinity binding, highly specific Cbl transporters emerged through gene duplication and mutations at the binding pocket.

Description: This work was funded by the Fundação para a Ciência e a Tecnologia (FCT) project PEst-C/MAR/LA0015/2013. PhD grant SFRH/BD/84238/2012 awarded to M.L-.M. Postdoctoral grant SFRH/BPD/72519/2010 awarded to R.R.

Description: Author Posting. © The Author(s), 2016. 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 205 (2016): 785-795, doi:10.1093/gji/ggw036.

Description: An L-configured, three-component short period seismic array was deployed on the Ross Ice Shelf, Antarctica during November 2014. Polarization analysis of ambient noise data from these stations shows linearly polarized waves for frequency bands between 0.2 and 2 Hz. A spectral peak at about 1.6 Hz is interpreted as the resonance frequency of the water column and is used to estimate the water layer thickness below the ice shelf. The frequency band from 4 to 18 Hz is dominated by Rayleigh and Love waves propagating from the north that, based on daily temporal variations, we conclude were generated by field camp activity. Frequency–slowness plots were calculated using beamforming. Resulting Love and Rayleigh wave dispersion curves were inverted for the shear wave velocity profile within the firn and ice to ∼150 m depth. The derived density profile allows estimation of the pore close-off depth and the firn–air content thickness. Separate inversions of Rayleigh and Love wave dispersion curves give different shear wave velocity profiles within the firn. We attribute this difference to an effective anisotropy due to fine layering. The layered structure of firn, ice, water and the seafloor results in a characteristic dispersion curve below 7 Hz. Forward modelling the observed Rayleigh wave dispersion curves using representative firn, ice, water and sediment structures indicates that Rayleigh waves are observed when wavelengths are long enough to span the distance from the ice shelf surface to the seafloor. The forward modelling shows that analysis of seismic data from an ice shelf provides the possibility of resolving ice shelf thickness, water column thickness and the physical properties of the ice shelf and underlying seafloor using passive-source seismic data.

Description: PDB, AD and PG were supported by NSF Grant PLR 1246151. RAS was supported by NSF Grant PLR-1246416. DAW, RA and AN were supported under NSF Grants PLR-1142518, 1141916 and 1142126, respectively. PDB also received support from the California Department of Parks and Recreation, Division of Boating and Waterways under contract 11-106-107.

Description: Author Posting. © The Author(s), 2019. This is the author's version of the work. It is posted here by permission of Oxford University Press for personal use, not for redistribution. The definitive version was published in ICES Journal of Marine Science 76(4), (2010): 781-786, doi:10.1093/icesjms/fsy194.

Description: Whales are federally protected by the Marine Mammal Protection Act; endangered species, such as the North Atlantic right whale, receive additional protection under the Endangered Species Act. However, their regulations have failed to satisfy conservation and animal welfare concerns. From 1990 to 2011 the North Atlantic right whale (Eubalaena glacialis, NARW) population grew at a mean of 2.8% annually. However, population trends reversed since 2011; the species is in decline, with only ∼100 reproductively active females remaining. This failure is driven by vessel collisions and increasingly fatal and serious entanglement in fixed fishing gear, whose rope strength has increased substantially. Chronic entanglement, drag, and associated morbidity have been linked to poor fecundity. Genuine solutions involve designating areas to be avoided and speed restrictions for ships and removing fishing trap ropes from the water column. A trap fishing closure for NARW habitat in the Cape Cod Bay (U.S.) area has been in place seasonally since 2015. 2017 mortalities in Eastern Canada elicited substantive management changes whereby the 2018 presence of NARW in active trap fishing areas resulted in an effective closure. To avoid these costly closures, the traditional trap fishery model of rope end lines attached to surface marker buoys has to be modified so that traps are marked virtually, and retrieved with gear that does not remain in the water column except during trap retrieval. Consumer demand for genuinely whale-safe products will augment and encourage the necessary regulatory changes so that trap fisheries conserve target and nontarget species.

Description: I thank Mark Baumgartner, Scott Kraus, Tim Werner, Amy Knowlton, Heather Pettis, Scott Landry, Stormy Mayo, Fred Penney, and Beth Casoni for discussions on this topic and Natalie Renier for drawing Figures 5 and 6. Funding was provided by the Woods Hole Oceanographic Institution Independent Research and Development Program.

Description: 2020-01-10

Description: Author Posting. © The Author(s), 2019. This is the author's version of the work. It is posted here by permission of Oxford University Press for personal use, not for redistribution. The definitive version was published in Geophysical Journal International, 219(1), (2019): 464-478, doi:10.1093/gji/ggz315.

Description: The electromagnetic (EM) field generated by ocean tidal flow is readily detectable in both satellite magnetic field data, and in ocean-bottom measurements of electric and magnetic fields. The availability of accurate charts of tidal currents, constrained by assimilation of modern satellite altimetry data, opens the possibility of using tidal EM fields as a source to image mantle electrical resistivity beneath the ocean basins, as highlighted by the recent success in defining the globally averaged lithosphere–asthenosphere boundary (LAB) with satellite data. In fact, seafloor EM data would be expected to provide better constraints on the structure of resistive oceanic lithosphere, since the toroidal magnetic mode, which can constrain resistive features, is a significant component of the tidal EM field within the ocean, but is absent above the surface (in particular in satellite data). Here we consider this issue in more detail, using a combination of simplified theoretical analysis and 1-D and 3-D numerical modelling to provide a thorough discussion of the sensitivity of satellite and seafloor data to subsurface electrical structure. As part of this effort, and as a step toward 3-D inversion of seafloor tidal data, we have developed a new flexible 3-D spherical-coordinate finite difference scheme for both global and regional scale modelling, with higher resolution models nested in larger scale solutions. We use the new 3-D model, together with Monte Carlo simulations of errors in tidal current estimates, to provide a quantitative assessment of errors in the computed tidal EM signal caused by uncertainty in the tidal source. Over the open ocean this component of error is below 0.01 nT in Bz at satellite height and 0.05 nT in Bx on the seafloor, well below typical signal levels. However, as coastlines are approached error levels can increase substantially. Both analytical and 3-D modelling demonstrate that the seafloor magnetic field is most sensitive to the lithospheric resistance (the product of resistivity and thickness), and is more weakly influenced (primarily in the phase) by resistivity of the underlying asthenosphere. Satellite data, which contain only the poloidal magnetic mode, are more sensitive to the conductive asthenosphere, but have little sensitivity to lithospheric resistance. For both seafloor and satellite data’s changes due to plausible variations in Earth parameters are well above error levels associated with source uncertainty, at least in the ocean interior. Although the 3-D modelling results are qualitatively consistent with theoretical analysis, the presence of coastlines and bathymetric variations generates a complex response, confirming that quantitative interpretation of ocean tidal EM fields will require a 3-D treatment. As an illustration of the nested 3-D scheme, seafloor data at five magnetic and seven electric stations in the northeastern Pacific (41○N, 165○W) are fit with trial-and-error forward modelling of a local domain. The simulation results indicate that the lithospheric resistance is roughly 7 × 108 Ωm2. The phase of the seafloor data in this region are inconsistent with a sharp transition between the resistive lithosphere and conductive asthenosphere.

Description: This work has been supported by National Natural Science Foundation of China grants 41804072 and 41574104, and NSF grant EAR-1447109. Special thanks to Dr Benjamin Murphy who provided the conductivity-depth profile for 1-D earth model, Dr Min Ding who provided valuable discussion about the oceanic lithosphere and Dr Jeffery Love who provided comments on the stylistics of the manuscript.

Description: © The Author(s), 2017. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Genome Biology and Evolution 9 (2017): 659-676, doi:10.1093/gbe/evx023.

Description: Understanding and predicting the fate of populations in changing environments require knowledge about the mechanisms that support phenotypic plasticity and the adaptive value and evolutionary fate of genetic variation within populations. Atlantic killifish (Fundulus heteroclitus) exhibit extensive phenotypic plasticity that supports large population sizes in highly fluctuating estuarine environments. Populations have also evolved diverse local adaptations. To yield insights into the genomic variation that supports their adaptability, we sequenced a reference genome and 48 additional whole genomes from a wild population. Evolution of genes associated with cell cycle regulation and apoptosis is accelerated along the killifish lineage, which is likely tied to adaptations for life in highly variable estuarine environments. Genome-wide standing genetic variation, including nucleotide diversity and copy number variation, is extremely high. The highest diversity genes are those associated with immune function and olfaction, whereas genes under greatest evolutionary constraint are those associated with neurological, developmental, and cytoskeletal functions. Reduced genetic variation is detected for tight junction proteins, which in killifish regulate paracellular permeability that supports their extreme physiological flexibility. Low-diversity genes engage in more regulatory interactions than high-diversity genes, consistent with the influence of pleiotropic constraint on molecular evolution. High genetic variation is crucial for continued persistence of species given the pace of contemporary environmental change. Killifish populations harbor among the highest levels of nucleotide diversity yet reported for a vertebrate species, and thus may serve as a useful model system for studying evolutionary potential in variable and changing environments.

Description: This work was primarily supported by a grant from the National Science Foundation (collaborative research grants DEB-1265282, DEB-1120512, DEB-1120013, DEB-1120263, DEB-1120333, DEB-1120398 to J.K.C., D.L.C., M.E.H., S.I.K., M.F.O., J.R.S., W.W., and A.W.). Further support was provided by the National Institute of Environmental Health Sciences (1R01ES021934-01 to A.W., P42ES7373 to T.H.H., P42ES007381 to M.E.H., and R01ES019324 to J.R.S.), the National Institute of General Medical Sciences (P20GM103423 and P20GM104318 to B.L.K.), and the National Science Foundation (DBI-0640462 and XSEDE-MCB100147 to D.G.).

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

Description: Recent studies have demonstrated that some hormones are present in baleen powder from bowhead (Balaena mysticetus) and North Atlantic right (Eubalaena glacialis) whales. To test the potential generalizability of this technique for studies of stress and reproduction in large whales, we sought to determine whether all major classes of steroid and thyroid hormones are detectable in baleen, and whether these hormones are detectable in other mysticetes. Powdered baleen samples were recovered from single specimens of North Atlantic right, bowhead, blue (Balaenoptera [B.]musculus), sei (B. borealis), minke (B. acutorostrata), fin (B. physalus), humpback (Megaptera novaeangliae) and gray (Eschrichtius robustus) whales. Hormones were extracted with a methanol vortex method, after which we tested all species with commercial enzyme immunoassays (EIAs, Arbor Assays) for progesterone, testosterone, 17β-estradiol, cortisol, corticosterone, aldosterone, thyroxine and tri-iodothyronine, representing a wide array of steroid and thyroid hormones of interest for whale physiology research. In total, 64 parallelism tests (8 species × 8 hormones) were evaluated to verify good binding affinity of the assay antibodies to hormones in baleen. We also tested assay accuracy, although available sample volume limited this test to progesterone, testosterone and cortisol. All tested hormones were detectable in baleen powder of all species, and all assays passed parallelism and accuracy tests. Although only single individuals were tested, the consistent detectability of all hormones in all species indicates that baleen hormone analysis is likely applicable to a broad range of mysticetes, and that the EIA kits tested here perform well with baleen extract. Quantification of hormones in baleen may be a suitable technique with which to explore questions that have historically been difficult to address in large whales, including pregnancy and inter-calving interval, age of sexual maturation, timing and duration of seasonal reproductive cycles, adrenal physiology and metabolic rate.

Description: This work was supported by (1) the Center for Bioengineering Innovation at Northern Arizona University and (2) the New England Aquarium.

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), 2019. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Aslan, C., Beckman, N. G., Rogers, H. S., Bronstein, J., Zurell, D., Hartig, F., Shea, K., Pejchar, L., Neubert, M., Poulsen, J., HilleRisLambers, J., Miriti, M., Loiselle, B., Effiom, E., Zambrano, J., Schupp, G., Pufal, G., Johnson, J., Bullock, J. M., Brodie, J., Bruna, E., Cantrell, R. S., Decker, R., Fricke, E., Gurski, K., Hastings, A., Kogan, O., Razafindratsima, O., Sandor, M., Schreiber, S., Snell, R., Strickland, C., & Zhou, Y. Employing plant functional groups to advance seed dispersal ecology and conservation. AoB Plants, 11(2), (2019):plz006, doi:10.1093/aobpla/plz006.

Description: Seed dispersal enables plants to reach hospitable germination sites and escape natural enemies. Understanding when and how much seed dispersal matters to plant fitness is critical for understanding plant population and community dynamics. At the same time, the complexity of factors that determine if a seed will be successfully dispersed and subsequently develop into a reproductive plant is daunting. Quantifying all factors that may influence seed dispersal effectiveness for any potential seed-vector relationship would require an unrealistically large amount of time, materials and financial resources. On the other hand, being able to make dispersal predictions is critical for predicting whether single species and entire ecosystems will be resilient to global change. Building on current frameworks, we here posit that seed dispersal ecology should adopt plant functional groups as analytical units to reduce this complexity to manageable levels. Functional groups can be used to distinguish, for their constituent species, whether it matters (i) if seeds are dispersed, (ii) into what context they are dispersed and (iii) what vectors disperse them. To avoid overgeneralization, we propose that the utility of these functional groups may be assessed by generating predictions based on the groups and then testing those predictions against species-specific data. We suggest that data collection and analysis can then be guided by robust functional group definitions. Generalizing across similar species in this way could help us to better understand the population and community dynamics of plants and tackle the complexity of seed dispersal as well as its disruption.

Description: Ideas for this manuscript initiated during the Seed Dispersal Workshop held in May 2016 at the Socio-Environmental Synthesis Center in Annapolis, MD and supported by the US National Science Foundation Grant DEB-1548194 to N.G.B. and the National Socio‐Environmental Synthesis Center under the US National Science Foundation Grant DBI‐1052875. D.Z. received funding from the Swiss National Science Foundation (SNF, grant: PZ00P3_168136/1) and from the German Science Foundation (DFG, grant: ZU 361/1- 1). Contributions by the authors C.A. led the development of the concepts, writing, and revising of the manuscript with input from N.G.B. and H.S.R. All authors contributed to the development of concepts and are listed in order of contribution and alphabetical order within each level of contribution.

Description: © The Author(s), 2019. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Alexander, H., Johnson, L. K., & Brown, C. T.. Keeping it light: (re)analyzing community-wide datasets without major infrastructure. Gigascience, 8(2),(2019): giy159, doi:10.1093/gigascience/giy159.

Description: DNA sequencing technology has revolutionized the field of biology, shifting biology from a data-limited to data-rich state. Central to the interpretation of sequencing data are the computational tools and approaches that convert raw data into biologically meaningful information. Both the tools and the generation of data are actively evolving, yet the practice of re-analysis of previously generated data with new tools is not commonplace. Re-analysis of existing data provides an affordable means of generating new information and will likely become more routine within biology, yet necessitates a new set of considerations for best practices and resource development. Here, we discuss several practices that we believe to be broadly applicable when re-analyzing data, especially when done by small research groups.

Description: Funding was provided by the Gordon and Betty Moore Foundation (award GBMF4551 to C.T.B.).

Description: © International Council for the Exploration of the Sea, 2016. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in ICES Journal of Marine Science 73 (2016): 1839-1850, doi: 10.1093/icesjms/fsw086.

Description: For terrestrial and marine benthic ecologists, landscape ecology provides a framework to address issues of complexity, patchiness, and scale—providing theory and context for ecosystem based management in a changing climate. Marine pelagic ecosystems are likewise changing in response to warming, changing chemistry, and resource exploitation. However, unlike spatial landscapes that migrate slowly with time, pelagic seascapes are embedded in a turbulent, advective ocean. Adaptations from landscape ecology to marine pelagic ecosystem management must consider the nature and scale of biophysical interactions associated with organisms ranging from microbes to whales, a hierarchical organization shaped by physical processes, and our limited capacity to observe and monitor these phenomena across global oceans. High frequency, multiscale, and synoptic characterization of the 4-D variability of seascapes are now available through improved classification methods, a maturing array of satellite remote sensing products, advances in autonomous sampling of multiple levels of biological complexity, and emergence of observational networks. Merging of oceanographic and ecological paradigms will be necessary to observe, manage, and conserve species embedded in a dynamic seascape mosaic, where the boundaries, extent, and location of features change with time.

Description: This work was supported by NASA grant NNX14AP62A “National Marine Sanctuaries as Sentinel Sites for a Demonstration Marine Biodiversity Observation Network (MBON)” funded under the National Ocean Partnership Program (NOPP RFP NOAA-NOS-IOOS-2014-2003803 in partnership between NOAA, BOEM, and NASA), the NOAA Integrated Ocean Observing System (IOOS) Program Office, and the LenFest Ocean Program.

Description: Author Posting. © The Authors, 2019. 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 218(3), (2019): 2165-2178, doi: 10.1093/gji/ggz280.

Description: A multitaper estimator is proposed that accommodates time-series containing gaps without using any form of interpolation. In contrast with prior missing-data multitaper estimators that force standard Slepian sequences to be zero at gaps, the proposed missing-data Slepian sequences are defined only where data are present. The missing-data Slepian sequences are frequency independent, as are the eigenvalues that define the energy concentration within the resolution bandwidth, when the process bandwidth is [−1/2,1/2) for unit sampling and the sampling scheme comprises integer multiples of unity. As a consequence, one need only compute the ensuing missing-data Slepian sequences for a given sampling scheme once, and then the spectrum at an arbitrary set of frequencies can be computed using them. It is also shown that the resulting missing-data multitaper estimator can incorporate all of the optimality features (i.e. adaptive-weighting, F-test and reshaping) of the standard multitaper estimator, and can be applied to bivariate or multivariate situations in similar ways. Performance of the missing-data multitaper estimator is illustrated using length of day, seafloor pressure and Nile River low stand time-series.

Description: The length of day utilized in Section 3 are available from http://hpiers.obspm.fr. The pressure data used in Section 4 are available from https://doi.org/10.1029/2018JC014586. A Matlab function MDmwps.m to compute missing-data power spectra is available from the Mathworks file exchange website. The author thanks Jeff Park and editor F.J. Simons for thorough reviews. This work was supported by an Internal Research and Development award at WHOI, and by the Walter A. and Hope Noyes Smith Chair for Excellence in Oceanography.

Description: © The Author(s), 2019. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Johnson, L. K., Alexander, H., & Brown, C. T. Re-assembly, quality evaluation, and annotation of 678 microbial eukaryotic reference transcriptomes. Gigascience, 8(4), (2019): giy158, doi: 10.1093/gigascience/giy158.

Description: Background: De novo transcriptome assemblies are required prior to analyzing RNA sequencing data from a species without an existing reference genome or transcriptome. Despite the prevalence of transcriptomic studies, the effects of using different workflows, or “pipelines,” on the resulting assemblies are poorly understood. Here, a pipeline was programmatically automated and used to assemble and annotate raw transcriptomic short-read data collected as part of the Marine Microbial Eukaryotic Transcriptome Sequencing Project. The resulting transcriptome assemblies were evaluated and compared against assemblies that were previously generated with a different pipeline developed by the National Center for Genome Research. Results: New transcriptome assemblies contained the majority of previous contigs as well as new content. On average, 7.8% of the annotated contigs in the new assemblies were novel gene names not found in the previous assemblies. Taxonomic trends were observed in the assembly metrics. Assemblies from the Dinoflagellata showed a higher number of contigs and unique k-mers than transcriptomes from other phyla, while assemblies from Ciliophora had a lower percentage of open reading frames compared to other phyla. Conclusions: Given current bioinformatics approaches, there is no single “best” reference transcriptome for a particular set of raw data. As the optimum transcriptome is a moving target, improving (or not) with new tools and approaches, automated and programmable pipelines are invaluable for managing the computationally intensive tasks required for re-processing large sets of samples with revised pipelines and ensuring a common evaluation workflow is applied to all samples. Thus, re-assembling existing data with new tools using automated and programmable pipelines may yield more accurate identification of taxon-specific trends across samples in addition to novel and useful products for the community.

Description: Funding was provided by the Gordon and Betty Moore Foundation under award number GBMF4551 to C.T.B. Jetstream cloud platform was used with XSEDE allocation TG-BIO160028 [66, 67].

Description: Author Posting. © Crown Copyright, 2015. 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 204 (2016): 1-20, doi:10.1093/gji/ggv416.

Description: The Canada Basin and the southern Alpha-Mendeleev ridge complex underlie a significant proportion of the Arctic Ocean, but the geology of this undrilled and mostly ice-covered frontier is poorly known. New information is encoded in seismic wide-angle reflections and refractions recorded with expendable sonobuoys between 2007 and 2011. Velocity–depth samples within the sedimentary succession are extracted from published analyses for 142 of these records obtained at irregularly spaced stations across an area of 1.9E + 06 km2. The samples are modelled at regional, subregional and station-specific scales using an exponential function of inverse velocity versus depth with regionally representative parameters determined through numerical regression. With this approach, smooth, non-oscillatory velocity–depth profiles can be generated for any desired location in the study area, even where the measurement density is low. Practical application is demonstrated with a map of sedimentary thickness, derived from seismic reflection horizons interpreted in the time domain and depth converted using the velocity–depth profiles for each seismic trace. A thickness of 12–13 km is present beneath both the upper Mackenzie fan and the middle slope off of Alaska, but the sedimentary prism thins more gradually outboard of the latter region. Mapping of the observed-to-predicted velocities reveals coherent geospatial trends associated with five subregions: the Mackenzie fan; the continental slopes beyond the Mackenzie fan; the abyssal plain; the southwestern Canada Basin; and, the Alpha-Mendeleev magnetic domain. Comparison of the subregional velocity–depth models with published borehole data, and interpretation of the station-specific best-fitting model parameters, suggests that sandstone is not a predominant lithology in any of the five subregions. However, the bulk sand-to-shale ratio likely increases towards the Mackenzie fan, and the model for this subregion compares favourably with borehole data for Miocene turbidites in the eastern Gulf of Mexico. The station-specific results also indicate that Quaternary sediments coarsen towards the Beaufort-Mackenzie and Banks Island margins in a manner that is consistent with the variable history of Laurentide Ice Sheet advance documented for these margins. Lithological factors do not fully account for the elevated velocity–depth trends that are associated with the southwestern Canada Basin and the Alpha-Mendeleev magnetic domain. Accelerated porosity reduction due to elevated palaeo-heat flow is inferred for these regions, which may be related to the underlying crustal types or possibly volcanic intrusion of the sedimentary succession. Beyond exploring the variation of an important physical property in the Arctic Ocean basin, this study provides comparative reference for global studies of seismic velocity, burial history, sedimentary compaction, seismic inversion and overpressure prediction, particularly in mudrock-dominated successions.

Description: Author Posting. © The Author, 2015. 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 203 (2015): 893-895, doi:10.1093/gji/ggv324.

Description: The statistics of directional data on a sphere can be modelled either using the Fisher distribution that is conditioned on the magnitude being unity, in which case the sample space is confined to the unit sphere, or using the latitude–longitude marginal distribution derived from a trivariate Gaussian model that places no constraint on the magnitude. These two distributions are derived from first principles and compared. The Fisher distribution more closely approximates the uniform distribution on a sphere for a given small value of the concentration parameter, while the latitude–longitude marginal distribution is always slightly larger than the Fisher distribution at small off-axis angles for large values of the concentration parameter. Asymptotic analysis shows that the two distributions only become equivalent in the limit of large concentration parameter and very small off-axis angle.

Description: © The Author(s), 2015. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Genome Biology and Evolution 7 (2015): 3207-3225, doi:10.1093/gbe/evv210.

Description: High-throughput sequencing of reduced representation libraries obtained through digestion with restriction enzymes—generically known as restriction site associated DNA sequencing (RAD-seq)—is a common strategy to generate genome-wide genotypic and sequence data from eukaryotes. A critical design element of any RAD-seq study is knowledge of the approximate number of genetic markers that can be obtained for a taxon using different restriction enzymes, as this number determines the scope of a project, and ultimately defines its success. This number can only be directly determined if a reference genome sequence is available, or it can be estimated if the genome size and restriction recognition sequence probabilities are known. However, both scenarios are uncommon for nonmodel species. Here, we performed systematic in silico surveys of recognition sequences, for diverse and commonly used type II restriction enzymes across the eukaryotic tree of life. Our observations reveal that recognition sequence frequencies for a given restriction enzyme are strikingly variable among broad eukaryotic taxonomic groups, being largely determined by phylogenetic relatedness. We demonstrate that genome sizes can be predicted from cleavage frequency data obtained with restriction enzymes targeting “neutral” elements. Models based on genomic compositions are also effective tools to accurately calculate probabilities of recognition sequences across taxa, and can be applied to species for which reduced representation data are available (including transcriptomes and neutral RAD-seq data sets). The analytical pipeline developed in this study, PredRAD (https://github.com/phrh/PredRAD), and the resulting databases constitute valuable resources that will help guide the design of any study using RAD-seq or related methods.

Description: This research was supported by the Office of Ocean Exploration and Research of the National Oceanic and Atmospheric Administration (NA09OAR4320129 to T.S.); the Division of Ocean Sciences of the National Science Foundation (OCE-1131620 to T.S.); the Astrobiology Science and Technology for Exploring Planets program of the National Aeronautics and Space Administration (NNX09AB76G to T.S.); and the Academic Programs Office (Ocean Ventures Fund to S.H.), the Ocean Exploration Institute (Fellowship support to T.M.S.), and the Ocean Life Institute of the Woods Hole Oceanographic Institution (internal grant to T.M.S. and S.H.).

Description: © The Author(s), 2017. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Molecular Biology and Evolution 34 (2017): 1890-1901, doi:10.1093/molbev/msx125.

Description: The highly conserved ADAR enzymes, found in all multicellular metazoans, catalyze the editing of mRNA transcripts by the deamination of adenosines to inosines. This type of editing has two general outcomes: site specific editing, which frequently leads to recoding, and clustered editing, which is usually found in transcribed genomic repeats. Here, for the first time, we looked for both editing of isolated sites and clustered, non-specific sites in a basal metazoan, the coral Acropora millepora during spawning event, in order to reveal its editing pattern. We found that the coral editome resembles the mammalian one: it contains more than 500,000 sites, virtually all of which are clustered in non-coding regions that are enriched for predicted dsRNA structures. RNA editing levels were increased during spawning and increased further still in newly released gametes. This may suggest that editing plays a role in introducing variability in coral gametes.

Description: This work was supported by the Australian Research Council (to PK), the European Research Council (grant 311257), the I-CORE Program of the Planning and Budgeting Committee in Israel (grants 41/11 and 1796/12), and the Israel Science Foundation (1380/14).

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: Author Posting. © The Author(s), 2015. 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 203 (2015): 1-21, doi:10.1093/gji/ggv251.

Description: We examine along-axis variations in melt content of the axial magma lens (AML) beneath the fast-spreading East Pacific Rise (EPR) using an amplitude variation with angle of incidence (AVA) crossplotting method applied to multichannel seismic data acquired in 2008. The AVA crossplotting method, which has been developed for and, so far, applied for hydrocarbon prospection in sediments, is for the first time applied to a hardrock environment. We focus our analysis on 2-D data collected along the EPR axis from 9°29.8′N to 9°58.4′N, a region which encompasses the sites of two well-documented submarine volcanic eruptions (1991–1992 and 2005–2006). AVA crossplotting is performed for a ∼53 km length of the EPR spanning nine individual AML segments (ranging in length from ∼3.2 to 8.5 km) previously identified from the geometry of the AML and disruptions in continuity. Our detailed analyses conducted at 62.5 m interval show that within most of the analysed segments melt content varies at spatial scales much smaller (a few hundred of metres) than the length of the fine-scale AML segments, suggesting high heterogeneity in melt concentration. At the time of our survey, about 2 yr after the eruption, our results indicate that the three AML segments that directly underlie the 2005–2006 lava flow are on average mostly molten. However, detailed analysis at finer-scale intervals for these three segments reveals AML pockets (from 〉62.5 to 812.5 m long) with a low melt fraction. The longest such mushy section is centred beneath the main eruption site at ∼9°50.4′N, possibly reflecting a region of primary melt drainage during the 2005–2006 event. The complex geometry of fluid flow pathways within the crust above the AML and the different response times of fluid flow and venting to eruption and magma reservoir replenishment may contribute to the poor spatial correlation between incidence of hydrothermal vents and presence of highly molten AML. The presented results are an important step forward in our ability to resolve small-scale characteristics of the AML and recommend the AVA crossplotting as a tool for examining mid-ocean ridge magma-systems elsewhere.

Description: This research was supported by NSF awards OCE0327872 to J.C.M. and S.M.C., OCE-0327885 to J.P.C., and OCE0624401 to M.R.N.

Description: © The Author(s), 2016. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Nucleic Acids Research 44 (2016): e157, doi:10.1093/nar/gkw738.

Description: Site-directed RNA editing (SDRE) is a strategy to precisely alter genetic information within mRNAs. By linking the catalytic domain of the RNA editing enzyme ADAR to an antisense guide RNA, specific adenosines can be converted to inosines, biological mimics for guanosine. Previously, we showed that a genetically encoded iteration of SDRE could target adenosines expressed in human cells, but not efficiently. Here we developed a reporter assay to quantify editing, and used it to improve our strategy. By enhancing the linkage between ADAR's catalytic domain and the guide RNA, and by introducing a mutation in the catalytic domain, the efficiency of converting a UAG premature termination codon (PTC) to tryptophan (UGG) was improved from ∼11% to ∼70%. Other PTCs were edited, but less efficiently. Numerous off-target edits were identified in the targeted mRNA, but not in randomly selected endogenous messages. Off-target edits could be eliminated by reducing the amount of guide RNA with a reduction in on-target editing. The catalytic rate of SDRE was compared with those for human ADARs on various substrates and found to be within an order of magnitude of most. These data underscore the promise of site-directed RNA editing as a therapeutic or experimental tool.

Description: National Institutes of Health [1R0111223855, 1R01NS64259]; Cystic Fibrosis Foundation Therapeutics [Rosent14XXO]; Infrastructural support was provided by the National Institutes of Health [NIGMS 1P20GM103642, NIMHD 8G12-MD007600]; National Science Foundation [DBI 0115825, DBI 1337284]; Department of Defense [52680-RT-ISP].

Description: Author Posting. © The Authors, 2016. 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 208 (2017): 1026-1042, doi:10.1093/gji/ggw435.

Description: In recent years, marine controlled source electromagnetics (CSEM) has found increasing use in hydrocarbon exploration due to its ability to detect thin resistive zones beneath the seafloor. It is the purpose of this paper to evaluate the physics of CSEM for an ocean whose electrical thickness is comparable to or much thinner than that of the overburden using the in-line configuration through examination of the elliptically polarized seafloor electric field, the time-averaged energy flow depicted by the real part of the complex Poynting vector, energy dissipation through Joule heating and the Fréchet derivatives of the seafloor field with respect to the subseafloor conductivity that is assumed to be isotropic. The deep water (ocean layer electrically much thicker than the overburden) seafloor EM response for a model containing a resistive reservoir layer has a greater amplitude and reduced phase as a function of offset compared to that for a half-space, or a stronger and faster response. For an ocean whose electrical thickness is comparable to or much smaller than that of the overburden, the electric field displays a greater amplitude and reduced phase at small offsets, shifting to a stronger amplitude and increased phase at intermediate offsets and a weaker amplitude and enhanced phase at long offsets, or a stronger and faster response that first changes to stronger and slower, and then transitions to weaker and slower. These transitions can be understood by visualizing the energy flow throughout the structure caused by the competing influences of the dipole source and guided energy flow in the reservoir layer, and the air interaction caused by coupling of the entire subseafloor resistivity structure with the sea surface. A stronger and faster response occurs when guided energy flow is dominant, while a weaker and slower response occurs when the air interaction is dominant. However, at intermediate offsets for some models, the air interaction can partially or fully reverse the direction of energy flux in the reservoir layer toward rather than away from the source, resulting in a stronger and slower response. The Fréchet derivatives are dominated by preferential sensitivity to the reservoir layer conductivity for all water depths except at high frequencies, but also display a shift with offset from the galvanic to the inductive mode in the underburden and overburden due to the interplay of guided energy flow and the air interaction. This means that the sensitivity to the horizontal conductivity is almost as strong as to the vertical component in the shallow parts of the subsurface, and in fact is stronger than the vertical sensitivity deeper down. However, the sensitivity to horizontal conductivity is still weak compared to the vertical component within thin resistive regions. The horizontal sensitivity is gradually decreased when the water becomes deep. These observations in part explain the success of shallow towed CSEM using only measurements of the in-line component of the electric field.

Description: © The Author(s), 2017. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Bioscience 67 (2017): 760–768, doi:10.1093/biosci/bix059.

Description: As the sampling frequency and resolution of Earth observation imagery increase, there are growing opportunities for novel applications in population monitoring. New methods are required to apply established analytical approaches to data collected from new observation platforms (e.g., satellites and unmanned aerial vehicles). Here, we present a method that estimates regional seasonal abundances for an understudied and growing population of gray seals (Halichoerus grypus) in southeastern Massachusetts, using opportunistic observations in Google Earth imagery. Abundance estimates are derived from digital aerial survey counts by adapting established correction-based analyses with telemetry behavioral observation to quantify survey biases. The result is a first regional understanding of gray seal abundance in the northeast US through opportunistic Earth observation imagery and repurposed animal telemetry data. As species observation data from Earth observation imagery become more ubiquitous, such methods provide a robust, adaptable, and cost-effective solution to monitoring animal colonies and understanding species abundances.

Description: We would like to thank generous support from International Fund for Animal Welfare, the Bureau of Ocean Energy, and the Oak Foundation for funding support for the telemetry devices.

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: 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: Author Posting. © The Authors, 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 211 (2017): 1046–1061, doi:10.1093/gji/ggx360.

Description: In recent years, marine controlled source electromagnetics (CSEM) has found increasing use in hydrocarbon exploration due to its ability to detect thin resistive zones beneath the seafloor. It is the purpose of this paper to evaluate the physics of CSEM for an ocean whose electrical thickness is comparable to or much thinner than that of the overburden using the in-line configuration through examination of the elliptically-polarized seafloor electric field, the time-averaged energy flow depicted by the real part of the complex Poynting vector, energy dissipation through Joule heating and the Fréchet derivatives of the seafloor field with respect to the sub-seafloor conductivity that is assumed to be transversely anisotropic, with a vertical-to-horizontal resistivity ratio of 3:1. For an ocean whose electrical thickness is comparable to that of the overburden, the seafloor electromagnetic response for a model containing a resistive reservoir layer has a greater amplitude and reduced phase as a function of offset compared to that for a halfspace, or a stronger and faster response, and displays little to no evidence for the air interaction. For an ocean whose electrical thickness is much smaller than that of the overburden, the electric field displays a greater amplitude and reduced phase at small offsets, shifting to a stronger amplitude and increased phase at intermediate offsets, and a weaker amplitude and enhanced phase at long offsets, or a stronger and faster response that first changes to stronger and slower, and then transitions to weaker and slower. By comparison to the isotropic case with the same horizontal conductivity, transverse anisotropy stretches the Poynting vector and the electric field response from a thin resistive layer to much longer offsets. These phenomena can be understood by visualizing the energy flow throughout the structure caused by the competing influences of the dipole source and guided energy flow in the reservoir layer, and the air interaction caused by coupling of the entire sub-seafloor resistivity structure with the sea surface. The Fréchet derivatives are dominated by preferential sensitivity to the vertical conductivity in the reservoir layer and overburden at short offsets. The horizontal conductivity Fréchet derivatives are weaker than to comparable to the vertical derivatives at long offsets in the substrate. This means that the sensitivity to the horizontal conductivity is present in the shallow parts of the subsurface. In the presence of transverse anisotropy, it is necessary to go to higher frequencies to sense the horizontal conductivity in the overburden as compared to an isotropic model with the same horizontal conductivity. These observations in part explain the success of shallow towed CSEM using only measurements of the in-line component of the electric field.

Description: This work was supported at WHOI by an Independent Research and Development award, and by the Walter A. and Hope Noyes Smith Chair for Excellence in Oceanography.

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: 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: 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: Repeated and independent emergence of trait divergence that matches habitat differences is a sign of parallel evolution by natural selection. Yet, the molecular underpinnings that are targeted by adaptive evolution often remain elusive. We investigate this question by combining genome-wide analyses of copy number variants (CNVs), single nucleotide polymorphisms (SNPs), and gene expression across four pairs of lake and river populations of the three-spined stickleback (Gasterosteus aculeatus). We tested whether CNVs that span entire genes and SNPs occurring in putative cis-regulatory regions contribute to gene expression differences between sticklebacks from lake and river origins. We found 135 gene CNVs that showed a significant positive association between gene copy number and gene expression, suggesting that CNVs result in dosage effects that can fuel phenotypic variation and serve as substrates for habitat-specific selection. Copy number differentiation between lake and river sticklebacks also contributed to expression differences of two immune-related genes in immune tissues, cathepsin A and GIMAP7. In addition, we identified SNPs in cis-regulatory regions (eSNPs) associated with the expression of 1,865 genes, including one eSNP upstream of a carboxypeptidase gene where both the SNP alleles differentiated and the gene was differentially expressed between lake and river populations. Our study highlights two types of mutations as important sources of genetic variation involved in the evolution of gene expression and in potentially facilitating repeated adaptation to novel environments.

Description: Cichlid fishes provide textbook examples of explosive phenotypic diversification and sympatric speciation, thereby making them ideal systems for studying the molecular mechanisms underlying rapid lineage divergence. Despite the fact that gene regulation provides a critical link between diversification in gene function and speciation, many genomic regulatory mechanisms such as microRNAs (miRNAs) have received little attention in these rapidly diversifying groups. Therefore, we investigated the posttranscriptional regulatory role of miRNAs in the repeated sympatric divergence of Midas cichlids (Amphilophus spp.) from Nicaraguan crater lakes. Using miRNA and mRNA sequencing of embryos from five Midas species, we first identified miRNA binding sites in mRNAs and highlighted the presences of a surprising number of novel miRNAs in these adaptively radiating species. Then, through analyses of expression levels, we identified putative miRNA/gene target pairs with negatively correlated expression level that were consistent with the role of miRNA in downregulating mRNA. Furthermore, we determined that several miRNA/gene pairs show convergent expression patterns associated with the repeated benthic/limnetic sympatric species divergence implicating these miRNAs as potential molecular mechanisms underlying replicated sympatric divergence. Finally, as these candidate miRNA/gene pairs may play a central role in phenotypic diversification in these cichlids, we characterized the expression domains of selected miRNAs and their target genes via in situ hybridization, providing further evidence that miRNA regulation likely plays a role in the Midas cichlid adaptive radiation. These results provide support for the hypothesis that extremely quickly evolving miRNA regulation can contribute to rapid evolutionary divergence even in the presence of gene flow.

Description: The Arctic Limnocalanus macrurus is a prominent representative of large copepods which performs several essential functions in freshwater and marine ecosystems. Being a cold stenotherm species, its distribution is primarily confined to deeper water layers. Based on the long-term observations from one of the largest spatially confined natural populations of this species in the Baltic Sea, we detected profound long-term variability of L. macrurus during 1958–2016: high abundances before the 1980s, then nearly disappearance in the 1990s and recovery in the 2000s. The main environmental parameters explaining the interannual variability of L. macrurus in spring were herring spawning stock biomass in preceding year, winter severity, and bottom water temperature in preceding summer. The effect of winter severity and water temperature was also non-linear. The sliding window correlation analysis pointed to a non-stationary relationship between the abundance of L. macrurus and the key variables. Given the observed pronounced seasonality in the population structure of L. macrurus (young stages dominated in the beginning of the year and only adults were left in the population in summer and autumn) we identified the dynamics of key environmental variables to understand this species under different ecosystem configurations and different combinations of drivers of change.

Description: The Swan Islands Transform Fault (SITF) marks the southern boundary of the Cayman Trough and the ocean–continent transition of the North American–Caribbean Plate boundary offshore Honduras. The CAYSEIS experiment acquired a 180-km-long seismic refraction and gravity profile across this transform margin, ∼70 km to the west of the Mid-Cayman Spreading Centre (MCSC). This profile shows the crustal structure across a transform fault system that juxtaposes Mesozoic-age continental crust to the south against the ∼10-Myr-old ultraslow spread oceanic crust to the north. Ocean-bottom seismographs were deployed along-profile, and inverse and forward traveltime modelling, supported by gravity analysis, reveals ∼23-km-thick continental crust that has been thinned over a distance of ∼70 km to ∼10 km-thick at the SITF, juxtaposed against ∼4-km-thick oceanic crust. This thinning is primarily accommodated within the lower crust. Since Moho reflections are not widely observed, the 7.0 km s−1 velocity contour is used to define the Moho along-profile. The apparent lack of reflections to the north of the SITF suggests that the Moho is more likely a transition zone between crust and mantle. Where the profile traverses bathymetric highs in the off-axis oceanic crust, higher P-wave velocity is observed at shallow crustal depths. S-wave arrival modelling also reveals elevated velocities at shallow depths, except for crust adjacent to the SITF that would have occupied the inside corner high of the ridge-transform intersection when on axis. We use a Vp/Vs ratio of 1.9 to mark where lithologies of the lower crust and uppermost mantle may be exhumed, and also to locate the upper-to-lower crustal transition, identify relict oceanic core complexes and regions of magmatically formed crust. An elevated Vp/Vs ratio suggests not only that serpentinized peridotite may be exposed at the seafloor in places, but also that seawater has been able to flow deep into the crust and upper mantle over 20–30-km-wide regions which may explain the lack of a distinct Moho. The SITF has higher velocities at shallower depths than observed in the oceanic crust to the north and, at the seabed, it is a relatively wide feature. However, the velocity–depth model subseabed suggests a fault zone no wider than ∼5–10 km, that is mirrored by a narrow seabed depression ∼7500 m deep. Gravity modelling shows that the SITF is also underlain, at 〉2 km subseabed, by a ∼20-km-wide region of density 〉3000 kg m−3 that may reflect a broad region of metamorphism. The residual mantle Bouguer anomaly across the survey region, when compared with the bathymetry, suggests that the transform may also have a component of left-lateral trans-tensional displacement that accounts for its apparently broad seabed appearance, and that the focus of magma supply may currently be displaced to the north of the MCSC segment centre. Our results suggest that Swan Islands margin development caused thinning of the adjacent continental crust, and that the adjacent oceanic crust formed in a cool ridge setting, either as a result of reduced mantle upwelling and/or due to fracture enhanced fluid flow.

Description: Standard seismic acquisition and processing require appropriate source-receiver offsets. P-cable technology represents the opposite, namely, very short source-receiver offsets at the price of increased spatial and lateral resolution with a high-frequency source. To use this advantage, a processing flow excluding offset information is required. This aim can be achieved with a processing tuned to diffractions because point diffractions scatter the same information in offset and midpoint direction. Usually, diffractions are small amplitude events and a careful diffraction separation is required as a first step. We suggest the strategy to use a multiparameter stacking operator, e.g, common-reflection surface, and stack along the midpoint direction. The obtained kinematic wavefront attributes are used to calculate time-migration velocities. A diffractivity map serves as filter to refine the velocities. This strategy is applied to a 3D P-cable data set to obtain a time-migrated image.

Description: This review shows that the presence of seagrass microbial community is critical for the development of seagrasses; from seed germination, through to phytohormone production and enhanced nutrient availability, and defence against pathogens and saprophytes. The tight seagrass-bacterial relationship highlighted in this review supports the existence of a seagrass holobiont and adds to the growing evidence for the importance of marine eukaryotic microorganisms in sustaining vital ecosystems. Incorporating a micro-scale view on seagrass ecosystems substantially expands our understanding of ecosystem functioning and may have significant implications for future seagrass management and mitigation against human disturbance.

Description: This article has been accepted for publication in Geophysical Journal International ©: The Authors 2019. Published by Oxford University Press on behalf of the Royal Astronomical Society. All rights reserved. Uploaded in accordance with the publisher's self-archiving policy.

Description: The Amatrice–Norcia–Visso sequence is characterized by complex behaviour that is somewhat atypical of main-shock–aftershock sequences, as there were multiple large main shocks that continued for months. In this study we focus on the Amatrice sequence (main shock 2016 August 24, Mw = 5.97) to evaluate the apparent stress values and magnitude-dependent scaling in order to improve our knowledge of processes that control small and large earthquakes within this active region of Italy. Apparent stress is proportional to the ratio of radiated seismic energy and seismic moment, and as such, these stress parameters play an important role in hazard prediction as they have a strong effect on the observed and predicted ground shaking. We analyse 83 events of the sequence from 2016 August 24 to October 16, within a radius of 20 km from the main shock and with an Mw ranging between 5.97 and 2.72. Taking advantage of the averaging nature of coda waves, we analyse coda-envelope-based spectral ratios between neighbouring event pairs.We use equations proposed byWalter et al. to consider stable, low-frequency and high-frequency spectral ratio levelswhich provide measures of the corner frequency and apparent stress ratios of the events within the sequence. The results demonstrate non-self-similar behaviour within the sequence, suggesting a change in dynamics between the largest events and the smaller aftershocks. The apparent stress and corner frequency estimates are compared to those obtained by Malagnini and Munaf`o who utilized hundreds of direct S-wave spectral ratio measurements to obtain their results. Although our analysis is based only on 83 events, our results are in very good agreement, demonstrating once more that the use of coda waves is very stable and provides lower variance measures than those using direct waves. A comparison with recent Central Apennines source scaling models derived from various seismic sequences (1997–1998 Colfiorito, 2002 San Giuliano di Puglia, 2009 L’Aquila) shows that the Amatrice sequence source scaling in this study is well represented by the models proposed by Pacor et al. and Malagnini and Mayeda.

Description: Published

Description: 446-455

Description: 3T. Sorgente sismica

Description: JCR Journal

Description: We present the analysis of rotational and translational ground motions from earthquakes recorded during October–November 2016, in association with the Central Italy seismic sequence. We use co-located measurements of the vertical ground rotation rate from a large ring laser gyroscope and the three components of ground velocity from a broad-band seismometer. Both instruments are positioned in a deep underground environment, within the Gran Sasso National Laboratories of the Istituto Nazionale di Fisica Nucleare. We collected dozens of events spanning the 3.5–5.9 magnitude range and epicentral distances between 30 and 70 km. This data set constitutes an unprecedented observation of the vertical rotational motions associated with an intense seismic sequence at local distance. Under the plane-wave approximation we process the data set in order to get an experimental estimation of the events backazimuth. Peak values of rotation rate (PRR) and horizontal acceleration (PGA) are markedly correlated, according to a scaling constant which is consistent with previous measurements from different earthquake sequences. We used a prediction model in use for Italy to calculate the expected PGA at the recording site, obtaining consequently predictions for PRR. Within the modelling uncertainties, predicted rotations are consistent with the observed ones, suggesting the possibility of establishing specific attenuation models for ground rotations, like the scaling of peak velocity and peak acceleration in empirical ground-motion prediction relationships. In a second step, after identifying the direction of the incoming wavefield, we extract phase-velocity data using the spectral ratio of the translational and rotational components. This analysis is performed over time windows associated with the P-coda, S-coda and Lg phase. Results are consistent with independent estimates of shear wave velocities in the shallow crust of the Central Apennines

Description: Published

Description: 705-715

Description: 4T. Sismicità dell'Italia

Description: JCR Journal

Description: Ubiquitous SAR11 Alphaproteobacteria numerically dominate marine planktonic communities. Because they are excruciatingly difficult to cultivate, there is comparatively little known about their physiology and metabolic responses to long- and short-term environmental changes. As surface oceans take up anthropogenic, atmospheric CO2, the consequential process of ocean acidification could affect the global biogeochemical significance of SAR11. Shipping accidents or inadvertent release of chemicals from industrial plants can have strong short-term local effects on oceanic SAR11. This study investigated the effect of 2.5-fold acidification of seawater on the metabolism of SAR11 and other heterotrophic bacterioplankton along a natural temperature gradient crossing the North Atlantic Ocean, Norwegian and Greenland Seas. Uptake rates of the amino acid leucine by SAR11 cells as well as other bacterioplankton remained similar to controls despite an instant ∼50% increase in leucine bioavailability upon acidification. This high physiological resilience to acidification even without acclimation, suggests that open ocean dominant bacterioplankton are able to cope even with sudden and therefore more likely with long-term acidification effects.

Description: Climate change influences the Baltic Sea ecosystem via its effects on oceanography and biogeochemistry. Sea surface temperature has been projected to increase by 2 to 4 °C until 2100 due to global warming; the changes will be more significant in the northern areas and less so in the south. The warming up will also diminish the annual sea ice cover by 57% to 71%, and ice season will be one to three months shorter than in the early 21st century, depending on latitude. A significant decrease in sea surface salinity has been projected because of an increase in rainfall and decrease of saline inflows into the Baltic Sea. The increasing surface flow has, in turn, been projected to increase leaching of nutrients from the soil to the watershed and eventually into the Baltic Sea. Also, acidification of the seawater and sea-level rise have been predicted. Increasing seawater temperature speeds up metabolic processes and increases growth rates of many secondary producers. Species associated with sea ice, from salt brine microbes to seals, will suffer. Due to the specific salinity tolerances, species’ geographical ranges may shift by tens or hundreds of kilometres with decreasing salinity. A decrease in pH will slow down calcification of bivalve shells, and higher temperatures also alleviate establishment of non-indigenous species originating from more southern sea areas. Many uncertainties still remain in predicting the couplings between atmosphere, oceanography and ecosystem. Especially projections of many oceanographic parameters, such as wind speeds and directions, the mean salinity level, and density stratification, are still ambiguous. Also, the effects of simultaneous changes in multiple environmental factors on species with variable preferences to temperature, salinity, and nutrient conditions are difficult to project. There is, however, enough evidence to claim that due to increasing runoff of nutrients from land and warming up of water, primary production and sedimentation of organic matter will increase; this will probably enhance anoxia and release of phosphorus from sediments. Such changes may keep the Baltic Sea in an eutrophicated state for a long time, unless strong measures to decrease nutrient runoff from land are taken. Changes in the pelagic and benthic communities are anticipated. Benthic communities will change from marine to relatively more euryhaline communities and will suffer from hypoxic events. The projected temperature increase and salinity decline will contribute to maintain the pelagic ecosystem of the Central Baltic and the Gulf of Finland in a state dominated by cyanobacteria, flagellates, small-sized zooplankton and sprat, instead of diatoms, large marine copepods, herring, and cod. Effects vary from area to area, however. In particular the Bothnian Sea, where hypoxia is less common and rivers carry a lot of dissolved organic carbon, primary production will probably not increase as much as in the other basins. The coupled oceanography-biogeochemistry ecosystem models have greatly advanced our understanding of the effects of climate change on marine ecosystems. Also, studies on climate associated “regime shifts” and cascading effects from top predators to plankton have been fundamental for understanding of the response of the Baltic Sea ecosystem to anthropogenic and climatic stress. In the future, modeling efforts should be focusing on coupling of biogeochemical processes and lower trophic levels to the top predators. Also, fine resolution species distribution models should be developed and combined with 3-D modelling, to describe how the species and communities are responding to climate-induced changes in environmental variables.

Description: We present high-resolution resistivity imaging of gas hydrate pipe-like structures, as derived from marine controlled-source electromagnetic (CSEM) inversions that combine towed and ocean-bottom electric field receiver data, acquired from the Nyegga region, offshore Norway. Two-dimensional CSEM inversions applied to the towed receiver data detected four new prominent vertical resistive features that are likely gas hydrate structures, located in proximity to a major gas hydrate pipe-like structure, known as the CNE03 pockmark. The resistivity model resulting from the CSEM data inversion resolved the CNE03 hydrate structure in high resolution, as inferred by comparison to seismically constrained inversions. Our results indicate that shallow gas hydrate vertical features can be delineated effectively by inverting both ocean-bottom and towed receiver CSEM data simultaneously. The approach applied here can be utilised to map and monitor seafloor mineralisation, freshwater reservoirs, CO2 sequestration sites and near-surface geothermal systems.

Description: Pseudovibrio is a marine bacterial genus members of which are predominantly isolated from sessile marine animals, and particularly sponges. It has been hypothesised that Pseudovibrio spp. form mutualistic relationships with their hosts. Here, we studied Pseudovibrio phylogeny and genetic adaptations that may play a role in host colonization by comparative genomics of 31 Pseudovibrio strains, including 25 sponge isolates. All genomes were highly similar in terms of encoded core metabolic pathways, albeit with substantial differences in overall gene content. Based on gene composition, Pseudovibrio spp. clustered by geographic region, indicating geographic speciation. Furthermore, the fact that isolates from the Mediterranean Sea clustered by sponge species suggested host-specific adaptation or colonization. Genome analyses suggest that Pseudovibrio hongkongensis UST20140214-015BT is only distantly related to other Pseudovibrio spp., thereby challenging its status as typical Pseudovibrio member. All Pseudovibrio genomes were found to encode numerous proteins with SEL1 and tetratricopeptide repeats, which have been suggested to play a role in host colonization. For evasion of the host immune system, Pseudovibrio spp. may depend on type III, IV and VI secretion systems that can inject effector molecules into eukaryotic cells. Furthermore, Pseudovibrio genomes carry on average seven secondary metabolite biosynthesis clusters, reinforcing the role of Pseudovibrio spp. as potential producers of novel bioactive compounds. Tropodithietic acid, bacteriocin and terpene biosynthesis clusters were highly conserved within the genus, suggesting an essential role in survival e.g. through growth inhibition of bacterial competitors. Taken together, these results support the hypothesis that Pseudovibrio spp. have mutualistic relations with sponges.

Description: Marine prey and predators will respond to future climate through physiological and behavioral adjustments. However, our understanding of how such direct effects may shift the outcome of predator–prey interactions is still limited. Here, we investigate the effects of ocean warming and acidification on foraging behavior and biomass of a common prey (shrimps, Palaemon spp.) tested in large mesocosms harboring natural resources and habitats. Acidification did not alter foraging behavior in prey. Under warming, however, prey showed riskier behavior by foraging more actively and for longer time periods, even in the presence of a live predator. No effects of longer-term exposure to climate stressors were detected on prey biomass. Our findings suggest that ocean warming may increase the availability of some prey to predators via a behavioral pathway (i.e., increased risk-taking by prey), likely by elevating metabolic demand of prey species.

Description: Ocean acidification (OA) is increasing due to anthropogenic CO2 emissions and poses a threat to marine species and communities worldwide. To better project the effects of acidification on organisms’ health and persistence, an understanding is needed of the 1) mechanisms underlying developmental and physiological tolerance and 2) potential populations have for rapid evolutionary adaptation. This is especially challenging in nonmodel species where targeted assays of metabolism and stress physiology may not be available or economical for large-scale assessments of genetic constraints. We used mRNA sequencing and a quantitative genetics breeding design to study mechanisms underlying genetic variability and tolerance to decreased seawater pH (-0.4 pH units) in larvae of the sea urchin Strongylocentrotus droebachiensis. We used a gene ontology-based approach to integrate expression profiles into indirect measures of cellular and biochemical traits underlying variation in larval performance (i.e., growth rates). Molecular responses to OA were complex, involving changes to several functions such as growth rates, cell division, metabolism, and immune activities. Surprisingly, the magnitude of pH effects on molecular traits tended to be small relative to variation attributable to segregating functional genetic variation in this species. We discuss how the application of transcriptomics and quantitative genetics approaches across diverse species can enrich our understanding of the biological impacts of climate change.

Description: To understand how ocean acidification (OA) influences sediment microbial communities, naturally CO2-rich sites are increasingly being used as OA analogues. However, the characterization of these naturally CO2-rich sites is often limited to OA-related variables, neglecting additional environmental variables that may confound OA effects. Here, we used an extensive array of sediment and bottom water parameters to evaluate pH effects on sediment microbial communities at hydrothermal CO2seeps in Papua New Guinea. The geochemical composition of the sediment pore water showed variations in the hydrothermal signature at seep sites with comparable pH, allowing the identification of sites that may better represent future OA scenarios. At these sites, we detected a 60% shift in the microbial community composition compared with reference sites, mostly related to increases inChloroflexisequences. pH was among the factors significantly, yet not mainly, explaining changes in microbial community composition. pH variation may therefore often not be the primary cause of microbial changes when sampling is done along complex environmental gradients. Thus, we recommend an ecosystem approach when assessing OA effects on sediment microbial communities under natural conditions. This will enable a more reliable quantification of OA effects via a reduction of potential confounding effects

Description: In this pilot study, we describe a high-pressure incubation system allowing multiple subsampling of a pressurized culture without decompression. The system was tested using one piezophilic (Photobacterium profundum), one piezotolerant (Colwellia maris) bacterial strain and a decompressed sample from the Mediterranean deep sea (3044 m) determining bacterial community composition, protein production (BPP) and cell multiplication rates (BCM) up to 27 MPa. The results showed elevation of BPP at high pressure was by a factor of 1.5 ± 1.4 and 3.9 ± 2.3 for P. profundum and C. maris, respectively, compared to ambient-pressure treatments and by a factor of 6.9 ± 3.8 fold in the field samples. In P. profundum and C. maris, BCM at high pressure was elevated (3.1 ± 1.5 and 2.9 ± 1.7 fold, respectively) compared to the ambient-pressure treatments. After 3 days of incubation at 27 MPa, the natural bacterial deep-sea community was dominated by one phylum of the genus Exiguobacterium, indicating the rapid selection of piezotolerant bacteria. In future studies, our novel incubation system could be part of an isopiestic pressure chain, allowing more accurate measurement of bacterial activity rates which is important both for modeling and for predicting the efficiency of the oceanic carbon pump.

Description: The ability of microorganisms to withstand long periods with extremely low energy input has gained increasing scientific attention in recent years. Starvation experiments in the laboratory have shown that a phylogenetically wide range of microorganisms evolve fitness-enhancing genetic traits within weeks of incubation under low-energy stress. Studies on natural environments that are cut off from new energy supplies over geologic time scales, such as deeply buried sediments, suggest that similar adaptations might mediate survival under energy limitation in the environment. Yet, the extent to which laboratory-based evidence of starvation survival in pure or mixed cultures can be extrapolated to sustained microbial ecosystems in nature remains unclear. In this review, we discuss past investigations on microbial energy requirements and adaptations to energy limitation, identify gaps in our current knowledge, and outline possible future foci of research on life under extreme energy limitation.

Description: Polynoidae contains ~900 species within 18 subfamilies, some of them restricted to the deep sea. Macellicephalinae is the most diverse among these deep-sea subfamilies. In the abyssal Equatorial Pacific Ocean, the biodiversity of benthic communities is at stake in the Clarion-Clipperton Fracture Zone (CCFZ) owing to increased industrial interest in polymetallic nodules. The records of polychaetes in this region are scarce. Data gathered during the JPI Oceans cruise SO239 made a significant contribution to fill this gap, with five different localities sampled between 4000 and 5000 m depth. Benthic samples collected using an epibenthic sledge or a remotely operated vehicle resulted in a large collection of polynoids. The aims of this study are as follows: (1) to describe new species of deep-sea polynoids using morphology and molecular data (COI, 16S and 18S); and (2) to evaluate the monophyly of Macellicephalinae. Based on molecular and morphological phylogenetic analyses, ten subfamilies are synonymized with Macellicephalinae in order to create a homogeneous clade determined by the absence of lateral antennae. Within this clade, the Anantennata clade was well supported, being determined by the absence of a median antenna. Furthermore, 17 new species and four new genera are described, highlighting the high diversity hidden in the deep. A taxonomic key for the 37 valid genera of the subfamily Macellicephalinae is provided.

Description: This review article aims to provide an overview and insight into the most relevant aspects of wind energy development and current state-of-the-art. The industry is in a very mature stage, so it seems to be the right time to take stock of the relevant areas of wind energy use for power generation. For this review, the authors considered the essential aspects of the development of wind energy technology: research, modeling, and prediction of wind speed as an energy source, the technology development of the plants divided into the mechanical and electrical systems and the plant control, and finally the optimal plant operation including the maintenance strategies. The focus is on the development in Europe, with a partial focus on Germany. The authors are employees of the Fraunhofer Institutes, Institute for Energy Economics and Energy Systems Technology and Institute for Wind Energy Systems, who have contributed to the development of this technology for decades.

Description: Metabolites give us a window into the chemistry of microbes and are split into two subclasses: primary and secondary. Primary metabolites are required for life whereas secondary metabolites have historically been classified as those appearing after exponential growth and are not necessarily needed for survival. Many microbial species are estimated to produce hundreds of metabolites and can be affected by differing nutrients. Using various analytical techniques, metabolites can be directly detected in order to elucidate their biological significance. Currently, a single experiment can produce anywhere from megabytes to terabytes of data. This big data has motivated scientists to develop informatics tools to help target specific metabolites or sets of metabolites. Broadly, it is imperative to identify clear biological questions before embarking on a study of metabolites (metabolomics). For instance, studying the effect of a transposon insertion on phenazine biosynthesis in Pseudomonas is a very different from asking what molecules are present in a specific banana-derived strain of Pseudomonas. This review is meant to serve as a primer for a ‘choose your own adventure’ approach for microbiologists with limited mass spectrometry expertise, with a strong focus on liquid chromatography mass spectrometry based workflows developed or optimized within the past five years.

Description: We extend the range of validity of the artis 3D radiative transfer code up to hundreds of days after explosion, when Type Ia supernovae (SNe Ia) are in their nebular phase. To achieve this, we add a non-local thermodynamic equilibrium population and ionization solver, a new multifrequency radiation field model, and a new atomic data set with forbidden transitions. We treat collisions with non-thermal leptons resulting from nuclear decays to account for their contribution to excitation, ionization, and heating. We validate our method with a variety of tests including comparing our synthetic nebular spectra for the well-known one-dimensional W7 model with the results of other studies. As an illustrative application of the code, we present synthetic nebular spectra for the detonation of a sub-Chandrasekhar white dwarf (WD) in which the possible effects of gravitational settling of 22Ne prior to explosion have been explored. Specifically, we compare synthetic nebular spectra for a 1.06 M⊙ WD model obtained when 5.5 Gyr of very efficient settling is assumed to a similar model without settling. We find that this degree of 22Ne settling has only a modest effect on the resulting nebular spectra due to increased 58Ni abundance. Due to the high ionization in sub-Chandrasekhar models, the nebular [Ni ii] emission remains negligible, while the [Ni iii] line strengths are increased and the overall ionization balance is slightly lowered in the model with 22Ne settling. In common with previous studies of sub-Chandrasekhar models at nebular epochs, these models overproduce [Fe iii] emission relative to [Fe ii] in comparison to observations of normal SNe Ia.

Description: We analyse the anisotropic clustering of the Sloan Digital Sky Survey-IV Extended Baryon Oscillation Spectroscopic Survey (eBOSS) Luminous Red Galaxy Data Release 14 (DR14) sample combined with Baryon Oscillation Spectroscopic Survey (BOSS) CMASS sample of galaxies in the redshift range 0.6 〈 z 〈 1.0, which consists of 80 118 galaxies from eBOSS and 46 439 galaxies from the BOSS-CMASS sample. The eBOSS-CMASS Luminous Red Galaxy sample has a sky coverage of 1844 deg2, with an effective volume of 0.9 Gpc3. The analysis was made in configuration space using a Legendre multipole expansion. The Redshift Space Distortion signal is modelled as a combination of the Convolution Lagrangian Perturbation Model and the Gaussian Streaming Model. We constrain the logarithmic growth of structure times the amplitude of dark matter density fluctuations, f(zeff)σ8(zeff)? = 0.454? ±? 0.134, and the Alcock-Paczynski dilation scales which constraints the angular diameter distance $D_A(z_{mathrm{ eff}})=1466.5 pm 133.2 (r_s/r_s^{ m fid})$ and $H(z_{ m eff})=105.8 pm 15.7 (r_s^{ m fid}/r_s) mathrm{km, s^{-1}, Mpc^{-1}}$, where rs is the sound horizon at the end of the baryon drag epoch and $r_s^{ m fid}$ is its value in the fiducial cosmology at an effective redshift zeff? = 0.72. These results are in full agreement with the current Λ-Cold Dark Matter (Λ-CDM) cosmological model inferred from Planck measurements. This study is the first eBOSS LRG full-shape analysis i.e. including Redshift Space Distortions simultaneously with the Alcock-Paczynski effect and the Baryon Acoustic Oscillation scale.

Description: Heterogeneity is a fundamental feature of complex phenotypes. So far, genomic screenings have profiled thousands of samples providing insights into the transcriptome of the cell. However, disentangling the heterogeneity of these transcriptomic Big Data to identify defective biological processes remains challenging. Here we present GSECA, a method exploiting the bimodal behavior of RNA-sequencing gene expression profiles to identify altered gene sets in heterogeneous patient cohorts. Using simulated and experimental RNA-sequencing data sets, we show that GSECA provides higher performances than other available algorithms in detecting truly altered biological processes in large cohorts. Applied to 5941 samples from 14 different cancer types, GSECA correctly identified the alteration of the PI3K/AKT signaling pathway driven by the somatic loss of PTEN and verified the emerging role of PTEN in modulating immune-related processes. In particular, we showed that, in prostate cancer, PTEN loss appears to establish an immunosuppressive tumor microenvironment through the activation of STAT3, and low PTEN expression levels have a detrimental impact on patient disease-free survival. GSECA is available at https://github.com/matteocereda/GSECA.

Description: In a consanguineous Pakistani family with two affected individuals, a homozygous variant Gly399Val in the eighth transmembrane domain of the taurine transporter SLC6A6 was identified resulting in a hypomorph transporting capacity of ~15% compared with normal. Three-dimensional modeling of this variant has indicated that it likely causes displacement of the Tyr138 (TM3) side chain, important for transport of taurine. The affected individuals presented with rapidly progressive childhood retinal degeneration, cardiomyopathy and almost undetectable plasma taurine levels. Oral taurine supplementation of 100 mg/kg/day resulted in maintenance of normal blood taurine levels. Following approval by the ethics committee, a long-term supplementation treatment was introduced. Remarkably, after 24-months, the cardiomyopathy was corrected in both affected siblings, and in the 6-years-old, the retinal degeneration was arrested, and the vision was clinically improved. Similar therapeutic approaches could be employed in Mendelian phenotypes caused by the dysfunction of the hundreds of other molecular transporters.

Description: A fluent economical collaboration between countries is a major need. European flows of trade and people are supported by efficient connections between main localities from a geographic region, in many cases overriding national borders. This paper introduces three traveling salesmen problem instances based on freely available geographic coordinates of the main cities of France, Portugal and Spain. These instances are unified, generating other four larger instances: three with all pairs of countries and one instance with the settlements from all the three countries. The study includes an analysis of quality of solutions for a version of branch & cut algorithm and some hybrid heuristics including the Lin–Kernighan algorithm. $Bormathring{u}vka$, Quick$Bormathring{u}vka$ and Greedy algorithms are also used in the hybrid approaches in order to obtain a potential beneficent initial solution for the Lin–Kernighan algorithm. Concorde solver, nowadays state-of-the-art exact software, together with the already mentioned algorithms, is used to test and furthermore analyze the new TSP instances. Some results are represented using online services such as Google Maps, showing the potential integration of the Concorde’s optimum results into commercial routing applications. The very good results provided by the Lin–Kernighan method allow its usage for real medium-sized routing instances.

Description: Specifically tailored industrial control systems (ICSs) attacks are becoming increasingly sophisticated, accentuating the need of ICS cyber security. The nature of these systems makes traditional IT security measures not suitable, requiring expressly developed security countermeasures. Within the past decades, research has been focused in network-based intrusion detection systems. With the appearance of software-defined networks (SDNs), new opportunities and challenges have shown up in the research community. This paper describes the potential benefits of using SDNs in industrial networks with security purposes and presents the set up and results of a pilot experiment carried out in a scaled physical implementation. The experimental set up consists in the detection of ICMP flood and packet payload alteration based on signature comparison. Results point to the potential viability of the technology for intrusion detection and the need of researching in architectural scalability.

Description: DNA damage is the driving force for mutation and genomic instability, which can both lead to cell death or carcinogenesis. DNA double strand breaks are detected and processed in part by the Mre11–Rad50–Nbs1 protein complex. Although the Mre11–Rad50–Nbs1 complex is essential, several spontaneous mutations have been noted in various cancers. One of these mutations, within a conserved motif of Rad50, resulted in an outlier curative response in a clinical trial. We show through biochemical and biophysical characterization that this cancer-associated mutation and a second mutation to the adjacent residue, previously described in a breast cancer patient, both have gain-of-function Rad50 ATP hydrolysis activity that results not from faster association of the ATP-bound form but faster dissociation leading to less stable Rad50 dimer. This disruption impairs the regulatory functions of the protein complex leading to a loss of exonuclease activity from Mre11. Interestingly, these two mutations affect Rad50 structure and dynamics quite differently. These studies describe the relationship between function, structure, and molecular motions in improperly regulated Rad50, which reveal the underlying biophysical mechanism for how these two cancer-associated mutations affect the cell.

Description: In this article, we study the effect of airports on local economic performance that arises from better access to domestic markets in the context of China’s recent airport network expansion. We measure access through the changes in network closeness centrality implied by the contraction in potential journey times between counties within China. Our key finding is that better access—primarily due to landside distance reductions to airports—increased manufacturing productivity. The analysis is carried out on a panel of counties built from micro data on industrial firms, administrative records and census data. To mitigate endogeneity issues, we focus on a subsample of ‘incidentally’ affected counties, whose location midway between existing and new airports implies that they were neither explicitly targeted for development nor directly affected by airport operations.

Description: Solution nuclear magnetic resonance (NMR) experiments allow RNA dynamics to be determined in an aqueous environment. However, when a limited number of peaks are assigned, it is difficult to obtain structural information. We here show a protocol based on the combination of experimental data (Nuclear Overhauser Effect, NOE) and molecular dynamics simulations with enhanced sampling methods. This protocol allows to (a) obtain a maximum entropy ensemble compatible with NMR restraints and (b) obtain a minimal set of metastable conformations compatible with the experimental data (maximum parsimony). The method is applied to a hairpin of 29 nt from an inverted SINEB2, which is part of the SINEUP family and has been shown to enhance protein translation. A clustering procedure is introduced where the annotation of base-base interactions and glycosidic bond angles is used as a metric. By reweighting the contributions of the clusters, minimal sets of four conformations could be found which are compatible with the experimental data. A motif search on the structural database showed that some identified low-population states are present in experimental structures of other RNA transcripts. The introduced method can be applied to characterize RNA dynamics in systems where a limited amount of NMR information is available.

Description: Early life adverse conditions can have major consequences on an individual’s life history traits. Oxidative stress has been hypothesized to be one main mechanism underlying the negative consequences of early life adverse conditions. To test this hypothesis, we restricted the food availability of Seba’s short-tailed bat (Carollia perspicillata) mothers of unweaned pups for 10 days, followed by ad libitum provisioning. We also had a control, unrestricted group. We explored the morphological consequences of dietary restriction during early life by measuring growth rate. We also measured four markers of blood oxidative balance during growth. We assessed the level of cortisol, and its inactive form cortisone, in the hair of the pups at the end of growth. Finally, we monitored survival during the first year. Food restriction triggered a slowdown in growth followed by catch-up growth when ad libitum feeding was restored which did not lead to full compensation in size or mass compared to control individuals. We found that higher growth rate was associated with elevated oxidative damage, suggesting an oxidative cost to growth. However, we found no clear evidence for physiological costs specific to the catch-up growth. Survival after a year was not impacted by the treatment, the oxidative balance or the level of glucocorticoids at the end of growth. In conclusion, our results show that individuals were able to efficiently mitigate the short-term consequences of adverse early life conditions. However, consequences might arise in the long-term, and could impact reproductive success or lifespan.

Description: The flat-spectrum radio quasar Ton 599 attained its highest ever γ-ray flux state during the first week of 2017 November. Observations of the source by the Swift satellite during this period made it possible to generate a simultaneous high flux state broad-band spectral energy distribution (SED). The high flux state activity of Ton 599 is modelled in this work for the first time. We modelled one high flux state and one quiescent state of the source in order to characterize the evolution of SEDs covering the entire dynamic range of γ-ray flux observed by Fermi-LAT. An attempt was made to model the 2017 November state of the source using an external Compton (EC) model in the leptonic scenario. We reproduce the broad-band flaring state SED using a two-component leptonic emission model. We considered one component as an EC+synchrotron self-Compton (SSC) component and the other as pure SSC, lying further down in the jet. The EC+SSC component was located outside the broad-line region (BLR). It mainly reproduces the GeV emission by an EC process with a dusty torus (DT) photon field providing seed photons. We reproduce the broad-band emission from Ton 599 satisfactorily during its peculiar flaring state with a leptonic two-component model. Besides this, we compare the model parameters of a quiescent-state SED with the available average state model parameters in the literature.

Description: Seismic diffractions are ideal carriers of information on small-scale, discontinuous objects and can therefore be used to detect these geologic objects. However, recognizing diffractions is difficult because specular reflection with strong energy masks the weak diffraction. In this study, we propose a diffraction separation and imaging method based on a Mahalanobis-based and phase-based attenuation function used to modify the Kirchhoff migration formula in the full-azimuth dip-angle domain. In this domain, reflections are restricted to within the first Fresnel zone and are distributed in the vicinity of the stationary point, while diffractions are located across a wide range of azimuth and dip angles. Synthetic and field data applications suggest that this new method can effectively separate and image diffractions. The results also demonstrate the efficiency of the new method in clarifying subsurface small-scale objects, which can provide finer information about these structures for seismic interpretation.

Description: We provide new very simple and compact expressions for the efficient calculation of gravitational lens optical scalars for Kerr space–time, which are exact along any null geodesic. These new results are obtained recurring to well-known results on geodesic motion that exploit obvious and hidden symmetries of Kerr space–time and contrast with the rather long and cumbersome expressions previously reported in the literature, providing a helpful improvement for the sake of an efficient integration of the geodesic deviation equation on Kerr geometry. We also introduce a prescription for the observer frame that captures a new notion of centre of the black hole, which can be used for any position of the observer, including those near the black hole. We compare the efficient calculation of weak lens optical scalars with the exact equations, finding an excellent agreement.

Description: Due to complex lithology, strong heterogeneity, low porosity and permeability; resistivity logging faces great challenges in oil saturation prediction of tight conglomerate reservoirs. First, 10 typical core samples were selected to measure and analyse the porosity, permeability, nuclear magnetic resonance (NMR) T2 spectrum and mercury injection capillary pressure (MICP) curve. Second, an empirical method was proposed for reconstructing the NMR T2 spectrum under completely watered conditions using MICP curve based on the ‘three-piece’ power function. The parameters of different models were calibrated via experimental data analysis, respectively. The 180 core experimental data from an MICP curve were used as the input database. Porosity and permeability were regarded as the MICP data selection criteria to apply this model in formation evaluation. The comparison results show good application effects. Finally, to reflect oil saturation, the ratio of T2 geometric means of NMR T2 spectra under oil-bearing and completely watered conditions was proposed. Then, the quantitative relation between oil saturation and the proposed ratio was established via experimental data from the sealed cores, which established a quantitative prediction on oil saturation of tight conglomerate reservoirs. This showed a good application effect. The average relative error and the root mean square error (RMSE) of the predicted oil saturation and sealed coring measurement were around 10 and 3%, respectively. As the proposed method is only influenced by the wettability of reservoir and viscosity of oil, it is not only appropriate for the studied area, but also for other water-wet reservoirs containing light oil. It is important for identifying oil layers, calculating oil saturation and improving log interpretation accuracy in tight conglomerate reservoirs.

Description: IRS 13E is an enigmatic compact group of massive stars located in projection only 3.6 arcsec away from Sgr A*. This group has been suggested to be bounded by an intermediate-mass black hole (IMBH). We present a multiwavelength study of the group and its interplay with the environment. Based on Chandra observations, we find the X-ray spectrum of IRS 13E can be well characterized by an optically thin thermal plasma. The emission peaks between two strongly mass-losing Wolf–Rayet stars of the group. These properties can be reasonably well reproduced by simulated colliding winds of these two stars. However, this scenario underpredicts the X-ray intensity in outer regions. The residual emission likely results from the ram-pressure confinement of the IRS 13E group wind by the ambient medium and is apparently associated with a shell-like warm gas structure seen in Pa α and in ALMA observations. These latter observations also show strongly peaked thermal emission with unusually large velocity spread between the two stars. These results indicate that the group is colliding with the bar of the dense cool gas mini-spiral around Sgr A*. The extended X-ray morphology of IRS 13E and its association with the bar further suggest that the group is physically much farther away than the projected distance from Sgr A*. The presence of an IMBH, while favourable to keep the stars bound together, is not necessary to explain the observed stellar and gas properties of IRS 13E.

Description: Westerlund 1 (Wd1) is potentially the largest star cluster in the Galaxy. That designation critically depends upon the distance to the cluster, yet the cluster is highly obscured, making luminosity-based distance estimates difficult. Using Gaia Data Release 2 (DR2) parallaxes and Bayesian inference, we infer a parallax of $0.35^{+0.07}_{-0.06}$ mas corresponding to a distance of $2.6^{+0.6}_{-0.4}$ kpc. To leverage the combined statistics of all stars in the direction of Wd1, we derive the Bayesian model for a cluster of stars hidden among Galactic field stars; this model includes the parallax zero-point. Previous estimates for the distance to Wd1 ranged from 1.0 to 5.5 kpc, although values around 5 kpc have usually been adopted. The Gaia DR2 parallaxes reduce the uncertainty from a factor of 3 to 18 per cent and rules out the most often quoted value of 5 kpc with 99 per cent confidence. This new distance allows for more accurate mass and age determinations for the stars in Wd1. For example, the previously inferred initial mass at the main-sequence turn-off was around 40 M⊙; the new Gaia DR2 distance shifts this down to about 22 M⊙. This has important implications for our understanding of the late stages of stellar evolution, including the initial mass of the magnetar and the LBV in Wd1. Similarly, the new distance suggests that the total cluster mass is about four times lower than previously calculated.

Description: The Scale-Invariant Vacuum (SIV) theory is based on Weyl’s Integrable Geometry, endowed with a gauge scalar field. The main difference between MOND and the SIV theory is that the first considers a global dilatation invariance of space and time, where the scale factor λ is a constant, while the second opens the likely possibility that λ is a function of time. The key equations of the SIV framework are used here to study the relationship between the Newtonian gravitational acceleration due to baryonic matter gbar and the observed kinematical acceleration gobs. The relationship is applied to galactic systems of the same age where the radial acceleration relation (RAR), between the gobs and gbar accelerations, can be compared with observational data. The SIV theory shows an excellent agreement with observations and with MOND for baryonic gravities gbar > 10−11.5 m s−2. Below this value, SIV still fully agrees with the observations, as well as with the horizontal asymptote of the RAR for dwarf spheroidals, while this is not the case for MOND. These results support the view that there is no need for dark matter and that the RAR and related dynamical properties of galaxies can be interpreted by a modification of gravitation.

Description: A number of scenarios for the formation of multiple populations in globular clusters (GCs) predict that second generation (2G) stars form in a compact and dense subsystem embedded in a more extended first-generation (1G) system. If these scenarios are accurate, a consequence of the denser 2G formation environment is that 2G binaries should be more significantly affected by stellar interactions and disrupted at a larger rate than 1G binaries. The fractions and properties of binary stars can thus provide a dynamical fingerprint of the formation epoch of multiple-population GCs and their subsequent dynamical evolution. We investigate the connection between binaries and multiple populations in five GCs, NGC 288, NGC 6121 (M 4), NGC 6352, NGC 6362, and NGC 6838 (M 71). To do this, we introduce a new method based on the comparison of Hubble Space Telescope observations of binaries in the F275W, F336W, F438W, F606W, and F814W filters with a large number of simulated binaries. In the inner regions probed by our data, we do not find large differences between the local 1G and the 2G binary incidences in four of the studied clusters, the only exception being M 4 where the 1G binary incidence is about three times larger than the 2G incidence. The results found are in general agreement with the results of simulations predicting significant differences in the global 1G and 2G incidences and in the local values in the clusters’ outer regions but similar incidences in the inner regions. The significant difference found in M 4 is consistent with simulations with a larger fraction of wider binaries. Our analysis also provides the first evidence of mixed (1G–2G) binaries, a population predicted by numerical simulations to form in a cluster’s inner regions as a result of stellar encounters during which one component of a binary is replaced by a star of a different population.

Description: Galactic bars are unstable to a vertical buckling instability which heats the disc and in some cases forms a boxy/peanut shaped bulge. We analyse the buckling instability as an application of classical Euler buckling followed by non-linear gravitational Landau damping in the collisionless system. We find that the buckling instability is dictated by the kinematic properties and geometry of the bar. The analytical result is compared to simulations of isolated galaxies containing the disc and dark matter components. Our results demonstrate that violent buckling does not destroy bars while a less energetic buckling can dissolve the bar. The discs that undergo gentle buckling remain stable to bar formation which may explain the observed bar fraction in the local Universe. Our results align with the results from recent surveys.

Description: Galactic outflows are thought to eject baryons back out to the circumgalactic medium. Studies based on metal absorption lines (Mg? ii in particular) in the spectra of background quasars indicate that the gas is ejected anisotropically, with galactic winds likely leaving the host in a bi-conical flow perpendicular to the galaxy disc. In this paper, we present a detailed analysis of an outflow from a z? = 0.7 ‘green-valley’ galaxy [log (M*/M⊙)? = 9.8; $mbox{SFR}=0.5, mathrm{M}_{odot }, mathrm{yr}^{-1}$] probed by two background sources from the MusE GAs FLOw and Wind (MEGAFLOW) survey. Thanks to a fortuitous configuration with a background quasar (SDSSJ1358? + 1145) and a bright background galaxy at z? = 1.4, both at impact parameters of $approx! 15, hbox{kpc}$, we can – for the first time – probe both the receding and approaching components of a putative galactic outflow around a distant galaxy. We measure a significant velocity shift between the Mg? ii absorption from the two sightlines ($84pm 17, hbox{km~s$^{-1}$}$), which is consistent with the expectation from our simple fiducial wind model, possibly combined with an extended disc contribution.

Description: Low-temperature plastic rheology of calcite plays a significant role in the dynamics of Earth's crust. However, it is technically challenging to study plastic rheology at low temperatures because of the high confining pressures required to inhibit fracturing. Micromechanical tests, such as nanoindentation and micropillar compression, can provide insight into plastic rheology under these conditions because, due to the small scale, plastic deformation can be achieved at low temperatures without the need for secondary confinement. In this study, nanoindentation and micropillar compression experiments were performed on oriented grains within a polycrystalline sample of Carrara marble at temperatures ranging from 23 to 175 °C, using a nanoindenter. Indentation hardness is acquired directly from nanoindentation experiments. These data are then used to calculate yield stress as a function of temperature using numerical approaches that model the stress state under the indenter. Indentation data are complemented by uniaxial micropillar compression experiments. Cylindrical micropillars ∼1 and ∼3 μm in diameter were fabricated using a focused ion beam-based micromachining technique. Yield stress in micropillar experiments is determined directly from the applied load and micropillar dimensions. Mechanical data are fit to constitutive flow laws for low-temperature plasticity and compared to extrapolations of similar flow laws from high-temperature experiments. This study also considered the effects of crystallographic orientation on yield stress in calcite. Although there is a clear orientation dependence to plastic yielding, this effect is relatively small in comparison to the influence of temperature.

Description: Across the boreal forest in North America, the black spruce (Picea mariana) cover type is ecologically and economically important, occupying roughly 10 percent of Minnesota’s, USA 17.4 million acres (7.0 million hectares) of forestland. Traditionally managed through clearcut regeneration harvests, alternative silvicultural systems are being increasingly used in Canada. Here, we examine the 10- and 57-year effects of six silvicultural treatments (clearcut strips, clearcut patches, thinning, group selection, single-tree selection, shelterwood) on stand structure and dynamics in lowland black spruce. Treatments were installed in 1948 in northern Minnesota, and remeasured and re-treated 10 years later. A subset of the clearcut strips, clearcut patches, and shelterwood treatments were remeasured in 2017. After 10 years, diameter growth of residual stems varied by treatment, with the shelterwood experiencing the greatest growth, and basal area increased in all but the shelterwood treatment. Over the long term, the shelterwood exhibited larger diameters and heights and greater crown ratios, basal area, structural complexity, and compositional diversity than the clearcuts. Our results suggest that managers may consider using a shelterwood instead of traditional large clearcuts to achieve increased structural and compositional diversity, particularly when eastern spruce dwarf mistletoe (Arceuthobium pusillum) does not necessitate a traditional clearcut.

Description: Summary Multiple databases provide valuable information about curated pathways and other resources that can be used to build and analyze networks. OmniPath combines 61 (and continuously growing) network resources into a comprehensive collection, with over 120 000 interactions. We present here the OmniPath App, a Cytoscape plugin to flexibly import data from OmniPath via a simple and intuitive interface. Thus, it makes possible to directly access the large body of high-quality knowledge provided by OmniPath within Cytoscape for inspection and further use with other tools. Availability and implementation The OmniPath App has been developed for Cytoscape 3 in the Java programing language. The latest source code and the plugin can be found at: https://github.com/saezlab/Omnipath_Cytoscape and http://apps.cytoscape.org/apps/omnipath, respectively. Supplementary information Supplementary data are available at Bioinformatics online.

Description: The distinctive security features of the Spanish electronic national identity card, known as Documento Nacional de Identidad electrónico, allow us to propose the usage of this cryptographic smart card in an authentication framework that can be used during the registration and login phases of internet services where the validation of the user’s age and real identity are key elements, as it is the case for example of the so-called social networks. Using this mechanism with NFC-capable devices, the identity and age of the potential user can be determined, allowing or denying the access to the service based on that information.

Description: We use the unsupervised and supervised neural network methods together to predict lithology of a gas hydrate reservoir from downhole data in the Krishna–Godavari (KG) offshore basin, India. In this study, we successfully identify the host litho-units of gas hydrate and show its effects in the identification of lithology using neural network techniques, which is not reported earlier. We use well log data acquired at three holes (10A, 03A and 04A) in 2006 during the first expedition of the Indian National Gas Hydrate Program (NGHP-01). Five different logging while drilling data (e.g. density, neutron porosity, gamma ray, resistivity and sonic) are considered for the mapping of lithology and gas hydrate. In the presence of gas hydrate, the resistivity and sonic velocity of the host sediments increase significantly, whereas density, neutron porosity and gamma ray are negligibly affected. Therefore, we calculate resistivity and sonic velocity for water-saturated sediment (without gas hydrate) theoretically to remove the effects of gas hydrate. At first, we apply the seven unsupervised classification methods (i.e. elbow, dendrogram, K-means, 3-D clustering, principal component analysis, Devies–Bouldin index and self-organizing map) to the data with gas hydrate (e.g. observed) and without gas hydrate (i.e. water-saturated/theoretical) to assess the data dimensionality and the number of clusters/litho-units. Each of the unsupervised schemes has its own pros and cons, and may provide different number of cluster/litho-units; sometimes, it is difficult to interpret from only one method. However, all seven methods provide same number of clusters in our study. Then, we apply the supervised classification method (i.e. Bayesian neural networks optimized by hybrid Monte Carlo searching technique) to the training data to refine the defined litho-units and map them with depth. Our approach identifies four types of litho-units and illustrates that the lithology in this area is dominated by clay (∼64 per cent) with some amount of silty clay, silt and minor sand. Gas hydrate is found in clay, silty clay and silt and not in sand. Results also show that, if gas hydrate is not considered as a separate unit, it is distributed as lithology in its hosts (i.e. clay, silty clay and silt). The method is very stable up to ∼15 per cent of random noise added to the data and results are well matched with the analysis of recovered core data. Identified lithologies at three wells correlate very well with seismic section crossing the wells. Very low permeability (〈0.1 mD) estimated at three wells also indicates the clay-dominated lithology in our study area.

Description: Avian brood parasites depend entirely on their hosts to raise their nestlings until independence. Thus, parasite females should select suitable host nests for egg laying according to traits that enhance offspring survival. The availability of nests of certain characteristics influencing the survival of parasitic offspring is, however, temporally dynamic and, thus, patterns of host selection should be evaluated considering characteristics of available host nests the day of parasitism. This allows detecting possible seasonal changes and, therefore, a more realistic picture of host selection by brood parasites. In this paper, we adopt such a new approach and consider daily availability of magpie (Pica pica) host nests at different breeding stage that were or were not parasitized by the great spotted cuckoo (Clamator glandarius). Theory predicts that cuckoos should select host nests at the laying stage. Accordingly, we detected that cuckoos preferred to parasitize magpie nests at the laying stage but, mainly, those that already harbored one or two cuckoo eggs, which may seem counterintuitive. We also showed that patterns of host selection by cuckoos varied during the breeding season, which implies that brood parasite–host interaction is dynamic depending on phenology. These patterns are hidden when not considering the temporally dynamic nature of the availability of host nests of characteristics of interest. We discuss the importance of such patterns and considering diary hosts nests availability for detecting them.

Description: A review, recently published in this journal by Fang (2019), showed that methods trained for the prediction of protein stability changes upon mutation have a very critical bias: they neglect that a protein variation (A- 〉 B) and its reverse (B- 〉 A) must have the opposite value of the free energy difference (ΔΔGAB = − ΔΔGBA). In this letter, we complement the Fang’s paper presenting a more general view of the problem. In particular, a machine learning-based method, published in 2015 (INPS), addressed the bias issue directly. We include the analysis of the missing method, showing that INPS is nearly insensitive to the addressed problem.

Description: Life-history theory predicts that animals adjust their resource allocation to somatic growth or to reproduction to maximize fitness. Resource allocation in Daphnia is known to respond to quantitative food limitation as well as to kairomones released from predators. Here we investigated in a full-factorial design how kairomone from larvae of Chaoborus flavicans, a gape-limited predator, and food quantity (0.5 mg C/L versus 1.5 mg C/L) affect the fatty acid allocation of D. pulex. Low food diminished somatic growth, clutch size and clutch biomass and increased neckteeth formation in response to the kairomone. Low food further led to increased fatty acid amounts per individual egg as well as to increased fatty acid content in eggs and to increased relative fatty acid allocation to reproduction. The latter effect was suppressed by kairomone of Chaoborus, whereas on high food the provision of eggs was further enhanced. We also found that more eicosapentaenoic acid was retained in the body of mothers in the presence of the predator at low food concentrations. These findings indicate that under food limitation and in the presence of kairomone from Chaoborus larvae, Daphnia switches from allocation into current reproduction to investment into future reproductive events.

Description: I aim to identify the presence of spatial interactions among local governments by exploiting a novel strategy. Specifically, I take advantage of the political cycle of Italian municipalities over the period of 2001–2011 to isolate the effect of the spending decisions of one municipality on neighboring municipalities. The results of this analysis point to the presence of strategic interactions between neighboring municipalities and indicate that such fiscal behavior is more pronounced during electoral years compared to non-electoral ones, when municipalities are governed by coalitions backed by a small majority, and in cities guided by a mayor who can run for re-election. Taken together, these results suggest that the observed spatial dependence in spending decisions seems to be consistent with the yardstick competition hypothesis.