Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “BLAST”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 991 records · Page 55
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Parasail: SIMD C library for global, semi-global, and local pairwise sequence alignments

Sequence alignment algorithms are a key component of many bioinformatics applications. Though various fast Smith-Waterman local sequence alignment implementations have been developed for x86 CPUs, most are embedded into larger database search tools. In addition, fast implementations of Needleman-Wunsch global sequence alignment and its semi-global variants are not as widespread. This article presents the first software library for local, global, and semi-global pairwise intra-sequence alignments and improves the performance of previous intra-sequence implementations. As a result, a faster intra-sequence pairwise alignment implementation is described and benchmarked. Using a 375 residue query sequence a speed of 136 billion cell updates per second (GCUPS) was achieved on a dual Intel Xeon E5-2670 12-core processor system, the highest reported for an implementation based on Farrar’s ’striped’ approach. When using only a single thread, parasail was 1.7 times faster than Rognes’s SWIPE. For many score matrices, parasail is faster than BLAST. The software library is designed for 64 bit Linux, OS X, or Windows on processors with SSE2, SSE41, or AVX2. Source code is available from https://github.com/jeffdaily/parasail under the Battelle BSD-style license. In conclusion, applications that require optimal alignment scores could benefit from the improved performance. For the first time, SIMD global, semi-global, and local alignments are available in a stand-alone C library.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

An FPGA-based hardware accelerator supporting sensitive sequence homology filtering with profile hidden Markov models

Abstract Background Sequence alignment lies at the heart of genome sequence annotation. While the BLAST suite of alignment tools has long held an important role in alignment-based sequence database search, greater sensitivity is achieved through the use of profile hidden Markov models (pHMMs). Here, we describe an FPGA hardware accelerator, called HAVAC, that targets a key bottleneck step (SSV) in the analysis pipeline of the popular pHMM alignment tool, HMMER. Results The HAVAC kernel calculates the SSV matrix at 1739 GCUPS on a $$\sim$$ ∼ $3000 Xilinx Alveo U50 FPGA accelerator card, $$\sim$$ ∼ 227× faster than the optimized SSV implementation in nhmmer . Accounting for PCI-e data transfer data processing, HAVAC is 65× faster than nhmmer’s SSV with one thread and 35× faster than nhmmer with four threads, and uses $$\sim$$ ∼ 31% the energy of a traditional high end Intel CPU. Conclusions HAVAC demonstrates the potential offered by FPGA hardware accelerators to produce dramatic speed gains in sequence annotation and related bioinformatics applications. Because these computations are performed on a co-processor, the host CPU remains free to simultaneously compute other aspects of the analysis pipeline.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2015

Antifungal wound penetration of amphotericin and voriconazole in combat-related injuries: case report

Background: Survivors of combat trauma can have long and challenging recoveries, which may be complicated by infection. Invasive fungal infections are a rare but serious complication with limited treatment options. Currently, aggressive surgical debridement is the standard of care, with antifungal agents used adjunctively with uncertain efficacy. Anecdotal evidence suggests that antifungal agents may be ineffective in the absence of surgical debridement, and studies have yet to correlate antifungal concentrations in plasma and wounds. Case presentation: Here we report the systemic pharmacokinetics and wound effluent antifungal concentrations of five wounds from two male patients, aged 28 and 30 years old who sustained combat-related blast injuries in southern Afghanistan, with proven or possible invasive fungal infection. Our data demonstrate that while voriconazole sufficiently penetrated the wound resulting in detectable effluent levels, free amphotericin B (unbound to plasma) was not present in wound effluent despite sufficient concentrations in circulating plasma. In addition, considerable between-patient and within-patient variability was observed in antifungal pharmacokinetic parameters. Conclusion: These data highlight the need for further studies evaluating wound penetration of commonly used antifungals and the role for therapeutic drug monitoring in providing optimal care for critically ill and injured war fighters.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2018

The orphan nuclear receptor EAR-2 (NR2F6) inhibits hematopoietic cell differentiation and induces myeloid dysplasia in vivo

