Search NASA⌕ Search

SEARCH · Search NASA

Results for “Lagrangian integration”

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

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

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Conservative High-Order Time Integration for Lagrangian Hydrodynamics

In tis study, we develop novel time integration methods for the compressible Euler equations in the Lagrangian frame that are of arbitrary high order and exactly preserve the mass, momentum, and total energy of the system. The equations are considered in nonconservative form, that is, common for staggered grid hydrodynamics (SGH) methods; namely, the evolved quantities are mass, momentum, and internal energy. A general family of time integration schemes is formulated, and practical pairs for orders three and four are derived. Numerical results on standard hydrodynamics benchmarks confirm the high-order convergence on smooth problems and the exact numerical preservation of all physically conserved quantities.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Numerical integrators for Lagrangian oceanography

Abstract. A common task in Lagrangian oceanography is to calculate a large number of drifter trajectories from a velocity field precalculated with an ocean model. Mathematically, this is simply numerical integration of an ordinary differential equation (ODE), for which a wide range of different methods exist. However, the discrete nature of the modelled ocean currents requires interpolation of the velocity field in both space and time, and the choice of interpolation scheme has implications for the accuracy and efficiency of the different numerical ODE methods. We investigate trajectory calculation in modelled ocean currents with 800 m, 4 km, and 20 km horizontal resolution, in combination with linear, cubic and quintic spline interpolation. We use fixed-step Runge–Kutta integrators of orders 1–4, as well as three variable-step Runge–Kutta methods (Bogacki–Shampine 3(2), Dormand–Prince 5(4) and 8(7)). Additionally, we design and test modified special-purpose variants of the three variable-step integrators, which are better able to handle discontinuous derivatives in an interpolated velocity field. Our results show that the optimal choice of ODE integrator depends on the resolution of the ocean model, the degree of interpolation, and the desired accuracy. For cubic interpolation, the commonly used Dormand–Prince 5(4) is rarely the most efficient choice. We find that in many cases, our special-purpose integrators can improve accuracy by many orders of magnitude over their standard counterparts, with no increase in computational effort. Equivalently, the special-purpose integrators can provide the same accuracy as standard methods at a reduced computational cost. The best results are seen for coarser resolutions (4 and 20 km), thus the special-purpose integrators are particularly advantageous for research using regional to global ocean models to compute large numbers of trajectories. Our results are also applicable to trajectory computations on data from atmospheric models.

58 GEOSCIENCES↗
DOE OSTI2025

Back Trajectories during SAIL

This data set contains air mass back trajectories generated during the Surface Atmosphere Integrated Field Laboratory (SAIL) campaign using the Hybrid Single-Particle Lagrangian Integrated Trajectory model (HYSPLIT, Stein et al. 2015, Rolph et al. 2017). For each hour from January 1, 2021 to June 30, 2023, a 96-hour back trajectory was initiated at the SAIL sampling site, starting 100 m above ground level. The calculations used the GDAS meteorological data set with model vertical velocity. The files within this dataset show location and meteorological parameters of the back trajectory analysis. There are 21 columns within each data file. The first 2 columns pertain to back trajectory parameters: the trajectory number (which is always 1) and the meteorological grid number. The next 9 columns pertain to spatiotemporal information: the year, month, day, hour (0-23), and minute of the trajectory, the previous time of the trajectory (in negative hours), and the airmass location (latitude, longitude, and height in m AGL). The rest of the parameters within each file pertains to surface and meteorological characteristics of the airmass: pressure (in hPa), potential temperature (in K), temperature (in K), precipitation rate (in mm/hr), mixing depth (in m), relative humidity (in % relative to liquid water), specific humidity (in g/kg), H2O mixing ratio (in g/kg), underlying terrain height (in meters above sea level), and incoming radiation (in W/m^2).

back trajectory↗
DOE OSTI2024

Estimated impacts of forest restoration scenarios on smoke exposures among outdoor agricultural workers in California

