Dynamic measurement of impurity ion transport in runaway electron plateaus in DIII-D
Explore the source record for details and available documents.
Hollmann, E. M.↗
SEARCH · Search NASA
Results for “dynamic parameter 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.
Stomatal response of three Arctic plant species, Seward Peninsula, Alaska, 2022
Stomatal response to changing irradiance measurements on three common tussock tundra shrub species. Data were collected in 2022 from two sites on the Seward Peninsula, Kougarok Mile 64 and Teller Mile 27. Measurements were made with LI-COR LI-6800 gas exchange systems using a slow curve light response method. Light curve observations enable the estimation of stomatal slope and intercept parameters, g1 and g0. The data package files include a data file, metadata files and the complete instrument output for all measurements in .csv format, and a pdf describing the experimental protocol. See the related data package NGA315 for foliar trait data (leaf mass per area, leaf nitrogen concentration) and NGA209 for co-located thaw depth and soil moisture. The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research.The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska.Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).
54 ENVIRONMENTAL SCIENCES↗
Ground-Support Algorithms for Simulation, Processing, and Calibration of Barnes Static Earth Sensor Measurements: Applications to Tropical Rainfall Measuring Mission Observatory
New algorithms are described covering the simulation, processing, and calibration of penetration angles of the Barnes static Earth sensor assembly (SESA) as implemented in the Goddard Space Flight Center Flight Dynamics Division ground support system for the Tropical Rainfall Measuring Mission (TRMM) Observatory. The new treatment involves a detailed analysis of the measurements by individual quadrants. It is shown that, to a good approximation, individual quadrant misalignments can be treated simply as penetration angle biases. Simple formulas suitable for real-time applications are introduced for computing quadrant-dependent effects. The simulator generates penetration angles by solving a quadratic equation with coefficients uniquely determined by the spacecraft's position and the quadrant's orientation in GeoCentric Inertial (GCI) coordinates. Measurement processing for attitude determination is based on linearized equations obtained by expanding the coefficients of the aforementioned quadratic equation as a Taylor series in both the Earth oblateness coefficient (alpha approx. 1/150) and the angle between the pointing axis and the geodetic nadir vector. A simple formula relating a measured value of the penetration angle to the deviation of the Earth-pointed axis from the geodetic nadir vector is derived. It is shown that even near the very edge of the quadrant's Field Of View (FOV), attitude errors resulting from quadratic effects are a few hundredths of a degree, which is small compared to the attitude determination accuracy requirement (0.18 degree, 3 sigma) of TRMM. Calibration of SESA measurements is complicated by a first-order filtering used in the TRMM onboard algorithm to compute penetration angles from raw voltages. A simple calibration scheme is introduced where these complications are avoided by treating penetration angles as the primary raw measurements, which are adjusted using biases and scale factors. In addition to three misalignment parameters, the calibration state vector contains only two average penetration angle biases (one per each pair of opposite quadrants) since, because of the very narrow sensor FOV (+/- 2.6 degrees), differences between biases of the penetration angles measured by opposite quadrants cannot be distinguished from roll and pitch sensor misalignments. After calibration, the estimated misalignments and average penetration angle biases are converted to the four penetration angle biases and to the yaw misalignment angle. The resultant biases and the estimated scale factors are finally used to update the coefficients necessary for onboard computations of penetration angles from measured voltages.
Natanson, G. A.↗
Assessing the Effect of Spaceflight on the Propensity for Astronauts to Develop Disk Herniation
BACKGROUND: A previous study [1] reported that the instantaneous risk of developing a Herniated Nucleus Pulposus (HNP) was higher in astronauts who had flown at least one mission, as compared with those in the corps who had not yet flown. However, the study only analyzed time to HNP after the first mission (if any) and did not account for the possible effects of multiple missions. While many HNP's occurred well into astronauts' careers or in some cases years after retirement, the higher incidence of HNPs relatively soon after completion of space missions appears to indicate that spaceflight may lead to an increased risk of HNP. The purpose of this study was to support the Human System Risk Board assessment of back pain, evaluate the risk of injury due to dynamic loads, and update the previous dataset which contained events up to December 31, 2006. METHODS: Data was queried from the electronic medical record and provided by the Lifetime Surveillance of Astronaut Health. The data included all 330 United States astronauts from 1959 through February 2014. Cases were confirmed by Magnetic Resonance Imaging, Computerized Tomography, Myelography, operative findings, or through clinical confirmation with a neurologist or neurosurgeon. In this analysis, astronauts who had an HNP at selection into the corps or had an HNP diagnosis prior to their first flight were excluded. The statistical challenges in using the available data to separate effects of spaceflight from those associated with general astronaut training and lifestyle on propensity to develop HNPs are many. The primary outcome is reported date of first HNP (if any), which at best is only an approximation to the actual time of occurrence. To properly analyze this data with a survival analysis model, one must also know the "exposure" time - i.e. how long each astronaut has been at risk for developing an HNP. If an HNP is reported soon after a mission, is it mission caused or general? If the former, exposure time should be counted from the time of landing (assuming the risk of HNP occurring during a mission is zero). If the latter, exposure time should be counted from the time of selection; however we can't directly know which one to use. In our analysis we take both of these possibilities into account with a competing risks model, wherein two distinct stochastic processes are going on: TG = time to HNP (general) and TS = time to HNP (spaceflight). Under this type of model, whichever of these occurs first is what we observe; in other words we don't observe TG or TS, only min(TG, TS). Here, we parameterized the model in terms of separate Weibull hazard functions for each process and estimated all parameters using maximum likelihood. In addition, we allowed for a "cured fraction" - i.e. the possibility that some astronauts may never develop an HNP. RESULTS: Results will include a depiction of the competing hazard functions as well as a probability curve for the relative likelihood that an HNP reported at a given time after a mission is actually mission caused. Other factors, such as dwell time in microgravity and vehicle landing environment will be explored. An overall assessment as to whether spaceflight truly exacerbates HNP risk will be made.
Feiveson, A. H.↗
NASA Tech Briefs, November 1995
The contents include: 1) Mission Accomplished; 2) Resource Report: Marshall Space Flight Center; 3) NASA 1995 Software of the Year Award; 4) Microbolometers Based on Epitaxial YBa2Cu3O(sub 7-x) Thin Films; 5) Garnet Random-Access Memory; 6) Fabrication of SNS Weak Links on SOS Substrates; 7) High-Voltage MOSFET Switching Circuit; 8) Asymmetric Switching for a PWM H-Bridge Power Circuit; 9) Better Ohmic Contacts for InP Semiconductor Devices; 10) Low-Bandgap Thermovoltaic Materials and Devices; 11) Digital Frequency-Differencing Circuit; 12) Imaging Magnetometer; 13) Computer-Assisted Monitoring of a Complex System; 14) Buffered Telemetry Demodulator; 15) Compact Multifunction Inspection Head; 16) Optical Detection of Fractures in Ceramic Diaphragms; 17) Eddy-Current Detection of Cracks in Reinforced Carbon/Carbon; 18) Apparent Thermal Conductivity of Multilayer Insulation; 19) Optimizing Misch-Metal Compositions in Metal Hydride Anodes; 20) Device for Sampling Surface Contamination; 21) Probabilistic Failure Assessment for Fatigue; 22) Probabilistic Fatigue and Flaw-Propagation Analysis; 23) Windows Program for Driving the TDU-850 Printer; 24) Subband/Transform MATLAB Functions for Processing Images; 25) Computing Equilibrium Chemical Compositions; 26) Program Processes Thermocouple Readings; 27) ICAN-Second-Generation Integrated Composite Analyzer; 28) Integrated Composite Analyzer with Damping Capabilities; 29) Computing Efficiency of Transfer of Microwave Power; 30) Program Calculates Power Demands of Electronic Designs; 31) Cost-Estimation Program; 32) Program Estimates Areas Required by Electronic Designs; 33) Program to Balance Mapped Turbopump Assemblies; 34) BiblioTech; 35) Controlling Mirror Tilt With a Bimorph Actuator; 36) Burst-Disk Device Simulates Effect of Pyrotechnic Device; 37) Bearing-Mounting Concept Accommodates Thermal Expansion; 38) Parallel-Plate Acoustic Absorbers for Hot Environments; 39) Adjustable-Length Strut Withstands Large Cyclic Loads; 40) Tool Indicates Contact Angles in Bearing Raceways; 41) Gravity Slides With Magnetic Braking; 42) High-Torque, Lightweight, Pneumatically Driven Wrench for Small Spaces; 43) Device for Testing Compatibility of an O-Ring; 44) Magnetic Heat Pump Containing Flow Diverters; 45) Variable-Tilt Helicopter Rotor Mast; 46) "Beach-Ball" Robotic Rovers; 47) Apparatus Would Measure Temperatures of Ball Bearings; 48) Flexible Borescope for Inspecting Ducts; 49) Texturing Copper To Reduce Secondary Emission of Electrons; 50) Automated Laser Cutting in Three Dimensions; 51) Algorithm Helps Monitor Engine Operation; 52) Flexible Revision of Data-Processing Communications; 53) Software for Managing the Use of Land; 54) Thermal Strap Increases Cryocooling Efficiency; 55) Reversible Nut With Engagement Indication; 56) Control Algorithms for Kinematically Redundant Manipulators; 57) Computed Hydrogen-Flow Splits in a Rocket Engine; 58) Pressure and Thermal Modeling of Rocket Launches; 59) Field of View of a Spacecraft Antenna: Analysis and Software; 60) Digital Controller for Laser-Beam-Steering Subsystem; 61) More About Beam-Steering Subsystem for Laser Communication; 62) Digital Controller for Laser-Beam-Steering Subsystem: Part 2; 63) Interface Circuit Board for Space-Shuttle Communications; 64) Automated Planning of Spacecraft Telecommunications; 65) Artifacts of Spectral Analysis of Instrument Readings; 66) Neural-Network Controller for Vibration Suppression; 67) Adaptive Finite-Element Computation in Fracture Mechanics; 68) Attitude Control for the Cassini Spacecraft; 69) Analytical Model for Fluid Dynamics in a Microgravity Environment; 70) Study of Rocket-Engine Joints Bonded by NVCU/NARloy-Z; 71) Improved Silicon Nitride for Advanced Heat Engines; 72) Parameters for Welding Aluminum/Lithium Alloys; 73) Lightweight Composite Intertank Structure; 74) Foil Patches Seal Small Vacuum Leaks; 75) Data Base on Cables and Connectors; 76) Effect of Clock Mode on Radiation Hardnessf an ADC; and 77) Fault-Tolerant Control for a Robotic Inspection System.
Source record↗
Molecular scale insight of pore morphology relation with mechanical properties of amorphous silica using ReaxFF
Porous materials are typically heterogeneous and they contain large variations of micro-/nano-pore structures, causing complicated behaviors. In continuum models, most mechanical properties of porous materials are estimated based on porosity, while the variations of micro/nano structures are ignored. That could be problematic as the microscopic heterogeneity may affect the mechanical response of porous materials. Thus, understanding micro/nano heterogeneity impact has been the focus in many scientific and engineering subjects. In the present study, we investigated the effect of nanopore structure (including pore shape and orientation) as well as porosity on mechanical properties of amorphous silica (a-SiO 2 ). The pore sizes in our simulations are comparable to the corresponding ones observed in a-SiO 2 based materials. We found that the existing of nanopores strongly influences Young’s modulus (E) and critical energy release rate (G IC ). These properties decrease with increasing porosity. Importantly, the impact of nanopores was characterized by structural parameters of porous materials. In addition to dependency on porosity, Young’s modulus also was found to vary as a function of potential energy per atom, which highly depends on nanopore shape. Furthermore, critical energy release rate was found to increase with increasing ligament length (also known as pore wall thickness). Finally, the results highlighted the importance of nanopore structures, which must be taken into account when studying fracture mechanisms in porous materials. Based on our findings, it was proposed that mechanical properties of porous materials can be controlled by nano-engineering pore structures.
36 MATERIALS SCIENCE↗
Evaluation of 18 Satellite- and Model-Based Soil Moisture Products Using in Situ Measurements From 826 Sensors
Information about the spatiotemporal variability of soil moisture is critical for many purposes, including monitoring of hydrologic extremes, irrigation scheduling, and prediction of agricultural yields. We evaluated the temporal dynamics of 18 state-of-the-art (quasi-)global near-surface soil moisture products, including six based on satellite retrievals, six based on models without satellite data assimilation (referred to hereafter as “open-loop” models), and six based on models that assimilate satellite soil moisture or brightness temperature data. Seven of the products are introduced for the first time in this study: one multi-sensor merged satellite product called MeMo (Merged soil Moisture) and six estimates from the HBV (Hydrologiska Byråns Vattenbalansavdelning) model with three precipitation inputs (ERA5, IMERG, and MSWEP) with and without assimilation of SMAPL3E satellite retrievals, respectively. As reference, we used in situ soil moisture measurements between 2015 and 2019 at 5 cm depth from 826 sensors, located primarily in the USA and Europe. The 3-hourly Pearson correlation (R) was chosen as the primary performance metric. We found that application of the Soil Wetness Index (SWI) smoothing filter resulted in improved performance for all satellite products. The best-to-worst performance ranking of the four single-sensor satellite products was SMAPL3E SWI , SMOS SWI , AMSR2 SWI , and ASCAT SWI , with the L-band-based SMAPL3E SWI (median R of 0.72) outperforming the others at 50 % of the sites. Among the two multi-sensor satellite products (MeMo and ESA-CCI SWI ), MeMo performed better on average (median R of 0.72 versus 0.67), probably due to the inclusion of SMAPL3E SWI . The best-to-worst performance ranking of the six open-loop models was HBV-MSWEP, HBV-ERA5, ERA5-Land, HBV-IMERG, VIC-PGF, and GLDAS-Noah. This ranking largely reflects the quality of the precipitation forcing. HBV-MSWEP (median R of 0.78) performed best not just among the open-loop models but among all products. The calibration of HBV improved the median R by +0.12 on average compared to random parameters, highlighting the importance of model calibration. The best-to-worst performance ranking of the six models with satellite data assimilation was HBV-MSWEP+SMAPL3E, HBV-ERA5+SMAPL3E, GLEAM, SMAPL4, HBV-IMERG+SMAPL3E, and ERA5. The assimilation of SMAPL3E retrievals into HBV-IMERG improved the median R by +0.06, suggesting that data assimilation yields significant benefits at the global scale.
Hylke E. Beck↗
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↗
Interactive Soil Dust Aerosol Model in the GISS GCM: Sensitivity of the Soil Dust Cycle to Radiative Properties of Soil Dust Aerosols - Part 1
The sensitivity of the soil dust aerosol cycle to the radiative forcing by soil dust aerosols is studied. Four experiments with the NASA/GISS atmospheric general circulation model, which includes a soil dust aerosol model, are compared, all using a prescribed climatological sea surface temperature as lower boundary condition. In one experiment, dust is included as dynamic tracer only (without interacting with radiation), whereas dust interacts with radiation in the other simulations. Although the single scattering albedo of dust particles is prescribed to be globally uniform in the experiments with radiatively active dust, a different single scattering albedo is used in those experiments to estimate whether regional variations in dust optical properties, corresponding to variations in mineralogical composition among different source regions, are important for the soil dust cycle and the climate state. On a global scale, the radiative forcing by dust generally causes a reduction in the atmospheric dust load corresponding to a decreased dust source flux. That is, there is a negative feedback in the climate system due to the radiative effect of dust. The dust source flux and its changes were analyzed in more detail for the main dust source regions. This analysis shows that the reduction varies both with the season and with the single scattering albedo of the dust particles. By examining the correlation with the surface wind, it was found that the dust emission from the Saharan/Sahelian source region and from the Arabian peninsula, along with the sensitivity of the emission to the single scattering albedo of dust particles, are related to large scale circulation patterns, in particular to the trade winds during Northern Hemisphere winter and to the Indian monsoon circulation during summer. In the other regions, such relations to the large scale circulation were not found. There, the dependence of dust deflation to radiative forcing by dust particles is probably dominated by physical processes with short time scales. The experiments show that dust radiative forcing can lead to significant changes both in the soil dust cycle and in the climate state. To estimate dust concentration and radiative forcing by dust more accurately, dust size distributions and dust single scattering albedo in the model should be a function of the source region, because dust concentration and climate response to dust radiative forcing are sensitive to dust radiative parameters.
Perlwitz, Jan↗
Characterization of Partially Observed Epidemics - Application to COVID-19
This report documents a statistical method for the "real-time" characterization of partially observed epidemics. Observations consist of daily counts of symptomatic patients, diagnosed with the disease. Characterization, in this context, refers to estimation of epidemiological parameters that can be used to provide short-term forecasts of the ongoing epidemic, as well as to provide gross information for the time-dependent infection rate. The characterization problem is formulated as a Bayesian inverse problem, and is predicated on a model for the distribution of the incubation period. The model parameters are estimated as distributions using a Markov Chain Monte Carlo (MCMC) method, thus quantifying the uncertainty in the estimates. The method is applied to the COVID-19 pandemic of 2020, using data at the country, provincial (e.g., states) and regional (e.g. county) levels. The epidemiological model includes a stochastic component due to uncertainties in the incubation period. This model-form uncertainty is accommodated by a pseudo-marginal Metropolis-Hastings MCMC sampler, which produces posterior distributions that reflect this uncertainty. We approximate the discrepancy between the data and the epidemiological model using Gaussian and negative binomial error models; the latter was motivated by the over-dispersed count data. For small daily counts we find the performance of the calibrated models to be similar for the two error models. For large daily counts the negative-binomial approximation is numerically unstable unlike the Gaussian error model. Application of the model at the country level (for the United States, Germany, Italy, etc.) generally provided accurate forecasts, as the data consisted of large counts which suppressed the day-to-day variations in the observations. Further, the bulk of the data is sourced over the duration before the relaxation of the curbs on population mixing, and is not confounded by any discernible country-wide second wave of infections. At the state-level, where reporting was poor or which evinced few infections (e.g., New Mexico), the variance in the data posed some, though not insurmountable, difficulties, and forecasts were able to capture the data with large uncertainty bounds. The method was found to be sufficiently sensitive to discern the flattening of the infection and epidemic curve due to shelter-in-place orders after around 90% quantile for the incubation distribution (about 10 days for COVID-19). The proposed model was also used at a regional level to compare the forecasts for the central and north-west regions of New Mexico. Modeling the data for these regions illustrated different disease spread dynamics captured by the model. While in the central region the daily counts peaked in the late April, in the north-west region the ramp-up continued for approximately three more weeks.
59 BASIC BIOLOGICAL SCIENCES↗
Solar Dynamics Observatory On-Orbit Jitter Testing, Analysis, and Mitigation Plans
The Solar Dynamics Observatory (SDO) was designed to understand the Sun and the Sun s influence on Earth. SDO was launched on February 11, 2010 carrying three scientific instruments: the Atmospheric Imaging Assembly (AIA), the Helioseismic and Magnetic Imager (HMI), and the Extreme Ultraviolet Variability Experiment (EVE). Both AIA and HMI are sensitive to high frequency pointing perturbations and have sub-arcsecond level line-of-sight (LOS) jitter requirements. Extensive modeling and analysis efforts were directed in estimating the amount of jitter disturbing the science instruments. To verify the disturbance models and to validate the jitter performance prior to launch, many jitter-critical components and subassemblies were tested either by the mechanism vendors or at the NASA Goddard Space Flight Center (GSFC). Although detailed analysis and assembly level tests were performed to obtain good jitter predictions, there were still several sources of uncertainties in the system. The structural finite element model did not have all the modes correlated to test data at high frequencies (greater than 50 Hz). The performance of the instrument stabilization system was not known exactly but was expected to be close to the analytical model. A true disturbance-to-LOS observatory level test was not available due to the tight schedule of the flight spacecraft, the cost in time and manpower, difficulties in creating gravity negation systems, and risks of damaging flight hardware. To protect the observatory jitter performance against model uncertainties, the SDO jitter team devised several on-orbit jitter reduction plans in addition to reserve margins on analysis results. Since some of these plans severely restricted the capabilities of several spacecraft components (e.g. wheels and High Gain Antennas), the SDO team performed on-orbit jitter tests to determine which jitter reduction plans, if any, were necessary to satisfy science LOS jitter requirements. The SDO on-orbit jitter tests were designed to satisfy the following four objectives: 1. Determine the acceptable reaction wheel operational speed range during Science Mode. 2. Determine HGA algorithm jitter parameters (number of stagger steps and enable/disable no-steprequests). 3. Determine acceptable EVE instrument filter wheels spin rates. 4. Determine if AIA instrument filter wheels excite the first AIA telescope structural mode. This paper provides detailed information on the SDO wheel jitter test plan, shows on-orbit jitter measurements and how ground predictions compare to those measurements, and describes the final jitter mitigation plan executed on SDO.
Liu, Kuo-Chia (Alice)↗
Development of a Comparison Framework for Evaluating Environmental Contours of Extreme Sea States
Environmental contours of extreme sea states are often utilized for the purposes of reliability-based offshore design. Many methods have been proposed to estimate environmental contours of extreme sea states, including, but not limited to, the traditional inverse first-order reliability method (I-FORM) and subsequent modifications, copula methods, and Monte Carlo methods. These methods differ in terms of both the methodology selected for defining the joint distribution of sea state parameters and in the method used to construct the environmental contour from the joint distribution. It is often difficult to compare the results of proposed methods to determine which method should be used for a particular application or geographical region. The comparison of the predictions from various contour methods at a single site and across many sites is important to making environmental contours of extreme sea states useful in practice. The goal of this paper is to develop a comparison framework for evaluating methods for developing environmental contours of extreme sea states. This paper develops generalized metrics for comparing the performance of contour methods to one another across a collection of study sites, and applies these metrics and methods to develop conclusions about trends in the wave resource across geographic locations, as demonstrated for a pilot dataset. These proposed metrics and methods are intended to judge the environmental contours themselves relative to other contour methods, and are thus agnostic to a specific device, structure, or field of application. The metrics developed and applied in this paper include measures of predictive accuracy, physical validity, and aggregated temporal performance that can be used to both assess contour methods and provide recommendations for the use of certain methods in various geographical regions. The application and aggregation of the metrics proposed in this paper outline a comparison framework for environmental contour methods that can be applied to support design analysis workflows for offshore structures. This comparison framework could be extended in future work to include additional metrics of interest, potentially including those to address issues pertinent to a specific application area or analysis discipline, such as metrics related to structural response across contour methods or additional physics-based metrics based on wave dynamics.
54 ENVIRONMENTAL SCIENCES↗
Modeling Streamflow and Sediment Transport Under Current and Future Climates at the West Valley Site - 20231
Performance Assessment models often require inputs representing specific hydrologic processes such as streamflow and sediment flux. This paper details the development of a process-level hydrologic model and its use in a Probabilistic Performance Assessment (PPA) system-level model for the West Valley Site in western New York. For this PPA model, a watershed model was developed using SWAT (Soil and Water Assessment Tool) and used to evaluate water and sediment fluxes in nearby streams draining the Site. A particular challenge of probabilistic modeling over long periods of time is scaling parameters appropriately so that outcomes are not grossly over- or underestimated. This challenge is especially important when the spatial and temporal scales of process and system level models differ. In this instance, the SWAT watershed model and GoldSim PPA model have similar spatial scales of stream reaches and sub-catchments within a single watershed (approximately hundreds to tens of thousands of square meters, but the temporal scales differ considerably. The SWAT model calculates streamflow and sediment transport rates at the daily scale, while the PPA model requires long-term average values that are applied over hundreds to thousands of year-long periods. As such, distributions of daily climatic parameters such as average precipitation rate, temperature, wind speed, and relative humidity were developed to force a transient SWAT model that captures the dynamic behavior of the watershed in response to daily changes in the weather. Output from the model, which includes stream flow (m{sup 3}/y) and sediment transport rate (Mg/y) is then averaged over a 100-year period as an approximation for the long-term average values which can be applied in GoldSim. A single SWAT model was run 5000 times with varying climatic inputs and deterministic soil, land use, and elevation data, and the range of outputs were compiled into a distribution which can be implemented stochastically into the PPA Model. Additionally, this paper includes a high-level discussion of how changes in future climate affect weather input distributions, and a summary of how those changes in weather affect the outcomes from the SWAT model and the inputs into the PPA model. In this analysis, future climate distributions were developed based on climate research documenting the expected changes to precipitation and temperature in western New York. The 5,000 SWAT simulations were then repeated with these future climate distributions, and the resulting output was implemented into the PPA Model as hydrologic conditions under a future climate. In the current version of the PPA Model, hydrologic parameters change linearly from the 'current climate' conditions to the 'future climate' conditions over a period of 100 years, after which they remain constant. (authors)
12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
Field and Model Data Associated with the Manuscript “Drivers of Streamflow Intermittency in Humid Regions: 1. Evaluating Above- and Below-ground Controls of Flow Persistence in a Forested Catchment”
This package contains field data, modeling files, and scripts supporting the investigation of the drivers of streamflow intermittency in a forested catchment. It includes the field data collected from electrical resistivity tomography (ERT) surveys, ground penetrating radar (GPR), continuous self-potential (SP) monitoring, electromagnetic (EM) imaging, groundwater and stilling well. In addition, it contains the data and results of the coupled water- and electrical-flow model developed using the COMSOL Multiphysics and Advanced Terrestrial Simulator (ATS), as well as software files and Jupyter notebooks used to process the data and generate figures in the manuscript submitted for peer review. The data archive is organized in the following directories: 1) Climate Includes hourly precipitation and daily evapotranspiration time series (2024 – 2025) provided as CSV files, alongside a text file detailing dataset units. 2) Coupled_model Contains two subfolders: Synthetic and Field_Application subfolder. Synthetic subfolder contains the ATS XML input script (can be opened using any code editor) for the four synthetic hydrological cases tested (Connected and gaining, Connected and losing, Disconnected and losing, and dry stream). It also includes other experimental cases to test the influence of precipitation and concentration gradient. For each synthetic case, the flow model simulation is executed using the ATS XML scripts and the included Python script (generate_data_set.py) to convert ATS output to COMSOL-ready input. COMSOL Multiphysics template (.mph can be opened with the commercial software COMSOL and requires a license) is executed using the ATS output data to simulate the potential field. It also includes the Synthetic_model_plot.ipynb (can be opened using any code editor) to visualize the SP result and generate manuscript figures. The data subfolder contains mesh files to run both the ATS (.exo and .stl files can be viewed using Paraview; .h5 files can be opened using HDFView software and h5py Python package) and COMSOL models. Field_Application subfolder contains two subfolders: ES_MDA_inversion and Final_Model. ES_MDA_inversion contains the Python script (.py can be opened using any code editor) and SP observation data used to run the Ensemble Smoother with Multiple Data Assimilation (ES-MDA) inversion sequence to get the optimal model parameters. The Final_model subfolder contains the ATS XML input scripts, data files, output data for the two SP sites. The same workflow steps outlined for the Synthetic subfolder apply here. It also contains the Jupyter notebook (Plot_final_calib.ipynb) to visualize the results of the modeled SP, stream-groundwater exchange and moisture content. 3) Discharge Includes the electrical conductivity (EC) time series (provided as CSV files) from salt slug injections. It also includes the Jupyter notebook (Discharge_process.ipynyb) used to estimate discharge. All discharge measurements collated into rating_curve_processed.csv 4) EM Contains the CSV file of the EM data from the DUALEM-42, including spatial coordinates (x, y, z), apparent conductivity, and in-phase measurements at 2 m coil separations for horizontal coplanar (HCP) and perpendicular (PRP) geometries. 5) ERT Contains raw resistivity data (provided as CSV files), spatial location of each of the electrodes (provided as CSV files), and files used for the resistivity inversion (.resipy can be opened with the open-source ResIPy software). 6) GPR Includes GPR field datasets collected at 100 MHz and 250 MHz antenna frequencies, along with the processing/interpretation project file (GPR_process.gpz can be viewed using EKKO_Project 6, a commercial software by Sensors & Software that requires a license). 7) Slug_test Includes the slug test data at all the groundwater wells provided as CSV files, as well as the Jupyter notebook (Slug_test.ipynb) for calculating hydraulic conductivity. 8) SP Contains the SP data collected in field at the two SP sites (one in the perennial reach and the other in the intermittent reach), provided as DAT files. 9) Well_data Contains two subfolders: 1) Raw, which provides unprocessed pressure, electrical conductivity and temperature timeseries downloaded from the loggers in all the groundwater and stilling wells, and 2) Processed, which contains sorted, QA/QC timeseries data for each well. The data archive also contains data_process.ipynb, a Jupyter notebook used for field data analysis and generating figures (plotting well, SP, climate, and discharge data, as well as calculating head gradient at sites with nested groundwater wells). It also includes DTW.ipynb, a Jupyter notebook containing the code for the dynamic time warping (DTW) with sliding window to evaluate SP signal synchronicity.
ATS↗
Predicting secondary organic aerosol phase state and viscosity and its effect on multiphase chemistry in a regional-scale air quality model
Atmospheric aerosols are a significant public health hazard and have substantial impacts on the climate. Secondary organic aerosols (SOAs) have been shown to phase separate into a highly viscous organic outer layer surrounding an aqueous core. This phase separation can decrease the partitioning of semi-volatile and low-volatile species to the organic phase and alter the extent of acid-catalyzed reactions in the aqueous core. A new algorithm that can determine SOA phase separation based on their glass transition temperature (T g ), oxygen to carbon (O:C) ratio and organic mass to sulfate ratio, and meteorological conditions was implemented into the Community Multiscale Air Quality Modeling (CMAQ) system version 5.2.1 and was used to simulate the conditions in the continental United States for the summer of 2013. SOA formed at the ground/surface level was predicted to be phase separated with core–shell morphology, i.e., aqueous inorganic core surrounded by organic coating 65.4 % of the time during the 2013 Southern Oxidant and Aerosol Study (SOAS) on average in the isoprene-rich southeastern United States. Our estimate is in proximity to the previously reported ~70 % in literature. The phase states of organic coatings switched between semi-solid and liquid states, depending on the environmental conditions. The semi-solid shell occurring with lower aerosol liquid water content (western United States and at higher altitudes) has a viscosity that was predicted to be 10 2 –10 12 Pa s, which resulted in organic mass being decreased due to diffusion limitation. Organic aerosol was primarily liquid where aerosol liquid water was dominant (eastern United States and at the surface), with a viscosity <10 2 Pa s. Phase separation while in a liquid phase state, i.e., liquid–liquid phase separation (LLPS), also reduces reactive uptake rates relative to homogeneous internally mixed liquid morphology but was lower than aerosols with a thick viscous organic shell. The sensitivity cases performed with different phase-separation parameterization and dissolution rate of isoprene epoxydiol (IEPOX) into the particle phase in CMAQ can have varying impact on fine particulate matter (PM 2.5 ) organic mass, in terms of bias and error compared to field data collected during the 2013 SOAS. This highlights the need to better constrain the parameters that govern phase state and morphology of SOA, as well as expand mechanistic representation of multiphase chemistry for non-IEPOX SOA formation in models aided by novel experimental insights.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Subsurface Characterization and Machine Learning Predictions at Brady Hot Springs Results
Geothermal power plants typically show decreasing heat and power production rates over time. Mitigation strategies include optimizing the management of existing wells - increasing or decreasing the fluid flow rates across the wells - and drilling new wells at appropriate locations. The latter is expensive, time-consuming, and subject to many engineering constraints, but the former is a viable mechanism for periodic adjustment of the available fluid allocations. Data and supporting literature from a study describing a new approach combining reservoir modeling and machine learning to produce models that enable strategies for the mitigation of decreased heat and power production rates over time for geothermal power plants. The computational approach used enables translation of sets of potential flow rates for the active wells into reservoir-wide estimates of produced energy and discovery of optimal flow allocations among the studied sets. In our computational experiments, we utilize collections of simulations for a specific reservoir (which capture subsurface characterization and realize history matching) along with machine learning models that predict temperature and pressure timeseries for production wells. We evaluate this approach using an "open-source" reservoir we have constructed that captures many of the characteristics of Brady Hot Springs, a commercially operational geothermal field in Nevada, USA. Selected results from a reservoir model of Brady Hot Springs itself are presented to show successful application to an existing system. In both cases, energy predictions prove to be highly accurate: all observed prediction errors do not exceed 3.68% for temperatures and 4.75% for pressures. In a cumulative energy estimation, we observe prediction errors that are less than 4.04%. A typical reservoir simulation for Brady Hot Springs completes in approximately 4 hours, whereas our machine learning models yield accurate 20-year predictions for temperatures, pressures, and produced energy in 0.9 seconds. This paper aims to demonstrate how the models and techniques from our study can be applied to achieve rapid exploration of controlled parameters and optimization of other geothermal reservoirs. Includes a synthetic, yet realistic, model of a geothermal reservoir, referred to as open-source reservoir (OSR). OSR is a 10-well (4 injection wells and 6 production wells) system that resembles Brady Hot Springs (a commercially operational geothermal field in Nevada, USA) at a high level but has a number of sufficiently modified characteristics (which renders any possible similarity between specific characteristics like temperatures and pressures as purely random). We study OSR through CMG simulations with a wide range of flow allocation scenarios. Includes a dataset with 101 simulated scenarios that cover the period of time between 2020 and 2040 and a link to the published paper about this project, where we focus on the Machine Learning work for predicting OSR's energy production based on the simulation data, as well as a link to the GitHub repository where we have published the code we have developed (please refer to the repository's readme file to see instructions on how to run the code). Additional links are included to associated work led by the USGS to identify geologic factors associated with well productivity in geothermal fields. Below are the high-level steps for applying the same modeling + ML process to other geothermal reservoirs: 1. Develop a geologic model of the geothermal field. The location of faults, upflow zones, aquifers, etc. need to be accounted for as accurately as possible 2. The geologic model needs to be converted to a reservoir model that can be used in a reservoir simulator, such as, for instance, CMG STARS, TETRAD, or FALCON 3. Using native state modeling, the initial temperature and pressure distributions are evaluated, and they become the initial conditions for dynamic reservoir simulations 4....
15 GEOTHERMAL ENERGY↗
Decision Analyses Related to Groundwater Remediation - 17051
Traditional decision analyses of groundwater remediation scenarios frequently fail because the probability of adverse, unanticipated events occurring is often high. This is driven by many factors including that (1) the models developed to represent the flow and transport in contaminated aquifers are typically simpler than reality and do not account for all the important governing processes, (2) probability distributions are assigned to critical performance parameters (such as model inputs or characteristics of the engineered remediation systems) even though some of these parameters might be unknown or not be very well constrained. However, there is a novel decision analysis methodology and a computational tool based on Bayesian Information-Gap Decision Theory (BIG-DT) that are designed to mitigate the shortcomings of the models and probabilistic decision analyses by leveraging a non-probabilistic decision theory -- information-gap decision theory. The new decision analysis methodology considers possible models that have not been explicitly enumerated and does not require us to commit to a particular probability distribution for model and remediation-design parameters. Both the set of possible models and the set of possible probability distributions grow as the degree of uncertainty increases. The fundamental question that BIG-DT asks is 'How large can these sets be before a particular decision results in an undesirable outcome?'. The decision that allows these sets to be the largest is considered to be the best option. In this way, BIG-DT enables robust decision-support for groundwater remediation problems. Here, we apply BIG-DT to a representative groundwater remediation scenario where different options for hydraulic containment and pump and treat are being considered. BIG-DT requires many model runs and high-performance computing resources are needed when complex models are used. We demonstrate the BIG-DT analyses on a series of synthetic problems that are designed to be consistent with real-world problems such as Los Alamos National Laboratory (LANL) contamination sites. We also discuss a general framework for how the BIG-DT analyses will be carried out at the chromium groundwater contamination site at LANL. BIG-DT is implemented in Julia (a high-level, high-performance dynamic programming language for technical computing) and is part of the MADS framework (http://mads.lanl.gov and http://madsjulia.github.io/Mads.jl). (authors)
12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
On the mechanism of ionization oscillations in Hall thrusters
Low-frequency ionization oscillations involving plasma and neutral density (breathing modes) are the most violent perturbations in Hall thrusters for electric propulsion. Because of its simplicity, the zero-dimensional (0D) predator–prey model of two nonlinearly coupled ordinary differential equations for plasma and neutral density has often been used for the characterization of such oscillations and scaling estimates. In this work, we investigate the properties of its continuum analog, the one-dimensional (1D) system of two nonlinearly coupled equations in partial derivatives (PDEs) for plasma and neutral density. This is a more general model, of which the standard 0D predator–prey model is a special limit case. We show that the 1D model is stable and does not show any oscillations for the boundary conditions relevant to Hall thrusters and the uniform ion velocity. We then propose a reduced 1D model based on two coupled PDEs for plasma and neutral densities that is unstable and exhibit oscillations if the ion velocity profile with the near-the-anode back-flow (toward the anode) region is used. Comparisons of the reduced model with the predictions of the full model that takes into account the self-consistent plasma response show that the main properties of the breathing mode are well captured. In particular, it is shown that the frequency of the breathing mode oscillations is weakly dependent on the final ion velocity but shows a strong correlation with the width of the ion back-flow region.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