Background: In patients with myelodysplastic syndrome (MDS), bone marrow cells have an increased predisposition to apoptosis, yet MDS cells outcompete normal bone marrow (BM)–suggesting that factors regulating growth potential may be important in MDS. We previously identified v-Erb A related-2 (EAR-2, NR2F6) as a gene involved in control of growth ability. Methods: Bone marrow obtained from C57BL/6 mice was transfected with a retrovirus containing EAR-2-IRES-GFP. Ex vivo transduced cells were flow sorted. In some experiments cells were cultured in vitro, in other experiments cells were injected into lethally irradiated recipients, along with non-transduced bone marrow cells. Short-hairpin RNA silencing EAR-2 was also introduced into bone marrow cells cultured ex vivo. Results: Here, we show that EAR-2 inhibits maturation of normal BM in vitro and in vivo and that EAR-2 transplant chimeras demonstrate key features of MDS. Competitive repopulation of lethally irradiated murine hosts with EAR2-transduced BM cells resulted in increased engraftment and increased colony formation in serial replating experiments. Recipients of EAR-2-transduced grafts had hypercellular BM, erythroid dysplasia, abnormal localization of immature precursors and increased blasts; secondary transplantation resulted in acute leukemia. Animals were cytopenic, having reduced numbers of erythrocytes, monocytes and granulocytes. Suspension culture confirmed that EAR-2 inhibits granulocytic and monocytic differentiation, while knockdown induced granulocytic differentiation. We observed a reduction in the number of BFU-E and CFU-GM colonies and the size of erythroid and myeloid colonies. Serial replating of transduced hematopoietic colonies revealed extended replating potential in EAR-2-overexpressing BM, while knockdown reduced re-plating ability. EAR-2 functions by recruitment of histone deacetylases, and inhibition of differentiation in 32D cells is dependent on the DNA binding domain. Conclusions: This data suggest that NR2F6 inhibits maturation of normal BM in vitro and in vivo and that the NR2F6 transplant chimera system demonstrates key features of MDS, and could provide a mouse model for MDS.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

High quality draft genome sequences of Pseudomonas fulva DSM 17717 T , Pseudomonas parafulva DSM 17004 T and Pseudomonas cremoricolorata DSM 17059 T type strains

Pseudomonas has the highest number of species out of any genus of Gram-negative bacteria and is phylogenetically divided into several groups. The Pseudomonas putida phylogenetic branch includes at least 13 species of environmental and industrial interest, plant-associated bacteria, insect pathogens, and even some members that have been found in clinical specimens. In the context of the Genomic Encyclopedia of Bacteria and Archaea project, we present the permanent, high-quality draft genomes of the type strains of 3 taxonomically and ecologically closely related species in the Pseudomonas putida phylogenetic branch: Pseudomonas fulva DSM 17717 T , Pseudomonas parafulva DSM 17004 T and Pseudomonas cremoricolorata DSM 17059T. All three genomes are comparable in size (4.6-4.9Mb), with 4,119-4,459 protein-coding genes. Average nucleotide identity based on BLAST comparisons and digital genome-to-genome distance calculations are in good agreement with experimental DNA-DNA hybridization results. The genome sequences presented here will be very helpful in elucidating the taxonomy, phylogeny and evolution of the Pseudomonas putida species complex.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2012

Using Comparative Genomics for Inquiry-Based Learning to Dissect Virulence of Escherichia coli O157:H7 and Yersinia pestis