Abstract As wildfires continue to worsen across western United States, forest managers are increasingly employing prescribed burns as a way to reduce excess fuels and future wildfire risk. While the ecological benefits of these fuel treatments are clear, little is known about the smoke exposure tradeoffs of using prescribed burns to mitigate wildfires, particularly among at-risk populations. Outdoor agricultural workers are a population at increased risk of smoke exposure because of their time spent outside and the physical demands of their work. Here, we assess the smoke exposure impacts among outdoor agricultural workers resulting from the implementation of six forest management scenarios proposed for a landscape in the Central Sierra, California. We leverage emissions estimates from LANDIS-II to model daily PM 2.5 concentrations with the Hybrid Single Particle Lagrangian Integrated Trajectory Model (HYSPLIT) and link those to agricultural employment data from the Bureau of Labor Statistics. We find a u-shaped result, in that moderate amounts of prescribed burning result in the greatest reduction in total smoke exposure among outdoor agricultural workers, particularly during months of peak agricultural activity due to wildfire-specific smoke reductions. The reduction in total smoke exposure, relative to scenarios with the least amount of management, decreases as more prescribed burning is applied to the landscape due to the contributions of the fuel treatments themselves to overall smoke burden. The results of this analysis may contribute to preparedness efforts aimed at reducing smoke exposures among outdoor agricultural workers, while also informing forest management planning for this specific landscape.

Environmental Sciences & Ecology↗
DOE OSTI2025

Forecasting for ESCAPE: A Multi-Institution Hybrid Forecasting and Nowcasting Operation for Sea-Breeze Convection Supporting a Ground-Based and Airborne Field Campaign

The Experiment of Sea-Breeze Convection, Aerosols, Precipitation and Environment (ESCAPE) field project deployed two aircraft and ground-based assets in the vicinity of Houston, Texas, between 27 May and 2 July 2022, examining how meteorological conditions, dynamics, and aerosols control the initiation, early growth stage, and evolution of coastal convective clouds. To ensure that airborne- and ground-based assets were deployed appropriately, a forecasting and nowcasting team was formed. Daily forecasts guided real-time decision-making by assessing synoptic weather conditions, environmental aerosol, and a variety of atmospheric modeling data to assign a probability for meeting specific ESCAPE campaign objectives. During the research flights, a small team of forecasters provided “nowcasting” support by analyzing radar, satellite, and new model data in real time. The nowcasting team proved invaluable to the campaign operation, as sometimes changing environmental conditions affected, for example, the timing of convective initiation. In addition to the success of the forecasting and nowcasting teams, the ESCAPE campaign offered a unique “testbed” opportunity where in-person and virtual support both contributed to campaign objectives. The forecasting and nowcasting teams were each composed of new and experienced forecasters alike, where new forecasters were given invaluable experience that would otherwise be difficult to attain. Both teams received training on forecast models, map analysis, Hybrid Single-Particle Lagrangian Integrated Trajectory model (HYSPLIT), and thermodynamic sounding analysis before the beginning of the campaign. In this article, the ESCAPE forecasting and nowcasting teams reflect on these experiences, providing potentially useful advice for future field campaigns requiring forecasting and nowcasting support in a hybrid virtual/in-person framework.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Initial Atmospheric Transport of Particles

