Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spatial variability”

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 1,027 records · Page 57
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Technical Advances for Sustainable Remediation: In Situ Monitoring and Climate Resiliency - 18436

Recent advances in sustainable remediation techniques show promise as solutions for contaminant removal in soil and groundwater media. Although passive remediation (or natural attenuation) has gained momentum as an acceptable technique, there are still challenges in regulatory frameworks as well as in public acceptance. Leaving contaminants in subsurface requires increased burden of proof to show that plumes are stable and residual contaminants do not pose a significant health risk. Climate changes and frequent extreme climate events contribute another level of complexity and make it difficult to demonstrate system stability and predict future plume conditions. In this study, we show two coupled approaches to tackle this challenge: (1) in situ groundwater monitoring based on advanced sensor and data analytics technologies (i.e., Big Data), and (2) state-of-art modeling for optimizing monitoring and site closure strategies. In the new monitoring approach, in situ sensors are used to continuously measure the key variables that control contaminant plume mobility and its spatial and temporal distribution (such as pH, redox potential, electrical conductivity, and groundwater level). Coupled with a limited number of groundwater sampling, the data analytics methods - data mining and machine learning - allow us to estimate contaminant concentrations, and also to detect any significant changes associated with the plume mobility. In parallel, the state-of-art numerical flow and reactive transport model - ASCEM software under DOE Office of Environmental Management - provides a physical and mechanistic understanding of the system behavior, and predicts the long-term plume migration for optimizing the monitoring and closure strategies. In particular, we focus on the effects of extreme precipitation or drought on the site closure as well as monitoring strategies. We demonstrate this approach at the Savannah River Site (SRS) F-Area, where groundwater is contaminated by various radionuclides, including uranium, tritium and technetium. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Condition-Based Maintenance of a Circulating Water System of a Canadian Nuclear Power Plant using Machine Learning and Statistical Tools

Canada Deuterium Uranium pressurized-heavy-water reactors (PHWR) are a type of nuclear power plant that generate clean and reliable energy. The scope of this work is to automate data analysis methodologies to inform a condition-based maintenance strategy of a circulating water system (CWS) of a PHWR. The multiunit CWS provides a continuous supply of water to cool steam condensers, even during transient scenarios, thereby improving the thermal efficiency. This work aims to develop a machine learning (ML) based approach to detect anomalies in heterogeneous data of a CWS in a PHWR to help inform a predictive maintenance strategy. The heterogeneous data include textual and numeric time series data for a PHWR. Natural-language-processing (NLP)-based models are used to analyze textual data contained in work orders and operator logs and an event-timeseries correlation detection method is applied to assist anomalies diagnoses for CWS. An ML model Robust Linear Model (RLM) is also used to remove the seasonal variations in the system variable distributions based on distributions of environmental variables. A machine learning model, Density-Based Spatial Clustering of Applications with Noise (DBSCAN), trained on both original data and data without any seasonal variations will then be used to detect if an anomaly exists. Thus, by moving to an automated methodology to detect, classify, and forecast anomalies, the maintenance strategy would be based on component condition instead of a time-based schedule.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2020

GPS-Based Gamma Survey for Characterizing and Decommissioning NORM Sites - 20389