Genomics and bioinformatics are topics of increasing interest in undergraduate biological science curricula. Many existing exercises focus on gene annotation and analysis of a single genome. In this paper, we present two educational modules designed to enable students to learn and apply fundamental concepts in comparative genomics using examples related to bacterial pathogenesis. Students first examine alignments of genomes of Escherichia coli O157:H7 strains isolated from three food-poisoning outbreaks using the multiple-genome alignment tool Mauve. Students investigate conservation of virulence factors using the Mauve viewer and by browsing annotations available at the A Systematic Annotation Package for Community Analysis of Genomes database. In the second module, students use an alignment of five Yersinia pestis genomes to analyze single-nucleotide polymorphisms of three genes to classify strains into biovar groups. Students are then given sequences of bacterial DNA amplified from the teeth of corpses from the first and second pandemics of the bubonic plague and asked to classify these new samples. Learning-assessment results reveal student improvement in self-efficacy and content knowledge, as well as students' ability to use BLAST to identify genomic islands and conduct analyses of virulence factors from E. coli O157:H7 or Y. pestis. Each of these educational modules offers educators new ready-to-implement resources for integrating comparative genomic topics into their curricula.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

Systems-Wide Prediction of Enzyme Promiscuity Reveals a New Underground Alternative Route for Pyridoxal 5’-Phosphate Production in E. coli

Recent insights suggest that non-specific and/or promiscuous enzymes are common and active across life. Understanding the role of such enzymes is an important open question in biology. Here we develop a genome-wide method, PROPER, that uses a permissive PSI-BLAST approach to predict promiscuous activities of metabolic genes. Enzyme promiscuity is typically studied experimentally using multicopy suppression, in which over-expression of a promiscuous ‘replacer’ gene rescues lethality caused by inactivation of a ‘target’ gene. We use PROPER to predict multicopy suppression in Escherichia coli, achieving highly significant overlap with published cases (hypergeometric p = 4.4e-13). We then validate three novel predicted target-replacer gene pairs in new multicopy suppression experiments. We next go beyond PROPER and develop a network-based approach, GEM-PROPER, that integrates PROPER with genome-scale metabolic modeling to predict promiscuous replacements via alternative metabolic pathways. GEM-PROPER predicts a new indirect replacer (thiG) for an essential enzyme (pdxB) in production of pyridoxal 5’-phosphate (the active form of Vitamin B 6 ), which we validate experimentally via multicopy suppression. Here, we perform a structural analysis of thiG to determine its potential promiscuous active site, which we validate experimentally by inactivating the pertaining residues and showing a loss of replacer activity. Thus, this study is a successful example where a computational investigation leads to a network-based identification of an indirect promiscuous replacement of a key metabolic enzyme, which would have been extremely difficult to identify directly.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2013

Genome analysis and physiological comparison of Alicycliphilus denitrificans strains BC and K601T

The genomes of the Betaproteobacteria Alicycliphilus denitrificans strains BC and K601T have been sequenced to get insight into the physiology of the two strains. Strain BC degrades benzene with chlorate as electron acceptor. The cyclohexanol-degrading denitrifying strain K601T is not able to use chlorate as electron acceptor, while strain BC cannot degrade cyclohexanol. The 16S rRNA sequences of strains BC and K601T are identical and the fatty acid methyl ester patterns of the strains are similar. Basic Local Alignment Search Tool (BLAST) analysis of predicted open reading frames of both strains showed most hits with Acidovorax sp. JS42, a bacterium that degrades nitro-aromatics. The genomes include strain-specific plasmids (pAlide201 in strain K601T and pAlide01 and pAlide02 in strain BC). Key genes of chlorate reduction in strain BC were located on a 120 kb megaplasmid (pAlide01), which was absent in strain K601T. Genes involved in cyclohexanol degradation were only found in strain K601T. Benzene and toluene are degraded via oxygenase-mediated pathways in both strains. Genes involved in the meta-cleavage pathway of catechol are present in the genomes of both strains. Strain BC also contains all genes of the ortho-cleavage pathway. The large number of mono- and dioxygenase genes in the genomes suggests that the two strains have a broader substrate range than known thus far.

Oosterkamp, Margreet J.↗
DOE OSTI2019

Assessment of BOLD and GenBank – Their accuracy and reliability for the identification of biological materials

