Search NASA⌕ Search

SEARCH · Search NASA

Results for “Estimated On Time”

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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Disproportionate impacts of COVID-19 in a large US city

COVID-19 has disproportionately impacted individuals depending on where they live and work, and based on their race, ethnicity, and socioeconomic status. Studies have documented catastrophic disparities at critical points throughout the pandemic, but have not yet systematically tracked their severity through time. Using anonymized hospitalization data from March 11, 2020 to June 1, 2021 and fine-grain infection hospitalization rates, we estimate the time-varying burden of COVID-19 by age group and ZIP code in Austin, Texas. During this 15-month period, we estimate an overall 23.7% (95% CrI: 22.5–24.8%) infection rate and 29.4% (95% CrI: 28.0–31.0%) case reporting rate. Individuals over 65 were less likely to be infected than younger age groups (11.2% [95% CrI: 10.3–12.0%] vs 25.1% [95% CrI: 23.7–26.4%]), but more likely to be hospitalized (1,965 per 100,000 vs 376 per 100,000) and have their infections reported (53% [95% CrI: 49–57%] vs 28% [95% CrI: 27–30%]). We used a mixed effect poisson regression model to estimate disparities in infection and reporting rates as a function of social vulnerability. We compared ZIP codes ranking in the 75th percentile of vulnerability to those in the 25th percentile, and found that the more vulnerable communities had 2.5 (95% CrI: 2.0–3.0) times the infection rate and only 70% (95% CrI: 60%-82%) the reporting rate compared to the less vulnerable communities. Inequality persisted but declined significantly over the 15-month study period. Our results suggest that further public health efforts are needed to mitigate local COVID-19 disparities and that the CDC’s social vulnerability index may serve as a reliable predictor of risk on a local scale when surveillance data are limited.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2026

Preliminary proof-of-concept of real-time divertor heat flux control from infrared cameras with nitrogen injection in the DIII-D tokamak

In future tokamak reactors like ITER and the Fusion Pilot Plant (FPP), real-time feedback control of heat flux to the plasma-facing components (PFC) will be critical for steady-state operation. This work presents the first experimental demonstration of real-time divertor heat flux estimation with infrared thermography and feedback control with impurity seeding on the DIII-D tokamak. The flexible infrastructure of the Plasma Control System (PCS) on DIII-D makes this new capability possible. The PCS software runs on a gateway computer system, and five real-time compute nodes. An array of low latency streaming digitizers from D-TACQ Solutions connects to these real-time computers to collect and process data, and send commands to actuators during plasma discharges. This system handles the signal IO from the tokamak and allows the PCS to utilize the diagnostic data necessary to perform control in real-time. Feedback control on heat flux was accomplished by feeding infrared camera data from the “Infrared TV” (IRTV) camera to a custom-developed User Datagram Protocol (UDP) server. This server transmits infrared data to a newly developed PCS algorithm that estimates the heat flux to PFC. Here, a proportional integral derivative (PID) controller minimizes the error between a heat flux reference and the real-time estimate by injecting nitrogen gas into the divertor.

Algorithm development↗
DOE OSTI2020

Economic Valuation of Coccidioidomycosis (Valley Fever) Projections in the United States in Response to Climate Change

Coccidioidomycosis, or valley fever, is an infectious fungal disease currently endemic to the southwestern United States. Symptoms of valley fever range in severity from flu-like illness to severe morbidity and mortality. Warming temperatures and changes in precipitation patterns may cause the area of endemicity to expand northward throughout the western United States, putting more people at risk for contracting valley fever. This may increase the health and economic burdens from this disease. We developed an approach to describe the relationship between climate conditions and valley fever incidence using historical data and generated projections of future incidence in response to both climate change and population trends using the Climate Change Impacts and Risk Analysis (CIRA) framework developed by the U.S. Environmental Protection Agency. We also developed a method to estimate economic impacts of valley fever that is based on case counts. For our 2000–15 baseline time period, we estimated annual medical costs, lost income, and economic welfare losses for valley fever in the United States were 400,000 USD per case, and the annual average total cost was 3.9 billion USD per year. For a high greenhouse gas emission scenario and accounting for population growth, we found that total annual costs for valley fever may increase up to 164% by year 2050 and up to 380% by 2090. By the end of the twenty-first century, valley fever may cost 620,000 USD per case and the annual average total cost may reach $18.5 billion per year. This work contributes to the broader effort to monetize climate change–attributable damages in the United States.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Elastic Changepoint Detection for Globally-indexed Functional Time Series Data with Climate Applications

