Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “regression models”

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 379 records · Page 21

Plasma phosphorylated tau217 strongly associates with memory deficits in the Alzheimer’s disease spectrum

Abstract Plasma phosphorylated tau (p-tau) biomarkers open unprecedented opportunities for identifying carriers of Alzheimer’s disease pathophysiology in early disease stages using minimally invasive techniques. Plasma p-tau biomarkers are believed to reflect tau phosphorylation and secretion. However, it remains unclear to what extent the magnitude of plasma p-tau abnormalities reflects neuronal network disturbance in the form of cognitive impairment. To address this question, we included 103 cognitively unimpaired elderly and 40 cognitively impaired, amyloid-β-positive individuals from the TRIAD cohort, in addition to 336 cognitively unimpaired and 216 cognitively impaired, amyloid-β-positive older adults from the BioFINDER-2 cohort. Participants had tau PET scans, amyloid PET scans or amyloid CSF, p-tau217, p-tau181 and p-tau231 blood measures, structural T1-MRI and cognitive assessments. In this cross-sectional study, we used regression models and correlation analyses to assess the relationship between plasma biomarkers and cognitive scores. Furthermore, we applied receiver operating characteristic curves to assess cognitive impairment across plasma biomarkers. Finally, we categorized participants into amyloid (A), p-tau (T1) and tau PET (T2) positive (+) or negative (−) profiles and ran non-parametric comparisons to assess differences across cognitive domains. We found that plasma p-tau217 was more associated with cognitive performance than p-tau181 and p-tau231 and that this relationship was particularly strong for memory scores (TRIAD: βp-tau217 = −0.53, βp-tau181 = −0.35 and βp-tau231 = −0.24; BioFINDER-2: βp-tau217 = −0.52, βp-tau181 = −0.24 and βp-tau231 = −0.29). Associations in amyloid-β-positive participants resembled these results, but other cognitive scores also showed strong associations in cognitively impaired individuals. Moreover, plasma p-tau217 outperformed plasma p-tau181 and plasma p-tau231 in identifying memory impairment (area under the curve values for TRIAD: p-tau217 = 0.86, p-tau181 = 0.77 and p-tau231 = 0.75; and for BioFINDER-2: p-tau217 = 0.86, p-tau181 = 0.76 and p-tau231 = 0.81) and in identifying executive function impairment only in the BioFINDER-2 cohort (p-tau217 = 0.82, p-tau181 = 0.76 and p-tau231 = 0.76). Lastly, we showed that subtle memory deficits were present in A+T1+T2− participants for plasma p-tau217 (P = 0.007) and plasma p-tau181 (P = 0.01) in the TRIAD cohort and for all biomarkers across cognitive domains in A+T1+T2− and A+T1+T2− individuals (P < 0.001 in all) in the BioFINDER-2 cohort. The A+T1+T2− individuals showed cognitive deficits in both cohorts (P < 0.001 in all). Together, our results suggest that plasma p-tau217 stands out as a biomarker capable of identifying memory deficits attributable to Alzheimer’s disease and that memory impairment certainly occurs in amyloid-β- and plasma p-tau-positive individuals who have no significant amounts of tau in the neocortex.

Neurosciences & Neurology↗

Symptoms and Systemic Drug Reactions in Persons Receiving Weekly Rifapentine Plus Isoniazid (3HP) Treatment for Latent Tuberculosis Infection

