Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “joint state estimation”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS2019

Advanced Modular Power System Electronics Enclosure Thermal Testing

An analysis was set up to model the temperature of the advanced modular power system (AMPS) power distribution cards when installed within the electronics enclosure case. The analysis was used to determine the steady-state temperature distribution of the cards within the case. To verify the analysis, an experiment was set up and conducted to simulate the operation of the cards within the enclosure. Four tests were conducted. The tests varied the position of the cold plate and evaluated the use of a thermal compound to reduce the contact resistance between the joints within the thermal path between the cards and the cold plate. Three of the four cases examined showed very good agreement between the analysis and the experiment with a less than 1-percent variation in the predicated temperatures determined through the analysis and the experimentally derived temperatures. In the remaining case, the difference between the analysis and experiment was approximately 12 percent. Both the experiment and analysis showed that the modular power conditioning cards can be maintained within their desired maximum operating temperature range of 40 to 45 °C through thermal conduction to a cold plate when operating with their estimated maximum heat output of 16 W per card.

Colozza, Anthony J.↗
DOE OSTI2025

Inverse rendering of fusion plasmas: Inferring plasma composition from imaging systems

In this work, we develop a differentiable rendering pipeline for visualising plasma emission within tokamaks, and estimating the gradients of the emission and estimating other physical quantities. Unlike prior work, we are able to leverage arbitrary representations of plasma quantities and easily incorporate them into a non-linear optimisation framework. The efficiency of our method enables not only estimation of a physically plausible image of plasma, but also recovery of the neutral Deuterium distribution from imaging and midplane measurements alone. We demonstrate our method with three different levels of complexity showing first that a poloidal neutrals density distribution can be recovered from imaging alone, second that the distributions of neutral Deuterium, electron density and electron temperature can be recovered jointly, and finally, that this can be done in the presence of realistic imaging systems that incorporate sensor cropping and quantisation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2024

XLSSC 122 Caught in the Act of Growing Up: Spatially Resolved SZ Observations of a z=1.98 Galaxy Cluster

Context. How protoclusters evolved from sparse galaxy overdensities to mature galaxy clusters is still not well understood. In this context, detecting and characterizing the hot intracluster medium (ICM) at high redshifts (z ∼ 2) is key to understanding how the continuous accretion from the filamentary large-scale structure and the mergers along it impact the first phases of cluster formation. Aims. We study the dynamical state and morphology of the z = 1:98 galaxy cluster XLSSC 122 with high-resolution observations (≈5″) of the ICM through the Sunyaev–Zeldovich (SZ) effect. XLSSC 122 is the highest redshift optically confirmed galaxy cluster found in an unbiased, widefield survey. Methods. Via Bayesian forward modeling, we mapped the ICM on scales from the virial radius down to the core of the cluster. To constrain such a broad range of spatial scales, we employed a new technique that jointly forward-models parametric descriptions of the pressure distribution to interferometric ACA and ALMA observations and multiband imaging data from ACT. Results. We detect the SZ effect with 11σ significance in the ALMA+ACA observations and find a flattened inner pressure profile that is consistent with a noncool core classification with a significance of ≥3σ. In contrast to the previous works, we find better agreement between the SZ effect signal and the X-ray emission as well as the cluster member distribution. Further, XLSSC 122 exhibits an excess of SZ flux in the south of the cluster where no X-ray emission is detected. By reconstructing the interferometric observations and modeling in the uv-plane, we obtain a tentative detection of an infalling group or filamentary-like structure in the southeast that is believed to boost and heat up the ICM while the density of the gas is still low. In addition to characterizing the dynamical state of the cluster, we provide an improved SZ mass estimate M500,c = 1.66−0.20+0.23 × 1014 M⊙ . Conclusions. Altogether, the observations indicate that we see XLSSC 122 in a dynamic phase of cluster formation while a large reservoir of gas is already thermalized.

large-scale structure of Universe↗
DOE OSTI2025

Polynomial Chaos Surrogate Construction for Random Fields with Parametric Uncertainty