Changepoint detection is a vital tool in the application of climate data analysis. Numerous types of climate observation data are most properly represented by functional time series, implying a need for accurate changepoint detection methods applicable to functional time series data. Such data taken at a global scale often contain both spatial heterogeneity and dependence as well as phase (time) misalignment. In this report, we present methods which can detect spatially-dependent changepoints while allowing different estimates of change time and change strength depending on location. Additionally, we provide extensions to this spatially-predicted model which controls for phase variability among observations. Our methods provide the ability to detect a single change, or control for epidemic changes (where a “return-to-normal” change is more likely to be detected than the initial change). We showcase results analyzing the June 1991 eruption of Mt. Pinatubo, where our methods demonstrate the ability to accurately detect both single and epidemic changepoints even in the presence of strong seasonal variability. We find that our spatially-predicted model improves the detection of relevant changepoints versus methods which do not take spatial information into account, and we find that controlling for phase variability helps to control the false discovery rate during the detection process.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Estimating the Performance Loss Rate of Photovoltaic Systems Using Time Series Change Point Analysis

The accurate quantification of the performance loss rate of photovoltaic systems is critical for project economics. Following the current research activities in the photovoltaic performance and reliability field, this work presents a comparative assessment between common change point methods for performance loss rate estimation of fielded photovoltaic installations. An extensive testing campaign was thus performed to evaluate time series analysis approaches for performance loss rate evaluation of photovoltaic systems. Historical electrical data from eleven photovoltaic systems installed in Nicosia, Cyprus, and the locations’ meteorological measurements over a period of 8 years were used for this investigation. The application of change point detection algorithms on the constructed monthly photovoltaic performance ratio series revealed that the obtained trend might not always be linear. Specifically, thin film photovoltaic systems showed nonlinear behavior, while nonlinearities were also detected for some crystalline silicon photovoltaic systems. When applying several change point techniques, different numbers and locations of changes were detected, resulting in different performance loss rate values (varying by up to 0.85%/year even for the same number of change points). The results highlighted the importance of the application of nonlinear techniques and the need to extract a robust nonlinear model for detecting significant changes in time series data and estimating accurately the performance loss rate of photovoltaic installations.

photovoltaics↗
DOE OSTI2022

Sulfuric acid in the Amazon basin: measurements and evaluation of existing sulfuric acid proxies

Abstract. Sulfuric acid is a key contributor to new particle formation, though measurements of its gaseous concentrations are difficult to make. Several parameterizations to estimate sulfuric acid exist, all of which were constructed using measurements from the Northern Hemisphere. In this work, we report the first measurements of sulfuric acid from the Amazon basin. These measurements are consistent with concentrations measured in Hyytiälä, Finland, though, unlike Hyytiälä, there is no clear correlation of sulfuric acid with global radiation. There was a minimal difference in sulfuric acid observed between the wet and dry seasons in the Amazon basin. We also test the efficacy of existing proxies to estimate sulfuric acid in this region. Our results suggest that nighttime sulfuric acid production is due to both a stabilized Criegee intermediate pathway and oxidation of SO2 by OH, the latter of which is not currently accounted for in existing proxies. These results also illustrate the drawbacks of the common substitution of radiation for OH concentrations. None of the tested proxies effectively estimate sulfuric acid measurements at night. For estimates at all times of day, a recently published proxy based on data from the boreal forest should be used. If only daytime estimates are needed, several recent proxies that do not include the Criegee pathway are sufficient. More investigation of nighttime sulfuric acid production pathways is necessary to close the gap between measurements and estimates with existing proxies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Canopy tree mortality and crown exposure data from the Amacayacu Forest Dynamics Plot, Northwestern Amazon

