SEARCH · Engineering Papers

Results for “causal analysis”

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.

SourceResearch connectionFrom yearThrough year

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

At least 145 records · Page 8

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Physics-guided machine learning approaches to predict the ideal stability properties of fusion plasmas

One of the biggest challenges to achieve the goal of producing fusion energy in tokamak devices is the necessity of avoiding disruptions of the plasma current due to instabilities. The disruption event characterization and forecasting (DECAF) framework has been developed in this purpose, integrating physics models of many causal events that can lead to a disruption. Two different machine learning approaches are proposed to improve the ideal magnetohydrodynamic (MHD) no-wall limit component of the kinetic stability model included in DECAF. First, a random forest regressor (RFR), was adopted to reproduce the DCON computed change in plasma potential energy without wall effects for a large database of equilibria from the national spherical torus experiment (NSTX). This tree-based method provides an analysis of the importance of each input feature, giving an insight into the underlying physics phenomena. Secondly, a fully-connected neural network has been trained on sets of calculations with the DCON code, to get an improved closed form equation of the no-wall β limit as a function of the relevant plasma parameters indicated by the RFR. The neural network has been guided by physics theory of ideal MHD in its extension outside the domain of the NSTX experimental data. The estimated value has been incorporated into the DECAF kinetic stability model and tested against a set of experimentally stable and unstable discharges. Moreover, the neural network results were used to simulate a real-time stability assessment using only quantities available in real-time. Finally, the portability of the model was investigated, showing encouraging results by testing the NSTX-trained algorithm on the mega ampere spherical tokamak (MAST).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Development of a Unified Taxonomy for HVAC System Faults

Detecting and diagnosing HVAC faults is critical for maintaining building operation performance, reducing energy waste, and ensuring indoor comfort. An increasing deployment of commercial fault detection and diagnostics (FDD) software tools in commercial buildings in the past decade has significantly increased buildings’ operational reliability and reduced energy consumption. A massive amount of data has been generated by the FDD software tools. However, efficiently utilizing FDD data for ‘big data’ analytics, algorithm improvement, and other data-driven applications is challenging because the format and naming conventions of those data are very customized, unstructured, and hard to interpret. This paper presents the development of a unified taxonomy for HVAC faults. A taxonomy is an orderly classification of HVAC faults according to their characteristics and causal relations. The taxonomy includes fault categorization, physical hierarchy, fault library, relation model, and naming/tagging scheme. The taxonomy employs both a physical hierarchy of HVAC equipment and a cause-effect relationship model to reveal the root causes of faults in HVAC systems. A structured and standardized vocabulary library is developed to increase data representability and interpretability. The developed fault taxonomy can be used for HVAC system ‘big data’ analytics such as HVAC system fault prevalence analysis or the development of an HVAC FDD software standard. A common type of HVAC equipment-packaged rooftop unit (RTU) is used as an example to demonstrate the application of the developed fault taxonomy. Two RTU FDD software tools are used to show that after mapping FDD data according to the taxonomy, the meta-analysis of the multiple FDD reports is possible and efficient.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗

DOE OSTI2026

Feedback and Oscillations: Constructing Feedback Systems for Root Cause Analysis of Oscillations in Power Grids

Dynamic phenomena linked to inverter-based resources (IBRs) have gained global attention. Several IBR-induced dynamics have caused bulk power system-connected wind or solar power plants to trip, and some have even led to widespread outages. In addition, many oscillations have been observed involving IBR power plants. In 2023, the IEEE Power & Energy Society (PES) IBR Subsynchronous Oscillations (SSO) task force published a journal article, “Real-World Subsynchronous Oscillation Events in Power Grids With High Penetrations of Inverter-Based Resources,” in which 19 IBR oscillation events were examined for their causation. Earlier in 2020, another PES task force article, “Definition and Classification of Power System Stability-Revisited & Extended,” authored by prominent academics, introduced converter-driven stability as a new category of stability. The international power grid industry community also took action by publishing the CIGRE Green Book, Power System Dynamic Modelling and Analysis in Evolving Networks (led by Babak Badrzadeh and Zia Emin) in 2024. In August 2024, the Energy Systems Integration Group (ESIG) released a practical guide led by Nick Miller, “Diagnosis and Mitigation of Observed Oscillations in IBR-Dominant Power System: A Practical Guide.” The goal of the guide is to assist practicing engineers in making initial judgments and conducting detailed analyses about oscillations. Finally, when addressing the classification of stability and oscillations, the guide emphasizes a causality-based taxonomy for grouping, such as voltage control-induced oscillations, synchronization-induced oscillations, and frequency or active power control-induced oscillations.