A new version of the Initial Atmospheric Transport (IAT) model has been developed. This report fully characterizes the new model. IAT is designed to model heat release into the atmosphere that contains contaminated particles. The heat release forms a buoyant puff or plume. Buoyancy is quickly dissipated in 10's to 100's of seconds through drag and entrainment of ambient air. The final location of contaminant particles in the atmosphere after dissipation of abnormal heat is an important input to long term transportation codes such as the National Oceanic and Atmospheric Administration's (NOAA) Hybrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) code. Ultimately the settling of contaminants after long term transport provides a basis for estimating health consequences. In comparison to the old version, the new model includes explicit simulation of particles, newly formulated equations of motion for a rising puff, and a new plume model. The plume model releases many instances of a single puff along a single trajectory in a way that balances energy released over a longer period than would be appropriate for a single puff. The puff model has additions for vorticity, turbulence, virtual mass, drag, and entrainment. The particle model leads to direct interfacing between IAT and the long-term transport code HYSPLIT. Previously, particle locations were given a pattern based on IAT's puff trajectory with no particle simulation. Now, particles are explicitly simulated, and 3-D spatial distributions of particles are output based on final positions. Ballistic particles are thrown outward at high velocity from explosive events and tend to hit the ground before the end of the simulation. Non-ballistic particles tend to get pulled into the rising puff s ring of vorticity with some escaping to ambient conditions as buoyant energy dissipates. A new model verification and validation effort was made for the new plume model. The Wallops Island Solid Rocket Propellant (SRP) fire plume was extracted to a 3-D transient dataset. After calibration, the model performs better than previous single puff attempts. Even so, error levels in the process used to estimate location and ambient weather conditions approach 30%. This leads to the conclusion that unplanned photogrammetry evaluations of the IAT model are useful but not sufficiently accurate to quantify model parameters. Like previous attempts, the vertical motion entrainment parameters come out much too high in comparison to values for other studies in the literature. It is suspected that a new model of shear entrainment would alleviate this problem. Unfortunately, IAT already has more calibration parameters than the data comparisons thus far can resolve with statistical certainty. This underscores the need to verify and validate the IAT model using detailed CFD studies that allow ideal statistical comparison with the IAT model for hundreds of cases. Even so, validation by comparison to carefully planned experiments (i.e. with acceptable error levels) is also necessary because the dynamics of entrainment are at a sub-grid level for CFD. Though better validation and verification are most important, many new enhancements are still envisioned for future IAT development. These include particle rain out due to condensation and agglomeration of water, direct connection to the North American Mesoscale (NAM) weather dataset, dynamic weather conditions, including ambient particle entrainment into the puff, and detrainment at the puff to ambient mixing interface. In general, including detailed aerosol dynamics is the greatest coming challenge which will require an entire new sub-model in IAT. This is very important if contaminant particles change in such a way that they increase or decrease in their potential harm to the environment.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Connecting the Radiative Influences of Aerosol upon the Mass Flux Profiles of Shallow Cumuli across the Southeast Atlantic Ocean Basin and its Boundaries (Final Report)

The Atlantic Ocean covers approximately 25% of Earth’s surface and the atmosphere above it is home to a complex array of clouds and aerosols that have important influences on regional and global weather and climate. These influences must be accurately depicted in short range, medium range, seasonal, and climate forecast models. Conditions over the tropical Atlantic are particularly complex due to continental scale plumes of dust from the Sahara Desert and smoke from agricultural burning in Africa that drift across the Atlantic Ocean basin toward the Americas. These plumes are often found meandering in the lower atmosphere above shallow tropical clouds that form above the ocean surface, presumably mingling with these clouds on occasion due to convective mixing processes. Elevated dust and smoke particles absorb incoming sunlight and substantially warm the marine atmosphere in the layer in which they are present. This warming may alter the thermal stability of the marine atmosphere and may throttle or enhance the development of clouds, change their internal structure and the rate at which they precipitate. Alternatively, it may isolate the lower atmosphere from drier layers above enabling water vapor to accumulate near the ocean surface potentially leading to the development of deeper convection. Our study was organized around the principal concept of determining the impact of African smoke and dust plumes upon cloud development at ASI and understanding how a popular shallow convection parameterization used in models responds to the presence of this aerosol. We analyzed observations collected during the US Department of Energy (DOE) Layered Atlantic Smoke Interactions with Clouds (LASIC) campaign using the Atmospheric Radiation Measurement Program’s Mobile Facility #1 (AMF-1), which was deployed to Ascension Island (ASI) in Southeastern Atlantic for a one-year period. Ascension Island is immediately downwind from an African source of these plumes, but far enough removed to enable the lower atmosphere to have reacted to their presence. A main goal of our study was to compute radiative heating rate profiles over Ascension Island and along the trajectory from the biomass burning regions along coastal Africa to Ascension Island. To compute the radiative heating rate due to aerosols and clouds, we employed observed profiles of temperature, humidity, and clouds from LASIC alongside aerosol optical properties from the Modern-era Retrospective analysis for Research and Applications Version 2 (MERRA-2), as input for the Rapid Radiation Transfer Model (RRTM). Radiative heating was also assessed across the southeast Atlantic Ocean using an ensemble of back trajectories from the Hybrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) model. We were successful in this effort and the resulting publication is already being cited despite its relatively short lifetime in the literature. The second part of our study involved the process-level representation of the clouds observed over the Southeastern Atlantic in models. Our initial task, upon which the balance of this portion of the study depended, was to evaluate the representativeness of the clouds observed by AMF-1. What orographic influence did Ascension Island have on the measured cloud structure? We set out to answer this question by trying to separate observations that were clearly indicative of orographic forcing from those that were representative of open ocean. The siting of AMF-1, which was 340-m above sea-level on the slope of a steep escarpment, proved demonstrably problematic for cloud process measurements. Despite considerable effort using multiple approaches including artificial intelligence (AI), we were unable to successfully compensate for the orographic forcing present at the AMF-1 site on ASI and, hence, produce process-level summaries of the convective mass flux representative of open ocean over the Southeastern Atlantic. Even though the project has officially ended, we are making a final attempt using two new types of clustering algorithm (i.e., AI) on the recommendation of a former student, who is an expert in this area. Since we have only recently begun to test these new AI schemes, the recommendations from the second part of our project outlined below are based on our experience at the time of this report.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