Data on the mortality of 984 canopy trees, their crown exposure to light (relative to total crown area), growth deviations (relative to conspecifics), tree size, and species’ wood density collected between 2013 and 2019 in 18 ha of the Amacayacu Forest Dynamics Plot, Northwestern Amazon. This dataset contains a single CSV data file. Variable definitions: 1. Species: [character] species identification 2. Family: [character] family of the species 3. tag: [character] unique consecutive for the tree 4. status: [character] status of the tree in the third census (ALIVE or DEAD) 5. wsg: [numeric]: species’ wood density (g cm-3) 6. growth_r1: [numeric]: annual growth rate of the tree between the first and second census (cm y-1) 7. growth_r2: [numeric]: annual growth rate of the tree between the second and third census (cm y-1) 8. gt1: [numeric] modulus transformed growth rate with a lambda of 0.4 for the annual growth rate of the tree between the first and second census (cm y-1) 9. gt2: [numeric]: modulus transformed growth rate with a lambda of 0.4 for the annual growth rate of the tree between the second and third census (cm y-1) 10. rgr1: [numeric] relative growth rate between the first and second census (cm y-1) 11. rgr2: [numeric] relative growth rate between the second and third census (cm y-1) 12. sa_gt: [numeric] species-adjusted modulus transformed growth rate 13. sa_rgr: [numeric] species-adjusted relative growth rate 14. gr_n: [numeric] number of individuals of the species used to calculate the mean and species modulus transformed growth rate 15. dbh_flight: [numeric] diameter at breast height (1.3 m) estimated at the time of the drone flight (cm) 16. eca: [numeric] exposed crown area calculated as the area of the crown polygon delineated in the orthomosaic (m2) 17. total_ca: [numeric] total crown area estimated from a crown area model (m2) 18. rcel: [numeric] relative crown exposure to light (m2) 19. time_flight_census3: [numeric] time in years from the drone flight date to the third census for that tree (yr)

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Simulation of Impedance Changes with Aging in Lithium Titanate-based Cells Using Physics-Based Dimensionless Modeling

Quantifying aging effects in lithium-ion cells with chemistries that have a flat open circuit potential is challenging. We implement a physics-based electrochemical model to track changes in the electrochemical impedance response of lithium titanate-based cells. Frequency domain equations of a pseudo two-dimensional model are made dimensionless, and the corresponding non-dimensional parameters are estimated using a Levenberg-Marquardt routine. The model weighs the relative contributions of changes in diffusion, ionic conduction within the electrolyte phase against solid phase electronic conduction towards cell aging. Solid-phase diffusion, charge transfer resistance and double layer capacitance at the solid-liquid interface are accounted for in the particle impedance. The estimation routine tracks dimensionless parameters using accelerated cycling data from full cells over 1000 cycles. The model can be deployed within a short time for state estimation using physics-based models without requiring prior knowledge of the battery chemistry, format, or capacity.

25 ENERGY STORAGE↗
DOE OSTI2021

Game-Theoretic Approach for Grace-Period Policy in Supercomputers

Job scheduling at supercomputing facilities is important for achieving high utilization of these valuable resources while ensuring effective execution of jobs submitted by users. The jobs are scheduled according to their specified resource demands such as expected job completion times, and the available resources based on allocations. Jobs that overrun their allocated times are terminated, for example, after a grace-period. It is non-trivial and often very complex for users to accurately estimate the completion times of their jobs, and consequently they face a dilemma: underestimate the job time to have a higher priority and risk job termination due to overrun, or overestimate it to ensure its completion and risk its delayed execution. In this paper, we investigate whether providing grace-period can benefit facility performance by developing a game- theoretic model between a facility provider and multiple users for a simplified scheduling scenario based on job execution times. We present closed-form expressions for the provider’s and user’s best-response strategies to maximize their respective utility functions. We describe conditions under which offering a grace-period is advantageous to both facility provider and users by deriving the Nash equilibrium of the game.

He, Fei↗
DOE OSTI2022

Probabilistic estimation of depth-resolved profiles of soil thermal diffusivity from temperature time series