Taxonomic identification of biological materials can be achieved through DNA barcoding, where an unknown “barcode” sequence is compared to a reference database. In many disciplines, obtaining accurate taxonomic identifications can be imperative (e.g., evolutionary biology, food regulatory compliance, forensics). The Barcode of Life DataSystems (BOLD) and GenBank are the main public repositories of DNA barcode sequences. In this study, an assessment of the accuracy and reliability of sequences in these databases was performed. To achieve this, 1) curated reference materials for plants, macro-fungi and insects were obtained from national collections, 2) relevant barcode sequences (rbcL, matK, trnH-psbA, ITS and COI) from these reference samples were generated and used for searching against both databases, and 3) optimal search parameters were determined that ensure the best match to the known species in either database. While GenBank outperformed BOLD for species-level identification of insect taxa (53% and 35%, respectively), both databases performed comparably for plants and macro-fungi (~81% and ~57%, respectively). Results illustrated that using a multi-locus barcode approach increased identification success. This study outlines the utility of the BLAST search tool in GenBank and the BOLD identification engine for taxonomic identifications and identifies some precautions needed when using public sequence repositories in applied scientific disciplines.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2012

A P-loop Mutation in G[alpha] Subunits Prevents Transition to the Active State: Implications for G-protein Signaling in Fungal Pathogenesis

Heterotrimeric G-proteins are molecular switches integral to a panoply of different physiological responses that many organisms make to environmental cues. The switch from inactive to active G{alpha}{beta}{gamma} heterotrimer relies on nucleotide cycling by the G{alpha} subunit: exchange of GTP for GDP activates G{alpha}, whereas its intrinsic enzymatic activity catalyzes GTP hydrolysis to GDP and inorganic phosphate, thereby reverting G{alpha} to its inactive state. In several genetic studies of filamentous fungi, such as the rice blast fungus Magnaporthe oryzae, a G42R mutation in the phosphate-binding loop of G{alpha} subunits is assumed to be GTPase-deficient and thus constitutively active. Here, we demonstrate that G{alpha}(G42R) mutants are not GTPase deficient, but rather incapable of achieving the activated conformation. Two crystal structure models suggest that Arg-42 prevents a typical switch region conformational change upon G{alpha}{sub i1}(G42R) binding to GDP {center_dot} AlF{sub 4}{sup -} or GTP, but rotameric flexibility at this locus allows for unperturbed GTP hydrolysis. G{alpha}(G42R) mutants do not engage the active state-selective peptide KB-1753 nor RGS domains with high affinity, but instead favor interaction with G{beta}{gamma} and GoLoco motifs in any nucleotide state. The corresponding G{alpha}{sub q}(G48R) mutant is not constitutively active in cells and responds poorly to aluminum tetrafluoride activation. Comparative analyses of M. oryzae strains harboring either G42R or GTPase-deficient Q/L mutations in the G{alpha} subunits MagA or MagB illustrate functional differences in environmental cue processing and intracellular signaling outcomes between these two G{alpha} mutants, thus demonstrating the in vivo functional divergence of G42R and activating G-protein mutants.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

Reconstructing the backbone of the Saccharomycotina yeast phylogeny using genome-scale data

Understanding the phylogenetic relationships among the yeasts of the subphylum Saccharomycotina is a prerequisite for understanding the evolution of their metabolisms and ecological lifestyles. In the last two decades, the use of rDNA and multilocus data sets has greatly advanced our understanding of the yeast phylogeny, but many deep relationships remain unsupported. In contrast, phylogenomic analyses have involved relatively few taxa and lineages that were often selected with limited considerations for covering the breadth of yeast biodiversity. Here we used genome sequence data from 86 publicly available yeast genomes representing nine of the 11 known major lineages and 10 nonyeast fungal outgroups to generate a 1233-gene, 96-taxon data matrix. Species phylogenies reconstructed using two different methods (concatenation and coalescence) and two data matrices (amino acids or the first two codon positions) yielded identical and highly supported relationships between the nine major lineages. Aside from the lineage comprised by the family Pichiaceae, all other lineages were monophyletic. Most interrelationships among yeast species were robust across the two methods and data matrices. Furthermore, eight of the 93 internodes conflicted between analyses or data sets, including the placements of: the clade defined by species that have reassigned the CUG codon to encode serine, instead of leucine; the clade defined by a whole genome duplication; and the species Ascoidea rubescens . These phylogenomic analyses provide a robust roadmap for future comparative work across the yeast subphylum in the disciplines of taxonomy, molecular genetics, evolutionary biology, ecology, and biotechnology. To further this end, we have also provided a BLAST server to query the 86 Saccharomycotina genomes, which can be found at http://y1000plus.org/blast.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