Abstract Background Three months of weekly rifapentine plus isoniazid (3HP) therapy for latent tuberculosis infection (LTBI) is recommended worldwide. The development of symptoms and systemic drug reactions (SDRs) on 3HP have not been fully characterized. We aimed to determine the patterns of symptom development and identify SDRs and associated factors in patients taking 3HP. Methods We analyzed symptoms data in participants receiving 3HP in the Tuberculosis Trials Consortium's iAdhere study (Study 33). We examined the patterns of symptom reporting across participants from baseline and 4 monthly visits. Bivariate analyses and multivariable regression models were used to identify factors associated with SDRs. Risk ratios (RRs) and 95% confidence intervals (CIs) were calculated. Results Among 1002 participants receiving 3HP, 768 (77%) reported at least 1 symptom; 97% of these symptoms were grade 1 (79%) or grade 2 (18%). Most symptoms developed in the first month and resolved. A total of 111 (11%) participants had symptoms that met criteria for SDRs; however, 53 (48%) of these participants completed therapy. Factors associated with SDRs and discontinuation included female sex (RR: 2.05; 95% CI: 1.19–3.54), age ≥45 years (RR: 1.99; 95% CI: 1.19–3.31), and use of concomitant medications (RR: 2.26; 95% CI: 1.15–4.42). Conclusions Although most patients receiving 3HP reported symptoms, most were mild, occurred early, and resolved without stopping treatment. Among patients experiencing SDRs, nearly half were able to complete therapy. Patient and provider education should focus on differentiating severe reactions where 3HP should be stopped from minor symptoms that will resolve. Clinical Trials Registration. NCT01582711.

Immunology↗

Rates of Seasonal Fuel Loading Do Not Differ by Sex or Overwintering Strategy in Three Species of Bats

For temperate-zone mammals, seasonal changes in weather and food availability often govern energy allocation. In addition, energy allocation strategies usually differ between males and females. Bats are an interesting group in which to evaluate energetic trade-offs as they are highly mobile and lead energetically demanding lives in habitats across a variety of seasonally variable climates. We evaluated year-round changes in body mass and fuel load for three species of bats in northern California: Silver-haired Bat (Lasionycteris noctivagans), California Myotis (Myotis californicus), and Yuma Myotis (Myotis yumanensis). Silver-haired bats are considered migratory species with females likely migrating farther than males. The two species of myotis are considered residents. Body mass of all species peaked in late autumn and was at a minimum during spring. We calculated a fuel load index to normalize size difference between species and sexes. We used sex- and season-specific multiple linear regression models to evaluate rates of change in seasonal fuel loading. Rates of change in fuel load did not differ among species or sexes except for male silver-haired bats that increased fuel loads rapidly during summer. Interspecific comparisons provided valuable insights into the energy allocation and overwintering strategies of these species and are an important initial step toward understanding their ecology over the full annual cycle.

body condition↗

Comparison of sleep parameters from wrist-worn ActiGraph and Actiwatch devices

Abstract Sleep and physical activity, two important health behaviors, are often studied independently using different accelerometer types and body locations. Understanding whether accelerometers designed for monitoring each behavior can provide similar sleep parameter estimates may help determine whether one device can be used to measure both behaviors. Three hundred and thirty one adults (70.7 ± 13.7 years) from the Baltimore Longitudinal Study of Aging wore the ActiGraph GT9X Link and the Actiwatch 2 simultaneously on the non-dominant wrist for 7.0 ± 1.6 nights. Total sleep time (TST), wake after sleep onset (WASO), sleep efficiency, number of wake bouts, mean wake bout length, and sleep fragmentation index (SFI) were extracted from ActiGraph using the Cole–Kripke algorithm and from Actiwatch using the software default algorithm. These parameters were compared using paired t-tests, Bland–Altman plots, and Deming regression models. Stratified analyses were performed by age, sex, and body mass index (BMI). Compared to the Actiwatch, the ActiGraph estimated comparable TST and sleep efficiency, but fewer wake bouts, longer WASO, longer wake bout length, and higher SFI (all p < .001). Both devices estimated similar 1-min and 1% differences between participants for TST and SFI (β = 0.99, 95% CI: 0.95, 1.03, and 0.91, 1.13, respectively), but not for other parameters. These differences varied by age, sex, and/or BMI. The ActiGraph and the Actiwatch provide comparable absolute and relative estimates of TST, but not other parameters. The discrepancies could result from device differences in movement collection and/or sleep scoring algorithms. Further comparison and calibration is required before these devices can be used interchangeably.