Abstract. Improving the quantification of soil thermal and physical properties is key to achieving a better understanding and prediction of soil hydro-biogeochemical processes and their responses to changes in atmospheric forcing. Obtaining such information at numerous locations and/or over time with conventional soil sampling is challenging. The increasing availability of low-cost, vertically resolved temperature sensor arrays offers promise for improving the estimation of soil thermal properties from temperature time series, and the possible indirect estimation of physical properties. Still, the reliability and limitations of such an approach need to be assessed. In the present study, we develop a parameter estimation approach based on a combination of thermal modeling, sliding time windows, Bayesian inference, and Markov chain Monte Carlo simulation to estimate thermal diffusivity and its uncertainty over time, at numerous locations and at an unprecedented vertical spatial resolution (i.e., down to 5 to 10 cm vertical resolution) from soil temperature time series. We provide the necessary framework to assess under which environmental conditions (soil temperature gradient, fluctuations, and trend), temperature sensor characteristics (bias and level of noise), and deployment geometries (sensor number and position) soil thermal diffusivity can be reliably inferred. We validate the method with synthetic experiments and field studies. The synthetic experiments show that in the presence of median diurnal fluctuations ≥ 1.5 ∘C at 5 cm below the ground surface, temperature gradients > 2 ∘C m−1, and a sliding time window of at least 4 d the proposed method provides reliable depth-resolved thermal diffusivity estimates with percentage errors ≤ 10 % and posterior relative standard deviations ≤ 5 % up to 1 m depth. Reliable thermal diffusivity under such environmental conditions also requires temperature sensors to be spaced precisely (with accuracy to a few millimeters), with a level of noise ≤ 0.02 ∘C, and with a bias defined by a standard deviation ≤ 0.01 ∘C. Finally, the application of the developed approach to field data indicates significant repeatability in results and similarity with independent measurements, as well as promise in using a sliding time window to estimate temporal changes in soil thermal diffusivity, as needed to potentially capture changes in bulk density or water content.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Deep learning to estimate permeability using geophysical data

Time-lapse electrical resistivity tomography (ERT) is a popular geophysical method to estimate three-dimensional (3D) permeability fields from electrical potential difference measurements. Traditional inversion and data assimilation methods are used to ingest this ERT data into hydrogeophysical models to estimate permeability. Due to ill-posedness and the curse of dimensionality, existing inversion strategies provide poor estimates and low resolution of the 3D permeability field. Recent advances in deep learning provide us with powerful algorithms to overcome this challenge. This paper presents a deep learning (DL) framework to estimate the 3D subsurface permeability from time-lapse ERT data. To test the feasibility of the proposed framework, we train DL-enabled inverse models on simulation data. Each measurement in both synthetic and field data is standardized by removing the mean and scaling the time-series to unit variance. This pre-processing step is necessary to bring simulation data closer to field observations. Subsurface process models based on hydrogeophysics are used to generate this synthetic data. Training performed on limited simulation data resulted in the DL model over-fitting. An advanced data augmentation based on mixup is implemented to generate additional training samples to overcome this issue. This mixup technique creates weakly labeled (low-fidelity) samples from strongly labeled (high-fidelity) data. The weakly labeled training data is then used to develop DL-enabled inverse models and reduce over-fitting. As both time-lapse ERT (1133048 features/realization) and 3D permeability (585453 features/realization) data samples are from a high-dimensional space, principal component analysis (PCA) is employed to reduce dimensionality. Encoded ERT and encoded permeability are generated using the trained PCA estimators. A deep neural network is then trained to map the encoded ERT to encoded permeability. This mixup training and unsupervised learning allowed us to build a fast and reasonably accurate DL-based inverse model under limited simulation data. Results show that proposed weak supervised learning can capture salient spatial features in the 3D permeability field. Quantitatively, the average mean squared error (in terms of the natural log) on the strongly labeled training, validation, and test datasets is less than 0.5. The R 2 -score (global metric) is greater than 0.75, and the percent error in each cell (local metric) is less than 10%. Finally, an added benefit in terms of computational cost is that the proposed DL-based inverse model is at least O(10 4 ) times faster than running a forward model once it is trained. Data generation, DL model training, and hyperparameter tuning to identify optimal neural network architectures utilized high-performance computing resources while the DL inference is performed on a standard laptop. Approximately, O(10 5 ) processor hours are used for generating data and DL tuning and training. We acknowledge that the data generation and DL model development are expensive. But once a DL model is trained, it can be re-used for inversion rapidly for the given system, with set physics and domain. Note that traditional inversion may require multiple forward model simulations (e.g., in the order of 10 to 1000), which are very expensive. This computational savings ≈ O(10 5 ) – O(10 7 )) makes the proposed DL-based inverse model attractive for subsurface imaging and real-time ERT monitoring applications due to fast and yet reasonably accurate estimations of permeability field.

58 GEOSCIENCES↗
DOE OSTI2022