Gamma survey techniques are an especially powerful decommissioning tool at naturally occurring radioactive material (NORM) sites due to both the low cost to obtain data over a large spatial scale and the abundance of gamma emitters in the uranium and thorium decay series. Gamma surveys are executed by coupling a detector - most often a sodium iodide crystal - to a global positioning system (GPS), then reporting a location and gross gamma reading coincidentally to a data logger. Systems may be carried by workers or mounted to a car, all-terrain vehicle, or unmanned aerial system (UAS). The resulting data set provides a high-resolution but low precision map of the gamma radiation field over the area surveyed. Frequently this map is also correlated to soil concentrations of NORM radionuclides (most often, Ra-226) and/or exposure rate. Gamma survey parameters such as movement speed, transect spacing, and data logging frequency define the spatial resolution of the resulting surface, and can be optimized depending on the desired survey sensitivity. This paper examines gamma survey as a tool for decommissioning NORM sites and provides an overview of current gamma survey technology designed to improve the efficiency and effectiveness of the decommissioning process. Topics to be discussed in the paper include: - An overview of gamma survey systems, and the utility of different delivery vehicles depending on desired cost, desired spatial resolution, and site topography. - The influence of physical detector characteristics on detection sensitivity and survey planning. - The tradeoff between high-resolution and large spatial extent, but inherently uncertain data, and low-resolution, low spatial extent, but highly certain data, as well as the specific utility of each of these types of data during NORM facility decommissioning. - Confounding variables that may limit the utility of gamma survey at some sites (e.g., radon gas and spatial heterogeneity / hot spots), and methods to plan for and control these conditions. Results show that the confounding variables, such as radon and data output can greatly influence the overall data quality associated with the decommissioning process. In addition, the use of real-time and aerial survey platforms provides a method for ensuring proper spatial extent of the data. When applied thoughtfully, gamma survey is a powerful tool for detecting NORM radionuclides in the environment and a cost-effective technique for identifying areas requiring remediation. However, entities performing or using gamma survey as a decommissioning tool must be aware of both its advantages and its limitations before basing remediation or regulatory action on gamma survey results. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Quantification of spatial and seasonal trends in the atmosphere and construction of statistical models for infrasonic propagation

SUMMARY Infrasonic waves are influenced by variations in the density, pressure and temperature as well as the ambient winds. Modelling infrasonic propagation can be challenging due to the dynamic nature of the atmosphere as well as the sparseness of measurements which result in variability and notable uncertainty. A framework is presented to quantify spatial and seasonal trends in atmospheric structure via analysis of the effective sound speed profile and identification of temporal trends in the middle atmospheric waveguide produced by the circumpolar vortex winds. Seasonal definitions identifying typical atmospheric structures during the summer, winter and spring/fall transition periods are identified using atmospheric data from 2010 through 2020. Seasonal trend analysis is conducted for a number of locations across the contiguous United States to quantify spatial variations in atmospheric structure that impact infrasonic propagation. For each season and location, empirical orthogonal function analysis is used to reduce the historical archive of atmospheric data into a smaller representative set that can be analysed using numerical tools more efficiently. Infrasonic ray tracing and finite-frequency modal propagation analyses are applied to construct propagation path geometry and transmission loss statistics which are useful in localization and yield estimation for infrasonic sources, respectively. An example application is detailed in which transmission loss statistics are combined with an explosive source model and noise statistics to quantify the capability of a network to detect nearby sources.

58 GEOSCIENCES↗
DOE OSTI2017

Socio-economic exposure to natural disasters

Even though the correct assessment of risks is a key aspect of the risk management analysis, we argue that limited effort has been devoted in the assessment of comprehensive measures of economic exposure at very low scale. For this reason, we aim at providing a series of suitable methodologies to provide a complete and detailed list of the exposure of economic activities to natural disasters. We use Input-Output models to provide information about several socio-economic variables, such as population density, employment density, firms' turnover and capital stock, that can be seen as direct and indirect socio-economic exposure to natural disasters. We then provide an application to the Italian context. These measures can be easily incorporated into risk assessment models to provide a clear picture of the disaster risk for local areas. - Highlights: • Ex ante assessment of economic exposure to disasters at very low geographical scale • Assessment of the cost of natural disasters in ex-post perspective • IO model and spatial autocorrelation to get information on socio-economic variables • Indicators supporting risk assessment and risk management models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Characterizing the multisectoral impacts of future global hydrologic variability