User guide for the Geospatial Risk Analysis Tool (V.1.0.1)

Developed by the Pacific Northwest National Laboratory (PNNL), the Geospatial Risk Analysis Tool (GRAT) links to Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT), a complete system for computing simple air parcel trajectories, as well as complex transport, dispersion, chemical transformation, and deposition simulations. The HYSPLIT model was developed by the National Oceanic and Atmospheric Administration (NOAA) Air Resources Laboratory and the Australian Bureau of Meteorology Research Center in 1998. GRAT is used to batch process meteorological data into ARL format. While HYSPLIT GUI limits users to convert meteorological data up to 6 date time points per operation, GRAT enables users to convert a significantly large amount of data (For example, 10 years of meteorological data measured with an interval of 5 minutes) within one operation. Incorporated with the Potential source distribution function (PSDF), GRAT can also be used to determine the areas influenced by the emission of hazardous chemicals.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

An LSTM Deep Learning Network for ¿Background Radiation Prediction

Determination of appropriate background radiation is important in any measurement application. Environmental radiation monitors and monitors used to assess dose to individuals outside of controlled areas are particularly susceptible to changes in readings due to fluctuations in the environmental conditions. These fluctuations (e.g. radon progeny concentrations) lead to changes in the detector response in the actual radiation environment, and they need to be taken into account when extracting the net operational doses. Work has been ongoing to apply advances in Deep Learning and Artificial Intelligence to account for changes in detector responses based on environmental parameters; in particular, a Long-Short Term Memory (LSTM) Deep Learning architecture has been utilized to incorporate time-series data into a prediction model. In this work, we present the current status of the project to predict radiation measurements based on meteorological conditions and air packet trajectories extracted using the National Oceanic and Atmospheric Administration's (NOAA) Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT4) model.

Degtiarenko, Pavel↗
DOE OSTI2024

Utilizing HYSPLIT for Emergency Response Modeling at SRS

At SRS, emergency responders use a variety of tools to detect, track, and mitigate hazardous material releases into the atmosphere. Two models currently used at SRS are Puff-Plume and the Lagrangian Particle Dispersion Model (LPDM), a Gaussian and Lagrangian model, respectively. A decision has been made to replace LPDM with the more widely-supported Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) model for evaluating inhalation and ingestion doses following a release. HYSPLIT is designed to compute complex dispersion and deposition simulation To achieve the implementation of HYSPLIT, we have developed a preliminary UI framework that will allow HYSPLIT to be run on ATG computers without the need for active network connections, thus avoiding the loss of capabilities in the event of a network outage during an emergency.

Earley, Ian↗
DOE OSTI2026

ARM Trajectories Data Set Value-Added Product Report