Liu, Fangyu (ORCID:0000000315541478)↗

Association of human T-cell leukemia virus type 1 with prevalent rheumatoid arthritis among atomic bomb survivors: A cross-sectional study

Previous studies have suggested that human T-cell leukemia virus type 1 (HTLV-1) might act as a pathogen in rheumatoid arthritis (RA), but epidemiological evidence of an association is scarce. We measured anti-HTLV-1 antibodies among Nagasaki atomic bomb survivors to determine whether HTLV-1 is related to RA and whether radiation exposure is associated with HTLV-1 and RA prevalence. This is a cross-sectional study among atomic bomb survivors who participated in biennial health examinations from 2006 to 2010. Serum levels of anti-HTLV-1 antibodies were measured using a chemiluminescent enzyme immunoassay and confirmed by Western blotting. Association between HTLV-1 and RA was analyzed by a logistic regression model. Of 2091 participants (women 61.5%; median age, 73 years), 215 (10.3%) had anti-HTLV-1 antibodies. HTLV-1 prevalence was higher among women (13.1% vs 5.8%; P<.001). Twenty-two participants (1.1%) were diagnosed with RA. HTLV-1 prevalence among RA participants was significantly higher than that among non-RA participants (27.3% vs 10.1%; P=.020). After adjustment for age, sex, and hepatitis C virus infection, HTLV-1 was significantly associated with prevalent RA (odds ratio, 2.89; 95% confidence interval, 1.06, 7.03). There was no association between radiation dose and either the prevalence of HTLV-1 or RA. This study, among a well-defined group of atomic bomb survivors, suggests that HTLV-1 is associated with RA. Abbreviations: AHS = Adult Health Study, ATL/ATLL = adult T cell leukemia/lymphoma, bioDMARDs = biological disease modifying antirheumatic drugs, CCP = cyclic citrullinated peptide, CI = confidence interval, CLEIA = chemiluminescent enzyme immunoassay, CRP = C-reactive protein, csDMARDs = conventional synthetic disease modifying antirheumatic drugs, Gy = gray, HCV = hepatitis C virus, HTLV-1 = human T-cell leukemia virus type 1, IF = immunofluorescent, RA = rheumatoid arthritis, RERF = Radiation Effects Research Foundation, RF = rheumatoid factor, SS = Sjogren syndrome.

59 BASIC BIOLOGICAL SCIENCES↗

Learning from metastable symmetric-tilt grain boundaries using physics-based descriptors

Grain boundaries (GBs) govern critical properties of polycrystalline materials. Although significant advancements have been made in characterizing minimum energy and ordered GBs, real GBs are seldom found in such well-defined states. This diversity of atomic arrangements in metastable states makes it challenging to establish structure-property relationships with physical insights. Here, to address this challenge, we use data-driven methods to explore these relationships and examine the underlying physics. In this study, we utilize a large atomistic database (~5000) of minimum energy and metastable states of symmetric-tilt copper GBs, combined with physically motivated local atomic environment (LAE) descriptors [strain functional descriptors (SFDs)], to predict GB properties and gain physical insights. Our regression models exhibit robust predictive capabilities using only 19 descriptors, generalizing to atomic environments in nanocrystals. A significant highlight of our work is the integration of an unsupervised method with SFDs to elucidate LAEs at GBs and their role in determining properties. The model, trained on these minimum energy and metastable GBs using SFDs, predicts the properties of unseen nanocrystals with good accuracy. Our research underscores the role of a physics-based representation of LAEs and the efficacy of data-driven methods in establishing GB structure-property relationships.

36 MATERIALS SCIENCE↗

MICCO: An Enhanced Multi-GPU Scheduling Framework for Many-Body Correlation Functions