There is significant uncertainty in how global water supply will evolve in the future, due to uncertain climate, socioeconomic, and land use change drivers and variability of hydrologic processes. It is critical to characterize the potential impacts of uncertainty in future water supply given its importance for food and energy production. In this work, we introduce a framework that integrates stochastic hydrology and human-environmental systems to characterize uncertainty in future water supply and its multisector impacts. We develop a global stochastic watershed model and demonstrate that this model can generate a large ensemble of realizations of basin-scale runoff with global coverage that preserves the mean, variance, and spatial correlation of a historical benchmark. We couple this model with a well-known human-environmental systems model to explore the impacts of runoff variability on the water and agricultural sectors across spatial scales. We find that the impacts of future hydrologic variability vary across sectors and regions. Impacts are felt most strongly in the water and agricultural sectors for basins that are expected to have unsustainable water use in the future, such as the Indus River basin. For this basin, we find that the variability in future irrigation water withdrawals and irrigated cropland increase over time due to uncertainty in renewable water supply. We also use the Indus basin to show how our stochastic ensemble can be leveraged to explore the global multisector consequences of local extreme runoff conditions. This work introduces a novel technique to explore the propagation of future hydrologic variability across human and natural systems and spatial scales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Soil Moisture Buffers the Impact of Precipitation Variability on Ecosystem Productivity

Water availability governs ecosystem productivity, yet estimates of vegetation sensitivity to water can differ greatly depending on whether the sensitivity is examined spatially or temporally. In particular, the spatial sensitivity is often reported to be much stronger than temporal sensitivities, leading to highly uncertain projections of ecosystem responses to future climate change when using space-for-time substitution. The large difference between spatial and temporal sensitivities remains unexplained. Prior research, however, primarily relied on precipitation as the water availability proxy, whereas vegetation responds to soil moisture. Here, we combined satellite estimates of vegetation productivity with soil moisture data across water-limited ecosystems of the continental United States (CONUS) to identify a convergent sensitivity of productivity to water availability. Using precipitation, we show that temporal sensitivity is 66% lower than spatial sensitivity overall. Our analysis identified the cause of the difference to be primarily driven by the seasonal variability of water availability, rooting depth, and soil properties. When using soil moisture instead of precipitation, we observed widespread convergence in the spatial and temporal sensitivities—that is, the two sensitivities became much more similar in magnitude across all water-limited ecosystems within CONUS. These results show that overlooking soil hydrology can inflate perceived discrepancies between spatial and temporal vegetation sensitivities, leading to biased projections of ecosystem dynamics under future hydro-climatic change.

