Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “statistical change detection”

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 109 records · Page 6

Sensitive neutron transverse polarization analysis using a 3 He spin filter

Here we report an experimental implementation for neutron transverse polarization analysis that is capable of detecting a small angular change ($\ll$10 –3 rad) in neutron spin orientation. This approach is demonstrated for monochromatic beams, and we show that it could be extended to polychromatic neutron beams. Our approach employs a 3 He spin filter inside a solenoid with an analyzing direction perpendicular to the incident neutron polarization direction. The method was tested with polarized neutron beams and a spin rotator placed inside a μ -metal shield just upstream of the analyzer. No cryogenic superconducting shields or additional neutron spin manipulations are needed. With a counting detector, we experimentally show that the angular resolution δ θ = 1 / ( P n A N ) rad is only determined by the counting statistics for the total counts N and the product of the neutron polarization P n and the analyzing power A . With a high-flux neutron beam, 10 –6 rad angular sensitivity is feasible within a day. This simple, classical-quantum-limited transverse polarization analysis scheme may reduce the overall complexity of experimental implementation for applications requiring sensitive neutron polarimetry and improve the precision in fundamental science studies and polarized neutron imaging.

47 OTHER INSTRUMENTATION↗

Time-Based CAN IDS Paper Results Code

Modern vehicles are complex cyber-physical systems made of hundreds of electronic control units (ECUs) that communicate over controller area networks (CANs). This inherited complexity has expanded the CAN attack surface which is vulnerable to message injection attacks. These injections change the overall timing characteristics of messages on the bus, and thus, to detect these malicious messages, time-based intrusion detection systems (IDSs) have been proposed. However, time-based IDSs are usually trained and tested on low-fidelity datasets with unrealistic, labeled attacks. This makes difficult the task of evaluating, comparing, and validating IDSs. Here we detail and benchmark four time-based IDSs against the newly published ROAD dataset, the first open CAN IDS dataset with real (non-simulated) stealthy attacks with physically verified effects. We found that methods that perform hypothesis testing by explicitly estimating message timing distributions have lower performance than methods that seek anomalies in a distribution related statistic. In particular, these “distribution-agnostic” based methods outperform “distribution-based” methods by at least 55% in area under the precision-recall curve (AUC-PR). Our results expand the body of knowledge of CAN time-based IDSs by providing details of these methods and reporting their results when tested on datasets with real advanced attacks. Finally, we develop an after-market plug-in detector using lightweight hardware, which can be used to deploy the best performing IDS method on nearly any vehicle.