Calculation of many-body correlation functions is one of the critical kernels utilized in many scientific computing areas, especially in Lattice Quantum Chromodynamics (Lattice QCD). It is formalized as a sum of a large number of contraction terms each of which can be represented by a graph consisting of vertices describing quarks inside a hadron node and edges designating quark propagations at specific time intervals. Due to its computation- and memory-intensive nature, real-world physics systems (e.g., multi-meson or multi-baryon systems) explored by Lattice QCD prefer to leverage multi-GPUs. Different from general graph processing, many-body correlation function calculations show two specific features: a large number of computation-/data-intensive kernels and frequently repeated appearances of original and intermediate data. The former results in expensive memory operations such as tensor movements and evictions. The latter offers data reuse opportunities to mitigate the data-intensive nature of many-body correlation function calculations. However, existing graph-based multi-GPU schedulers cannot capture these data-centric features, thus resulting in a sub-optimal performance for many-body correlation function calculations. To address this issue, this paper presents a multi-GPU scheduling framework, MICCO, to accelerate contractions for correlation functions particularly by taking the data dimension (e.g., data reuse and data eviction) into account. This work first performs a comprehensive study on the interplay of data reuse and load balance, and designs two new concepts: local reuse pattern and reuse bound to study the opportunity of achieving the optimal trade-off between them. Based on this study, MICCO proposes a heuristic scheduling algorithm and a machine-learning-based regression model to generate the optimal setting of reuse bounds. Specifically, MICCO is integrated into a real-world Lattice QCD system, Redstar, for the first time running on multiple GPUs. The evaluation demonstrates MICCO outperforms other state-of-art works, achieving up to 2.25× speedup in synthesized datasets, and 1.49× speedup in real-world correlation functions.

Wang, Qihan↗

A Data-Driven Method for Estimating Behind-the-Meter Photovoltaic Generation in Hawaii

Due to the increasing penetration of distributed behind-the-meter photovoltaic (PV) systems and the installed utility revenue metering limited to monitoring only the net power import/export of the household, it is increasingly challenging for utilities to effectively plan and operate the grid. This paper proposes a methodology that estimates behind-the-meter PV generation using a selected subset of monitored PV systems. It is a data-driven approach, and the PV output is estimated utilizing a statistic regression model. A Minimum Redundancy Maximum Relevance (MRMR) algorithm is applied to preselect the optimal subset of the monitored PV systems. The performance of this approach is compared with a spatial interpolation method and a model-based approach. The proposed method is validated using high-resolution meter data recorded from 18 residential rooftop PV systems located on the island of Maui, Hawaii.

Data-driven modeling↗

Signatures of cysteine oxidation on muscle structural and contractile proteins are associated with physical performance and muscle function in older adults: Study of Muscle, Mobility and Aging (SOMMA)

Oxidative stress is considered a contributor to declining muscle function and mobility during aging; however, the underlying molecular mechanisms remain poorly described. We hypothesized that greater levels of cysteine (Cys) oxidation on muscle proteins are associated with decreased measures of mobility. Herein, we applied a novel redox proteomics approach to measure reversible protein Cys oxidation in vastus lateralis muscle biopsies collected from 56 subjects in the Study of Muscle, Mobility and Aging (SOMMA), a community-based cohort study of individuals aged 70 years and older. We tested whether levels of Cys oxidation on key muscle proteins involved in muscle structure and contraction were associated with muscle function (leg power and strength), walking speed, and fitness (VO 2 peak on cardiopulmonary exercise testing) using linear regression models adjusted for age, sex, and body weight. Higher oxidation levels of select nebulin Cys sites were associated with lower VO 2 peak, while greater oxidation of myomesin-1, myomesin-2, and nebulin Cys sites was associated with slower walking speed. Higher oxidation of Cys sites in key proteins such as myomesin-2, alpha-actinin-2, and skeletal muscle alpha-actin were associated with lower leg power and strength. We also observed an unexpected correlation (R = 0.48) between a higher oxidation level of eight Cys sites in alpha-actinin-3 and stronger leg power. Despite this observation, the results generally support the hypothesis that Cys oxidation of muscle proteins impairs muscle power and strength, walking speed, and cardiopulmonary fitness with aging.