Moment tensor analysis of very shallow sources

An issue for moment tensor (MT) inversion of shallow seismic sources is that some components of the Green’s functions have vanishing amplitudes at the free surface, which can result in bias in the MT solution. The effects of the free surface on the stability of the MT method become important as we continue to investigate and improve the capabilities of regional full MT inversion for source–type identification and discrimination. It is important to understand free–surface effects on discriminating shallow explosive sources for nuclear monitoring purposes. It may also be important in natural systems that have very shallow seismicity, such as volcanic and geothermal systems. We examine the effects of the free surface on the MT via synthetic testing and apply the MT–based discrimination method to three quarry blasts from the HUMMING ALBATROSS experiment. These shallow chemical explosions at ~10 m depth and recorded up to several kilometers distance represent rather severe source–station geometry in terms of free–surface effects. We show that the method is capable of recovering a predominantly explosive source mechanism, and the combined waveform and first–motion method enables the unique discrimination of these events. Furthermore, recovering the design yield using seismic moment estimates from MT inversion remains challenging, but we can begin to put error bounds on our moment estimates using the network sensitivity solution technique.

58 GEOSCIENCES↗
DOE OSTI2018

Depth Discrimination Using Rg-to-Sg Spectral Amplitude Ratios for Seismic Events in Utah Recorded at Local Distances

Short-period fundamental-mode Rayleigh waves (Rg) are commonly observed on seismograms of anthropogenic seismic events and shallow, naturally occurring tectonic earthquakes (TEs) recorded at local distances. In the Utah region, strong Rg waves traveling with an average group velocity of about 1.8 km/s are observed at ~1 Hz on waveforms from shallow events ( depth<10 km ) recorded at distances up to about 150 km. At these distances, Sg waves, which are direct shear waves traveling in the upper crust, are generally the dominant signals for TEs. Here in this study, we leverage the well-known notion that Rg amplitude decreases dramatically with increasing event depth to propose a new depth discriminant based on Rg-to-Sg spectral amplitude ratios. The approach is successfully used to discriminate shallow events (both earthquakes and anthropogenic events) from deeper TEs in the Utah region recorded at local distances ( <150 km ) by the University of Utah Seismographic Stations (UUSS) regional seismic network. Using Mood’s median test, we obtained probabilities of nearly zero that the median Rg-to-Sg spectral amplitude ratios are the same between shallow events on the one hand (including both shallow TEs and anthropogenic events), and deeper earthquakes on the other, suggesting that there is a statistically significant difference in the estimated Rg-to-Sg ratios between the two populations. We also observed consistent disparities between the different types of shallow events (e.g., mining blasts vs. mining-induced earthquakes), implying that it may be possible to separate the subpopulations that make up this group. Lastly, this suggests that using local distance Rg-to-Sg spectral amplitude ratios one can not only discriminate shallow events from deeper events but may also be able to discriminate among different populations of shallow events.

58 GEOSCIENCES↗
DOE OSTI2018

Discrimination of Anthropogenic Events and Tectonic Earthquakes in Utah Using a Quadratic Discriminant Function Approach with Local Distance Amplitude Ratios

