Search NASA⌕ Search

SEARCH · Search NASA

Results for “Model initialization”

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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Regional Seismic Travel Time v.3.2.0

The RSTT software computes predictions of the travel time of seismic energy traveling from a source to a receiver through 2.5D models of the seismic velocity distribution within the Earth. The two primary applications for the RSTT library are tomographic inversion studies and seismic event location calculations. In tomographic inversions studies, a seismologist begins with number of source-receiver travel time observations and an initial starting model of the velocity distribution within the Earth. A forward travel time calculator, such as the RSTT library, is used to compute predictions of each observed travel time and all of the residuals (observed minus predicted travel time) are calculated. The Earth model is then modified in some systematic way with the goal of minimizing the residuals. The Earth model obtained in this way is assumed to be a better model than the starting model if it has lower residuals. The other major application for the RSTT library is seismic event location. Given an Earth model, an initial estimate of the location of a seismic event, and some number of observations of seismic travel time thought to have originated from that event, location codes systematically modify the estimate of the location of the event with the goal of minimizing the difference between the observed and predicted travel times. The second application, seismic event location, is routinely implemented by the military as part of its effort to monitor the Earth for nuclear tests conducted by foreign countries.

BALLARD, SANFORD↗
NASA NTRS2010

Geometric Model of a Coronal Cavity

We observed a coronal cavity from August 8-18 2007 during a multi-instrument observing campaign organized under the auspices of the International Heliophysical Year (IHY). Here we present initial efforts to model the cavity with a geometrical streamer-cavity model. The model is based the white-light streamer mode] of Gibson et a]. (2003 ), which has been enhanced by the addition of a cavity and the capability to model EUV and X-ray emission. The cavity is modeled with an elliptical cross-section and Gaussian fall-off in length and width inside the streamer. Density and temperature can be varied in the streamer and cavity and constrained via comparison with data. Although this model is purely morphological, it allows for three-dimensional, multi-temperature analysis and characterization of the data, which can then provide constraints for future physical modeling. Initial comparisons to STEREO/EUVI images of the cavity and streamer show that the model can provide a good fit to the data. This work is part of the effort of the International Space Science Institute International Team on Prominence Cavities

Kucera, Therese A.↗
DOE OSTI2018

Effects of surface initial condition on aqueous corrosion of glass-A study with a model nuclear waste glass

Being a nonequilibrium material, the structure of glass varies with the sample history. Thus, the initial surface condition of a glass can vary with the preparation condition and have a large impact on its reactivity. This paper shows that the aqueous corrosion behavior of international simple glass (ISG) varies depending on the initial surface state. The ISG glass samples were prepared as-polished-only and polished-then-annealed and they were immersed in aqueous solution saturated with soluble SiO 2 at 30°C (modeling a mild condition) and at 90°C (modeling a severe condition). Molecular dynamics simulations were performed to obtain coordination numbers of network formers of ISG to assist oxygen speciation calculations. The surface structures of as-prepared and corroded ISG samples were analyzed using various imaging and spectroscopic techniques. Among these analyses, only the oxygen speciation with x-ray photoelectron spectroscopy showed discernable differences between two uncorroded surfaces with different preparation histories; all other methods could not differentiate the surface preparation history before aqueous corrosion. Such minor difference in chemical structures was found to have a profound impact on corrosion behaviors in the mild condition. In the harsh condition, the surface history dependence was not as drastic as the corrosion in the mild condition. The analysis results of the corroded surfaces suggested that the thickness and structure of the alteration layer formed on ISG in aqueous corrosion can vary with the initial surface state.

nuclear, defects, corrosion↗
DOE OSTI2022

The Application of a Genetic Algorithm to the Optimization of a Mesoscale Model for Emergency Response

