SEARCH · Search NASA

Results for “backward”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 343 records · Page 19

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Effects of Nonzero-frequency Fluctuations on Turbulence Spectral Observations

In situ observations of turbulence spectra in space plasmas are usually interpreted as wavenumber spectra, assuming that the fluctuation frequency is negligible in the plasma flow frame. We explore the effects of nonzero frequency in the plasma flow frame on turbulence spectral observations. The finite frequency can be caused by either propagating waves or nonlinear broadening of nonpropagating structures. We show that the observed frequency spectrum can be modified by the nonzero frequency of turbulent fluctuations in several ways. Specifically, (i) frequency broadening results in a minor modification to the observed spectrum, primarily acting as a smoothing kernel of the spectrum near the spectral break, while the asymptotic spectral index remains unchanged; (ii) wave propagation can affect the observed spectral index for anisotropic turbulence. The effect is significant at low frequencies and weaker at high frequencies, leading to a "concave" shape of the observed perpendicular spectrum; (iii) the Doppler shift for forward- and backward-propagating Elsasser modes can result in a nonzero cross helicity for critical-balanced turbulence since the effect of the Doppler shift favors outward-propagating waves systematically, resulting in an observed imbalance. These results may have important implications for the interpretation of solar wind flows observed by Parker Solar Probe.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Assimilating partial observation to enhance feedback control of stochastic dynamical systems

Here, in this paper, we present a novel methodology to tackle feedback optimal control problems in scenarios where the exact state of the controlled process is unknown. It integrates data assimilation techniques and optimal control solvers to manage partial observation of the state process, a common occurrence in practical scenarios. Traditional stochastic optimal control methods assume full state observation, which is often not feasible in real-world fluid dynamics control problems. Our approach underscores the significance of utilizing observational data to inform control policy design. Specifically, we introduce a kernel learning backward stochastic differential equation (SDE) filter to enhance data assimilation efficiency and propose a sample-wise stochastic optimization method within the stochastic maximum principle framework. We demonstrate the efficacy and accuracy of our method in the control of advection-diffusion-reaction flow problem and the Dubins airplane maneuvering problem with model uncertainty.

data driven↗

DOE OSTI2025

Scrambling Dynamics with Imperfections in a Solvable Model

We study how probes of quantum scrambling dynamics respond to two kinds of imperfections -- unequal forward and backward evolutions and decoherence -- in a solvable Brownian circuit model. We calculate a ``renormalized'' out-of-time-order correlator (ROTOC) in the model with $N$ qubits, and we show that the circuit-averaged ROTOC is controlled by an effective probability distribution in operator weight space which obeys a system of $N$ non-linear equations of motion. These equations can be easily solved numerically for large system sizes which are beyond the reach of exact methods. Moreover, for an operator initially concentrated on weight one $w_0=1$, we provide an exact solution to the equations in the thermodynamic limit of many qubits that is valid for all times, all non-vanishing perturbation strengths $p\gtrsim 1/\sqrt{N}$, and all decoherence strengths. We also show that a generic initial condition $w_0 >1$ leads to a metastable state that eventually collapses to the $w_0=1$ case after a lifetime $\sim \log(N/w_0)$. Our results highlight situations where it is still possible to extract the unperturbed chaos exponent even in the presence of imperfections, and we comment on the applications of our results to existing experiments with nuclear spins and to future scrambling experiments.

FOS: Physical sciences↗

DOE OSTI2021

A global study of hygroscopicity-driven light-scattering enhancement in the context of other in situ aerosol optical properties

The scattering and backscattering enhancement factors (ƒ(RH) and ƒ b (RH)) describe how aerosol particle light scattering and backscattering, respectively, change with relative humidity (RH). They are important parameters in estimating direct aerosol radiative forcing (DARF). In this study we use the dataset presented in Burgos et al. (2019) that compiles ƒ(RH) and ƒ b (RH) measurements at three wavelengths (i.e., 450, 550 and 700 nm) performed with tandem nephelometer systems at multiple sites around the world. We present an overview of ƒ(RH) and ƒ b (RH) based on both long-term and campaign observations from 23 sites representing a range of aerosol types. The scattering enhancement shows a strong variability from site to site, with no clear pattern with respect to the total scattering coefficient. In general, higher ƒ(RH) is observed at Arctic and marine sites, while lower values are found at urban and desert sites, although a consistent pattern as a function of site type is not observed. The backscattering enhancement ƒ b (RH) is consistently lower than ƒ(RH) at all sites, with the difference between ƒ(RH) and ƒ b (RH) increasing for aerosol with higher ƒ(RH). This is consistent with Mie theory, which predicts higher enhancement of the light scattering in the forward than in the backward direction as the particle takes up water. Our results show that the scattering enhancement is higher for PM 1 than PM 10 at most sites, which is also supported by theory due to the change in scattering efficiency with the size parameter that relates particle size and the wavelength of incident light. At marine-influenced sites this difference is enhanced when coarse particles (likely sea salt) predominate. For most sites, ƒ(RH) is observed to increase with increasing wavelength, except at sites with a known dust influence where the spectral dependence of ƒ(RH) is found to be low or even exhibit the opposite pattern. The impact of RH on aerosol properties used to calculate radiative forcing (e.g., single-scattering albedo, ω 0 , and backscattered fraction, b) is evaluated. The single-scattering albedo generally increases with RH, while b decreases. The net effect of aerosol hygroscopicity on radiative forcing efficiency (RFE) is an increase in the absolute forcing effect (negative sign) by a factor of up to 4 at RH = 90 % compared to dry conditions (RH < 40 %). Because of the scarcity of scattering enhancement measurements, an attempt was made to use other more commonly available aerosol parameters (i.e., ω 0 and scattering Ångström exponent, α sp ) to parameterize ƒ(RH). The majority of sites (75 %) showed a consistent trend with ω 0 (higher ƒ(RH = 85 %) for higher ω 0 ), while no clear pattern was observed between ƒ(RH = 85 %) and αsp. This suggests that aerosol ω 0 is more promising than α sp as a surrogate for the scattering enhancement factor, although neither parameter is ideal. Nonetheless, the qualitative relationship observed between ω 0 and ƒ(RH) could serve as a constraint on global model simulations.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Physicochemical characterization and source apportionment of Arctic ice-nucleating particles observed in Ny-Ålesund in autumn 2019