Moriano, Pablo [Oak Ridge National Lab. (ORNL), Oa↗

Robust Anthropogenic Signal Identified in the Seasonal Cycle of Tropospheric Temperature

Previous work identified an anthropogenic fingerprint pattern in T AC (x, t), the amplitude of the seasonal cycle of mid- to upper-tropospheric temperature (TMT), but did not explicitly consider whether fingerprint identification in satellite T AC (x, t) data could have been influenced by real-world multidecadal internal variability (MIV). Here we address this question here using large ensembles (LEs) performed with five climate models. LEs provide many different sequences of internal variability noise superimposed on an underlying forced signal. Despite differences in historical external forcings, climate sensitivity, and MIV properties of the five models, their T AC (x, t) fingerprints are similar and statistically identifiable in 239 of the 240 LE realizations of historical climate change. Comparing simulated and observed variability spectra reveals that consistent fingerprint identification is unlikely to be biased by model underestimates of observed MIV. Even in the presence of large (factor of 3–4) intermodel and inter-realization differences in the amplitude of MIV, the anthropogenic fingerprints of seasonal cycle changes are robustly identifiable in models and satellite data. This is primarily due to the fact that the distinctive, global-scale fingerprint patterns are spatially dissimilar to the smaller-scale patterns of internal T AC (x, t) variability associated with the Atlantic multidecadal oscillation and El Niño–Southern Oscillation. The robustness of the seasonal cycle detection and attribution results shown here, taken together with the evidence from idealized aquaplanet simulations, suggest that basic physical processes are dictating a common pattern of forced T AC (x, t) changes in observations and in the five LEs. The key processes involved include GHG-induced expansion of the tropics, lapse-rate changes, land surface drying, and sea ice decrease.

54 ENVIRONMENTAL SCIENCES↗

Evaluation of CMIP6 models in simulating the statistics of extreme precipitation over Eastern Africa

We report the Eastern Africa region experiences frequent extreme precipitation events that can cause destruction of property and environment, and loss of lives. Thus, there is a need to understand how these events may change in the future and how well the global climate models that are used to make projections can simulate precipitation extremes in this region before they can be used in downscaling or flood and drought impact assessment studies. In this work, we evaluated the ability of sixteen Coupled Model Intercomparison Project Phase 6 (CMIP6) models to simulate present-day precipitation extremes over the Eastern Africa region during the two rainy seasons (March–May and September–November). We used nine extreme precipitation indices (including seven (one) indices of wet (dry) extremes) defined by the Expert Team on Climate Change Detection and Indices. The CMIP6 models were evaluated against two gridded observation datasets: Global Precipitation Climatology Project One-Degree Daily Dataset and Tropical Rainfall Measuring Mission Multi-satellite Precipitation Analysis 3B42. Three model performance metrics (percentage bias, normalized root-mean-square error, and pattern correlation coefficient) were employed to further assess the strengths and weakness of the models. Our results show that the multi-model ensemble mean generally provides a better representation of observed precipitation and related extremes compared to individual models when considering all metrics and seasons. Several consistent biases are evident across CMIP6 models, which tend to overestimate the total-wet day precipitation and consecutive wet days, and underestimate very wet days and maximum 5-day precipitation in both seasons. Furthermore, no single model consistently performs best, model performance varies with the season and index under consideration and is generally independent of horizontal resolution.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Data Analysis of X-ray Panels Used in X-ray CT Scientific Instruments

This report is a follow-up study on M. Skeate’s study on System Drift Detection for Health Monitoring (Skeate, 2023). The following report highlights analyses of the panel’s behavior with extended use. This includes effects of burn-in on the panel from over-use and bad pixels (defined in Methodology.) Using data-forward statistical analyses across a sequence of scans, and given that the other conditions present in the data can be replicated, this study show that overtime use of panels does not affect bad pixel count. Also, I present conclusive evidence of a direct relationship between the regions of the panel that are exposed to radiation and burn-in damage to that region overtime. Additionally, this study also highlights the effect of the duty cycle on the dark current change in the panel.

42 ENGINEERING↗

Screening and qualification methodology for SiC end plug processing methods

Deployment of SiC-ceramic-based fuel cladding for light water reactors requires a hermetic end plug–to–cladding joint that can withstand neutron irradiation during normal operation and maintain integrity during design-basis accidents. Reactor experiments have shown that some SiC composite tubes with SiC end plugs can retain hermeticity after irradiation. However, achieving consistent joint performance under irradiation remains a key challenge. Resolving this issue is essential to enable integral irradiation testing and to demonstrate fuel integrity under commercial-reactor irradiation conditions. This report aims to: (1) provide guidance for designing radiation-tolerant end plug joints for SiC cladding; (2) demonstrate experimental methods to detect processing defects that are unstable under neutron irradiation at light-water-reactor-relevant temperatures and doses; and (3) outline a step-by-step approach for designing and conducting reactor experiments to screen joining methods. The resulting data will be used to improve joint processing and to define critical defect types and sizes that must be detected and eliminated through non-destructive evaluation for quality assurance. Based on prior irradiation experiments at the High Flux Isotope Reactor, differential swelling among the cladding, bonding layer, and end plug was identified as an underlying mechanism for irradiation-induced joint degradation. Accordingly, this effect must be considered in the design of radiation-tolerant joining techniques. In this work, miniature SiC end plug joint specimens irradiated during the previous project were analyzed using X-ray computed tomography to characterize the joint microstructure. Digital volume correlation of the tomography data quantified radiation-induced microstructural changes and enabled evaluation of defect-related risks. Finally, ongoing neutron irradiation efforts using larger specimen volumes are presented. These efforts aim to statistically assess joint performance and to build a microstructure–performance (e.g., leak-tightness) dataset to inform processing improvements and quality control.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DECADE+DES Y3 Weak Lensing Mass Map: A 13,000 deg$^2$ View of Cosmic Structure from 270 Million Galaxies

We present the largest galaxy weak lensing mass map of the late-time Universe, reconstructed from 270 million galaxies in the DECADE and DES Year 3 datasets, covering 13,000 square degrees. We validate the map through systematic tests against observational conditions (depth, seeing, etc.), finding the map is statistically consistent with no contamination. The large area covered by the mass map makes it a well-suited tool for cosmological analyses, cross-correlation studies and the identification of large-scale structure features. We demonstrate its potential by detecting cosmic filaments directly from the mass map for the first time and validating them through their association with galaxy clusters selected using the Sunyaev-Zeldovich effect from Planck and ACT DR6.

Gatti, M. [Chicago U., KICP] (ORCID:00000001613487↗

Learning new physics from data: A symmetrized approach

Thousands of person years have been invested in searches for new physics (NP), the majority of them motivated by theoretical considerations. Yet, no evidence of beyond the Standard Model physics has been found. This suggests that model-agnostic searches might be an important key to explore NP, and help discover unexpected phenomena which can inspire future theoretical developments. A possible strategy for such searches is identifying asymmetries between data samples that are expected to be symmetric within the Standard Model. We propose exploiting neural networks (NNs) to quickly fit and statistically test the differences between two samples. Our method is based on an earlier work, originally designed for inferring the deviations of an observed dataset from that of a much larger reference dataset. We present a symmetric formalism, generalizing the original one, avoiding fine-tuning of the NN parameters and any constraints on the relative sizes of the samples. Our formalism could be used to detect small symmetry violations, extending the discovery potential of current and future particle physics experiments.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Robust neural network-enhanced estimation of local primordial non-Gaussianity

When applied to the nonlinear matter distribution of the universe, neural networks have been shown to be very statistically sensitive probes of cosmological parameters, such as the linear perturbation amplitude σ 8 . However, when used as a “black box,” neural networks are not robust to baryonic uncertainty. We propose a robust architecture for constraining primordial non-Gaussianity f NL , by training a neural network to locally estimate σ 8 , and correlating these local estimates with the large-scale density field. We apply our method to N-body simulations, and show that σ⁡(f NL ) is 3.5 times better than the constraint obtained from a standard halo-based approach. Finally, we show that our method has the same robustness property as large-scale halo bias: baryonic physics can change the normalization of the estimated f NL , but cannot change whether f NL is detected.

79 ASTRONOMY AND ASTROPHYSICS↗

SDSS IV MaNGA: visual morphological and statistical characterization of the DR15 sample

ABSTRACT We present a detailed visual morphological classification for the 4614 MaNGA galaxies in SDSS Data Release 15, using image mosaics generated from a combination of r band (SDSS and deeper DESI Legacy Surveys) images and their digital post-processing. We distinguish 13 Hubble types and identify the presence of bars and bright tidal debris. After correcting the MaNGA sample for volume completeness, we calculate the morphological fractions, the bi-variate distribution of type and stellar mass M* – where we recognize a morphological transition ‘valley’ around S0a-Sa types – and the variations of the g − i colour and luminosity-weighted age over this distribution. We identified bars in 46.8 per cent of galaxies, present in all Hubble types later than S0. This fraction amounts to a factor ∼2 larger when compared with other works for samples in common. We detected 14 per cent of galaxies with tidal features, with the fraction changing with M* and morphology. For 355 galaxies, the classification was uncertain; they are visually faint, mostly of low/intermediate masses, low concentrations, and discy in nature. Our morphological classification agrees well with other works for samples in common, though some particular differences emerge, showing that our image procedures allow us to identify a wealth of added value information as compared to SDSS-based previous estimates. Based on our classification, we also propose an alternative criteria for the E–S0 separation, in the structural semimajor to semiminor axis versus bulge to total light ratio (b/a − B/T) and concentration versus semimajor to semiminor axis (C − b/a) space.

Vázquez-Mata, J. A. (ORCID:0000000186941204)↗

The case for digital twins in metal additive manufacturing

The digital twin (DT) is a relatively new concept that is finding increased acceptance in industry. A DT is generally considered as comprising a physical entity, its virtual replica, and two-way digital data communications in-between. Its primary purpose is to leverage the process intelligence captured within digital models—or usually their faster-solving surrogates—towards generating increased value from the physical entities. The surrogate models are created using machine learning based on data obtained from the field, experiments and digital models, which may be physics-based or statistics-based. Anomaly detection and correction, and diagnostic closed-loop process control are examples of how a process DT can be deployed. In the manufacturing industry, its use can achieve improvements in product quality and process productivity. Metal additive manufacturing (AM) stands to gain tremendously from the use of DTs. This is because the AM process is inherently chaotic, resulting in poor repeatability. However, a DT acting in a supervisory role can inject certainty into the process by actively keeping it within bounds through real-time control commands. Closed-loop feedforward control is achieved by observing the process through sensors that monitor critical parameters and, if there are any deviations from their respective optimal ranges, suitable corrective actions are triggered. The type of corrective action (e.g. a change in laser power or a modification to the scanning speed) and its magnitude are determined by interrogating the surrogate models. Because of their artificial intelligence (AI)-endowed predictive capabilities, which allow them to foresee a future state of the physical twin (e.g. the AM process), DTs proactively take context-sensitive preventative steps, whereas traditional closed-loop feedback control is usually reactive. Apart from assisting a build process in real-time, a DT can help with planning the build of a part by pinpointing the optimum processing window relevant to the desired outcome. Again, the surrogate models are consulted to obtain the required information. In this article, we explain how the application of DTs to the metal AM process can significantly widen its application space by making the process more repeatable (through quality assurance) and cheaper (by getting builds right the first time).

36 MATERIALS SCIENCE↗

The influence of environmental microseismicity on detection and interpretation of small-magnitude events in a polar glacier setting

Glacial environments exhibit temporally variable microseismicity. To investigate how microseismicity influences event detection, we implement two noise-adaptive digital power detectors to process seismic data from Taylor Glacier, Antarctica. We add scaled icequake waveforms to the original data stream, run detectors on the hybrid data stream to estimate reliable detection magnitudes and compare analytical magnitudes predicted from an ice crack source model. We find that detection capability is influenced by environmental microseismicity for seismic events with source size comparable to thermal penetration depths. When event counts and minimum detectable event sizes change in the same direction (i.e. increase in event counts and minimum detectable event size), we interpret measured seismicity changes as ‘true’ seismicity changes rather than as changes in detection. Generally, one detector (two degree of freedom (2dof)) outperforms the other: it identifies more events, a more prominent summertime diurnal signal and maintains a higher detection capability. We conclude that real physical processes are responsible for the summertime diurnal inter-detector difference. One detector (3dof) identifies this process as environmental microseismicity; the other detector (2dof) identifies it as elevated waveform activity. Our analysis provides an example for minimizing detection biases and estimating source sizes when interpreting temporal seismicity patterns to better infer glacial seismogenic processes.

54 ENVIRONMENTAL SCIENCES↗

First Limits on Light Dark Matter Interactions in a Low Threshold Two-Channel Athermal Phonon Detector from the TESSERACT Collaboration

We present results of a search for spin-independent dark matter-nucleus interactions in a 1 cm 2 by 1 mm thick (0.233 g) high-resolution silicon athermal phonon detector operated above ground. For interactions in the substrate, this detector achieves an rms baseline energy resolution of 361.5⁢(4) m⁢ eV (statistical error), the best for any athermal phonon detector to date. With an exposure of 0.233 g ×12 hours, we place the most stringent constraints on dark matter masses between 44 and 87 M⁢ eV/c 2 , with the lowest unexplored cross section of 4⁢ × 10 −32 c⁢m 2 at 87 M⁢ eV/c 2 . We employ a conservative salting technique to reach the lowest dark matter mass ever probed via direct detection experiment. This constraint is enabled by two-channel rejection of low energy backgrounds that are coupled to individual sensors.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Predictive analytics to direct clinical attention to complex patients with elevated suicide risk: enhancement of the Veterans Health Administration REACH VET model

Suicide is a major public health concern, particularly among Veterans. The U.S. Department of Veterans Affairs Veterans Health Administration (VHA) employs the Recovery Engagement and Coordination for Health–Veterans Enhanced Treatment (REACH VET) model to prioritise high-risk patients for targeted clinical attention. REACH VET 1.0 (RV 1.0) was developed on 2008–2011 data. To reflect changes in clinical practice and populations, VHA updated it to REACH VET 2.0 (RV 2.0). This study describes its development and validation. RV 2.0 used longitudinal data from 7,248,170 VHA patients (4,967 suicide deaths) in 2018–2019, with 650 time-varying demographic, clinical and area-level predictors derived from a 2-year lookback (2016–2019). An ensemble of Elastic-Net logistic regression models was trained on 2018 data and evaluated monthly at the population level in 2019, focusing on the top 0.1% intervention risk tier. Analyses assessed model discrimination, suicide detection, risk concentration, subgroup consistency (sex, age and race/ethnicity) and performance relative to RV 1.0 using the same percentile-based risk strata. RV 2.0 outperformed RV 1.0 across all risk strata, with better discrimination (C-statistic 0.76 vs 0.69) and consistent performance across demographic subgroups. Within the top 0.1% of predicted risk, RV 2.0 identified more deaths, higher suicide rates and greater mortality risk concentration both when averaged across the 12 monthly 2019 test sets (5.6 vs 3.6; 83.6 vs 53.7 per 100,000 person-years; 21.0 vs 14.1) and when annualised for 2019 (67 vs 43; 2.7% vs 1.7%; 1,003 vs 644 per 100,000 person-years; 26.7 vs 17.1). RV 2.0 improves suicide risk stratification among Veterans, demonstrating better performance and consistent prediction across subgroups and highlighting the need for regular model updates and evaluation.

Peluso, Alina [Oak Ridge National Laboratory (ORNL↗

Defining Lipidomic Responses to Coronavirus Infection

Highly pathogenic human coronavirus infection can cause a severe atypical, rapid onset pneumonia with a mortality rate of up to 10% for severe acute respiratory syndrome coronavirus 1 (SARS-CoV 1), 35% for Middle East respiratory syndrome coronavirus (MERS-CoV), or 1% for severe acute respiratory syndrome coronavirus 2 (SARS-CoV 2 causative agent of COVID 19). While medical countermeasures successfully controlled the worldwide SARS-CoV epidemic, the MERS-CoV epidemic is still ongoing and continues to be a concern for travelers in the Middle East and the multi-year SARS-CoV 2 pandemic underscore the importance of defining biomarkers that are diagnostic and/or predictive of severe disease outcomes for respiratory viruses. Systems biology approaches provide global snapshots of infection induced changes in host cells/tissues and provide extremely rich datasets for understanding host pathogen interactions. Metabolites, especially lipids, are critical for viral replication but less is understood about infection induced changes to lipids due to limits in lipid species detection and identification. To characterize how individual lipid species and proteins with lipid associated functions contribute to highly pathogenic human coronavirus replication and disease severity, existing datasets were probed to determine cell type specific lipid responses to MERS-CoV infection and verification studies were performed to determine if modification of the host lipid signature (by inhibiting enzymatic functions that produce specific lipid species) can perturb CoV replication in human lung cells. MERS-CoV infects human lung epithelial, endothelial and fibroblast cells. All three cell types were infected with MERS-CoV and samples collected to analyze lipids, proteins, metabolites, and transcripts from 12 to 48 hours post infection. Time matched mock-infected cells were collected in parallel for each cell type. Following sample and statistical analysis, functional enrichment and bioinformatic analysis was performed to identify differentially expressed lipids and proteins. Two lipid species were found to be significantly upregulated following MERS-CoV infection, ceramides, and triacylglycerol both of which are indicative of cells undergoing apoptotic cell death. In contrast, sphingomyelins (lipid molecules that can serve as a precursor for one pathway for ceramide synthesis) had significantly decreased expression in MERS-CoV infected cells. An inhibitor of acid sphingomyelinase (that converts sphingomyelin to ceramides and phosphorylcholine) reduced MERS-CoV replication suggesting that production of ceramides is key for successful viral replication and transmission. Acyl-CoA-synthetase 3 (ACSL3), the only protein whose function is lipid associated and had increased differential expression in our dataset, regulates the synthesis of triacylglycerol (increased expression). Inhibitors that directly block ACSL3 (Triacsin C) but not steps later in the triacylglycerol synthesis pathway (Etomoxir) inhibit MERS-CoV replication suggesting that triacylglycerol production and/or ACSL3 activity is also key for viral replication. ACSL3 expression is also upregulated in lung cancer cells and is predicted to promote continued cell viability which would also enhance viral replication. The differentially expressed lipid and lipid-associated protein species identified in our studies suggest that MERS-CoV infection is activating cellular death pathways to limit the number of cells that become infected but also stimulating the production of lipid-associated enzymes that can prolong host cell viability and the amount of time progeny virions can be produced and released. As the inhibitors that worked against MERS-CoV infection were also efficacious against SARS-CoV 2 infection, countermeasures that target these host pathways may provide novel ways to block highly pathogenic human coronavirus infection and/or prevent severe disease outcomes.

59 BASIC BIOLOGICAL SCIENCES↗

Unique & challenging aspects of plutonium metal standards exchange program for actinide measurements

The Los Alamos National Laboratory exchange program is the only program of its kind for the distribution of plutonium (Pu) standards materials with a range of impurity contents to multiple laboratories for destructive measurements of elemental concentration. This paper discusses statistical methods used to address challenges in Pu metal exchange data by way of two case studies. Challenges include how to evaluate a data set when a large fraction of the values are minimum detection limits (MDLs), and how to determine potential outliers with limited in-formation on the true spread of the data.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DNA Viral Diversity, Abundance, and Functional Potential Vary across Grassland Soils with a Range of Historical Moisture Regimes

Soil viruses are abundant, but the influence of the environment and climate on soil viruses remains poorly understood. Here, we addressed this gap by comparing the diversity, abundance, lifestyle, and metabolic potential of DNA viruses in three grassland soils with historical differences in average annual precipitation, low in eastern Washington (WA), high in Iowa (IA), and intermediate in Kansas (KS). Bioinformatics analyses were applied to identify a total of 2,631 viral contigs, including 14 complete viral genomes from three deep metagenomes (1 terabase [Tb] each) that were sequenced from bulk soil DNA. An additional three replicate metagenomes (~0.5 Tb each) were obtained from each location for statistical comparisons. Identified viruses were primarily bacteriophages targeting dominant bacterial taxa. Both viral and host diversity were higher in soil with lower precipitation. Viral abundance was also significantly higher in the arid WA location than in IA and KS. More lysogenic markers and fewer clustered regularly interspaced short palindromic repeats (CRISPR) spacer hits were found in WA, reflecting more lysogeny in historically drier soil. More putative auxiliary metabolic genes (AMGs) were also detected in WA than in the historically wetter locations. The AMGs occurring in 18 pathways could potentially contribute to carbon metabolism and energy acquisition in their hosts. Structural equation modeling (SEM) suggested that historical precipitation influenced viral life cycle and selection of AMGs. The observed and predicted relationships between soil viruses and various biotic and abiotic variables have value for predicting viral responses to environmental change.

59 BASIC BIOLOGICAL SCIENCES↗

Climate Change Signal in Atlantic Tropical Cyclones Today and Near Future

This manuscript discusses the challenges in detecting and attributing recently observed trends in the Atlantic tropical cyclone (TC) and the epistemic uncertainty we face in assessing future risk. We use synthetic storms downscaled from five CMIP5 models by the Columbia HAZard model (CHAZ), and directly simulated storms from high-resolution climate models. We examine three aspects of recent TC activity: the upward trend and multi-decadal oscillation of the annual frequency, the increase in storm wind intensity, and the decrease in forward speed. Some data sets suggest that these trends and oscillation are forced while others suggest that they can be explained by natural variability. Projections under warming climate scenarios also show a wide range of possibilities, especially for the annual frequencies, which increase or decrease depending on the choice of moisture variable used in the CHAZ model and on the choice of climate model. The uncertainties in the annual frequency lead to epistemic uncertainties in TC risk assessment. Here, we investigate the potential for reduction of these epistemic uncertainties through a statistical practice, namely likelihood analysis. We find that historical observations are more consistent with the simulations with increasing frequency than those with decreasing frequency, but we are not able to rule out the latter. We argue that the most rational way to treat epistemic uncertainty is to consider all outcomes contained in the results. In the context of risk assessment, since the results contain possible outcomes in which TC risk is increasing, this view implies that the risk is increasing.

54 ENVIRONMENTAL SCIENCES↗