Besides solving the equations of momentum, heat, and moisture transport on the model grid, mesoscale weather models must account for subgrid-scale processes that affect the resolved model variables. These are simulated with model parameterizations, which often rely on values preset by the user. Such “free” model parameters, along with others set to initialize the model, are often poorly constrained, requiring that a user select each from a range of plausible values. Finding the values to optimize any forecasting tool can be accomplished with a search algorithm, and one such process—the genetic algorithm (GA)—has become especially popular. As applied to modeling, GAs represent a Darwinian process: an ensemble of simulations is run with a different set of parameter values for each member, and the members subsequently judged to be most accurate are selected as “parents” who pass their parameters onto a new generation. At the U.S. Department of Energy’s Savannah River Site in South Carolina, we are applying a GA to the Regional Atmospheric Modeling System (RAMS) mesoscale weather model, which supplies input to a model to simulate the dispersion of an airborne contaminant as part of the site’s emergency response preparations. An ensemble of forecasts is run each day, weather data are used to “score” the individual members of the ensemble, and the parameters from the best members are used for the next day’s forecasts. As meteorological conditions change, the parameters change as well, maintaining a model configuration that is best adapted to atmospheric conditions. Significance Statement We wanted to develop a forecasting system by which a weather model is run over the Savannah River Site each day and repeatedly adjusted according to how well it performed the previous day. To run the model, a series of values (parameters) must be set to control how the model will calculate winds, temperatures, and other desired variables. Each day the model was run several times using different combinations of these parameters and later compared with observed meteorological conditions. Parameters that produced the most accurate forecasts were preferentially reused to create the forecasts for the next day. The process was tested for the summer of 2020 and exhibited lower errors than forecasts produced by the model using default values of the parameters.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Simulated CO 2 storage efficiency factors for saline formations of various lithologies and depositional environments using new experimental relative permeability data

Saline formations are attractive geologic reservoirs for permanent carbon dioxide (CO 2 ) storage. Here, the U.S. Department of Energy's National Energy Technology Laboratory (DOE-NETL) has worked to develop and refine methods and tools for the calculation of CO 2 storage potential in subsurface reservoirs. DOE-NETL's CO 2 -SCREEN provides an online tool for executing these storage methods. CO 2 storage efficiency terms are input parameters in DOE-NETL's methods and equations embedded in the CO 2 -SCREEN, which assesses pore space available for CO 2 storage. In this work, a modeling workflow was initiated to refine two CO 2 storage efficiency terms - volumetric displacement (E V ) and microscopic displacement (E d ). The models are based on new experimental relative permeability data that are specific to homogenous lithology and depositional environments of key subsurface saline formations targeted for CO 2 storage. In future work, heterogenous features will be added to this initial modeling effort to update efficiency factors as described in DOE-NETL's methods and CO 2 -SCREEN tool. E V accounts for the volume utilized in the reservoir under the areal plume, while E d accounts for saturation values in the plume to assess efficiency of CO 2 storage at the pore scale. The results of this work are significant in that prior values were based on a limited geologically non-specific relative permeability data set that were collected prior to 2009. Specifically, we applied numerical simulations using TOUGH3 models to update CO 2 storage efficiency values for supercritical CO 2 injection into brine-saturated reservoirs for three lithologies (clastics, limestone, dolomite) and six depositional environments (Marginal Marine, Strand Plain, Deltaic Complex Fluvial, Aeolian, Shallow Marine, and Reef) that have a high potential for geologic CO 2 storage. Experimental relative permeability data in cores from these environments were utilized in the models with corresponding rock type/sedimentary environment. Results of this study showed that dolomite followed by limestone generated higher ranges of storage efficiency compared to clastics. The updated values provided a tighter efficiency range for clastics, lower P 10 but higher P 90 range for limestone, and higher P 10 and P 90 for dolomite. In general, tighter reservoirs with relatively low permeability and porosity were associated with higher E V and E d , showing efficient reservoir and pore utilization in these scenarios. High reservoir pressure and temperature associated with increasing depth increased the E V , and high CO 2 injection rates resulted in increases in E V and E d , while the impact of permeability anisotropy was minimal after the 30-year injection period.

03 NATURAL GAS↗
NASA NTRSDate not supplied

Advancing Open Science in Atmospheric Research: Integrating Data Usability and Machine Learning