The U.S. Department of Energy Atmospheric Radiation Measurement (ARM) user facility’s ARM Trajectories Data Set (ARMTRAJ) Value-Added Product (VAP) provides trajectory data sets initialized at ARM deployment coordinates and configured using ARM data sets. The four trajectory data sets support aerosol, cloud, and planetary boundary-layer research. Trajectory calculations use the Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) model informed by the European Centre for Medium-Range Weather Forecasts (ECMWF) fifth-generation atmospheric reanalysis (ERA5) data set at its highest spatial resolution (~31 km). HYSPLIT also runs at multiple initial starting locations surrounding ARM deployments (in latitude/longitude and/or vertical coordinates), facilitating an ensemble for each sample in the data sets. The ensemble mean and variability reported in ARMTRAJ improve the fidelity and provide uncertainty estimates of trajectory coordinates, thermodynamic properties, and other output fields.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

ARM Trajectories Data Set Value-Added Product Report

The U.S. Department of Energy Atmospheric Radiation Measurement (ARM) user facility’s ARM Trajectories Data Set (ARMTRAJ) Value-Added Product (VAP) provides trajectory data sets initialized at ARM deployment coordinates and configured using ARM data sets. The six trajectory data sets support aerosol, cloud, planetary boundary layer, and related research (aerosol-cloud interactions, etc.), as well as studies using ARM Aerial Facility (AAF) and tethered balloon system (TBS) measurements. Trajectory calculations use the Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) model informed by the European Centre for Medium-Range Weather Forecasts (ECMWF) fifth-generation atmospheric reanalysis (ERA5) data set at its highest spatial resolution (~31 km). HYSPLIT also runs at multiple initial starting locations surrounding ARM deployments (in latitude/longitude and/or vertical coordinates), facilitating an ensemble for each sample in the data sets. The ensemble mean and variability reported in ARMTRAJ improve the fidelity and provide uncertainty estimates of trajectory coordinates, thermodynamic properties, and other output fields.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Beryllium Isotopes in Rainfall for Stratospheric-Tropospheric Exchange Dynamics Experiment (BIRDIE) (Field Campaign Report)