59 BASIC BIOLOGICAL SCIENCES↗

Crowd-based spatial risk assessment of urban flooding: Results from a municipal flood hotline in Detroit, MI

Climate change is increasing the frequency and intensity of extreme precipitation events, raising the risk of urban flood disasters. This study uses a crowd-sourced municipal call database to characterize the spatial distribution of flood risk in Detroit, MI. Call data including dates and addresses were obtained from the City of Detroit Department of Public Works for 2021. Calls were mapped and aggregated to census tract counts and merged with neighborhood-level data. Associations of predictors with flood calls were tested using spatial regression models. Flooding calls were located throughout the city but were concentrated in specific areas. Multivariate models of census tract level call counts indicated that increased poverty and Black, immigrant, and older residents were positively associated with flood calls, while increased elevation was associated with protective effects. Longer distances from waste water interceptors were associated with higher risk for calls. Crowd-sourced flood hotline call data can be used for effective spatial flood risk assessment. Though flooding occurs throughout the city of Detroit, infrastructural, neighborhood, and household factors influence flooding extent. Limitations included the self-reported nature of calls. Future modeling efforts might include input from local stakeholders to improve spatial risk assessment.

54 ENVIRONMENTAL SCIENCES↗

Source:sink imbalance detected with leaf and canopy-level spectroscopy in a field-grown crop

The finely-tuned balance between sources and sinks determines plant resource partitioning and regulates growth and development. Understanding and measuring metabolic indicators of source or sink limitation forms a vital part of global efforts to increase crop yield for future food security. Here, we measured metabolic profiles of Cucurbita pepo (zucchini) grown in the field under carbon sink limitation and control conditions. We demonstrate that these profiles can be measured non-destructively using hyperspectral reflectance at both leaf and canopy scales. Total non-structural carbohydrates (TNC) increased 82% in sink-limited plants; leaf mass per unit area (LMA) increased 38%; free amino acids increased 22%. Partial least-squares regression models link these measured functional traits with reflectance data, enabling high-throughput estimation of traits comprising the sink limitation response. Leaf- and canopy-scale models for TNC had R 2 values of 0.93 and 0.64 and %RMSE of 13% and 38% respectively. For LMA, R 2 values were 0.91 and 0.60 and %RMSE 7% and 14%; for free amino acids, R 2 was 0.53 and 0.21 with %RMSE 20% and 26%. Remote sensing can enable accurate, rapid detection of sink limitation in the field at the leaf and canopy scale, greatly expanding our ability to understand and measure metabolic responses to stress.

54 ENVIRONMENTAL SCIENCES↗

A Bayesian nonparametric analysis for zero-inflated multivariate count data with application to microbiome study

High-throughput sequencing technology has enabled researchers to profile microbial communities from a variety of environments, but analysis of multivariate taxon count data remains challenging. Here, we develop a Bayesian nonparametric (BNP) regression model with zero inflation to analyse multivariate count data from microbiome studies. A BNP approach flexibly models microbial associations with covariates, such as environmental factors and clinical characteristics. The model produces estimates for probability distributions which relate microbial diversity and differential abundance to covariates, and facilitates community comparisons beyond those provided by simple statistical tests. We compare the model to simpler models and popular alternatives in simulation studies, showing, in addition to these additional community-level insights, it yields superior parameter estimates and model fit in various settings. The model's utility is demonstrated by applying it to a chronic wound microbiome data set and a Human Microbiome Project data set, where it is used to compare microbial communities present in different environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Numerical Investigation of a Central Fuel Property Hypothesis Under Boosted Spark-Ignition Conditions