Abstract. Ice-nucleating particles (INPs) initiate primary ice formation in Arctic mixed-phase clouds (MPCs), altering cloud radiative properties and modulating precipitation. For atmospheric INPs, the complexity of their spatiotemporal variations, heterogeneous sources, and evolution via intricate atmospheric interactions challenge the understanding of their impact on microphysical processes in Arctic MPCs and induce an uncertain representation in climate models. In this work, we performed a comprehensive analysis of atmospheric aerosols at the Arctic coastal site in Ny-Ålesund (Svalbard, Norway) from October to November 2019, including their ice nucleation ability, physicochemical properties, and potential sources. Overall, INP concentrations (NINP) during the observation season were approximately up to 3 orders of magnitude lower compared to the global average, with several samples showing degradation of NINP after heat treatment, implying the presence of proteinaceous INPs. Particle fluorescence was substantially associated with INP concentrations at warmer ice nucleation temperatures, indicating that in the far-reaching Arctic, aerosols of biogenic origin throughout the snow- and ice-free season may serve as important INP sources. In addition, case studies revealed the links between elevated NINP and heat lability, fluorescence, high wind speeds originating from the ocean, augmented concentration of coarse-mode particles, and abundant organics. Backward trajectory analysis demonstrated a potential connection between high-latitude dust sources and high INP concentrations, while prolonged air mass history over the ice pack was identified for most scant INP cases. The combination of the above analyses demonstrates that the abundance, physicochemical properties, and potential sources of INPs in the Arctic are highly variable despite its remote location.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Long-range transported continental aerosol in the eastern North Atlantic: three multiday event regimes influence cloud condensation nuclei

Abstract. The eastern North Atlantic (ENA) is a region dominated by pristine marine environment and subtropical marine boundary layer clouds. Under unperturbed atmospheric conditions, the regional aerosol regime in the ENA varies seasonally due to different seasonal surface-ocean biogenic emissions, removal processes, and meteorological regimes. However, during periods when the marine boundary layer aerosol in the ENA is impacted by particles transported from continental sources, aerosol properties within the marine boundary layer change significantly, affecting the concentration of cloud condensation nuclei (CCN). Here, we investigate the impact of long-range transported continental aerosol on the regional aerosol regime in the ENA using data collected at the U.S. Department of Energy's (DOE) Atmospheric Radiation Measurement (ARM) user facility on Graciosa Island in 2017 during the Aerosol and Cloud Experiments in the Eastern North Atlantic (ACE-ENA) campaign. We develop an algorithm that integrates number concentrations of particles with optical particle dry diameter (Dp) between 100 and 1000 nm, single scattering albedo, and black carbon concentration to identify multiday events (with duration >24 consecutive hours) of long-range continental aerosol transport in the ENA. In 2017, we detected nine multiday events of long-range transported particles that correspond to ∼ 7.5 % of the year. For each event, we perform HYSPLIT 10 d backward trajectories analysis, and we evaluate CALIPSO aerosol products to assess, respectively, the origins and compositions of aerosol particles arriving at the ENA site. Subsequently, we group the events into three categories, (1) mixture of dust and marine aerosols, (2) mixture of marine and polluted continental aerosols from industrialized areas, and (3) biomass burning aerosol from North America and Canada, and we evaluate their influence on aerosol population and cloud condensation nuclei in terms of potential activation fraction and concentrations at supersaturation of 0.1 % and 0.2 %. The arrival of plumes dominated by the mixture of dust and marine aerosol in the ENA in the winter caused significant increases in baseline Ntot. Simultaneously, the baseline particle size modes and CCN potential activation fraction remained almost unvaried, while cloud condensation nuclei concentrations increased proportionally to Ntot. Events dominated by a mixture of marine and polluted continental aerosols in spring, fall, and winter led to a statistically significant increase in baseline Ntot, a shift towards larger particular sizes, a higher CCN potential activation fractions, and cloud condensation nuclei concentrations of >170 % and up to 240 % higher than during baseline regime. Finally, the transported aerosol plumes characterized by elevated concentration of biomass burning aerosol from continental wildfires detected in the summertime did not statistically contribute to increase baseline aerosol particle concentrations in the ENA. However, particle diameters were larger than under baseline conditions, and CCN potential activation fractions were >75 % higher. Consequentially, cloud concentration nuclei concentrations increased by ∼ 115 % during the period affected by the biomass burning events. Our results suggest that, through the year, multiday events of long-range continental aerosol transport periodically affect the ENA and represent a significant source of CCN in the marine boundary layer. Based on our analysis, in 2017, the multiday aerosol plume transport dominated by a mixture of dust and marine aerosol, a mixture of marine and polluted continental aerosols, and biomass burning aerosols caused increases in the NCCN baseline regime of, respectively, 6.6 %, 8 %, and 7.4 % at SS 0.1 % (and, respectively, 6.5 %, 8.2 %, and 7.3 % at SS 0.2 %) in the ENA.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Measurement report: A comparison of ground-level ice-nucleating-particle abundance and aerosol properties during autumn at contrasting marine and terrestrial locations