The Beryllium Isotopes in Rainfall for Stratospheric-Tropospheric Exchange Dynamics Experiment or BIRDIE was a precipitation collection campaign designed to measure cosmogenic isotopes of beryllium in rainfall and other forms of precipitation across the U.S. Department of Energy Atmospheric Radiation Measurement (ARM) user facility’s Southern Great Plains (SGP) atmospheric observatory. Precipitation samples were collected over three weeks (May 21-June 10, 2022) at 11 SGP sites across a 120-km 2 domain in northern and central Oklahoma. The field campaign targeted summertime convection rain events. We aimed to capture any spatial and/or temporal variability in the concentrations of beryllium-10 (10Be) and beryllium-7 (7Be) across the ARM SGP site during these rainfall events. The ratio of 10Be to 7Be (10/7Be) has been shown to be a powerful tracer for transport processes in the stratosphere-troposphere exchange (STE) system (e.g., Koch and Rind 1998). Any observed variability is analyzed to see if it could be linked to STE processes, including overshooting of deep convection clouds that penetrate the lower stratosphere. Six wet weather system events were collected within our three-week span. Five out of six of these rainfall events are classified as moderately to strongly convective. Those events are studied in more detail. The five systems also follow each other over a short period of 11 days. In our rainfall samples we observed spatial and event variability in the concentrations of beryllium-10 and -7 that resulted in variability of 10/7Be ratios. Each weather event appeared to exhibit a unique 10/7Be pattern across the domain that was not directly linked to storm total precipitation amount. To better understand the causes for this spatial variability, cloud characteristics in the domain were looked at to assess cloud type, cloud height, and depth information. We also ran Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) model back trajectories to identify possible air mass and moisture origins across our events. Detailed analysis is ongoing. The purpose of this technical report is to provide background on using cosmogenic isotopes as atmospheric tracers, describe the experimental design, and present the beryllium concentration and ratio data. Lastly, we note that BIRDIE coincided with two other deep convection field campaigns in the area: The National Aeronautics and Space Administration (NASA)’s Dynamics and Chemistry of the Summer Stratosphere (DCOTSS) experiment (https://dcotss.org/) centered in Salina, Kansas, and the Targeted Observation by Radars and (uncrewed aerial systems) UAS of Supercells (TORUS) experiment (https://catalog.eol.ucar.edu/torus_2022) across the U.S. Central Plains. Additionally, a Lawrence Livermore National Laboratory ( LLNL) team at DCOTSS attempted to measure 10Be/7Be ratios on aerosols collected in situ in the lower and mid-stratosphere over Kansas using high-altitude balloons. Precipitation samples were also collected on the ground in Salina, Kansas.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Radiative heating rate profiles over the southeast Atlantic Ocean during the 2016 and 2017 biomass burning seasons

Marine boundary layer clouds, including the transition from stratocumulus to cumulus, are poorly represented in numerical weather prediction and general circulation models. Further uncertainties in the cloud structure arise in the presence of biomass burning carbonaceous aerosol, as is the case over the southeast Atlantic Ocean, where biomass burning aerosol is transported from the African continent. As the aerosol plume progresses across the southeast Atlantic Ocean, radiative heating within the aerosol layer has the potential to alter the thermodynamic environment and therefore the cloud structure; however, limited work has been done to quantify this along the trajectory of the aerosol plume in the region. The deployment of the first Atmospheric Radiation Measurement (ARM) Mobile Facility (AMF1) in support of the Layered Atlantic Smoke Interactions with Clouds field campaign provided a unique opportunity to collect observations of cloud and aerosol properties during two consecutive biomass burning seasons during July through October of 2016 and 2017 over Ascension Island (7.96° S, 14.35° W). Using observed profiles of temperature, humidity, and clouds from the field campaign alongside aerosol optical properties from Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2), as input for the Rapid Radiation Transfer Model (RRTM), profiles of the radiative heating rate due to aerosols and clouds were computed. Radiative heating is also assessed across the southeast Atlantic Ocean using an ensemble of back trajectories from the Hybrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) model. Idealized experiments using the RRTM with and without aerosols and a range of values for the single-scattering albedo (SSA) demonstrate that shortwave (SW) heating within the aerosol layer above Ascension Island can locally range between 2 and 8 K d -1 depending on the aerosol optical properties, though impacts of the aerosol can be felt elsewhere in the atmospheric column. When considered under clear conditions, the aerosol has a cooling effect at the TOA, and based on the observed cloud properties at Ascension Island, the cloud albedo is not large enough to overcome this. Shortwave radiative heating due to biomass burning aerosol is not balanced by additional longwave (LW) cooling, and the net radiative impact results in a stabilization of the lower troposphere. However, these results are extremely sensitive to the single-scattering albedo assumptions in models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Observations of supermicron-sized aerosols originating from biomass burning in southern Central Africa

During the 3 years of the ObseRvations of Aerosols above CLouds and their intEractionS (ORACLES) campaign, the NASA Orion P-3 was equipped with a 2D stereo (2D-S) probe that imaged particles with maximum dimension (D) ranging from 10 < D < 1280 µm. The 2D-S recorded supermicron-sized aerosol particles (SAPs) outside of clouds within biomass burning plumes during flights over the southeastern Atlantic off Africa's coast. Numerous SAPs with 10 < D < 1520 µm were observed in 2017 and 2018 at altitudes between 1230 and 4000 m, 1000 km from the coastline, mostly between 7–11° S. No SAPs were observed in 2016 as flights were conducted further south and further from the coastline. Number concentrations of refractory black carbon (rBC) measured by a single particle soot photometer ranged from 200 to 1200 cm –3 when SAPs were observed. Transmission electron microscopy images of submicron particulates, collected on Holey carbon grid filters, revealed particles with potassium salts, black carbon (BC), and organics. Energy-dispersive X-ray spectroscopy spectra also detected potassium, a tracer for biomass burning. These measurements provided evidence that the submicron particles originated from biomass burning. NOAA Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) 3 d back trajectories show a source in northern Angola for times when large SAPs were observed. Fire Information for Resource Management System (FIRMS) Moderate Resolution Imaging Spectroradiometer (MODIS) 6 active fire maps showed extensive biomass burning at these locations. Given the back trajectories, the high number concentrations of rBC, and the presence of elemental tracers indicative of biomass burning, it is hypothesized that the SAPs imaged by the 2D-S are examples of BC aerosol, ash, or unburned plant material.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Exploring the elevated water vapor signal associated with the free tropospheric biomass burning plume over the southeast Atlantic Ocean