In the dynamic realm of atmospheric sciences, the convergence of data science methodologies and open data marks a transformative era, driving research advancements and nurturing aspiring scientists. This abstract highlights two pivotal projects that epitomize open science principles, aligning seamlessly with the session's objective of interdisciplinary synergy and the cultivation of emerging talent. As a NASA-certified data center, our foremost endeavor focuses on enhancing the visibility and traceability of NASA datasets within atmospheric science research. This initiative not only elevates these datasets' prominence but also establishes a robust framework ensuring their credibility in scholarly discourse. By bridging the gap between data sources and research publications, this project serves as an educational catalyst, nurturing a new generation of scholars in open collaboration and dataset authenticity. Concurrently, our second project pioneers an early warning system for flooding events, utilizing machine learning algorithms to predict flooded fractions. Through multi-source data fusion and predictive modeling, this initiative goes beyond forecasting; it embodies the core of open science by enabling proactive risk mitigation strategies. This project not only advances atmospheric sciences but also fosters an environment where young scholars engage in practical, data-driven solutions. These intertwined projects exemplify the fusion of data science with open data solutions, ensuring both the usability of quality datasets and the cultivation of scientific knowledge among emerging scholars. By spotlighting these impactful use cases, our aim is to foster discussions emphasizing the importance of open collaboration, data integrity, and the nurturing of scientific talent in atmospheric sciences." "In the dynamic realm of atmospheric sciences, the convergence of data science methodologies and open data marks a transformative era, driving research advancements and nurturing aspiring scientists. This abstract highlights two pivotal projects that epitomize open science principles, aligning seamlessly with the session's objective of interdisciplinary synergy and the cultivation of emerging talent. As a NASA-certified data center, our foremost endeavor focuses on enhancing the visibility and traceability of NASA datasets within atmospheric science research. This initiative not only elevates these datasets' prominence but also establishes a robust framework ensuring their credibility in scholarly discourse. By bridging the gap between data sources and research publications, this project serves as an educational catalyst, nurturing a new generation of scholars in open collaboration and dataset authenticity. Concurrently, our second project pioneers an early warning system for flooding events, utilizing machine learning algorithms to predict flooded fractions. Through multi-source data fusion and predictive modeling, this initiative goes beyond forecasting; it embodies the core of open science by enabling proactive risk mitigation strategies. This project not only advances atmospheric sciences but also fosters an environment where young scholars engage in practical, data-driven solutions. These intertwined projects exemplify the fusion of data science with open data solutions, ensuring both the usability of quality datasets and the cultivation of scientific knowledge among emerging scholars. By spotlighting these impactful use cases, our aim is to foster discussions emphasizing the importance of open collaboration, data integrity, and the nurturing of scientific talent in atmospheric sciences.

Jennifer Wei↗
DOE OSTI2022

Experiments and modelling of adhesive failure initiation of an epoxy underfill: a study of electronic packaging survivability

Epoxy underfills can be implemented in electronic packaging to enhance solder joint reliability of surface mounted components. However, it is important for an engineer to have a failure criterion that can be used for failure predictions and redesign of electronic assemblies. For this study, data from epoxy bond failure in mock electronic part assemblies were correlated to finite element analyses to predict adhesive failure initiation. Experiments were performed to determine failure loads for various loading locations and nonlinear viscoelastic analyses were performed for the same loading locations to determine a maximum principal strain failure parameter. Predictions showed that a maximum principal strain failure parameter defined from one test could be used as an indicator of adhesive failure of an epoxy bond undergoing other modes of loading. Failure initiation predictions matched experimental data using a maximum principal strain failure parameter for an epoxy bond undergoing mixed modes of loading for both unfilled and alumina oxide filled 828DEA epoxy. Such experimental setup is deemed appropriate for future epoxy testing.

42 ENGINEERING↗
DOE OSTI2019

Collision system and beam energy dependence of anisotropic flow fluctuations

New measurements of elliptic flow from two- and four-particle correlations are used to investigate flow fluctuations in collisions of U+U at √ s NN = 193 GeV, Cu+Au at √ s NN = 200 GeV and Au+Au at several beam energies. These measurements highlight the dependence of these fluctuations on the event shape, system size and beam energy and indicate a dominant role for initial-state-driven fluctuations. These measurements could provide further constraints for initial-state models, as well as for precision extraction of the temperature-dependent specific shear viscosity η/s( T ).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