Abstract. Ice-nucleating particles (INPs) are an essential class of aerosols found worldwide that have far-reaching but poorly quantified climate feedback mechanisms through interaction with clouds and impacts on precipitation. These particles can have highly variable physicochemical properties in the atmosphere, and it is crucial to continuously monitor their long-term concentration relative to total ambient aerosol populations at a wide variety of sites to comprehensively understand aerosol–cloud interactions in the atmosphere. Hence, our study applied an in situ forced expansion cooling device to measure ambient INP concentrations and test its automated continuous measurements at atmospheric observatories, where complementary aerosol instruments are heavily equipped. Using collocated aerosol size, number, and composition measurements from these sites, we analyzed the correlation between sources and abundance of INPs in different environments. Toward this aim, we have measured ground-level INP concentrations at two contrasting sites, one in the Southern Great Plains (SGP) region of the United States with a substantial terrestrially influenced aerosol population and one in the Eastern North Atlantic Ocean (ENA) region with a primarily marine-influenced aerosol population. These measurements examined INPs mainly formed through immersion freezing and were performed at a ≤ 12 min resolution and with a wide range of heterogeneous freezing temperatures (Ts above −31 °C) for at least 45 d at each site. The associated INP data analysis was conducted in a consistent manner. We also explored the additional offline characterization of ambient aerosol particle samples from both locations in comparison to in situ data. From our ENA data, on average, INP abundance ranges from ≈ 1 to ≈ 20 L−1 (−30 °C ≤ T ≤ −20 °C) during October–November 2020. Backward air mass trajectories reveal a strong marine influence at ENA with 75.7 % of air masses originating over the Atlantic Ocean and 96.6 % of air masses traveling over open water, but analysis of particle chemistry suggests an additional INP source besides maritime aerosols (e.g., sea spray aerosols) at ENA. In contrast, 90.8 % of air masses at the SGP location originated from the North American continent, and 96.1 % of the time, these air masses traveled over land. As a result, organic-rich SGP aerosols from terrestrial sources exhibited notably high INP abundance from ≈ 1 to ≈ 100 L−1 (−30 °C ≤ T ≤ −15 °C) during October–November 2019. The probability density function of aerosol surface area-scaled immersion freezing efficiency (ice nucleation active surface site density; ns) was assessed for selected freezing temperatures. While the INP concentrations measured at SGP are higher than those of ENA, the ns(T) values of SGP (≈ 105 to ≈ 107 m−2 for −30 °C ≤ T ≤ −15 °C) are reciprocally lower than ENA for approximately 2 orders of magnitude (≈ 107 to ≈ 109 m−2 for −30 °C ≤ T ≤ −15 °C). The observed difference in ns(T) mainly stems from varied available aerosol surface areas, Saer, from two sites (Saer,SGP > Saer,ENA). INP parameterizations were developed as a function of examined freezing temperatures from SGP and ENA for our study periods.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Evaluating the feasibility of using downwind methods to quantify point source oil and gas emissions using continuously monitoring fence-line sensors