Engineering and applied science rely on computational experiments to rigorously study physical systems. The mathematical models used to probe these systems are highly complex, and sampling-intensive studies often require prohibitively many simulations for acceptable accuracy. Surrogate models provide a means of circumventing the high computational expense of sampling such complex models. In particular, polynomial chaos expansions (PCEs) have been successfully used for uncertainty quantification studies of deterministic models where the dominant source of uncertainty is parametric. We discuss an extension to conventional PCE surrogate modeling to enable surrogate construction for stochastic computational models that have intrinsic noise in addition to parametric uncertainty. We develop a PCE surrogate on a joint space of intrinsic and parametric uncertainty, enabled by Rosenblatt transformations, which are evaluated via kernel density estimation of the associated conditional cumulative distributions. Furthermore, we extend the construction to random field data via the Karhunen–Loève expansion. We then take advantage of closed-form solutions for computing PCE Sobol indices to perform a global sensitivity analysis of the model which quantifies the intrinsic noise contribution to the overall model output variance. Additionally, the resulting joint PCE is generative in the sense that it allows generating random realizations at any input parameter setting that are statistically approximately equivalent to realizations from the underlying stochastic model. The method is demonstrated on a chemical catalysis example model and a synthetic example controlled by a parameter that enables a switch from unimodal to bimodal response distributions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Tree-level carbon stock estimations across diverse species using multi-source remote sensing integration

Forests are critical carbon sinks, and remote sensing has been increasingly widely used for forest monitoring and biomass estimations. However, species-specific tree-level studies remain limited. In this study, we demonstrated the feasibility of integrating UAV-based LiDAR with high-resolution optical satellite imagery (0.5 m) to estimate biomass for individual trees across different species. The proposed method accurately estimated biomass for 53 trees (R² = 0.82, rRMSE = 0.44), with species-specific datasets, showing an average 25.2% increase in R² and a 14.8% reduction in rRMSE. A novel vegetation index combining forest structure parameters with vegetation indices (VIs) was developed using high-resolution multispectral satellite data (3 m) to explore its relationship with individual tree biomass. Combining forest structural parameters with VIs further improved estimation accuracy, achieving an R²of 0.89 and an rRMSE of 0.34. Species-specific datasets show an 11.6% increase in R²compared to methods without VIs, and a 22.2% improvement over methods using only VIs. SHapley Additive exPlanations (SHAP) analysis shows that the volume feature played a key role in model performance and remained stable throughout the training process. Altogether, the proposed approach enhances individual tree biomass and carbon sink estimations, showing great potential for large-scale precise forest carbon monitoring using multi-source remote sensing data.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2019

Three CA-Related Covariance Issues and Their Solutions

This presentation addresses three issues that arise in the use of DoD-produced satellite state estimate covariances in the conjunction assessment process: the realism of the provided covariances, how to address correlated error between two satellites' covariance matrices, and how to proceed when a furnished covariance is non-positive-definite. To address the first, DoD has implemented a set of two "consider parameters" with which to alter/expand the covariance to account for atmospheric density forecast error and satellite frontal area uncertainty, the two largest sources of unmodeled position prediction error for LEO orbits; these values are governed by a satellite's orbital parameters and ballistic coefficient histories, as well as the current and predicted space weather situation. This approach substantially improves the realism of the covariance by accounting for expected prediction errors that are not part of the fit process that generates the covariance. The second issue of covariance correlation is a complicated one, but the main shared error source—global atmospheric density error that is common to both satellites can be characterized and, through sensitivity vectors, quantified at the conjunction's time of closest approach and removed from the joint covariance, which is used in the CA probability of collision (Pc) calculation. Finally, while according to the orbit determination theory a non-positive-definite (NPD) covariance is not possible, numerical truncation and covariance interpolation can conspire to produce NPD results, which render the covariance unusable for certain CA risk assessment calculations, such as Monte Carlo Pc determination. Three similar techniques are profiled and the simplest of them recommended as a reasonable remediation technique when NPD covariances are encountered by CA practitioners.

Dynamic Consider Parameter↗
NASA NTRS2020

Assimilation of Vegetation Optical Depth Retrievals from Passive Microwave Radiometry