A Model for the Signal Initiation Complex Between Arrestin-3 and the Src Family Kinase Fgr

Arrestins regulate a wide range of signaling events, most notably when bound to active G protein-coupled receptors (GPCRs). Among the known effectors recruited by GPCR-bound arrestins are Src family kinases, which regulate cellular growth and proliferation. Here, we focus on arrestin-3 interactions with Fgr kinase, a member of the Src family. Previous reports demonstrated that Fgr exhibits high constitutive activity, but can be further activated by both arrestin-dependent and arrestin-independent pathways. We report that arrestin-3 modulates Fgr activity with a hallmark bell-shaped concentration-dependence, consistent with a role as a signaling scaffold. Here, we further demonstrate using NMR spectroscopy that a polyproline motif within arrestin-3 interacts directly with the SH3 domain of Fgr. To provide a framework for this interaction, we determined the crystal structure of the Fgr SH3 domain at 1.9 Å resolution and developed a model for the GPCR-arrestin-3-Fgr complex that is supported by mutagenesis. This model suggests that Fgr interacts with arrestin-3 at multiple sites and is consistent with the locations of disease-associated Fgr mutations. Collectively, these studies provide a structural framework for arrestin-dependent activation of Fgr.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

The CYBER security – Competency Health and Maturity Progression (CYBER-CHAMP) model: Extending the National Initiative for Cybersecurity Education (NICE) Framework Across Organizational Security

Problem Statement: There is a pervasive talent deficit in the cybersecurity industry that prevents employers from being able to fill their open positions efficiently. A holistic approach to security is required to ensure organizations have adequate prevention and response capabilities in case of a cyberattack. Specifically, industrial control systems (ICS’s) and their operational technology (OT) components have become a constant target for cyberattacks. Research Questions: It is proposed that the NICE Framework should be extended in the following areas: 1) Include guidance regarding the job roles and competencies for both IT and OT professionals. 2) Offer step-by-step solutions, based on the work role mappings from the NICE Framework, to increase cybersecurity through employee training and education. 3) Provide a streamlined, lifecycle approach to building a cybersecurity program. Contribution: The CYBER security – Competency Health and Maturity Progression (CYBER-CHAMP©) model provides a customized solution for businesses to understand their education gaps in organizational security and target areas for improvement. Rationale: The Framework for Improving Critical Infrastructure Cybersecurity v1.1 addresses ICS but does not offer a measurement of cybersecurity maturity or clear methods to ascertain an organization’s current risk profile. In Phases 1 and 5 of the model, measurements are provided to help an organization build their current and target risk profiles. The NICE framework provides a structure for planning an IT cybersecurity workforce, but the OT aspects of cybersecurity are only briefly discussed. The model uses Phases 2-3 to examine the competencies of an organization’s workforce, which includes both IT and OT roles. Current frameworks do not offer next steps to increase an organization’s cybersecurity. During Phase 4, employees’ roles are mapped to training, education, and/or certifications from common vendors. Investigative Approach: The model provides measurements and metrics for both an organization’s status and continual improvement. This improvement methodology includes guidance for creating an overall strategic plan for security improvement via products designed to increase an organization’s operational readiness through workforce competency health. Lessons Learned: Depending on who was participating, there were contradicting answers given in Phase 1 due to different security cultures in the organization. This revelation has influenced the steps listed in the User’s Guide, where Phase 1’s first recommended step is to assemble a team that champions the facilitation and implementation of the model in the organization. During Phase 2, the discovery was made that organizations may be missing roles that are necessary to perform critical cybersecurity functions. By understanding the functional roles and competencies needed, they can contract or hire cybersecurity help to fill these gaps. Implications: Using the model, organizations can discuss quantitative measures for improvement as a business case for advancing their security program. Future research can validate and extend the present theory and model to a variety of environments. It is of interest to investigate additional security roles and knowledge domains that are used to build standardized cybersecurity curriculum.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

A bespoke model of Arctic river basins based on hillslope delineation: Model Archive