Fan, Lingling [Univ. of South Florida, Tampa, FL (↗

DOE OSTI2025

AUTOMATIC GENERATION OF EVENT TREES AND FAULT TREES: A MODEL-BASED APPROACH

In the past few decades, increasing complexity in modern engineering systems has been driven by the integration of a large number of components and by the fact that the system operations involve many disciplines (e.g., thermal-hydraulics, plant operations, cyber-security). Current safety/reliability modeling approaches to such systems are labor intensive, difficult to learn, and rely heavily on simplistic Boolean logic to depict failure propagation and accident progression. While these methods serve well for simple systems (i.e., linear causal systems with limited small inter- and intra-system interactions), their results are difficult to verify when modeling complex systems (typically performed through the extensive use of modeling assumptions). The development of new methods is addressed to meet these challenges through a model-based system engineering (MBSE) lens. Under MBSE philosophy, every aspect of the system (form or function) is represented by a model that completely characterizes its architecture or behavior. MBSE approach greatly improves the management of design, analysis and verification of complex systems. An integration of Dynamic Probabilistic Risk Assessment (DPRA) methods with MBSE models is proposed to perform safety/reliability analyses of engineering systems. In particular, MBSE representation of the system (performed using Systems Modeling Language [SysML]) is coupled with DPRA methods to automatically generate event trees and fault trees.

97 - MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Collagen Turnover in Relation to Risk Factors and Hemodynamics in Human Intracranial Aneurysms

Determinants for molecular and structural instability, that is, impending growth or rupture, of intracranial aneurysms (IAs) remain uncertain. To elucidate this, we endeavored to estimate the actual turnover rates of the main molecular constituent in human IA (collagen) on the basis of radiocarbon ( 14 C) birth dating in relation to IA hemodynamics. Collagen turnover rates in excised human IA samples were calculated using mathematical modeling of 14 C birth dating data of collagen in relation to risk factors and histological markers for collagen maturity/turnover in selected IA. Hemodynamics were simulated using image-based computational fluid dynamics. Correlation, logistic regression, and receiver operating characteristic analyses were performed. Collagen turnover rates were estimated in 46 IA (43 patients); computational fluid dynamics could be performed in 20 IA (20 patients). The mean collagen turnover rate (γ) constituted 126% (±1% error) per year. For patients with arterial hypertension, γ was greater than 2600% annually, whereas γ was distinctly lower with 32% (±1% error) per year for patients without risk factors, such as smoking and hypertension. There was a distinct association between histological presence of rather immature collagen in human IA and the presence of modifiable risk factors. Spatial-temporal averaged wall shear stress predicted rapid collagen turnover (odds ratio, 1.6 [95% CI, 1.0–2.7]). Receiver operating characteristic analysis demonstrated a good test accuracy (area under the curve, 0.798 [95% CI, 0.598–0.998]) for average wall shear stress with a threshold ≥4.9 Pa for rapid collagen turnover. Our data indicate that turnover rates and stability of collagen in human IA are strongly associated with the presence of modifiable risk factors and aneurysmal hemodynamics. Finally, these findings underline the importance of strict risk factor modification in patients with unruptured IA. Future should include more detailed risk factor data to establish a more causal understanding of hemodynamics and the rupture risk of individual IA.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2022

Divertor plasma detachment: roles of plasma momentum, energy, and particle balances

Abstract Contrary to statements from some recent papers, the bifurcation of the scrape-off layer and divertor plasma parameters, as described within a 1D model, does exist. The analysis of both numerical simulations and experimental data shows that the real physical reasons that can cause a strong reduction of the plasma particle flux to the target j d are impurity radiation loss and volumetric plasma recombination. Fitting the parameter describing the ‘plasma momentum removal’, based on the results of 2D numerical simulations, then using this fit in subsequent analytic expressions for j d , and finally interpreting the evolution of j d as the result of elastic plasma-neutral interactions break the causality principle and lead to a wrong physical picture of divertor plasma detachment.

Physics↗

DOE OSTI2025

Virulence and Genetic Diversity of Puccinia spp., Causal Agents of Rust on Switchgrass (Panicum virgatum L.) in the USA

Switchgrass (Panicum virgatum L.) is an important cellulosic biofuel grass native to North America. Rust, caused by Puccinia spp. is the most predominant disease of switchgrass and has the potential to impact biomass conversion. In this study, virulence patterns were determined on a set of 38 switchgrass genotypes for 14 single-spore rust isolates from 14 field samples collected in seven states. Single nucleotide polymorphism (SNP) variation was also assessed in 720 sequenced cloned amplicons representing 654 base pairs of the elongation factor 1-α gene from the field samples. Five major haplotypes were identified differing by 11 out of the 39 SNP positions identified. STRUCTURE, Principal Coordinate Analysis, and phylogenetic analyses divided the rust population into two genetic clusters. Virginia and Georgia had the highest and lowest rust genetic diversity, respectively. Only nine accessions showed a differential disease response between the 14 isolates, allowing the identification of eight races, differing by 1–3 virulence factors. Overall, the results suggested clonal reproduction of the pathogen and a North–South differentiation via local adaptation. However, similar haplotypes and races were also recovered from several states, suggesting migration events, and highlighting the need to further investigate the switchgrass rust population structure and evolution in the USA.

Bahri, Bochra A. (ORCID:0000000159055880)↗

DOE OSTI2024

Predicting nutrition and environmental factors associated with female reproductive disorders using a knowledge graph and random forests

Female reproductive disorders (FRDs) are common health conditions that may present with significant symptoms. Diet and environment are potential areas for FRD interventions. We utilized a knowledge graph (KG) method to predict factors associated with common FRDs (for example, endometriosis, ovarian cyst, and uterine fibroids). We harmonized survey data from the Personalized Environment and Genes Study (PEGS) on internal and external environmental exposures and health conditions with biomedical ontology content. We merged the harmonized data and ontologies with supplemental nutrient and agricultural chemical data to create a KG. We analyzed the KG by embedding edges and applying a random forest for edge prediction to identify variables potentially associated with FRDs. We also conducted logistic regression analysis for comparison. Across 9765 PEGS respondents, the KG analysis resulted in 8535 significant or suggestive predicted links between FRDs and chemicals, phenotypes, and diseases. Amongst these links, 32 were exact matches when compared with the logistic regression results, including comorbidities, medications, foods, and occupational exposures. Mechanistic underpinnings of predicted links documented in the literature may support some of our findings. Our KG methods are useful for predicting possible associations in large, survey-based datasets with added information on directionality and magnitude of effect from logistic regression. These results should not be construed as causal but can support hypothesis generation. This investigation enabled the generation of hypotheses on a variety of potential links between FRDs and exposures. Future investigations should prospectively evaluate the variables hypothesized to impact FRDs.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2022

Effects of vessel traffic and ocean noise on gray whale stress hormones

Human use of marinescapes is rapidly increasing, especially in populated nearshore regions where recreational vessel traffic can be dense. Marine animals can have a physiological response to such elevated human activity that can impact individual health and population dynamics. To understand the physiological impacts of vessel traffic on baleen whales, we investigated the adrenal stress response of gray whales (Eschrichtius robustus) to variable vessel traffic levels through an assessment of fecal glucocorticoid metabolite (fGC) concentrations. This analysis was conducted at the individual level, at multiple temporal scales (1–7 days), and accounted for factors that may confound fGC: sex, age, nutritional status, and reproductive state. Data were collected in Oregon, USA, from June to October of 2016–2018. Results indicate significant correlations between fGC, month, and vessel counts from the day prior to fecal sample collection. Furthermore, we show a significant positive correlation between vessel traffic and underwater ambient noise levels, which indicates that noise produced by vessel traffic may be a causal factor for the increased fGC. This study increases knowledge of gray whale physiological response to vessel traffic and may inform management decisions regarding regulations of vessel traffic activities and thresholds near critical whale habitats.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2022

Identification of driver genes for critical forms of COVID-19 in a deeply phenotyped young patient cohort

The drivers of critical coronavirus disease 2019 (COVID-19) remain unknown. Given major confounding factors such as age and comorbidities, true mediators of this condition have remained elusive. We used a multi-omics analysis combined with artificial intelligence in a young patient cohort where major comorbidities were excluded at the onset. The cohort included 47 “critical” (in the intensive care unit under mechanical ventilation) and 25 “non-critical” (in a non-critical care ward) patients with COVID-19 and 22 healthy individuals. The analyses included whole-genome sequencing, whole-blood RNA sequencing, plasma and blood mononuclear cell proteomics, cytokine profiling, and high-throughput immunophenotyping. An ensemble of machine learning, deep learning, quantum annealing, and structural causal modeling were used. Patients with critical COVID-19 were characterized by exacerbated inflammation, perturbed lymphoid and myeloid compartments, increased coagulation, and viral cell biology. Among differentially expressed genes, we observed up-regulation of the metalloprotease ADAM9. This gene signature was validated in a second independent cohort of 81 critical and 73 recovered patients with COVID-19 and was further confirmed at the transcriptional and protein level and by proteolytic activity. Ex vivo ADAM9 inhibition decreased severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) uptake and replication in human lung epithelial cells. In conclusion, within a young, otherwise healthy, cohort of individuals with COVID-19, we provide the landscape of biological perturbations in vivo where a unique gene signature differentiated critical from non-critical patients. We further identified ADAM9 as a driver of disease severity and a candidate therapeutic target.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Natural Language Processing-Enhanced Nuclear Industry Operating Experience Data Analysis: Aggregation and Interpretation of Multi-Report Analysis Results

Industry-wide operating experience is a critical source of raw data for reliability and risk model parameter estimations for nuclear power plants. A large portion of operating experience data are failure events stored as reports that contain unstructured data, such as narratives. In current practice, a failure report is usually reviewed and manually coded by analysts. The coding is based on extracting several event characteristics such as system name, component type, sub-part type, failure mode, and failure cause. Event narratives are mostly used to help understand events and extract their characteristics. In this line of research, we aim to maximize the usage of event narratives by leveraging natural language processing (NLP) methods to automatically convert an event narrative to a causal graph. This research has promise to improve physical understanding of failure initiation and propagation and to facilitate use of non-failure data (e.g., near-misses and degradations) to complement the limited data pool of failures. In our previous work, we developed an NLP tool and applied it to analyze a number of licensee event reports submitted by U.S. nuclear power plants to the Nuclear Regulatory Commission. In this paper, we will report our recent research progress in aggregating the results of multiple reports, developing network model(s), and drawing statistical insights.

99 GENERAL AND MISCELLANEOUS↗

DOE OSTI2024

Implementation of a realistic artificial data generator for crash data generation

In this paper, a framework is outlined to generate realistic artificial data (RAD) as a tool for comparing different models developed for safety analysis. The primary focus of transportation safety analysis is on identifying and quantifying the influence of factors contributing to traffic crash occurrence and its consequences. The current framework of comparing model structures using only observed data has limitations. With observed data, it is not possible to know how well the models mimic the true relationship between the dependent and independent variables. Further, real datasets do not allow researchers to evaluate the model performance for different levels of complexity of the dataset. RAD offers an innovative framework to address these limitations. Hence, we propose a RAD generation framework embedded with heterogeneous causal structures that generates crash data by considering crash occurrence as a trip level event impacted by trip level factors, demographics, roadway and vehicle attributes. Within our RAD generator we employ three specific modules: (a) disaggregate trip information generation, (b) crash data generation and (c) crash data aggregation. For disaggregate trip information generation, we employ a daily activity-travel realization for an urban region generated from an established activity-based model for the Chicago region. We use this data of more than 2 million daily trips to generate a subset of trips with crash data. For trips with crashes crash location, crash type, driver/vehicle characteristics, and crash severity. The daily RAD generation process is repeated for generating crash records at yearly or multi-year resolution. In conclusion, the crash databases generated can be employed to compare frequency models, severity models, crash type and various other dimensions by facility type - possibly establishing a universal benchmarking system for alternative model frameworks in safety literature.

42 ENGINEERING↗

DOE OSTI2023

Operation Optimization Using Reinforcement Learning with Integrated Artificial Reasoning Framework

In large and complex systems, operational decision-making requires a systematic analysis with a vast amount of data from both process parameters and component status monitoring. In this paper, we present an integrated artificial reasoning approach for system state transition models that can help operational decision-making with explainable and traceable reasoning. The integrated artificial reasoning framework is a physics-based approach of defining the system structure in a Bayesian network, so we leveraged it in a Markov decision process (MDP) for finding optimal operational solutions. In our proposed framework, the MDP is implemented on a dynamic Bayesian network (DBN), which represents causalities in a system. The multilevel flow modeling was utilized in order to extract these causalities in a more efficient and objective manner. Since multilevel flow modeling is based on the fundamental energy and mass conservation laws, the target system is decomposed into several mass, energy, and information structures, which serve as the basis for a DBN. The MDP consists of the processes of finding a solution for the Bellman equation, which can be derived from the conditional probability equations of the constructed DBN. System operators can capture stochastic system dynamics as multiple subsystem state transitions based on their physical relations and uncertainties coming from the component degradation process or random failures. We analyzed a simplified example system to illustrate finding an optimal operational policy with this approach.

Kim, Junyung↗

DOE OSTI2023

Structure-informed clustering for population stratification in association studies

Background: Identifying variants associated with complex traits is a challenging task in genetic association studies due to linkage disequilibrium (LD) between genetic variants and population stratification, unrelated to the disease risk. Existing methods of population structure correction use principal component analysis or linear mixed models with a random effect when modeling associations between a trait of interest and genetic markers. However, due to stringent significance thresholds and latent interactions between the markers, these methods often fail to detect genuinely associated variants. Results: To overcome this, we propose CluStrat, which corrects for complex arbitrarily structured populations while leveraging the linkage disequilibrium induced distances between genetic markers. It performs an agglomerative hierarchical clustering using the Mahalanobis distance covariance matrix of the markers. In simulation studies, we show that our method outperforms existing methods in detecting true causal variants. Applying CluStrat on WTCCC2 and UK Biobank cohorts, we found biologically relevant associations in Schizophrenia and Myocardial Infarction. CluStrat was also able to correct for population structure in polygenic adaptation of height in Europeans. Conclusions: CluStrat highlights the advantages of biologically relevant distance metrics, such as the Mahalanobis distance, which captures the cryptic interactions within populations in the presence of LD better than the Euclidean distance.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2021

Earth as a Proxy Exoplanet: Deconstructing and Reconstructing Spectrophotometric Light Curves

Point-source spectrophotometric (single-point) light curves of Earth-like planets contain a surprising amount of information about the spatial features of those worlds. Spatially resolving these light curves is important for assessing time-varying surface features and the existence of an atmosphere, which in turn is critical to life on Earth and significant for determining habitability on exoplanets. Given that Earth is the only celestial body confirmed to harbor life, treating it as a proxy exoplanet by analyzing time-resolved spectral images provides a benchmark in the search for habitable exoplanets. The Earth Polychromatic Imaging Camera (EPIC) on the Deep Space Climate Observatory (DSCOVR) provides such an opportunity, with observations of ∼5000 full-disk sunlit Earth images each year at 10 wavelengths with high temporal frequency. We disk-integrate these spectral images to create single-point light curves and decompose them into principal components (PCs). Using machine-learning techniques to relate the PCs to six preselected spatial features, we find that the first and fourth PCs of the single-point light curves, contributing ∼83.23% of the light-curve variability, contain information about low and high clouds, respectively. Surface information relevant to the contrast between land and ocean reflectance is contained in the second PC, while individual land subtypes are not easily distinguishable (<0.1% total light-curve variation). We build an Earth model by systematically altering the spatial features to derive causal relationships to the PCs. This model can serve as a baseline for analyzing Earth-like exoplanets and guide wavelength selection and sampling strategies for future observations.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2021

Classification and computation of extreme events in turbulent combustion

In the design of practical combustion systems, ensuring safety and reliability is an important requirement. For instance, reliably avoiding lean blowout, flame flashback or inlet unstart is critical for ensuring safe operation. Currently, the science of predicting such events is based on prior experience, limited modeling or diagnostic tools and purely statistical approaches. Even though computational and experimental tools for studying combustion devices have vastly advanced in the last three decades, the analysis of such failure events has not been pursued widely. While the use of data for model development and calibration is being widely accepted, the extension to failure events introduces numerous challenges. In particular, the focus here is on so-called data-poor problems, where the cost of generating data is extremely high and is not easily amenable to existing computational and experimental approaches. Data-poor problems are particularly relevant when related to extreme events (also called anomalous events) that can lead to catastrophic failure of the system. It is argued that transient events that describe such failure can have different causal mechanisms. To develop the scientific inference process, a classification of such problems is used to determine specific modeling paths as well as computational tools needed. Research opportunities in the emerging field of extreme event prediction are highlighted in order to identify critical and immediate needs.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Spatio-Temporal Assessment of Heavy-Duty Truck Incident and Inspection Data

Vehicular incidents, especially those involving tractor trailers, are increasing in number every year. These events are extremely costly for fleets, in terms of damage or loss of property, loss of efficiency, and certainly in terms of loss of life. Although the U.S. Department of Transportation (DOT) is responsible for performing inspections, and fleet managers are encouraged to maintain their fleet and participate in regular inspections, it is uncertain whether these inspections are occurring at a frequency that is necessary to prevent incidents. The Federal Motor Carrier Safety Administration (FMCSA) of the DOT manages and maintains the Motor Carrier Management Information System (MCMIS) dataset, which contains all incident and inspection data regarding commercial vehicles in the U.S. The purpose of this preliminary analysis was to explore the MCMIS dataset through spatiotemporal analyses, to uncover findings that may hint at potential improvements in the DOT inspection process and highlight location-specific trends in the dataset. These analyses are novel, as previous research using the MCMIS dataset only examined the data at the state or county level, not at a national scale. The results from the analyses pinpointed specific major metropolitan areas, namely Harris County (Houston), Texas, and three of the New York boroughs (Kings, Queens, and the Bronx), which were found to have increasing incident rates during the study period (2016–2020). An overview of potential causal factors contributing to this increase are provided as well as an overview of the inspection process, and suggestions for improvement relative to the highlighted locations in Texas and New York are also provided. Ultimately, it is suggested that the incorporation of advanced technology and automation may prove beneficial in reducing the occurrence of events that lead to incidents and may also help in the inspection process.

99 GENERAL AND MISCELLANEOUS↗

DOE OSTI2023

Causes and evolution of winter polynyas north of Greenland

During the 42-year period (1979–2020) of satellite measurements, four major winter (December–March) polynyas have been observed north of Greenland: one in December 1986 and three in the last decade, i.e., February of 2011, 2017, and 2018. The 2018 polynya was unparalleled in its magnitude and duration compared to the three previous events. Given the apparent recent increase in the occurrence of these extreme events, this study aims to examine their evolution and causality, in terms of forced versus natural variability. The limited weather station and remotely sensed sea ice data are analyzed combining with output from the fully coupled Regional Arctic System Model (RASM), including one hindcast and two ensemble simulations. We found that neither the accompanying anomalous warm surface air intrusion nor the ocean below had an impact (i.e., no significant ice melting) on the evolution of the observed winter open-water episodes in the region. Instead, the extreme atmospheric wind forcing resulted in greater sea ice deformation and transport offshore, accounting for the majority of sea ice loss in all four polynyas. Our analysis suggests that strong southerly winds (i.e., northward wind with speeds greater than 10 m s –1 ) blowing persistently over the study region for at least 2 d or more were required over the study region to mechanically redistribute some of the thickest Arctic sea ice out of the region and thus to create open-water areas (i.e., a latent heat polynya). To assess the role of internal variability versus external forcing of such events, we carried out and examined results from the two RASM ensembles dynamically downscaled with output from the Community Earth System Model (CESM) Decadal Prediction Large Ensemble (DPLE) simulations. Out of 100 winters in each of the two ensembles (initialized 30 years apart: one in December 1985 and another in December 2015), 17 and 16 winter polynyas were produced north of Greenland, respectively. The frequency of polynya occurrence had no apparent sensitivity to the initial sea ice thickness in the study area pointing to internal variability of atmospheric forcing as a dominant cause of winter polynyas north of Greenland. We assert that dynamical downscaling using a high-resolution regional climate model offers a robust tool for process-level examination in space and time, synthesis with limited observations, and probabilistic forecasts of Arctic events, such as the ones being investigated here and elsewhere.

54 ENVIRONMENTAL SCIENCES↗