Load Disaggregation (Modeling Individual Appliance Power Generation & Consumption in Real-Time)

We developed a machine learning-based load disaggregation method to estimate the real-time output of individual appliances from the whole-house measurements. We first learn the important features associated with each type of appliances using the ground truth consumption data of individual appliances potentially available for a small set of houses equipped with submeters. The learned features are then be used to estimate the power generation/consumption of the appliances from the whole-house consumption. This developed load disaggregation software includes two steps. The first step is to identify the on/off status of different appliances using a classification method, and the second step is to estimate the appliance output using a regression method.

Hao, Jun↗
DOE OSTI2022

Recursive Gaussian Process over graphs for Integrating Multi-timescale Measurements in Low-Observable Distribution Systems

The transition to a smarter grid is empowered by enhanced sensor deployments and smart metering infrastructure in the distribution system. Measurements from these sensors and meters can be used for many applications, including distribution system state estimation (DSSE). However, these measurements are typically sampled at different rates and could be intermittent due to losses during the aggregation process. These multi timescale measurements should be reconciled in real-time to perform accurate grid monitoring. This paper tackles this problem by formulating a recursive multi-task Gaussian process (RGP-G) approach that sequentially aggregates sensor measurements. Specifically, we formulate a recursive multi-task GP with and without network connectivity information to reconcile the multi time-scale measurements in distribution systems. Here, the proposed framework is capable of aggregating the multi-time scale measurements batch-wise or in real-time. Following the aggregation of the multi time-scale measurements, the spatial states of the consistent time-series are estimated using matrix completion based DSSE approach. Simulation results on IEEE 37 and IEEE 123 bus test systems illustrate the efficiency of the proposed methods from the standpoint of both multi time-scale data aggregation and DSSE.

42 ENGINEERING↗
DOE OSTI2022

Chemical abundances in Seyfert galaxies – Ⅸ. Helium abundance estimates

For the first time, the helium abundance relative to hydrogen (He/H), which relied on direct measurements of the electron temperature, has been derived in the narrow line regions (NLRs) from a local sample of Seyfert 2 nuclei. In view of this, optical emission line intensities [3000 < λ(Å) < 7000] of 65 local Seyfert 2 nuclei (z < 0.2), taken from Sloan Digital Sky Survey Data Release 15 and additional compilation from the literature, were considered. Here we used photoionization model grid to derive an Ionization Correction Factor (ICF) for the neutral helium. The application of this ICF indicates that the NLRs of Seyfert 2 present a neutral helium fraction of ~50 percent in relation to the total helium abundance. We find that Seyfert 2 nuclei present helium abundance ranging from 0.60 to 2.50 times the solar value, while ~85 percent of the sample present oversolar abundance values. The derived (He/H)–(O/H) abundance relation from the Seyfert 2 is stepper than that of star-forming regions (SFs) and this difference could be due to excess of helium injected into the interstellar medium by the winds of Wolf–Rayet stars. From a regression to zero metallicity, by using Seyfert 2 estimates combined with SFs estimates, we obtained a primordial helium mass fraction Y p = 0.2441 ± 0.0037, a value in good agreement with the one inferred from the temperature fluctuations of the cosmic microwave background by the Planck Collaboration, i.e. $Y_{p}^{Planck}$ =0.2471 ± 0.0003.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Recent Development of Frequency Estimation Methods for Future Smart Grid

The frequency estimated by the Phasor Measurement Unit (PMU) is a critical index of power system status and supports many smart grid applications. The future smart grid features high penetration of renewables and more fast-moving power electronics inverters but raises challenges to the reliable frequency estimation. This article presents three methods to address these challenges. First, an enhanced zero-crossing algorithm was developed to track the fast-changing frequency in system dynamics. Second, we propose a technology that can tolerate the system transient and suppress the outliers. Third, an algorithm was developed to export high time-resolution frequency estimations with minimum computational effort. All of the proposed methods are realized in hardware and compared with classical frequency estimation methods. The testing results indicate that the proposed methods have excellent performance. They can be used in future PMUs and provide reliable and high time resolution data for smart grid applications.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

An Online Approach to Solve the Dynamic Vehicle Routing Problem with Stochastic Trip Requests for Paratransit Services