This dataset is a model archive of the paper A bespoke model of Arctic river basins based on hillslope delineation (in prep), which introduces a watershed decomposition and parameterization method for large scale permafrost hydrology simulation. With this dataset, this study aims to address the research question: whether a computationally efficient hillslope-based modeling framework can reliably simulate discharge at Arctic river-basin scales. This dataset contains model input and output data for five modeling scenarios at a study site located in the Sagavanirktok River basin. The five modeling scenarios include three modeling cases under temperate conditions using full 3D, decomposed 3D, and decomposed 2D modeling strategies; and two modeling cases under actual Arctic conditions with permafrost using full 3D and decomposed 2D modeling strategies. Simulations were performed using the Advanced Terrestrial Simulator (ATS, v1.6 for three temperate scenarios and v1.5 for two Arctic scenarios), a physics-rich integrated surface–subsurface hydrologic model with cryo-hydrology features. For the three temperate models, simulations were conducted for the period of 10/01/1993 - 09/30/2002; and for the two Arctic models, simulations were conducted for the period of 01/01/1994 - 12/31/2002. To facilitate reproducibility of simulations, all datasets are organized hierarchically. The dataset contains: (1) Mesh files (.exo) for full 3D model, decomposed 3D models, and decomposed 2D models, located in huc/190604020802_gauge15906000/mesh/. Mesh files can be visualized through Paraview or read by Python. (2) Climate forcings (.h5) for full 3D model and decomposed 3D/2D models are located in huc/190604020802_gauge15906000/daymet_onePiece/, and huc/190604020802_gauge15906000/vp_pr_revised_daymet_1980_2006_with_wind/ separately. Accessible by Python. (3) Raw measured gage discharge (.csv) from USGS, located in huc/190604020802_gauge15906000/gaged_basin15906000_discharge_usgs/. Accessible by Python. (4) Delineated subdomain raster (.tif) and shape files (.shp), and the final parameterized results (.npy) for decomposed models, located in huc/190604020802_gauge15906000/data_preprocessed-meshing. Accessible by Python. (5) Temperate models are located in nonpermaf_huc190604020802_gauge15906000/, which includes three cases: decomposed 2D models (inside model_0*-hillslope_*), decomposed 3D models (inside model_1*-subcatchment_*), and full 3D model (inside model_2*-onepiece_*). Two step spin-up results (checkpoint_final.h5) are located in model_*1-*_spinup_steadystate and model_*2-*_spinup_cycle, separately, which are used to initialize real transient models. The input files (.xml) and output results (.dat) of the real transient models are located in model_*3-*_transient/. Especially, for two example hillslope models (ID=-11 and 11), additional h5py files are included in model_03-hillslope_transient/hillslope-11/, model_03-hillslope_transient/hillslope11, model_13-subcatchment_transient/subcatchment-11/, model_13-subcatchment_transient/subcatchment/11, respectively, which are used to plot the saturation figure (Figure 5) in the manuscript. Accessible by Python. (6) Arctic models are located in huc190604020802_gauge15906000/, which includes two cases: decomposed 2D models (inside model_04-hillslope_transient), and full 3D model (inside model_05-onepiece_transient_mannp1_ra). Three step spin-up results (checkpoint_final.h5) are located in model_01-column_freezeup/, model_02-column_spinup/, model_03-hillslope_spinup/, respectively, which are used to initialize real 2D transient hillslope models. The input files (.xml) and output results (.dat) of transient 2D hillslope models are located in model_04-hillslope_transient/. The input files (.xml) and output results (.dat) of the full 3D transient model is located in model_05-onepiece_transient_mannp1_ra/. The full 3D transient model is initialized by model_02-column_spinup/. Accessible by Python. (7) The MOSART routed discharge results (.csv) under Arctic conditions is located in huc190604020802_gauge15906000/MOSART/. Accessible by Python. (8) All Python codes (.py) used to parameterize full 3D model to decomposed 2D models are located in script/. These codes fit with watershed workflow (a watershed delineation tool) v1.4 under the branch gaob/v1.4 from https://github.com/gaobhub/watershed-workflow.git.

EARTH SCIENCE > CRYOSPHERE↗
NASA NTRS2001

Analysis of EIT/LASCO Observations Using Available MHD Models: Investigation of CME Initiation Propagation and Geoeffectiveness