In the present work, a central fuel property hypothesis (CFPH), which states that fuel properties are sufficient to provide an indication of a fuel’s performance irrespective of its chemical composition, was numerically investigated. In particular, the objective of the study was to determine whether Research Octane Number (RON) and Motor Octane Number (MON), as fuel properties, are sufficient to describe a fuel’s knock-limited performance under boosted spark-ignition (SI) conditions within the framework of CFPH. To this end, four TPRF-bioblendstock surrogates having different compositions but matched RON (=98) and MON (=90), were first generated using a non-linear regression model based on artificial neural network (ANN). Additionally, three unconventional bioblendstocks were included in the analysis: di-isobutylene (DIB), isobutanol, and Anisole. Skeletal reaction mechanisms were generated for the TPRF-DIB, TPRF-isobutanol, and TPRF-anisole blends from a detailed kinetic mechanism. Thereafter, numerical simulations were performed for the fuel surrogates using the skeletal mechanisms and a virtual cooperative fuel research (CFR) engine model, under a representative boosted operating condition. In the computational fluid dynamics (CFD) model, the G-equation approach was employed to track the turbulent flame front and the well-stirred reactor model combined with the multi-zone binning strategy was used to capture auto-ignition in the end-gas. In addition, laminar flame speed (LFS) was tabulated for each blend as a function of pressure, temperature, and equivalence ratio a priori, and the lookup tables were used to prescribe laminar flame speed as an input to the G-equation model. Parametric spark timing sweeps were performed for each fuel blend to determine the corresponding knock-limited spark advance (KLSA) and 50% burn point (CA50) at the respective KLSA timing. It was observed that despite same RON, MON, and engine operating conditions, the TPRF-anisole blend exhibited markedly different knock-limited performance from the other three blends. This deviation from the octane index (OI) expectation was shown to be caused by differences in laminar flame speed. However, it was found that relatively large fuel-specific differences in LFS (>20%) would have to be present to cause any appreciable deviation from the OI framework. Otherwise, RON and MON would still be robust enough to predict a fuel’s knock-limited performance.

42 ENGINEERING↗

Factors affecting powerhouse passage of spring migrant smolts at federally operated hydroelectric dams of the Snake and Columbia rivers

From 2008 to 2018, acoustic telemetry studies were conducted to evaluate dam passage survival of spring migrant Chinook salmon and steelhead smolts at seven of the eight federally operated dams on the lower Snake and Columbia rivers. Data from over 87 000 dam passage events were evaluated using regression modeling to identify the effect of spill operations, environmental conditions, and fish characteristics on powerhouse passage probability. In general, powerhouse passage was positively correlated with discharge, negatively correlated with forebay temperature and fish size, and higher for fish that passed the dam at night and for those that approached from the powerhouse side of the river, suggesting powerhouse passage is largely a function of smolt activity level and swimming ability. As such, spilling large volumes of water to reduce powerhouse passage is likely to be most effective during times of reduced activity and swimming ability (e.g., at night, high flows, and cold temperatures). This information can be used to develop dam- and time-specific spill operations that optimize smolt passage, power generation, and other competing demands, such as adult passage.

60 APPLIED LIFE SCIENCES↗

Quantifying Uncertainty in HPC Job Queue Time Predictions

High Performance Computing (HPC) has developed at an unprecedented pace in recent decades. This growth has demanded corresponding development in the area of HPC Operational Data Analytics (ODA), which encompasses a wide range of data analysis techniques, ML/AI efforts, tools, and visualizations. Published studies in ODA offer a variety of practical ways to inform HPC users, administrators, procurement managers, and other stakeholders. Uncertainty analysis, however, is rare in the related published literature. For instance, we identify only 1 out of 14 existing studies focused on job queue time prediction that investigates the uncertainty aspect of their proposed predictions. We recognize the utmost importance uncertainty quantification can have in such predictive analytics solutions, with consequences in how users interpret information they receive, and attempt to bridge this gap. With the goal of improving access to such insights, we develop a process for determining upper and lower bounds of the predicted queue times of a regression model at a specified confidence level. Our current research is focused on the uncertainty in predicting job queue times, yet our approach may be employed in predicting other metrics.