Vegetation optical depth (VOD) retrievals from passive microwave sensors provide analog estimates of above-ground canopy biomass. This study presents the development and analysis of assimilating VOD retrievals from X-, C-, and L-band passive microwave instruments within the Noah-MP land surface model over the Continental U.S. The results from this study demonstrate that the assimilation of VOD retrievals have a significant beneficial impact on the simulation of evapotranspiration and GPP, particularly over the agricultural areas of the U.S. The improvements in the water and carbon fluxes from the assimilation of VOD from X- and C-band sensors are found to be comparable to those obtained from the assimilation of vegetation indices from optical sensors.The study also quantifies the relative and joint impact of assimilating surface soil moisture and VOD from the Soil Moisture Active Passive (SMAP) mission. The utility of soil moisture assimilation for improving ET is more significant over water-limited regions, whereas VOD DA is more impactful over areas where soil moisture is not the primary controlling factor on ET. The results also indicate that the information on moisture and vegetation states from SMAP can be simultaneously exploited through the joint assimilation of surface soil moisture and VOD. Since passive microwave-based VOD retrievals are available in nearly all weather conditions, their use within data assimilation systems offers the ability to extend and improve the utility obtained from the use of optical/infrared based vegetation retrievals

Sujay V Kumar↗
DOE OSTI2025

Constraints on Local Primordial Non-Gaussianity with 3D Velocity Reconstruction from the Kinetic Sunyaev-Zeldovich Effect

The cosmic velocity field is an unbiased probe of the total matter distribution but is challenging to measure directly at intermediate and high redshifts. The large-scale velocity field imprints a signal in the cosmic microwave background (CMB) through the kinetic Sunyaev-Zeldovich (kSZ) effect. We perform the first 3D reconstruction of the large-scale velocity field from the kSZ effect by applying a quadratic estimator to CMB temperature maps and the 3D positions of galaxies. We do so by combining CMB data from the fifth data release of the Atacama Cosmology Telescope (in combination with Planck) and a spectroscopic galaxy sample from the Sloan Digital Sky Survey. We then measure the galaxy-velocity cross-power spectrum and detect the presence of the kSZ signal at a signal-to-noise ratio of 7.2⁢𝜎. Using this galaxy-velocity cross-correlation alone, we constrain the amplitude of local primordial non-Gaussianity finding 𝑓 NL =−9⁢0$^{+210}_{−350}$. In conclusion, this pathfinder measurement sets the stage for joint galaxy-CMB kSZ constraints to significantly enhance the 𝑓 NL information obtained from galaxy surveys through sample variance cancellation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Code and Solution Verification Assessment of the CTF Thermal Hydraulic Subchannel Code

CTF is a thermal-hydraulics subchannel code jointly developed by Oak Ridge National Laboratory and North Carolina State University. Over the past seven years, the Consortium for Advanced Simulation of Light Water Reactors (CASL) has made a significant investment in developing CTF so it can be used to model light water reactors, including nominal operating conditions, departure from nucleate boiling analysis, and transients ranging from loss of flow to reactivity insertion accidents. In addition to implementing new modeling capabilities and developing the user input and output interface, extensive work has been performed to improve the code’s quality assurance program, resulting in a development process that conforms with NQA-1 requirements. The CASL program follows the Predictive Capability Maturity Model (PCMM) approach for assessing code quality, which emphasizes performing code verification(ensuring the code converges to the correct answer) and solution verification (ensuring the code converges for the intended application). Code and solution verification are used to identify uncertainty errors introduced by numerical approximations in the code and are important for demonstrating that the model is coded without error, which is an important aspect of the Best Estimate plus Uncertainty method. This paper presents a comprehensive overview of the code and solution verification testing that has been performed on CTF. A top-down approach is taken in which the intended CTF applications are presented, followed by the code features required for their modeling. These features are then linked to the applicable code and solution verification tests that demonstrate proper functioning. Past testing efforts are summarized, and new tests are added to help close gaps in the presented test matrix. Rather than performing “one-off” exercises, these tests are added to the automated CTF regression test suite to ensure continual code quality.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRS2008

Close Approach Prediction Analysis of the Earth Science Constellation with the Fengyun-1C Debris