The Sun's activity drives the variability of geospace (i.e., near-earth environment). Observations show that the ejection of plasma from the sun, called coronal mass ejections (CMEs), are the major cause of geomagnetic storms. This global-scale solar dynamical feature of coronal mass ejection was discovered almost three decades ago by the use of space-borne coronagraphs (OSO-7, Skylab/ATM and P78-1). Significant progress has been made in understanding the physical nature of the CMEs. Observations show that these global-scale CMEs have size in the order of a solar radius (approximately 6.7 x 10(exp 5) km) near the sun, and each event involves a mass of about 10(exp 15) g and an energy comparable to that of a large flare on the order of 10(exp 32) ergs. The radial propagation speeds of CMEs have a wide range from tens to thousands of kilometers per second. Thus, the transit time to near earth's environment [i.e., 1 AU (astronomical unit)] can be as fast as 40 hours to 100 hours. The typical transit time for geoeffective events is approximately 60-80 h. This paper consists of two parts: 1) A summary of the observed CMEs from Skylab to the present SOHO will be presented. Special attention will be made to SOHO/ LASCO/ EIT observations and their characteristics leading to a geoeffectiv a CME 2) The chronological development of theory and models to interpret the physical nature of this fascinating phenomenon will be reviewed. Finally, an example will be presented to illustrate the geoeffectiveness of the CMEs by using both observation and model.

Wu, S. T.↗
DOE OSTI2020

Evaluation of the Quasi‐Biennial Oscillation in global climate models for the SPARC QBO‐initiative

Abstract Quasi‐biennial oscillations (QBOs) in thirteen atmospheric general circulation models forced with both observed and annually repeating sea surface temperatures (SSTs) are evaluated. In most models the QBO period is close to, but shorter than, the observed period of 28 months. Amplitudes are within ±20 % of the observed QBO amplitude at 10 hPa, but typically about half of that observed at lower altitudes (50 and 70 hPa). For almost all models, the oscillation's amplitude profile shows an overall upward shift compared to reanalysis and its meridional extent is too narrow. Asymmetry in the duration of eastward and westward phases is reasonably well captured, though not all models replicate the observed slowing of the descending westward shear. Westward phases are generally too weak, and most models have an eastward time mean wind bias throughout the depth of the QBO. The intercycle period variability is realistic and in some models is enhanced in the experiment with observed SSTs compared to the experiment with repeated annual cycle SSTs. Mean periods are also sensitive to this difference between SSTs, but only when parametrized non‐orographic gravity wave (NOGW) sources are coupled to tropospheric parameters and not prescribed with a fixed value. Overall, however, modelled QBOs are very similar whether or not the prescribed SSTs vary interannually. A portrait of the overall ensemble performance is provided by a normalized grading of QBO metrics. To simulate a QBO, all but one model used parametrized NOGWs, which provided the majority of the total wave forcing at altitudes above 70 hPa in most models. Hence the representation of NOGWs either explicitly or through parametrization is still a major uncertainty underlying QBO simulation in these present‐day experiments.

Bushell, A. C.↗
NASA NTRS1997

Examination of Soil Moisture Retrieval Using SIR-C Radar Data and a Distributed Hydrological Model

A major objective of soil moisture-related hydrological-research during NASA's SIR-C/X-SAR mission was to determine and compare soil moisture patterns within humid watersheds using SAR data, ground-based measurements, and hydrologic modeling. Currently available soil moisture-inversion methods using active microwave data are only accurate when applied to bare and slightly vegetated surfaces. Moreover, as the surface dries down, the number of pixels that can provide estimated soil moisture by these radar inversion methods decreases, leading to less accuracy and, confidence in the retrieved soil moisture fields at the watershed scale. The impact of these errors in microwave- derived soil moisture on hydrological modeling of vegetated watersheds has yet to be addressed. In this study a coupled water and energy balance model operating within a topographic framework is used to predict surface soil moisture for both bare and vegetated areas. In the first model run, the hydrological model is initialized using a standard baseflow approach, while in the second model run, soil moisture values derived from SIR-C radar data are used for initialization. The results, which compare favorably with ground measurements, demonstrate the utility of combining radar-derived surface soil moisture information with basin-scale hydrological modeling.