Most of the commonly used seismic discrimination approaches are designed for teleseismic and regional data. To monitor for the smallest events, some of these discriminants have been adapted for local distances (< 200 km), with mixed level of success. For this, we take advantage of the variety of seismic sources, including nontraditionally studied anthropogenic sources and the existence of a dense regional seismic network in the Utah region to evaluate amplitude ratio seismic discrimination at local distances. First, we explored phase-amplitude Pg-to-Sg ratios for multiple frequency bands to classify events in a dataset that comprises populations of single-shot surface explosions, shallow and deep ripple-fired mining blasts, mining-induced events (MIEs), and tectonic earthquakes. We achieved a success rate of about 59%–83%. Then, for the same dataset, we combined the Pg-to-Sg phase-amplitude ratios with Sg-to-Rg spectral amplitude ratios in a multivariate quadratic discriminant function (QDF) approach. For two-category pairwise classification, seven of ten population pairs show misclassification rates of about 20% or less, with five pairs showing rates of about 10% or less. The approach performs best for the pair involving the populations of single-shot explosions and MIEs. By combining both Pg-to-Sg and Rg-to-Sg ratios in the multivariate QDFs, we are able to achieve an average improvement of about 4%–14% in misclassification rates compared with Pg-to-Sg ratios alone. When all five event populations are considered simultaneously, as expected, the potential for misclassification increases, and our QDF approach using both Pg-to-Sg and Rg-to-Sg ratios achieves an average success rate of about 74% compared with the rate of about 86% for two-category pairwise classification.

58 GEOSCIENCES↗
DOE OSTI2019

Classification of Local Seismic Events in the Utah Region: A Comparison of Amplitude Ratio Methods with a Spectrogram-Based Machine Learning Approach

Abstract The capability to discriminate low-magnitude earthquakes from low-yield anthropogenic sources, both detectable only at local distances, is of increasing interest to the event monitoring community. Here, we used a dataset of seismic events in Utah recorded during a 14-day period (1–14 January 2011) by the University of Utah Seismic Stations network to perform a comparative study of event classification at local scale using amplitude ratio (AR) methods and a machine learning (ML) approach. The event catalog consists of 7377 events with magnitudes M C ranging from –2 and lower up to 5.8. Events were subdivided into six populations based on location and source type: tectonic earthquakes (TEs), mining-induced events (MIEs), and mining blasts from four known mines (WMB, SMB, LMB, and CQB). The AR approach jointly exploits Pg-to-Sg phase ARs and Rg-to-Sg spectral ARs in multivariate quadratic discriminant functions and was able to classify 370 events with high signal quality from the three groups with sufficient size (TE, MIE, and SMB). For that subset of the events, the method achieved success rates between about 80% and 90%. The ML technique used trained convolutional neural network (CNN) models to classify the populations. The CNN approach was able to classify the subset of events with accuracies between about 91% and 98%. Because the neural network approach does not have a minimum signal quality requirement, we applied it to the entire event catalog, including the abundant extremely low-magnitude events, and achieved accuracies of about 94%–100%. Furthermore, we compare the AR and ML methodologies using a broad set of criteria and conclude that a major advantage to ML methods is their robustness to low signal-to-noise ratio data, allowing them to classify significantly smaller events.

58 GEOSCIENCES↗
DOE OSTI2020

Discrimination of Small Earthquakes and Buried Single-Fired Chemical Explosions at Local Distances (<150 km) in the Western United States from Comparison of Local Magnitude ( M L ) and Coda Duration Magnitude ( M C )

Seismologists distinguish underground nuclear explosions from more commonly occurring earthquakes using moment tensor inversion, high-frequency P/S amplitude ratios, m b :M s comparisons, and P–pP differential travel times. These methods are generally successful for large seismic events (M >3-4) well recorded at regional-to-teleseismic distances (>150 km); however, it is unclear whether they can be modified to work for small events (M <3) well recorded only at local distances (<150 km). Here, we evaluate a recently proposed, local-distance seismic source discriminant—the difference between local magnitude (M L ) and coda duration magnitude (M C )—using seismograms of earthquakes and buried, single-fired chemical explosions recorded in three regions of the western United States. The quantity M L –M C was previously found to be sensitive to source depth, effectively discriminating mine blasts, induced earthquakes, and very shallow tectonic earthquakes from deeper crustal earthquakes. In this study, we report the first evaluation of M L –M C as a depth discriminant using data from buried, single-fired explosions that, unlike the seismic sources studied earlier, are good analogs for underground nuclear explosions. We find that even when using generic, uncalibrated methods of assigning magnitudes, M L –M C separates single-fired explosions and earthquakes. The area under the receiver operating characteristic curve is 0.92 for 19 explosions and 14 earthquakes in Washington, 0.90 for 22 explosions and 90 earthquakes in Wyoming, and 0.99 for three explosions and 149 earthquakes in Nevada. Finally, M L :M C comparisons have the potential to enhance discrimination based on high-frequency P/S amplitudes ratios—which perform less well at local than regional distances—because the two metrics have complementary sensitivities.