Routine satellite operations for the Earth Science Constellation (ESC) include collision risk assessment between members of the constellation and other orbiting space objects. Each day, close approach predictions are generated by a U.S. Department of Defense Joint Space Operations Center Orbital Safety Analyst using the high accuracy Space Object Catalog maintained by the Air Force's 1" Space Control Squadron. Prediction results and other ancillary data such as state vector information are sent to NASAJGoddard Space Flight Center's (GSFC's) Collision Risk Assessment analysis team for review. Collision analysis is performed and the GSFC team works with the ESC member missions to develop risk reduction strategies as necessary. This paper presents various close approach statistics for the ESC. The ESC missions have been affected by debris from the recent anti-satellite test which destroyed the Chinese Fengyun- 1 C satellite. The paper also presents the percentage of close approach events induced by the Fengyun-1C debris, and presents analysis results which predict the future effects on the ESC caused by this event. Specifically, the Fengyun-1C debris is propagated for twenty years using high-performance computing technology and close approach predictions are generated for the ESC. The percent increase in the total number of conjunction events is considered to be an estimate of the collision risk due to the Fengyun-1C break- UP.

Duncan, Matthew↗
NASA NTRS2011

Global Assimilation of Multi-Sensor Snow Observations for Improved Characterization of Snow Processes

Snow conditions on the land surface are recognized to be key components of the global hydrological cycle as they play a critical role in the determination of local and regional climate. In many mid-latitude and high-latitude regions, the seasonal water storage and associated spring snowmelt dominate the local hydrology. The contribution to the runoff and moisture conditions from snow is vital in supporting agriculture and in determining water resources management practices. Consequently, accurate characterization of snow properties becomes important for both end-use applications and weather and climate research. Recently a joint effort between the u.S. Air Force and NASA has enabled a blended, multi-sensor snow product known as the AFWA NASA Snow Algorithm (ANSA). This global snow dataset has been generated by utilizing the Earth Observation System (EOS) Moderate Resolution Imaging Spectroradiometer (MODIS) and Advanced Microwave Scanning Radiometer for EOS (AMSR-E) datasets. ANSA product includes estimates of snow cover extent, snow water equivalent (SWE) and SWE-derived snow depth fields. The MODIS-based products enable snow cover mappings under cloud-free conditions whereas the passive microwave data from AMSR-E provides measurements under cloudy conditions. These remotely-sensed snow observations are further augmented with the information from ground-based snow measurements through data fusion techniques. The resulting ANSA products are employed in the NASA Land Information System (LIS) data assimilation framework, which provides a comprehensive environment for integrating community land surface models, ground and satellite-based observations, and ensemble-based data assimilation tools. LIS incorporates the multisensor ANSA snow retrievals with the land surface model estimates to generate spatially and temporally continuous estimates of snow states, through data assimilation. A suite of experiments to assimilate ANSA snow cover, SWE and snow depth estimates with different land surface models in LIS are conducted and the resulting estimates of snow conditions are evaluated against a number of in-situ observational datasets, over several regions of the world. These evaluations are used to compare and contrast the advantages and disadvantages of these multi-sensor snow observations.

Kumar, Sujay↗
NASA NTRS1997

Anthropogenic Sulfate, Clouds, and Climate Forcing

This research work is a joint effort between research groups at the Battelle Pacific Northwest Laboratory, Virginia Tech University, Georgia Institute of Technology, Brookhaven National Laboratory, and Texas A&M University. It has been jointly sponsored by the National Aeronautics and Space Administration, the U.S. Department of Energy, and the U.S. Environmental Protection Agency. In this research, a detailed tropospheric aerosol-chemistry model that predicts oxidant concentrations as well as concentrations of sulfur dioxide and sulfate aerosols has been coupled to a general circulation model that distinguishes between cloud water mass and cloud droplet number. The coupled model system has been first validated and then used to estimate the radiative impact of anthropogenic sulfur emissions. Both the direct radiative impact of the aerosols and their indirect impact through their influence on cloud droplet number are represented by distinguishing between sulfuric acid vapor and fresh and aged sulfate aerosols, and by parameterizing cloud droplet nucleation in terms of vertical velocity and the number concentration of aged sulfur aerosols. Natural sulfate aerosols, dust, and carbonaceous and nitrate aerosols and their influence on the radiative impact of anthropogenic sulfate aerosols, through competition as cloud condensation nuclei, will also be simulated. Parallel simulations with and without anthropogenic sulfur emissions are performed for a global domain. The objectives of the research are: To couple a state-of-the-art tropospheric aerosol-chemistry model with a global climate model. To use field and satellite measurements to evaluate the treatment of tropospheric chemistry and aerosol physics in the coupled model. To use the coupled model to simulate the radiative (and ultimately climatic) impacts of anthropogenic sulfur emissions.

Ghan, Steven J.↗
DOE OSTI2025

Source Characterization at Low Yields from Physics Experiment One A

The detonation of a well-recorded coupled explosion at the well-characterized Nevada National Security Site, as part of Physics Experiment One, has provided us with an opportunity to compare and contrast methods and results of yield estimation and other event characterization methods. Yield estimates are made using several methods, including local/regional magnitudes, waveform envelopes, acoustic amplitudes, joint seismoacoustic yield, and moment tensors. We include in our analysis a discussion and comparison of applicability, calibrations, assumptions, and uncertainties involved in each. With the exception of magnitude-based methods, all the techniques were more or less successful in recovering the known yield. It is, however, useful to be cognizant of requirements and calibrations, such as station coverage, emplacement conditions, propagation characterization, and uncertainties, to apply the methods to other areas with less complete information. We demonstrate that we can be successful at monitoring low-yield explosions using a variety of methods assuming that the signal-to-noise ratio is adequate at recording stations, and that the material properties are taken into account. Depth estimates are also important, particularly if one allows the possibility of a surface explosion.

Pasyanos, Michael E. [Lawrence Livermore National ↗
NASA NTRS2012

TRMM Latent Heating Retrieval and Comparisons with Field Campaigns and Large-Scale Analyses

Rainfall production is a fundamental process within the Earth's hydrological cycle because it represents both a principal forcing term in surface water budgets, and its energetics corollary, latent heating (LH), is one of the principal sources of atmospheric diabatic heating. Latent heat release itself is a consequence of phase changes between the vapor, liquid, and frozen states of water. The vertical distribution of LH has a strong influence on the atmosphere, controlling large-scale tropical circulations, exciting and modulating tropical waves, maintaining the intensities of tropical cyclones, and even providing the energetics of midlatitude cyclones and other mobile midlatitude weather systems. Moreover, the processes associated with LH result in significant non-linear changes in atmospheric radiation through the creation, dissipation and modulation of clouds and precipitation. Yanai et al. (1973) utilized the meteorological data collected from a sounding network to present a pioneering work on thermodynamic budgets, which are referred to as the apparent heat source (Q1) and apparent moisture sink (Q2). Yanai's paper motivated the development of satellite-based LH algorithms and provided a theoretical background for imposing large-scale advective forcing into cloud-resolving models (CRMs). These CRM-simulated LH and Q1 data have been used to generate the look-up tables used in LH algorithms. This paper examines the retrieval, validation, and application of LH estimates based on rain rate quantities acquired from the Tropical Rainfall Measuring Mission satellite (TRMM). TRMM was launched in November 1997 as a joint enterprise between the American and Japanese space agencies -- with overriding goals of providing accurate four-dimensional estimates of rainfall and LH over the global Tropics and subtropics equatorward of 35o. Other literature has acknowledged the achievement of the first goal of obtaining an accurate rainfall climatology. This paper describes the second major goal of obtaining credible LH estimates as well as their applications within TRMM's zone of coverage, the standard TRMM LH products, and areas for further improvement.

latent heating↗
DOE OSTI2021

Disaggregating Future Retail Electricity Rate Growth [Slides]

Recent Berkeley Lab research found that modest retail rate increases over the past 10 years were mostly driven by large increases in capital expenditures (CapEx) that were offset in part by substantial wholesale price reductions. Decision-makers are increasingly concerned about the potential future rate impacts of a number of policies and industry trends that support rapid decarbonization, electrification, and grid modernization. Using historical FERC Form 1 data and the existing literature on policies and industry trends that are likely to affect utility-incurred costs and retail sales, Berkeley Lab researchers developed ranges of forecasted growth rates for cost-related rate drivers (i.e., fuel and purchased power; transmission, distribution, generation, and other categories of both non-fuel operations & maintenance and CapEx) and non-cost related rate drivers (i.e., retail sales, peak demand, and customers). These were then used as inputs to a pro-forma utility financial model (FINDER) that estimated the growth in retail electric rates between 2020 and 2030 for a prototypical vertically-integrated investor-owned utility in the United States. The analysis produced the following results: 1. Assuming average growth rates in all rate drivers, future retail rate growth is driven by sizable increases in all CapEx costs, where fuel and purchased power costs are replaced by generation CapEx as the largest rate component between 2020 and 2030. 2. Growth in sales/peak demand/customers, generation CapEx costs, and fuel and purchased power (FPP) costs, in isolation, produce the most uncertainty in rate growth. Specifically, a 1% increase in the compound annual growth rate (CAGR) of retail sales, coincident peak demand (CP), and customers (Sales-CP-Cust) results in a 0.88-0.93% decrease in the CAGR of rates, in isolation. However, a 1% increase in the CAGR of generation CapEx budgets results in a 0.07-0.14% increase in the CAGR of rates, while a 1% increase in the CAGR of FPP costs causes a 0.10-0.14% increase in the CAGR of rates, all else being equal. 3. Taking into account the correlation and variability of the growth in all rate drivers jointly, generation CapEx is expected to be both the largest and most uncertain rate component by 2030 (20-25% share of the retail rate). Transmission and distribution CapEx, along with fuel and purchased power costs are each expected to comprise between 12% and 17% of retail rates.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
NASA NTRS2012

Validation of Atmosphere/Ionosphere Signals Associated with Major Earthquakes by Multi-Instrument Space-Borne and Ground Observations

The latest catastrophic earthquake in Japan (March 2011) has renewed interest in the important question of the existence of pre-earthquake anomalous signals related to strong earthquakes. Recent studies have shown that there were precursory atmospheric/ionospheric signals observed in space associated with major earthquakes. The critical question, still widely debated in the scientific community, is whether such ionospheric/atmospheric signals systematically precede large earthquakes. To address this problem we have started to investigate anomalous ionospheric / atmospheric signals occurring prior to large earthquakes. We are studying the Earth's atmospheric electromagnetic environment by developing a multisensor model for monitoring the signals related to active tectonic faulting and earthquake processes. The integrated satellite and terrestrial framework (ISTF) is our method for validation and is based on a joint analysis of several physical and environmental parameters (thermal infrared radiation, electron concentration in the ionosphere, lineament analysis, radon/ion activities, air temperature and seismicity) that were found to be associated with earthquakes. A physical link between these parameters and earthquake processes has been provided by the recent version of Lithosphere-Atmosphere-Ionosphere Coupling (LAIC) model. Our experimental measurements have supported the new theoretical estimates of LAIC hypothesis for an increase in the surface latent heat flux, integrated variability of outgoing long wave radiation (OLR) and anomalous variations of the total electron content (TEC) registered over the epicenters. Some of the major earthquakes are accompanied by an intensification of gas migration to the surface, thermodynamic and hydrodynamic processes of transformation of latent heat into thermal energy and with vertical transport of charged aerosols in the lower atmosphere. These processes lead to the generation of external electric currents in specific regions of the atmosphere and the modifications, by dc electric fields, in the ionosphere-atmosphere electric circuit. We retrospectively analyzed temporal and spatial variations of four different physical parameters (gas/radon counting rate, lineaments change, long-wave radiation transitions and ionospheric electron density/plasma variations) characterizing the state of the lithosphere/atmosphere coupling several days before the onset of the earthquakes. Validation processes consist in two phases: A. Case studies for seven recent major earthquakes: Japan (M9.0, 2011), China (M7.9, 2008), Italy (M6.3, 2009), Samoa (M7, 2009), Haiti (M7.0, 2010) and, Chile (M8.8, 2010) and B. A continuous retrospective analysis was preformed over two different regions with high seismicity- Taiwan and Japan for 2003-2009. Satellite, ground surface, and troposphere data were obtained from Terra/ASTER, Aqua/AIRS, POES and ionospheric variations from DEMETER and COSMIC-I data. Radon and GPS/TEC were obtaining from monitoring sites in Taiwan, Japan and Italy and from global ionosphere maps (GIM) respectively. Our analysis of ground and satellite data during the occurrence of 7 global earthquakes has shown the presence of anomalies in the atmosphere. Our results for Tohoku M9.0 earthquake show that on March 7th, 2011 (4 days before the main shock and 1 day before the M7.2 foreshock of March 8, 2011) a rapid increase of emitted infrared radiation was observed by the satellite data and an anomaly was developed near the epicenter. The GPS/TEC data indicate an increase and variation in electron density reaching a maximum value on March 8. From March 3 to 11 a large increase in electron concentration was recorded at all four Japanese ground-based ionosondes, which returned to normal after the main earthquake. Similar approach for analyzing atmospheric and ionospheric parameters has been applied for China (M7.9, 2008), Italy (M6.3, 2009), Samoa (M7, 2009), Haiti (M7.0, 2010) and Chile (M8.8, 2010) eahquakes. Results have revealed the presence of related variations of these parameters implying their connection with the earthquake process. The second phase (B) of this validation included 102 major earthquakes (M>5.9) in Taiwan and Japan. We have found anomalous behavior before all of these events with no false negatives. False alarm ratio for false positives is less then 10% and has been calculated for the same month of the earthquake occurrence for the entire period of analysis (2003-2009). The commonalities for detecting atmospheric/ionospheric anomalies are: i.) Regularly appearance over regions of maximum stress (i.e., along plate boundaries); ii.) Anomaly existence over land and sea; and iii) association with M>5.9 earthquakes not deeper than 100km. Due to their long duration over the same region these anomalies are not consistent with a meteorological origin. Our initial results from the ISTF validation of multi-instrument space-borne and ground observations show a systematic appearance of atmospheric anomalies near the epicentral area, one to seven (average) days prior to the largest earthquakes, and suggest that it could be explained by a coupling process between the observed physical parameters and the pre-earthquake preparation processes.

Ouzounov, Dimitar↗
DOE OSTI2021

Validation of the TOPAS Monte Carlo toolkit for HDR brachytherapy simulations

The goal of this work is to validate the user-friendly Geant4-based Monte Carlo toolkit TOol for PArticle Simulation (TOPAS) for brachytherapy applications. Brachytherapy simulations performed with TOPAS were systematically compared with published TG-186 reference data. The photon emission energy spectrum, the air-kerma strength, and the dose-rate constant of the model-based dose calculation algorithm (MBDCA)-WG generic Ir-192 source were extracted. For dose calculations, a track-length estimator was implemented. Additionally, the four Joint AAPM/ESTRO/ABG MBDCA-WG test cases were evaluated through histograms of the local and global dose difference volumes. A prostate, a palliative lung, and a breast case were simulated. For each case, the dose ratio map, the histogram of the global dose difference volume, and cumulative dose-volume histograms were calculated. The air-kerma strength was (9.772 ± 0.001) × 10 -8 U Bq -1 (within 0.3% of the reference value). The dose-rate constant was 1.1107 ± 0.0005 cGy h -1 U -1 (within 0.01% of the reference value). For all cases, at least 96.9% of voxels had a local dose difference within [-1%, 1%] and at least 99.9% of voxels had a global dose difference within [-0.1%, 0.1%]. The implemented track-length estimator scorer was more efficient than the default analog dose scorer by a factor of 237. For all clinical cases, at least 97.5% of voxels had a global dose difference within [-1%, 1%]. Dose-volume histograms were consistent with the reference data. TOPAS was validated for high-dose-rate brachytherapy simulations following the TG-186 recommended approach for MBDCAs. Built on top of Geant4, TOPAS provides broad access to a state-of-the-art Monte Carlo code for brachytherapy simulations.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Efficient Active Learning for Gaussian Process Classification by Error Reduction

Active learning sequentially selects the best instance for labeling by optimizing an acquisition function to enhance data/label efficiency. The selection can be either from a discrete instance set (pool-based scenario) or a continuous instance space (query synthesis scenario). In this work, we study both active learning scenarios for Gaussian Process Classification (GPC). The existing active learning strategies that maximize the Estimated Error Reduction (EER) aim at reducing the classification error after training with the new acquired instance in a onestep-look-ahead manner. The computation of EER-based acquisition functions is typically prohibitive as it requires retraining the GPC with every new query. Moreover, as the EER is not smooth, it can not be combined with gradient-based optimization techniques to efficiently explore the continuous instance space for query synthesis. To overcome these critical limitations, we develop computationally efficient algorithms for EER-based active learning with GPC. Further, we derive the joint predictive distribution of label pairs as a one-dimensional integral, as a result of which the computation of the acquisition function avoids retraining the GPC for each query, remarkably reducing the computational overhead. We also derive the gradient chain rule to efficiently calculate the gradient of the acquisition function, which leads to the first query synthesis active learning algorithm implementing EER-based strategies. Our experiments clearly demonstrate the computational efficiency of the proposed algorithms. We also benchmark our algorithms on both synthetic and real-world datasets, which show superior performance in terms of sampling efficiency compared to the existing state-of-the-art algorithms.

97 MATHEMATICS AND COMPUTING↗