Hsu, A. Y.↗
DOE OSTI2021

An initial exploration of Bayesian model calibration for estimating the composition of rocks and soils on Mars

The Mars Curiosity rover carries an instrument, ChemCam, designed to measure the composition of surface rocks and soil using laser-induced breakdown spectroscopy (LIBS). The measured spectra from this instrument must be analyzed to identify the component elements in the target sample, as well as their relative proportions. This process, which we call disaggregation, is complicated by so-called matrix effects, which describe nonlinear changes in the relative heights of emission lines as an unknown function of composition due to atomic interactions within the LIBS plasma. In this work, we explore the use of the plasma physics code ATOMIC, developed at Los Alamos National Laboratory, for the disaggregation task. ATOMIC has recently been used to model LIBS spectra and can robustly reproduce matrix effects from first principles. The ability of ATOMIC to predict LIBS spectra presents an exciting opportunity to perform disaggregation in a manner not yet tried in the LIBS community, namely via Bayesian model calibration. However, using it directly to solve our inverse problem is computationally intractable due to the large parameter space and the computation time required to produce a single output. Therefore, we also explore the use of emulators as a fast solution for this analysis. We discuss a proof of concept Gaussian process emulator for disaggregating two-element compounds of sodium and copper. The training and test datasets were simulated with ATOMIC using a Latin hypercube design. After testing the performance of the emulator, we successfully recover the composition of 25 test spectra with Bayesian model calibration.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Partitioning uncertainty in projections of Arctic sea ice

Abstract Improved knowledge of the contributing sources of uncertainty in projections of Arctic sea ice over the 21st century is essential for evaluating impacts of a changing Arctic environment. Here, we consider the role of internal variability, model structure and emissions scenario in projections of Arctic sea-ice area (SIA) by using six single model initial-condition large ensembles and a suite of models participating in Phase 5 of the Coupled Model Intercomparison Project. For projections of September Arctic SIA change, internal variability accounts for as much as 40%–60% of the total uncertainty in the next decade, while emissions scenario dominates uncertainty toward the end of the century. Model structure accounts for 60%–70% of the total uncertainty by mid-century and declines to 30% at the end of the 21st century in the summer months. For projections of wintertime Arctic SIA change, internal variability contributes as much as 50%–60% of the total uncertainty in the next decade and impacts total uncertainty at longer lead times when compared to the summertime. In winter, there exists a considerable scenario dependence of model uncertainty with relatively larger model uncertainty under strong forcing compared to weak forcing. At regional scales, the contribution of internal variability can vary widely and strongly depends on the calendar month and region. For wintertime SIA change in the Greenland-Iceland-Norwegian and Barents Seas, internal variability contributes 60%–70% to the total uncertainty over the coming decades and remains important much longer than in other regions. We further find that the relative contribution of internal variability to total uncertainty is state-dependent and increases as sea ice volume declines. These results demonstrate that internal variability is a significant source of uncertainty in projections of Arctic sea ice.

Bonan, David B. (ORCID:0000000338676009)↗
DOE OSTI2025

Climate models show colorado drying sooner and with greater certainty east of the continental divide

Many studies have examined the aridity of the Colorado River Basin and the possible impacts of climate change which could further strain already over-allocated water resources in the region. Fewer studies have examined the multiple Colorado Rocky Mountain headwater regions specifically. This is especially true of areas East of the Continental Divide, despite water originating there being critical to cities and agriculture in Eastern Colorado and further downstream. This paper explores and compares drying trends in the Eastern and Western Colorado Rocky Mountains using single-model initial-condition large ensembles from ten global climate models. The use of multiple models allows us to identify signals that are consistent across different physics parameterizations, model grids, and other model intricacies. The large ensembles also allow us to quantify the time of emergence of these climate change signals--that is, when did (or when will) the long term change due to anthropogenic greenhouse gasses exceed the internal variability of the climate system. Consistent with previous studies, we find evidence of drying on both sides of the Continental Divide. That drying is more pronounced, occurs sooner, and is more consistent across global climate models in the East, however, highlighting the region’s importance despite generally receiving less attention than the West.

Rugg, Allyson↗