The dependable reporting of methane (CH 4 ) emissions from point sources, such as fugitive leaks from oil and gas infrastructure, is important for profit maximization (retaining more hydrocarbons), evaluating climate impacts, assessing CH 4 fees for regulatory programs, and validating CH 4 intensity in differentiated gas programs. Currently, there are disagreements between emissions reported by different quantification techniques for the same sources. It has been suggested that downwind CH 4 quantification methods using CH 4 measurements on the fence line of production facilities could be used to generate emission estimates from oil and gas operations at the site level, but it is currently unclear how accurate the quantified emissions are. To investigate the accuracy of downwind methods, this study uses fence-line simulated data collected during controlled-release experiments as input for a non-standard closed-path eddy covariance (EC), the Gaussian plume inverse model (GPIM), and the backward Lagrangian stochastic (bLs) model in a range of atmospheric conditions. This study's EC attempt was unsuccessful due to data collection and instrumentation issues, resulting in invalid results characterized by underestimated emissions, large negative fluxes, and cospectra/ogives that deviated from their ideal shapes. Consequently, the EC results could not be compared with the GPIM and bLS model. The bLs model demonstrated the highest accuracy for single-release single-point emissions, though it exhibited greater uncertainty than GPIM under multi-release conditions. Across the GPIM and bLs model, the most reliable quantification was achieved with 15 min averaging and a narrow 5° wind sector range. Although EC was limited in this context, future studies should consider employing a standard EC system and further optimizing GPIM and bLs approaches – particularly for complex multi-source scenarios – to enhance quantification accuracy and reduce uncertainty.

03 NATURAL GAS↗

DOE OSTI2025

Observational ozone datasets over the global oceans and polar regions (version 2024)

Studying tropospheric ozone over the remote areas of the planet, such as the open oceans and the polar regions, is crucial to understand the role of ozone as a global climate forcer and regulator of atmospheric oxidative capacity. A focus on the pristine oceanic and polar regions complements the available land-based datasets and provides insights into key photochemical and depositional loss processes that control the concentrations and spatiotemporal variability in ozone as well as the physicochemical mechanisms driving these patterns. However, an assessment of the role of ozone over the oceanic and polar regions has been hampered by a lack of comprehensive observational datasets. Here, we present the first comprehensive collection of ozone data over the oceans and the polar regions. The overall dataset consists of 77 ship cruises/buoy-based observations and 48 aircraft-based campaigns. The dataset, consisting of more than 630 000 independent ozone measurement data points covering the period from 1977 to 2022 and an altitude range from the surface to 5000 m (with a focus on the lowest 2000 m), allows systematic analyses of the spatiotemporal distribution and long-term trends over the 11 defined ocean/polar regions. The datasets from ships, buoys, and aircraft are complemented by ozonesonde data from 29 launch sites or field campaigns and by 21 non-polar and 17 polar ground-based station datasets. The datasets contain information on how long the observed air masses were isolated from land, as estimated by backward trajectories from the individual observation points. To extract observations representative of oceanic conditions, we recommend using a subset of the data with an isolation time of 72 h or longer, from the analysis with coincident radon observations. These filtered oceanic and polar data showed typically flat diurnal cycles at high latitudes, whereas daytime decreases in ozone (11 %–16 %) were observed at lower latitudes. The ship/buoy- and aircraft-based datasets presented here will supplement the land-based ones in the TOAR-II (Tropospheric Ozone Assessment Report Phase II) database to provide a fully global assessment of tropospheric ozone. The described dataset is available at https://doi.org/10.17596/0004044 (Kanaya et al., 2025).