58 GEOSCIENCES↗
DOE OSTI2021

Local-Distance Seismic Event Relocation and Relative Magnitude Estimation, Applications to Mining Related Seismicity in the Powder River Basin, Wyoming

Recent efforts to characterize small ( M w < 3 ) seismic events at local distances have become more important because of the increased observation of human-triggered and induced seismicity and the need to advance nuclear explosion monitoring capabilities. The signals generated by low-magnitude seismic sources necessitate the use of nearby short-period observations, which are sensitive to local geological heterogeneity. Local to near-regional distance ( < 300 km ) surface and shear waves can dominate short-period observations from small, shallow seismic sources. In this work, we utilize these observations to estimate precise, relative locations and magnitudes of ~ 700 industrial mining events in Wyoming, using nearly 360,000 observations. The precise, relative location estimates (with formal location uncertainty estimates of less than 1 km) collapse a diffuse collection of mining events into discrete clusters associated with individual blasting operations. We also invert the cross-correlation amplitudes to estimate precise, relative moment magnitude estimates, which help validate and identify disparities in the event sizes reported by regional network catalogs. Joint use of multiple phases allows for the inclusion of more seismic events due to the increase in the number of observations. In some cases, using a single phase allowed us to relocate only 50% of the original reported seismic events within a cluster. Combining shear- and surface-wave phases increased the number of events to above 90% of the original events, allowing us to characterize a broader range of event sizes, source to station distances, and event distributions. This analysis takes a step toward making a fuller characterization of small industrial seismic events observed at local distances.

58 GEOSCIENCES↗
DOE OSTI2024

Inferring the Focal Depths of Small Earthquakes in Southern California Using Physics-Based Waveform Features

Determining the depths of small crustal earthquakes is challenging in many regions of the world, because most seismic networks are too sparse to resolve trade-offs between depth and origin time with conventional arrival-time methods. Precise and accurate depth estimation is important, because it can help seismologists discriminate between earthquakes and explosions, which is relevant to monitoring nuclear test ban treaties and producing earthquake catalogs that are uncontaminated by mining blasts. Here, we examine the depth sensitivity of several physics-based waveform features for ~8000 earthquakes in southern California that have well-resolved depths from arrival-time inversion. We focus on small earthquakes (2 < M L < 4) recorded at local distances (<150 km), for which depth estimation is especially challenging. We find that differential magnitudes (M w /M L –M c ) are positively correlated with focal depth, implying that coda wave excitation decreases with focal depth. We analyze a simple proxy for relative frequency content, Φ≡log 10 (M 0 )+3log 10 (f c ), and find that source spectra are preferentially enriched in high frequencies, or “blue-shifted,” as focal depth increases. Here, we also find that two spectral amplitude ratios Rg 0.5–2 Hz/Sg 0.5–8 Hz and Pg/Sg at 3–8 Hz decrease as focal depth increases. Using multilinear regression with these features as predictor variables, we develop models that can explain 11%–59% of the variance in depths within 10 subregions and 25% of the depth variance across southern California as a whole. We suggest that incorporating these features into a machine learning workflow could help resolve focal depths in regions that are poorly instrumented and lack large databases of well-located events. Some of the waveform features we evaluate in this study have previously been used as source discriminants, and our results imply that their effectiveness in discrimination is partially because explosions generally occur at shallower depths than earthquakes.

58 GEOSCIENCES↗