HPC↗

Did You Win the GPU Cloud Lottery? Benchmarking from TFLOPS to Tokens/$

Cloud GPUs are commonly assumed to deliver consistent performance for a given GPU model. This assumption does not always hold: cloud providers employ diverse system configurations and virtualization mechanisms, and GPUs themselves exhibit non-negligible manufacturing variability (the silicon lottery). In this work, we present a large-scale measurement study of GPU performance variability across 11 cloud providers, covering over 3,500 physical GPUs and 6,800 benchmark runs. Our hierarchical analysis shows that while execution-level variation stays below 9%, performance varies by up to 38% across devices and providers for the same GPU model. Regression analysis indicates that driver- and OS-related software factors contribute less than 1% of the variance; instead, silicon lottery effects dominate observed performance variation, and cloud providers further amplify them through persistent, systematic second-order effects.

Slynko, Platon [Silicon Data, New York, USA] (ORCI↗

The updated ITPA global H-mode confinement database: description and analysis

The multi-machine ITPA Global H-mode Confinement Database has been upgraded with new data from JET with the ITER-like wall and ASDEX Upgrade with the full tungsten wall. This paper describes the new database and presents results of regression analysis to estimate the global energy confinement scaling in H-mode plasmas using a standard power law. Various subsets of the database are considered, focusing on type of wall and divertor materials, confinement regime (all H-modes, ELMy H or ELM-free) and ITER-like constraints. Apart from ordinary least squares, two other, robust regression techniques are applied, which take into account uncertainty on all variables. Regression on data from individual devices shows that, generally, the confinement dependence on density and the power degradation are weakest in the fully metallic devices. Using the multi-machine scalings, predictions are made of the confinement time in a standard ELMy H-mode scenario in ITER. The uncertainty on the scaling parameters is discussed with a view to practically useful error bars on the parameters and predictions. One of the derived scalings for ELMy H-modes on an ITER-like subset is studied in particular and compared to the IPB98(y,2) confinement scaling in engineering and dimensionless form. Transformation of this new scaling from engineering variables to dimensionless quantities is shown to result in large error bars on the dimensionless scaling. Regression analysis in the space of dimensionless variables is therefore proposed as an alternative, yielding acceptable estimates for the dimensionless scaling. The new scaling, which is dimensionally correct within the uncertainties, suggests that some dependencies of confinement in the multi- machine database can be reconciled with parameter scans in individual devices. This includes vanishingly small dependence of confinement on line-averaged density and normalized plasma pressure (β), as well as a noticeable, positive dependence on effective atomic mass and plasma triangularity. Extrapolation of this scaling to ITER yields a somewhat lower confinement time compared to the IPB98(y, 2) prediction, possibly related to the considerably weaker dependence on major radius in the new scaling (slightly above linear). Further studies are needed to compare more flexible regression models with the power law used here. In addition, data from more devices concerning possible ‘hidden variables’ could help to determine their influence on confinement, while adding data in sparsely populated areas of the parameter space may contribute to further disentangling some of the global confinement dependencies in tokamak plasmas.

Database↗

cjohnson-LANL/GRL_Kilauea

The python routines are outlined in detail to perform the methods and results in the manuscript under review in the journal Geophysical Research Letters titled “Seismic features predict ground motions during repeating caldera collapse sequence” with LA-UR-23-33345. All routines are written in open source python and were applied to publicly available data sets. The codes formats the data into the appropriate structure required to train a boosted tree regression model. Other codes produce figure results.

Johnson, Christopher W↗