Kanaya, Yugo [Japan Agency for Marine-Earth Scienc↗

DOE OSTI2021

Testing the reliability of interpretable neural networks in geoscience using the Madden–Julian oscillation

Abstract. We test the reliability of two neural network interpretation techniques, backward optimization and layerwise relevance propagation, within geoscientific applications by applying them to a commonly studied geophysical phenomenon, the Madden–Julian oscillation. The Madden–Julian oscillation is a multi-scale pattern within the tropical atmosphere that has been extensively studied over the past decades, which makes it an ideal test case to ensure the interpretability methods can recover the current state of knowledge regarding its spatial structure. The neural networks can, indeed, reproduce the current state of knowledge and can also provide new insights into the seasonality of the Madden–Julian oscillation and its relationships with atmospheric state variables. The neural network identifies the phase of the Madden–Julian oscillation twice as accurately as a linear regression approach, which means that nonlinearities used by the neural network are important to the structure of the Madden–Julian oscillation. Interpretations of the neural network show that it accurately captures the spatial structures of the Madden–Julian oscillation, suggest that the nonlinearities of the Madden–Julian oscillation are manifested through the uniqueness of each event, and offer physically meaningful insights into its relationship with atmospheric state variables. We also use the interpretations to identify the seasonality of the Madden–Julian oscillation and find that the conventionally defined extended seasons should be shifted later by 1 month. More generally, this study suggests that neural networks can be reliably interpreted for geoscientific applications and may thereby serve as a dependable method for testing geoscientific hypotheses.

58 GEOSCIENCES↗

DOE OSTI2024

Development of the tangent linear and adjoint models of the global online chemical transport model MPAS-CO 2 v7.3

We describe the development of the tangent linear (TL) and adjoint models of the Model for Prediction Across Scales (MPAS)-CO 2 transport model, which is a global online chemical transport model developed upon the non-hydrostatic Model for Prediction Across Scales – Atmosphere (MPAS-A). The primary goal is to make the model system a valuable research tool for investigating atmospheric carbon transport and inverse modeling. First, we develop the TL code, encompassing all CO 2 transport processes within the MPAS-CO 2 forward model. Then, we construct the adjoint model using a combined strategy involving re-calculation and storage of the essential meteorological variables needed for CO 2 transport. This strategy allows the adjoint model to undertake a long-period integration with moderate memory demands. To ensure accuracy, the TL and adjoint models undergo vigorous verifications through a series of standard tests. The adjoint model, through backward-in-time integration, calculates the sensitivity of atmospheric CO 2 observations to surface CO 2 fluxes and the initial atmospheric CO 2 mixing ratio. To demonstrate the utility of the newly developed adjoint model, we conduct simulations for two types of atmospheric CO 2 observations, namely the tower-based in situ CO 2 mixing ratio and satellite-derived column-averaged CO 2 mixing ratio (X CO 2 ). A comparison between the sensitivity to surface flux calculated by the MPAS-CO 2 adjoint model with its counterpart from CarbonTracker–Lagrange (CT-L) reveals a spatial agreement but notable magnitude differences. These differences, particularly evident for X CO 2 , might be attributed to the two model systems' differences in the simulation configuration, spatial resolution, and treatment of vertical mixing processes. Moreover, this comparison highlights the substantial loss of information in the atmospheric CO 2 observations due to CT-L's spatial domain limitation. Furthermore, the adjoint sensitivity analysis demonstrates that the sensitivities to both surface flux and initial CO 2 conditions spread out throughout the entire Northern Hemisphere within a month. MPAS-CO 2 forward, TL, and adjoint models stand out for their calculation efficiency and variable-resolution capability, making them competitive in computational cost. In conclusion, the successful development of the MPAS-CO 2 TL and adjoint models, and their integration into the MPAS-CO 2 system, establish the possibility of using MPAS's unique features in atmospheric CO 2 transport sensitivity studies and in inverse modeling with advanced methods such as variational data assimilation.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Combined Mesonet and Tracker

Title: Combined Mesonet and Trackers (UNL Mobile Mesonets) Authors University of Nebraska PI: Adam Houston, UNL Professor (ahouston2@unl.edu) Mailing Address: 126 Bessey Hall P.O. Box 880340 Lincoln, NE 68588-0340 CoMeT Overview The University of Nebraska-Lincoln operates three Combined Mesonet and Tracker (CoMeTs). CoMeTs are Ford Explorers (model years 2015, 2017, and 2019) with forward-mounted suites of meteorological sensors and dual moonroofs, combining the capability of a mobile mesonet to collect near-surface observations with the capability of an unmanned aircraft systems (UAS) tracker vehicle, which enables an observer in the second row of seats to see the aircraft and maintain compliance with Federal Aviation Administration policies on UAS operation. The CoMeTs collect observations of slow temperature and humidity at ~2 m above ground level (AGL) using a Vaisala HMP155A, fast temperature at ~2 m AGL using a Campbell Scientific 109SS-L thermistor, pressure at ~2.5 m AGL using a Vaisala PTB210, wind speed and direction at ~3.25 m AGL using an R.M. Young 05103 propeller anemometer, and vehicle heading using a KVH Industries C-100 fluxgate compass (Barbieri et al. 2019). The HMP155A and 109SS-L thermistor are shielded and aspirated within a U-tube (Waugh and Frederickson 2010; Houston et al. 2016). This list of sensors is also included in the CoMeT data file metadata. Manufacturer specifications for these instruments are given in Table 1 of Hanft and Houston (2018). The reported measured quantities are summarized below.CoMeT-3 was funded through an equipment allocation included in the NSF TORUS award (AGS-1824649). Instrument Description The specific sensors included on each CoMeT are summarized in the table at the end of this section. In general each CoMeT collects observations of slow temperature and humidity at ~2 m above ground level (AGL) using a Vaisala HMP155, fast temperature at ~2 m AGL using a Campbell Scientific 109SS-L thermistor, pressure at ~2.5 m AGL using a Vaisala PTB210 barometer with a Gill pressure port, wind speed and direction at ~3.25 m AGL using an R.M. Young 05103 propeller anemometer, position using a Garmin 19x HVS receiver, and vehicle heading using a KVH Industries C-100 fluxgate compass. The HMP155 and 109SS are shielded and aspirated within a U-tube (Waugh and Frederickson 2010; Houston et al. 2016). Fast temperature and corrected RH measurements (using sensors housed within the U-tube) have a time constant of 10-12 s based on data collected across a temperature and RH shock during the CLOUD-MAP 2017 calibration/validation tests on June 26, 2017. Vehicle speed was < 10 kts for this test. CoMeT-1 CoMeT-2 CoMeT-3 Slow Temperature Slow RH Vaisala HMP155A-L20-PT Part #: 22280-7 Vaisala HMP155E Part #: E1AA11A0B1A1A0A Vaisala HMP155E Part #: E1AA11A0B1A1A0A Fast temperature Campbell Scientific 109SS-L20-PT Part #: 21448-3 Campbell Scientific 109SS-L12-PW Part #: 21448-109 Campbell Scientific 109SS-L12-PW Part #: 21448-150 Pressure Vaisala PTB-210 Part #: A1A1B Gill Pressure Port Part #: 61002 Vaisala PTB-210 Part #: A1A1B Gill Pressure Port Part #: 61002 Vaisala PTB-210 Part #: A1A1B Gill Pressure Port Part #: 61002 Wind RM Young 05103-L20-PT Part #: 18435-310 RM Young 05103-L20-PW Part #: 18435-244 RM Young 05103-L20-PW Part #: 18435-244 GPS Garmin GPS 19x HVS (NMEA 0183) Part #: 010-01010-00 Garmin GPS 19x HVS (NMEA 0183) Part #: 010-01010-00 Garmin GPS 19x HVS (NMEA 0183) Part #: 010-01010-00 Compass KVH C-100 Part #: 01-0177-15 KVH C-100 Part #: 01-0177-15 KVH C-100 Part #: 01-0177-15 Logger Campbell Scientific CR6-NA-XT-SW Part #: 28385-9 Campbell Scientific CR6-WIFI-XT-SW Part #: 28385-6 Campbell Scientific CR6-WIFI-XT-SW Part #: 28385-6 Data Collection and Real-Time Processing The reported measured quantities are summarized in the table below. Quantity Units Source Epoch time Seconds GPS Latitude and longitude Degrees GPS Altitude m GPS Pressure hPa PTB210 Temperature (fast) deg C 109SS-L Temperature (slow) deg C HMP155 RH (slow) % HMP155 Vehicle speed m/s GPS Vehicle heading deg C-100 and GPS In addition to the measured variables, several derived variables are calculated. Corrected/fast relative humidity (%) Relative humidity is adjusted to the fast temperature following Richardson et al. (1998) and Houston et al. (2016). Water vapor mixing ratio (g/kg) Dew point temperature (°C) Potential temperature (Kelvin) Virtual potential temperature (Kelvin) Equivalent potential temperature (Kelvin) Regular intercomparisons between all three CoMeTs were performed during TORUS 2019. Comparisons were also conducted between CoMeT-1 and CoMeT-2 during LAPSE-RATE (2018) on 14 July. In these intercomparisons, the vehicles were parked adjacent to each other aligned perpendicular to (and facing into) the wind. To minimize engine heating effects, intercomparisons were only conducted when the wind speed was >10 kts. Data Format Original data files for each deployment are saved as text files and then converted to NetCDF. NetCDF versions have units that are CF compliant and may not match the original units in the txt files. The naming convention for the NetCDF files is as follows: UNL.CoMeT3.{deployment date YYYYMMDD}.{start time of observation collection in UTC HHMM}.L2.{post-processing codes}.cdf example: UNL.CoMeT3.20190627.1931.L2.g1.f1.cdf Post-processing codes are included to track modifications to the raw data. These codes are closely connected to error flags associated with each record. Each letter corresponds to a particular instrument: g: GPS p: Barometer tf: Fast temperature ts: Slow temperature rh: Relative humidity f: Compass w: Wind monitor a: All instruments Each number corresponds to a particular post-processing action described more below. Measured and derived variables are included in the following table. Variable Heading Standard Name Units time Time seconds since 00:00:00, 01-01-1970 Alt Altitude meters lat Latitude degrees north lon Longitude degrees east fast_temp Air Temperature kelvin slow_temp Air Temperature kelvin pressure Air Pressure pascals logger_RH Relative Humidity percent calc_corr_RH Relative Humidity percent wind_speed Wind Speed meters per second wind_dir Wind From Direction degrees vehicle_dir Vehicle Direction degrees dewpoint Dew Point Temperature kelvin mixing_ratio Humidity Mixing Ratio g/g theta Air Potential Temperature kelvin theta_v Virtual Potential Temperature kelvin theta_e Equivalent Potential Temperature Kelvin error_flag The error_flag variable is a string that matches the post-processing codes listed above. All instruments will have an associated code, but will have a “0” if the datum is unchanged from the initial processed value. Error Codes The following table summarizes the error codes for data collected before 2020: Error Code Relevant CoMeT Description g1 1,2,3 Exact correction. GPS position and time reprocessed from raw data g2 1 As far as we can tell this is an exact correction to an error in the GPS time. During the correct time periods the time suddenly went backwards ~250s and stayed at this offset for 750s when it corrected itself. The offset was applied to the “time warp” period. p1 2 Approximate correction. Hole in the pressure tube connecting the pressure port to the barometer. Resulted in erroneously low air pressure measurements when the vehicle was in motion. Derived variables recalculated (dew point temperature [e depends on qv and p], water vapor mixing ratio, potential temperature, virtual potential temperature, equivalent potential temperature) a1 3 Exact correction. Missing data reprocessed from raw data a2 1 Bug fix to bias correction for ts1, ts2, and rh1: water vapor mixing ratio was off by a factor of 10 and virtual potential temperature was wrong because of this. f1 3 No correction, missing data. Fluxgate compass inoperable. Wind speed and direction calculated using GPS-derived vehicle heading instead. rh1 1 Approximate correction. Constant bias of +1.7% removed from relative humidity. Derived variables recalculated (corrected/fast relative humidity, dew point temperature, water vapor mixing ratio, virtual potential temperature, equivalent potential temperature) ts1 1 Approximate correction. Constant bias of +0.6 K removed from slow temperature. Derived variables recalculated (corrected/fast relative humidity, dew point temperature, water vapor mixing ratio, virtual potential temperature, equivalent potential temperature) ts2 1 Approximate correction. Constant bias of +1.0 K removed. Derived variables recalculated (corrected/fast relative humidity, dew point temperature, water vapor mixing ratio, virtual potential temperature, equivalent potential temperature) References Bolton, D., 1980: The Computation of Equivalent Potential Temperature. Mon. Wea. Rev., 108, 1046–1053, https://doi.org/10.1175/1520-0493(1980)108<1046:TCOEPT>2.0.CO;2. Hanft, W., and A. L. Houston, 2018: An Observational and Modeling Study of Mesoscale Air Masses with High Theta-E. Mon. Wea. Rev., 146, 2503–2524, https://doi.org/10.1175/MWR-D-17-0389.1.Wexler Houston, A. L., R. J. Laurence III, T. W. Nichols, S. Waugh, B. Argrow, and C. L. Ziegler, 2016: Intercomparison of unmanned aircraft-borne and mobile mesonet atmospheric sensors. Journal of Atmospheric and Oceanic Technology. 33, 1569-1582, doi: 10.1175/JTECH-D-15-0178.1. Lowe, P. R., 1977: An Approximating Polynomial for the Computation of Saturation Vapor Pressure. J. Applied Meteorology, 16, 100–103. Richardson, S. J., S. E. Frederickson, F. V. Brock, and J. A. Brotzge, 1998: Combination temperature and relative humidity probes: Avoiding large air temperature errors and associated relative humidity errors. Preprints, 10th Symp. On Meteorological Observations and Instrumentation, Phoenix, AZ, Amer. Meteor. Soc., 278–283. Waugh, S., and S. E. Frederickson, 2010: An improved aspirated temperature system for mobile meteorological observations, especially in severe weather. 25th Conf. on Severe Local Storms, Denver, CO, Amer. Meteor. Soc., P5.2. [Available online at https://ams.confex.com/ams/25SLS/techprogram/paper_176205.htm.]

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

A geospatial risk analysis graphical user interface for identifying hazardous chemical emission sources

Background: Performing back trajectory and forward trajectory using the Hybrid Single-Particle Lagrangian Integrated Trajectory Model (HYSPLIT) is a reliable approach for assessing particle transport after release among mid-field atmospheric models. HYSPLIT has an externally facing online interface that allows non-expert users to run the model trajectories without requiring extensive training or programming. However, the existing HYSPLIT interface is limited if simulations have a large amount of meteorological data and timesteps that are not coincident. The objective of this study is to design and develop a more robust tool to rapidly evaluate hazard transport conditions and to perform risk analysis, while still maintaining an intuitive and user-friendly interface. Methods: HYSPLIT calculates forward and backward trajectories of particles based on wind speed, wind direction, and the corresponding location, timestamp, and Pasquill stability classes of the regions of the atmosphere in terms of the wind speed, the amount of solar radiation, and the fractional cloud cover. The computed particle transport trajectories, combined with the online Proton Transfer Reaction-Mass Spectrometry (PTR-MS) data (https://figshare.com/articles/dataset/ARL_Data_from_PROS_station_at_Hanford_site/19993964), can be used to identify and quantify the sources and affected area of the hazardous chemicals’ emission using the potential source distribution function (PSDF). PSDF is an improved statistical function based on the well-known potential source contribution function (PSCF) in establishing the air pollutant source and receptor relationship. Performing this analysis requires a range of meteorological and pollutant concentration measurements to be statistically meaningful. The existing HYSPLIT graphical user interface (GUI) does not easily permit computations of trajectories of a dataset of meteorological data in high temporal frequency. To improve the performance of HYSPLIT computations from a large dataset and enhance risk analysis of the accidental release of material at risk, a geospatial risk analysis tool (GRAT-GUI) is created to allow large data sets to be processed instantaneously and to provide ease of visualization. Results: The GRAT-GUI is a native desktop-based application and can be run in any Windows 10 system without any internet access requirements, thus providing a secure way to process large meteorological datasets even on a standalone computer. GRAT-GUI has features to import, integrate, and convert meteorological data with various formats for hazardous chemical emission source identification and risk analysis as a self-explanatory user interface. The tool is available at https://figshare.com/articles/software/GRAT/19426742.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Metamaterial-Inspired Vacuum Electron Devices and Accelerators

Metamaterials (MTMs) are structured materials with subwavelength features that can be engineered to have some unique properties not found in nature, such as negative refractive index, reversed Doppler effect, and reversed Cherenkov radiation. Based on these novel MTMs, several important research groups have made great attempts to develop novel MTM-inspired vacuum electron devices (VEDs) and accelerators. Just as solid-state power devices are innovated by incessant emerging of new semiconductor materials, VEDs can also be inspired by MTMs to have very remarkable advantages, such as smaller size, higher power, higher efficiency, and/or larger gain relative to conventional VEDs, such as traveling-wave tubes, backward-wave oscillators, and klystrons. Similarly, relative to conventional accelerators, MTM-inspired accelerators have obvious advantages, such as smaller size and higher accelerating gradient. Furthermore, MTM-inspired devices have promising applications in areas such as radar, communication, electronic warfare, microwave heating, and imaging.

43 PARTICLE ACCELERATORS↗

DOE OSTI2020

Compressive sensing optical design and simulation tool

A system or method for designing optical systems includes generating a mathematical model of an imaging system; tracing a plurality of light rays from a source through the model of the imaging system in a forward direction and a backward direction incident upon an arbitrarily defined detector location; calculating the total light rays incident upon the arbitrarily defined detector location in response to the plurality of light rays; tracing a series of calibration images first through the computational or compressive imager mathematical model; measuring an output set of parameters based on a plurality of detected rays calibration images at the defined detector location; storing the output set of parameters in software, and determining an estimated system transfer function.

42 ENGINEERING↗

DOE OSTI2020

Simulation and Development of the Radial Time Projection Chamber For the Bonus12 Experiment in CLAS12

Knowledge of the structure of nucleons (i.e. protons and neutrons) is a central topic of interest to nuclear/particle physicists. Much more is known about the structure of the proton than the neutron due to the lack of high-density free neutron targets. The Barely Off-shell Nucleon Structure experiment (BONuS12) at Jefferson Lab (JLab) is a second generation experiment upgraded/optimized to advance our knowledge of the neutron's structure using the deep-inelastic scattering of electrons off deuterium. Typically, since deuterium is a nuclear target, corrections for off-shell and nuclear binding effects must be taken into account in order to extract results on the neutron. These corrections are model-dependent and therefore have limited our success in extracting neutron information using deuterium targets. In the BONuS12 experiment, 10.6 GeV electrons are scattered off of a deuterium target. By detecting the low momentum spectator proton at backward angles, the uncertainty due to final state interactions is minimized. The goal of the experiment is to measure the ratio of the neutron to proton structure functions ($F^n_2/F^p_2$) as the Bjroken scaling variable x approaches 1. The newly designed Radial Time Projection Chamber (RTPC) for BONuS12 detects the spectator proton in coincidence with the scattered electron, which is detected in the CEBAF Large Acceptance Spectrometer (CLAS12). This work presents the simulation and development of the new BONuS12 RTPC. The design, construction, and testing of the Drift-gas Monitoring Sysytem (DMS) for the BONuS12 experiment is also described. The results of the DMS operation as well as the first preliminary data from the BONuS12 experimental run are given. Because the BONuS12 data analysis depends on CLAS12 working effectively, an effort to verify the CLAS12 operation with the extraction of the inclusive deep inelastic cross section from the first experiment in CLAS12 (Run Group A) will be presented.

Dzbenski, Nathan↗

DOE OSTI2022

Neumann Series in MGS-GMRES and Inner-Outer Iterations: Preprint

A low-synchronization MGS-GMRES Krylov solver employing a truncated Neumann series for the inverse compact WY MGS correction matrix T is presented. A corollary to the backward stability result of Paige et al. [1] establishes that T = I - Lk is sufficient for convergence of GMRES when kLkp F = O("p)_p F (B), where the strictly lower triangular matrix L is defined by the inner products of Krylov vectors V T 1:k-2 vk-1. The preconditioner is the classical Ruge-Stuben AMG algorithm with compatible relaxation and inner-outer Gauss-Seidel smoother. This smoother may also be expressed as a truncated Neumann series. Drop tolerances are applied to the lower triangular matrices arising in the smoother in order to reduce the number of non-zeros and accelerate the time to solution. The number of small matrix elements are found to increase from fine to coarse levels and thus the effciency gains are greater for large problems with many levels in the V -cycle. The solver is applied to the pressure continuity equation for the incompressible Navier-Stokes equations. Unlike the inner-outer iteration, the solver convergence rate with the standard Gauss-Seidel smoother deteriorates with dropping. The solver compute time is reduced by up to 50% without a change in the convergence rate.

Gauss-Seidel smoother↗

DOE OSTI2021

IEEE P2988: Recommended Practice of Use and Functions of Virtual Synchronous Machines

Power systems are going through a paradigm shift from centralized generation to distributed generation. This brings unprecedented challenges to power systems stability, reliability, and security. Operating power electronic converters as Virtual Synchronous Machines (VSMs) can bring backward compatibility to power systems, unify and harmonize the integration and interaction of future power systems players, and improve power system stability, reliability, and security. This is vital for the large-scale adoption of distributed energy resources, the advancement of sustainability, and the development of a low-carbon economy. This IEEE 2988 Recommended Practice will define the fundamental principles, mandatory functions, and optional functions of a VSM, i.e., a power electronic converter that is operated to behave like conventional synchronous machines. It will not describe or specify power semiconductor devices, hardware topologies, or micro-controllers that are used to build a VSM. This talk will discuss the activities and development plan for IEEE 2988 with the aim to engaging broad stakeholders and widening participation.

Zhong, Qing-Chang↗