Wang, Huiqi [University of California, Berkeley, C↗
DOE OSTI2018

Domain-Adapted Convolutional Networks for Satellite Image Classification: A Large-Scale Interactive Learning Workflow

Satellite imagery often exhibits large spatial extent areas that encompass object classes with considerable variability. This often limits large-scale model generalization with machine learning algorithms. Notably, acquisition conditions, including dates, sensor position, lighting condition, and sensor types, often translate into class distribution shifts introducing complex nonlinear factors and hamper the potential impact of machine learning classifiers. Here, this article investigates the challenge of exploiting satellite images using convolutional neural networks (CNN) for settlement classification where the class distribution shifts are significant. We present a large-scale human settlement mapping workflow based-off multiple modules to adapt a pretrained CNN to address the negative impact of distribution shift on classification performance. To extend a locally trained classifier onto large spatial extents areas we introduce several submodules: First, a human-in-the-loop element for relabeling of misclassified target domain samples to generate representative examples for model adaptation; second, an efficient hashing module to minimize redundancy and noisy samples from the mass-selected examples; and third, a novel relevance ranking module to minimize the dominance of source example on the target domain. The workflow presents a novel and practical approach to achieve large-scale domain adaptation with binary classifiers that are based-off CNN features. Experimental evaluations are conducted on areas of interest that encompass various image characteristics, including multisensors, multitemporal, and multiangular conditions. Domain adaptation is assessed on source–target pairs through the transfer loss and transfer ratio metrics to illustrate the utility of the workflow.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

A Satellite-Derived Upper-Ocean Stratification Data Set for the Tropical North Atlantic With Potential Applications for Hurricane Intensity Prediction

Upper-ocean stratification strongly impacts vertical mixing and the heat flux between the ocean and atmosphere, especially under extreme conditions of tropical cyclones (TCs). Knowledge of prestorm stratification is important for accurate TC intensity prediction. In situ observations of the tropical ocean have significantly increased in the past decade. However, they are still too sparse to resolve ocean stratification variability in near-real time and on small spatial scales. Here, in this study, based on long-term observations and an ocean reanalysis data set from 2004–2017, we investigate the possibility of retrieving upper-ocean stratification from sea surface temperature (SST), sea surface salinity (SSS), and sea surface height (SSH) using a simple regression method. It is found that more than 90% of the mean seasonal cycle and about 30% to 80% of temperature and salinity stratification anomalies can be reconstructed using surface data from either observations or an ocean reanalysis. Simple regression can be used with satellite observations to create a high-resolution, near-real-time-gridded ocean stratification data set that successfully reproduces both the large and mesoscale variability of ocean stratification. When used in a simple expression for TC-induced SST cooling, the satellite-derived stratification shows improvements over an ocean analysis in terms of variance explained of SST cooling, offering promise as a near-real-time indicator of the ocean's impact on TC intensification.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

Unprecedented study of the broadband emission of Mrk 421 during flaring activity in March 2010

Because of its proximity, Mrk 421 is one of the best sources on which to study the nature of BL Lac objects. Its proximity allows us to characterize its broadband spectral energy distribution (SED). Here, the goal is to better understand the mechanisms responsible for the broadband emission and the temporal evolution of Mrk 421. These mechanisms may also apply to more distant blazars that cannot be studied with the same level of detail. A flare occurring in March 2010 was observed for 13 consecutive days (from MJD 55 265 to MJD 55 277) with unprecedented wavelength coverage from radio to very high energy (VHE; E> 100 GeV) γ-rays with MAGIC, VERITAS, Whipple , Fermi -LAT, MAXI, RXTE, Swift , GASP-WEBT, and several optical and radio telescopes. We modeled the day-scale SEDs with one-zone and two-zone synchrotron self-Compton (SSC) models, investigated the physical parameters, and evaluated whether the observed broadband SED variability can be associated with variations in the relativistic particle population. As a result, the activity of Mrk 421 initially was high and then slowly decreased during the 13-day period. The flux variability was remarkable at the X-ray and VHE bands, but it was minor or not significant at the other bands. The variability in optical polarization was also minor. These observations revealed an almost linear correlation between the X-ray flux at the 2–10 keV band and the VHE γ-ray flux above 200 GeV, consistent with the γ-rays being produced by inverse-Compton scattering in the Klein-Nishina regime in the framework of SSC models. The one-zone SSC model can describe the SED of each day for the 13 consecutive days reasonably well, which once more shows the success of this standard theoretical scenario to describe the SEDs of VHE BL Lacs such as Mrk 421. This flaring activity is also very well described by a two-zone SSC model, where one zone is responsible for the quiescent emission, while the other smaller zone, which is spatially separated from the first, contributes to the daily variable emission occurring at X-rays and VHE γ-rays. The second blob is assumed to have a smaller volume and a narrow electron energy distribution with 3 × 104<γ< 6 × 105, where γ is the Lorentz factor of the electrons. Such a two-zone scenario would naturally lead to the correlated variability at the X-ray and VHE bands without variability at the optical/UV band, as well as to shorter timescales for the variability at the X-ray and VHE bands with respect to the variability at the other bands. In conclusion, both the one-zone and the two-zone SSC models can describe the daily SEDs via the variation of only four or five model parameters, under the hypothesis that the variability is associated mostly with the underlying particle population. This shows that the particle acceleration and cooling mechanism that produces the radiating particles might be the main mechanism responsible for the broadband SED variations during the flaring episodes in blazars. The two-zone SSC model provides a better agreement with the observed SED at the narrow peaks of the low- and high-energy bumps during the highest activity, although the reported one-zone SSC model could be further improved by varying the parameters related to the emitting region itself (δ, B and R), in addition to the parameters related to the particle population.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

THE ISLANDS PROJECT. I. ANDROMEDA XVI, AN EXTREMELY LOW MASS GALAXY NOT QUENCHED BY REIONIZATION

Based on data aquired in 13 orbits of Hubble Space Telescope time, we present a detailed evolutionary history of the M31 dSph satellite Andromeda XVI, including its lifetime star formation history (SFH), the spatial distribution of its stellar populations, and the properties of its variable stars. And XVI is characterized by prolonged star formation activity from the oldest epochs until star formation was quenched ∼6 Gyr ago, and, notably, only half of the mass in stars of And XVI was in place 10 Gyr ago. And XVI appears to be a low-mass galaxy for which the early quenching by either reionization or starburst feedback seems highly unlikely, and thus it is most likely due to an environmental effect (e.g., an interaction), possibly connected to a late infall in the densest regions of the Local Group. Studying the SFH as a function of galactocentric radius, we detect a mild gradient in the SFH: the star formation activity between 6 and 8 Gyr ago is significantly stronger in the central regions than in the external regions, although the quenching age appears to be the same, within 1 Gyr. We also report the discovery of nine RR Lyrae (RRL) stars, eight of which belong to And XVI. The RRL stars allow a new estimate of the distance, (m − M){sub 0} = 23.72 ± 0.09 mag, which is marginally larger than previous estimates based on the tip of the red giant branch.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Second-generation downscaled earth system model data using generative machine learning

The second-generation Sup3rCC dataset provides high-resolution meteorological data generated through the downscaling of multiple earth system models (ESMs) from the Coupled Model Intercomparison Project Phase 6 (CMIP6). This downscaling is performed through application of a generative machine learning approach called Super-Resolution for Renewable Resource Data (sup3r). This dataset builds on the first-generation Sup3rCC data by applying improved bias correction methods and adding downscaled precipitation to the output variables. As with the first Sup3rCC version, the data still include temperature, wind speed and direction at multiple heights, pressure, three components of downwelling solar radiation, and relative humidity—all at 4-kilometer (km) hourly resolution over the contiguous United States. This is a 25x spatial enhancement and 24x temporal enhancement of the source 100-km daily-average ESM data. This extension of the Sup3rCC dataset includes data from six ESMs from two shared socioeconomic pathways (SSPs) totaling 400 years of data with multiple future projections of changing meteorological conditions. The scenario selection was based on a structured evaluation of historical ESM skill and comprehensive representation of possible trajectories of future climate change in temperature, humidity, precipitation, solar irradiance, and near-surface wind speeds. The inclusion of multiple future projections is intended to enable users to assess key drivers of un 36 certainty and variability. All data are double-bias corrected, resulting in a product that can be used out-of-the-box for energy system analysis with minimal historical bias. The potential applications of Sup3rCC data extend to various topics in renewable energy resource assessment, energy systems modeling, and grid resilience studies. High-resolution future meteorological projections are critical for evaluating the effects of changing meteorological conditions on renewable energy generation, energy demand, and for optimizing energy storage and grid infrastructure. The 4-km hourly resolution of the downscaled data enables understanding of spatial and temporal variability at the scales necessary for energy system operational planning. In addition, the dataset can support risk assessments by providing detailed information on possible future extreme weather events and long-term meteorological variability at scales relevant to energy infrastructure. By offering an enhanced representation of possible future meteorological conditions, the second-generation Sup3rCC dataset enables more precise modeling of energy resilience and adaptation strategies in response to changing meteorological conditions.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2018

ENSO's Changing Influence on Temperature, Precipitation, and Wildfire in a Warming Climate

On interannual to decadal time scales, the climate mode with many of the strongest societal impacts is the El Niño–Southern Oscillation (ENSO). However, quantifying ENSO's changes in a warming climate remains a formidable challenge, due to both the noise arising from internal variability and the complexity of air-sea feedbacks in the tropical Pacific Ocean. In this paper, we use large (≥30-member) ensembles of climate simulations to show that anthropogenic climate change can produce systematic increases in ENSO teleconnection strength over many land regions, driving increased interannual variability in regional temperature extremes and wildfire frequency. As the spatial character of this intensification exhibits strong land-ocean contrasts, a causal role for land-atmosphere feedbacks is suggested. The identified increase in variance occurs in multiple model ensembles, independent of changes in sea surface temperature variance. In conclusion, this suggests that in addition to changes in the overall likelihoods of heat and wildfire extremes, the variability in these events may also be a robust feature of future climate.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Genetic Identification of Lamprey Genera and Anadromous Ecotypes in Watersheds of the Northeastern Pacific Ocean

ABSTRACT Nonparasitic, nonmigratory Western Brook Lamprey (WBL; Lampetra ayresii ), and parasitic, anadromous Western River Lamprey (WRL; L. ayresii ) are sympatric lampreys that likely represent different life history variations of a single species. Novel genetic tools are critical for differentiating WBL and WRL, whose larvae preclude morphological identification (ID) and will enable comprehensive assessment of imperiled native lampreys of the Northeastern Pacific (including WBL, WRL, and Pacific Lamprey, Entosphenus tridentatus ). We developed 47 candidate single nucleotide polymorphism (SNP) markers using whole genome resequencing of WBL ( N = 24) and WRL ( N = 15) from Ksi Ts'oohl Ts'ap Creek (Nass River, British Columbia, Canada) which are likely ecotypes distinguished by few divergent SNPs across multiple chromosomes. We used five novel candidate SNPs to perform genetic ID of WBL and WRL ecotypes in collections of mixed native lampreys from lower Columbia River tributaries ( N = 1474), Ksi Ts'oohl Ts'ap Creek ( N = 352), and ocean phase WRL from the Georgia Basin (Salish Sea, British Columbia, Canada; N = 91). Two previously published SNPs were used to ID genera, Entosphenus versus Lampetra . Morphological ID utilized photographs collected from a subset of genotyped lampreys, and high concordance was demonstrated between ID methods for genera (99%) and Lampetra ecotypes (> 98%). We characterized spatial and temporal composition of lamprey genera and ecotypes surveyed across NE Pacific tributaries under the expectation these compositions would be similar across nearby sites and across years at the same site. Proportions of lamprey genera were highly variable within regions and across years; however, Lampetra ecotypic proportions were spatially and temporally stable. WRL were rare in lower Columbia tributaries (~1% average rate among Lampetra ) and common further north (> 40% of Lampetra ). Genetic ID methods are powerful monitoring tools that create the novel ability to ascertain genera and ecotypes regardless of life stage, while increasing the efficiency of surveys by eliminating time‐intensive morphological data collection.

Silver, G. S. [Columbia River Inter‐Tribal Fish Co↗
DOE OSTI2026

Short-term electricity load forecasting: Application-driven evaluation of machine learning models across spatial and temporal scales

As we transition towards a decarbonized economy, the integration of variable renewable energy resources and new demands (e.g., electric vehicles, heat pumps) into the electricity grid places unprecedented pressure on grid operators to effectively anticipate and manage peak load. In this context, machine learning algorithms are proving to be indispensable for accurate short-term load forecasting, a crucial task to address these challenges. This study benchmarks 6 machine learning algorithms, including three neural networks and three tree-based algorithms, across various levels of spatial aggregation and time horizons (1, 4, 8, 24, and 48 h). The central contribution of this work is the comparison and analysis of load forecasting models not only based on statistical metrics, but also based on a novel error metric, which evaluates the cost implications of forecast errors for power system stakeholders. Results show that tree-based models outperform neural networks, based on statistical metrics, and yield less skewed error distributions for most spatial scales. However, through the lens of the novel error metric, neural networks are the more competitive choice, especially for forecast horizons that exceed 8 h. The study concludes with actionable recommendations to grid operators and highlights the need for the development of error metrics that link forecasting accuracy to operational costs. To promote transparency and open science, the datasets and Python code are open-sourced via a supplementary repository.

Houben, Nikolaus↗
DOE OSTI2013

Phase-space finite elements in a least-squares solution of the transport equation

The linear Boltzmann transport equation is solved using a least-squares finite element approximation in the space, angular and energy phase-space variables. The method is applied to both neutral particle transport and also to charged particle transport in the presence of an electric field, where the angular and energy derivative terms are handled with the energy/angular finite elements approximation, in a manner analogous to the way the spatial streaming term is handled. For multi-dimensional problems, a novel approach is used for the angular finite elements: mapping the surface of a unit sphere to a two-dimensional planar region and using a meshing tool to generate a mesh. In this manner, much of the spatial finite-elements machinery can be easily adapted to handle the angular variable. The energy variable and the angular variable for one-dimensional problems make use of edge/beam elements, also building upon the spatial finite elements capabilities. The methods described here can make use of either continuous or discontinuous finite elements in space, angle and/or energy, with the use of continuous finite elements resulting in a smaller problem size and the use of discontinuous finite elements resulting in more accurate solutions for certain types of problems. The work described in this paper makes use of continuous finite elements, so that the resulting linear system is symmetric positive definite and can be solved with a highly efficient parallel preconditioned conjugate gradients algorithm. The phase-space finite elements capability has been built into the Sceptre code and applied to several test problems, including a simple one-dimensional problem with an analytic solution available, a two-dimensional problem with an isolated source term, showing how the method essentially eliminates ray effects encountered with discrete ordinates, and a simple one-dimensional charged-particle transport problem in the presence of an electric field. (authors)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Understanding the Drivers of Atlantic Multidecadal Variability using a Stochastic Model Hierarchy

The relative importance of ocean and atmospheric dynamics in generating Atlantic Multidecadal Variability (AMV) remains an open question. Comparisons between climate models with SLAB and fully-dynamic (FULL) ocean components are often used to explore this question, but cannot reveal how individual ocean processes generate these differences. We build a hierarchy of physically interpretable stochastic models to investigate the contribution of two upper-ocean processes to AMV: the role of seasonal variation and mixed-layer entrainment. This interpretability arises from the stochastic model’s simplified representation of sea surface temperature (SST), considering only the local upper ocean response to white-noise atmospheric forcing and its impact on surface heat exchange. We focus on understanding differences between SLAB and FULL non-eddy resolving pre-industrial control simulations of the Community Earth System Model 1 (CESM), and estimate the stochastic model parameters from each respective simulation. Despite its simplicity, the stochastic model reproduces temporal characteristics of SST variability in the SPG, including reemergence, seasonal-to-interannual persistence and power spectra. Furthermore, unrealistically persistent SST of the CESM-SLAB ocean simulation is reproduced in the equivalent stochastic model configuration where the mixed-layer depth (MLD) is constant. The stochastic model also reveals that vertical entrainment primarily damps SST variability, thus explaining why SLAB exhibits larger SST variance than FULL. Here, the stochastic model driven by temporally stochastic, spatially coherent forcing patterns reproduces the canonical AMV pattern. However, the amplitude of low-frequency variability remains underestimated, suggesting a role for ocean dynamics beyond entrainment.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

PROBING THE IONIZATION STATES OF POLYCYCLIC AROMATIC HYDROCARBONS VIA THE 15–20 μm EMISSION BANDS

We report new correlations between ratios of band intensities of the 15–20 μm emission bands of polycyclic aromatic hydrocarbons (PAHs) in a sample of 57 sources observed with the Spitzer/Infrared Spectrograph. This sample includes Large Magellanic Cloud point sources from the SAGE-Spec survey, nearby galaxies from the Spitzer Infrared Nearby Galaxies Survey survey, two Galactic interstellar medium cirrus sources, and the spectral maps of the Galactic reflection nebulae NGC 2023 and NGC 7023. We find that the 16.4, 17.4, and 17.8 μm band intensities are inter-correlated in all environments. In NGC 2023 and NGC 7023 these bands also correlate with the 11.0 and 12.7 μm band intensities. The 15.8 μm band correlates only with the 15–18 μm plateau and the 11.2 μm emission. We examine the spatial morphology of these bands and introduce radial cuts. We find that these bands can be spatially organized into three sets: the 12.7, 16.4, and 17.8 μm bands; the 11.2, 15.8 μm bands and the 15–18 μm plateau; and the 11.0 and 17.4 μm bands. We also find that the spatial distribution of the 12.7, 16.4, and 17.8 μm bands can be reconstructed by averaging the spatial distributions of the cationic 11.0 μm and neutral 11.2 μm bands. We conclude that the 17.4 μm band is dominated by cations, the 15.8 μm band by neutral species, and the 12.7, 16.4, and 17.8 μm bands by a combination of the two. These results highlight the importance of PAH ionization for spatially differentiating sub-populations by their 15–20 μm emission variability.

79 ASTRONOMY AND ASTROPHYSICS↗