In southern Africa, widespread agricultural fires produce substantial biomass burning (BB) emissions over the region. The seasonal smoke plumes associated with these emissions are then advected westward over the persistent stratocumulus cloud deck in the southeast Atlantic (SEA) Ocean, resulting in aerosol effects which vary with time and location. Much work has focused on the effects of these aerosol plumes, but previous studies have also described an elevated free tropospheric water vapor signal over the SEA. Water vapor influences climate in its own right, and it is especially important to consider atmospheric water vapor when quantifying aerosol-cloud interactions and aerosol radiative effects. Here we present airborne observations made during the NASA ORACLES (ObseRvations of Aerosols above CLouds and their intEractionS) campaign over the SEA Ocean. In observations collected from multiple independent instruments on the NASA P-3 aircraft (from near-surface to 6-7 km), we observe a strongly linear correlation between pollution indicators (carbon monoxide (CO) and aerosol loading) and atmospheric water vapor content, seen at all altitudes above the boundary layer. The focus of the current study is on the especially strong correlation observed during the ORACLES-2016 deployment (out of Walvis Bay, Namibia), but a similar relationship is also observed in the August 2017 and October 2018 ORACLES deployments. Using reanalyses from the European Centre for Medium-Range Weather Forecasts (ECMWF) and Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2), and specialized WRF-Chem simulations, we trace the plume-vapor relationship to an initial humid, smoky continental source region, where it mixes with clean, dry upper tropospheric air and then is subjected to conditions of strong westward advection, namely the southern African easterly jet (AEJ-S). Our analysis indicates that air masses likely left the continent with the same relationship between water vapor and carbon monoxide as was observed by aircraft. This linear relationship developed over the continent due to daytime convection within a deep continental boundary layer (up to ~5-6 km) and mixing with higher-altitude air, which resulted in fairly consistent vertical gradients in CO and water vapor, decreasing with altitude and varying in time, but this water vapor does not originate as a product of the BB combustion itself. Due to a combination of conditions and mixing between the smoky, moist continental boundary layer and the dry and fairly clean upper-troposphere air above (~6 km), the smoky, humid air is transported by strong zonal winds and then advected over the SEA (to the ORACLES flight region) following largely isentropic trajectories. Hybrid Single-Particle Lagrangian Integrated Trajectory model (HYSPLIT) back trajectories support this interpretation. This work thus gives insights into the conditions and processes which cause water vapor to covary with plume strength. Better understanding of this relationship, including how it varies spatially and temporally, is important to accurately quantify direct, semi-direct, and indirect aerosol effects over this region.

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

Twistor Actions for Integrable Systems

Many integrable systems can be reformulated as holomorphic vector bundles on twistor space. This is a powerful organizing principle in the theory of integrable systems. One shortcoming is that it is formulated at the level of the equations of motion. From this perspective, it is mysterious that integrable systems have Lagrangians. In this paper, we study a Chern-Simons action on twistor space and use it to derive the Lagrangians of some integrable sigma models. Our focus is on examples that come from dimensionally reduced gravity and supergravity. The dimensional reduction of general relativity to two spacetime dimensions is an integrable coset sigma model coupled to a dilaton and 2d gravity. The dimensional reduction of supergravity to two spacetime dimensions is an integrable coset sigma model coupled to matter fermions, a dilaton, and 2d supergravity. We derive Lax operators and Lagrangians for these 2d integrable systems using the Chern-Simons theory on twistor space. In the supergravity example, we use an extended setup in which twistor Chern-Simons theory is coupled to a pair of matter fermions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