Many transit agencies operating paratransit and microtransit services have to respond to trip requests that arrive in real-time, which entails solving hard combinatorial and sequential decision-making problems under uncertainty. To avoid decisions that lead to significant inefficiency in the long term, vehicles should be allocated to requests by optimizing a non-myopic utility function or by batching requests together and optimizing a myopic utility function. While the former approach is typically offline, the latter can be performed online. We point out two major issues with such approaches when applied to paratransit services in practice. First, it is difficult to batch paratransit requests together as they are temporally sparse. Second, the environment in which transit agencies operate changes dynamically (e.g., traffic conditions can change over time), causing the estimates that are learned offline to become stale. To address these challenges, we propose a fully online approach to solve the dynamic vehicle routing problem (DVRP) with time windows and stochastic trip requests that is robust to changing environmental dynamics by construction. We focus on scenarios where requests are relatively sparse—our problem is motivated by applications to paratransit services. We formulate DVRP as a Markov decision process and use Monte Carlo tree search to evaluate actions for any given state. Accounting for stochastic requests while optimizing a non-myopic utility function is computationally challenging; indeed, the action space for such a problem is intractably large in practice. To tackle the large action space, we leverage the structure of the problem to design heuristics that can sample promising actions for the tree search. Our experiments using real-world data from our partner agency show that the proposed approach outperforms existing state-of-the-art approaches both in terms of performance and robustness.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

Observational Constraints on the Cloud Feedback Pattern Effect

Abstract Model evidence for the “pattern effect” assumes that global climate models (GCMs) faithfully simulate how clouds respond to varying sea surface temperature (SST) patterns and associated meteorological perturbations. We exploit time-invariant satellite-based estimates of the sensitivity of marine low clouds to meteorological perturbations to estimate how these clouds responded to time-varying SST patterns and meteorology between 1870 and 2014. GCMs and reanalyses provide estimates of the historical meteorological changes. Observations suggest that increasing estimated inversion strength (EIS) between 1980 and 2014 produced a negative low cloud feedback, opposite to the positive feedback expected from increasing CO 2 . This indicates that the processes responsible for marine cloud changes from 1980 to the near present are distinct from those associated with an increase in CO 2 . We also observationally constrain the difference between the historical near-global marine low cloud feedback, , and that arising from increasing CO 2 , . We find that this cloud feedback pattern effect depends strongly on time period and reanalysis dataset, and that varying changes in EIS and SST with warming explain much of its variability. Between 1980 and 2014, we estimate that (90% confidence) assuming meteorological changes from the Multiple Reanalysis Ensemble, implying a total pattern effect (that arising from all climate feedbacks) of 1.86 ± 0.45 W m −2 K −1 . This observational evidence corroborates previous quantitative estimates of the pattern effect, which heretofore relied largely upon GCM-based cloud changes. However, disparate historical meteorological changes across individual reanalyses contribute to considerable uncertainty in its magnitude.

58 GEOSCIENCES↗
DOE OSTI2022

Earthquake Phase Association Using a Bayesian Gaussian Mixture Model

Earthquake phase association algorithms aggregate picked seismic phases from a network of seismometers into individual seismic events and play an important role in earthquake monitoring and research. Dense seismic networks and improved phase picking methods produce massive seismic phase datasets, particularly for earthquake swarms and aftershocks occurring closely in time and space, making phase association a challenging problem. Here, we present a new association method, the Gaussian Mixture Model Association (GaMMA), that combines the Gaussian mixture model with earthquake location, origin time, and magnitude estimation. We treat earthquake phase association as an unsupervised clustering problem in a probabilistic framework, where each earthquake corresponds to a cluster of P and S phases with a hyperbolic moveout of arrival times and a decay of amplitude with distance. We use the multivariate Gaussian distribution to model the collection of phase picks of an event; and the mean of the multivariate Gaussian distribution is given by the predicted arrival time and amplitude from the causative event. We carry out the pick assignment to each earthquake and determine earthquake source parameters (i.e., earthquake location, origin time, and magnitude) under the maximum likelihood criterion using the Expectation-Maximization algorithm. The GaMMA method does not require typical association steps of other algorithms, such as grid-search or supervised training. The results for both synthetic tests and for the 2019 Ridgecrest earthquake sequence show that GaMMA effectively associates phases from a temporally and spatially dense earthquake sequence while producing useful estimates of earthquake location and magnitude.

58 GEOSCIENCES↗