Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “local sensitivity factor”

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 217 records · Page 12

Dissipation and Bathymetric Sensitivities in an Unstructured Mesh Global Tidal Model

Abstract The mechanisms and geographic distribution of global tidal dissipation in barotropic tidal models are examined using a high resolution unstructured mesh finite element model. Mesh resolution varies between 2 and 25 km and is especially focused on inner shelves and steep bathymetric gradients. Tidal response sensitivities to bathymetric changes are examined to put into context response sensitivities to frictional processes. We confirm that the Ronne Ice Shelf dramatically affects Atlantic tides but also find that bathymetry in the Hudson Bay system is a critical control. We follow a sequential frictional parameter optimization process and use TPXO9 data‐assimilated tidal elevations as a reference solution. From simulated velocities and depths, dissipation within the global model is estimated and allows us to pinpoint dissipation at high resolution. Boundary layer dissipation is extremely focused with 1.4% of the ocean accounting for 90% of the total. Internal tide friction is much more distributed with 16.7% of the ocean accounting for 90% of the total. Often highly regional dissipation can impact basin‐scale and even ocean wide tides. Optimized boundary layer friction parameters correlate very well with the physical characteristics of the locality with high friction factors associated with energetic tidal regions, deep ocean island chains, and ice covered areas. Global complex M 2 tide errors are 1.94 cm in deep waters. Total global boundary layer and internal tide dissipation are estimated, respectively, at 1.83 and 1.49 TW. This continues the trend in the literature toward attributing more dissipation to internal tides.

54 ENVIRONMENTAL SCIENCES↗

An Algorithm for Efficient Maximum Likelihood Estimation and Confidence Interval Determination in Nonlinear Estimation Problems

An algorithm for maximum likelihood (ML) estimation is developed with an efficient method for approximating the sensitivities. The algorithm was developed for airplane parameter estimation problems but is well suited for most nonlinear, multivariable, dynamic systems. The ML algorithm relies on a new optimization method referred to as a modified Newton-Raphson with estimated sensitivities (MNRES). MNRES determines sensitivities by using slope information from local surface approximations of each output variable in parameter space. The fitted surface allows sensitivity information to be updated at each iteration with a significant reduction in computational effort. MNRES determines the sensitivities with less computational effort than using either a finite-difference method or integrating the analytically determined sensitivity equations. MNRES eliminates the need to derive sensitivity equations for each new model, thus eliminating algorithm reformulation with each new model and providing flexibility to use model equations in any format that is convenient. A random search technique for determining the confidence limits of ML parameter estimates is applied to nonlinear estimation problems for airplanes. The confidence intervals obtained by the search are compared with Cramer-Rao (CR) bounds at the same confidence level. It is observed that the degree of nonlinearity in the estimation problem is an important factor in the relationship between CR bounds and the error bounds determined by the search technique. The CR bounds were found to be close to the bounds determined by the search when the degree of nonlinearity was small. Beale's measure of nonlinearity is developed in this study for airplane identification problems; it is used to empirically correct confidence levels for the parameter confidence limits. The primary utility of the measure, however, was found to be in predicting the degree of agreement between Cramer-Rao bounds and search estimates.

Murphy, Patrick Charles↗

Revisiting the Effectiveness Of HCHO/NO2 Ratios for Inferring Ozone Sensitivity to Its Precursors Using High Resolution Airborne Remote Sensing Observations in A High Ozone Episode During the KORUS-AQ Campaign

The nonlinear chemical processes involved in ozone production (P(O3)) have necessitated using proxy indicators to convey information about the primary dependence of P(O3) on volatile organic compounds (VOCs) or nitrogen oxides (NOx). In particular, the ratio of remotely sensed columns of formaldehyde (HCHO) to nitrogen dioxide (NO2) has been widely used for studying O3 sensitivity. Previous studies found that the errors in retrievals and the incoherent relationship between the column and the near-surface concentrations are a barrier in applying the ratio in a robust way. In addition to these obstacles, we provide calculational-observational evidence, using an ensemble of 0-D photochemical box models constrained by DC-8 aircraft measurements on an ozone event during the Korea-United States Air Quality (KORUS-AQ) campaign over Seoul, to demonstrate the chemical feedback of NO2 on the formation of HCHO is a controlling factor for the transition line between NOx-sensitive and NOx-saturated regimes. A fixed value (~2.7) of the ratio of the chemical loss of NOx (LNOx) to the chemical loss of HO2+RO2 (LROx) perceptibly differentiates the regimes. Following this value, data points with a ratio of HCHO/NO2 less than 1 can be safely classified as NOx-saturated regime, whereas points with ratios between 1 and 4 fall into one or the other regime. We attribute this mainly to the HCHO-NO2 chemical relationship causing the transition line to occur at larger (smaller) HCHO/NO2 ratios in VOC-rich (VOC-poor) environments. We then redefine the transition line to LNOx/LROx~2.7 that accounts for the HCHO-NO2 chemical relationship leading to HCHO = 3.7 × (NO2 – 1.14 × 1016 molec.cm-2). Although the revised formula is locally calibrated (i.e., requires for readjustment for other regions), its mathematical format removes the need for having a wide range of thresholds used in HCHO/NO2 ratios that is a result of the chemical feedback. Therefore, to be able to properly take the chemical feedback into consideration, the use of HCHO = a × (NO2 – b) formula should be preferred to the ratio in future works. We then use the Geostationary Trace gas and Aerosol Sensor Optimization (GeoTASO) airborne instrument to study O3 sensitivity in Seoul. The unprecedented spatial (250 × 250 m2) and temporal (~every 2 h) resolutions of HCHO and NO2 observations form the sensor enhance our understanding of P(O3) in Seoul; rather than providing a crude label for the entire city, more in-depth variabilities in chemical regimes are observed that should be able to inform mitigation strategies correspondingly.

KORUS-AQ; formaldehyde; NOx; VOC; OMI; GeoTASO; oz↗

Environmental controls on observed spatial variability of soil pore water geochemistry in small headwater catchments underlain with permafrost

Abstract. Soil pore water (SPW) chemistry can vary substantially across multiple scales in Arctic permafrost landscapes. The magnitude of these variations and their relationship to scale are critical considerations for understanding current controls on geochemical cycling and for predicting future changes. These aspects are especially important for Arctic change modeling where accurate representation of sub-grid variability may be necessary to predict watershed-scale behaviors. Our research goal is to characterize intra- and inter-watershed soil water geochemical variations at two contrasting locations in the Seward Peninsula of Alaska, USA. We then attempt to identify the key factors controlling concentrations of important pore water solutes in these systems. The SPW geochemistry of 18 locations spanning two small Arctic catchments was examined for spatial variability and its dominant environmental controls. The primary environmental controls considered were vegetation, soil moisture and/or redox condition, water–soil interactions and hydrologic transport, and mineral solubility. The sampling locations varied in terms of vegetation type and canopy height, presence or absence of near-surface permafrost, soil moisture, and hillslope position. Vegetation was found to have a significant impact on SPW NO3- concentrations, associated with the localized presence of nitrogen-fixing alders and mineralization and nitrification of leaf litter from tall willow shrubs. The elevated NO3- concentrations were, however, frequently equipoised by increased microbial denitrification in regions with sufficient moisture to support it. Vegetation also had an observable impact on soil-moisture-sensitive constituents, but the effect was less significant. The redox conditions in both catchments were generally limited by Fe reduction, seemingly well-buffered by a cache of amorphous Fe hydroxides, with the most reducing conditions found at sampling locations with the highest soil moisture content. Non-redox-sensitive cations were affected by a wide variety of water–soil interactions that affect mineral solubility and transport. Identification of the dominant controls on current SPW hydrogeochemistry allows for qualitative prediction of future geochemical trends in small Arctic catchments that are likely to experience warming and permafrost thaw. As source areas for geochemical fluxes to the broader Arctic hydrologic system, geochemical processes occurring in these environments are particularly important to understand and predict with regards to such environmental changes.

Conroy, Nathan Alec (ORCID:0000000305973373)↗

Insights into Constraining Rate Coefficients in Fuel Oxidation Mechanisms Using Genetic Algorithm Optimization

Accurate fuel oxidation mechanisms can enable predictive capabilities that aid in advancing combustion technologies. High-level computational kinetics can yield reasonable rate coefficients with uncertainties, in some cases, below a factor of 2. Computed rate coefficients can be constrained further by optimizing against experimental data. Here, we explore the application of genetic algorithm (GA) optimization to constrain computed rate coefficients in complex fuel oxidation mechanisms in conjunction with temperature-dependent species mole fractions from jet-stirred reactor (JSR) measurements. Cyclohexane is a model candidate for understanding the reactivity of cyclic fuels. In this work, we optimize the rate coefficients of the most recent literature cyclohexane mechanism, which incorporates theoretically computed rate coefficients for the reaction networks stemming from the first and second O 2 addition pathways, against the experimental results of two separate literature JSR studies. Optimization consistency is evaluated by carrying out three GA optimizations: fitting to the temperature-dependent species mole fractions in each JSR experiment separately and simultaneously fitting the species mole fractions in both experiments. Local sensitivity analyses are used to identify five influential low-temperature oxidation reactions for optimization. Although the three optimizations do not yield identical rate coefficients, the direction of change in all five rate coefficients is consistent among the three optimizations. Performance of the models from the three optimizations is assessed against literature ignition delay times with differences in the level of agreement observed among the different optimizations. In conclusion, comparisons are made with our recent optimization work of a cyclopentane oxidation master-equation model against time-resolved species concentrations, and insights and improvements of the strategy for constraining rate coefficients using GA optimization are discussed.

33 ADVANCED PROPULSION SYSTEMS↗

Modeling Molecular Hydrogen in Low-metallicity Galaxies

We use a suite of hydrodynamics simulations of the interstellar medium (ISM) within a galactic disk, which includes radiative transfer, a nonequilibrium model of molecular hydrogen, and a realistic model for star formation and feedback, to study the structure of the ISM and H 2 abundance as a function of local ISM properties. We show that the star formation rate and structure of the ISM are sensitive to the metallicity of the gas with a progressively smoother density distribution with decreasing metallicity. In addition to the well-known trend of the H I–H 2 transition shifting to higher densities with decreasing metallicity, the maximum achieved molecular fraction in the ISM drops drastically at Z ≲ 0.2 Z ⊙ as the formation time of H 2 becomes much longer than a typical lifetime of dense regions of the ISM. We present accurate fitting formulae for both volumetric and projected ƒ H 2 measured on different scales as a function of gas metallicity, UV radiation field, and gas density. We show that when the formulae are applied to the patches in the simulated galaxy, the overall molecular gas mass is reproduced to better than a factor of ≲1.5 across the entire range of metallicities and scales. We also show that the presented fit is considerably more accurate than any of the previous ƒ H 2 models and fitting formulae in the low-metallicity regime. The fit can thus be used for modeling molecular gas in low-resolution simulations and semi-analytic models of galaxy formation in the dwarf and high-redshift regimes.

79 ASTRONOMY AND ASTROPHYSICS↗

Galaxy Strategy for Ligo-Virgo Gravitational Wave Counterpart Searches

In this work we continue a line of inquiry begun in Kanner et al. which detailed a strategy for utilizing telescopes with narrow fields of view, such as the Swift X-Ray Telescope (XRT), to localize gravity wave (GW) triggers from LIGO (Laser Interferometer Gravitational-Wave Observatory) / Virgo. If one considers the brightest galaxies that produce 50 percent of the light, then the number of galaxies inside typical GW error boxes will be several tens. We have found that this result applies both in the early years of Advanced LIGO when the range is small and the error boxes large, and in the later years when the error boxes will be small and the range large. This strategy has the beneficial property of reducing the number of telescope pointings by a factor 10 to 100 compared with tiling the entire error box. Additional galaxy count reduction will come from a GW rapid distance estimate which will restrict the radial slice in search volume. Combining the bright galaxy strategy with a convolution based on anticipated GW localizations, we find that the searches can be restricted to about 18 plus or minus 5 galaxies for 2015, about 23 plus or minus 4 for 2017, and about 11 plus or minus for 2020. This assumes a distance localization at the putative neutron star-neutron star (NS-NS) merger range mu for each target year, and these totals are integrated out to the range. Integrating out to the horizon would roughly double the totals. For localizations with r (rotation) greatly less than mu the totals would decrease. The galaxy strategy we present in this work will enable numerous sensitive optical and X-ray telescopes with small fields of view to participate meaningfully in searches wherein the prospects for rapidly fading afterglow place a premium on a fast response time.

galaxies: statistics aEuro" gamma-ray burst: gener↗

Incorporating climate change into risk-informed resilience planning

In response to the development of portfolio-wide Climate Action Plans by federal agencies, federal sites are working to incorporate the impacts of climate change into their resilience assessments. However, it can be challenging to incorporate climate change scenarios into resilience assessments given the uncertainty inherent in climate change modeling. Incorporating these factors into a resilience plan requires an understanding of what the different climate scenarios mean, as well as how to estimate potential impacts of climate change on hazard occurrence on a regional, or even local, scale under different scenarios. We discuss approaches to incorporating this data into risk-informed resilience assessment processes, such as those implemented in the Department of Energy’s (DOE) Federal Energy Management Program’s (FEMP) Technical Resilience Navigator (TRN) and the Sustainability Performance Division’s (SPD) Vulnerability Assessment and Resilience Plan (VARP) Risk Assessment Tool. We also describe the climate scenarios and the availability of hazard data for site resilience planning, based on modeling included in the International Panel on Climate Change (IPCC), the National Climate Assessment (NCA), and state-level reports. We present examples from the TRN risk assessment and the VARP Risk Assessment Tool to illustrate how sensitivity analysis can be used to incorporate climate change projections into the resilience planning process.

Rabinowitz, Hannah S.↗

Volcanism, global catastrophe and mass mortality

The effects of very large volcanic eruptions are well documented in many studies, mostly based on observations made on three historic eruptions, Laki 1783; Tambora 1815 and Krakatau 1883. Such eruptions have effects that are catastrophic locally and measurable globally, but it is not clear that even the largest volcanic eruptions have had global catastrophic effects, nor caused mass extinctions. Two different types of volcanic eruption were considered as likely to have the most serious widespread effects: large silicic explosive eruptions producing hundreds or thousands of cubic kilometers of pyroclastic materials, and effusive basaltic eruptions producing of approximately 100 cubic kilometers of lava. In both cases, the global effects are climatic, and attributable to production of stratospheric aerosols. Other possibilities need to be explored. Recent research on global change has emphasized the extreme sensitivity of the links between oceanic circulation, atmospheric circulation and climate. In particular, it was argued that the pattern of ocean current circulation (which strongly influences climate) is unstable; it may rapidly flip from one pattern to a different one, with global climatic consequences. If volcanism has been a factor in global environmental change and a cause of mass extinctions, it seems most likely that it has done so by providing a trigger to other processes, for example by driving oceanic circulation from one mode to another.

Francis, P. W.↗

Improving Photometry and Stellar Signal Preservation with Pixel-Level Systematic Error Correction

The Kepler Mission has demonstrated that excellent stellar photometric performance can be achieved using apertures constructed from optimally selected CCD pixels. The clever methods used to correct for systematic errors, while very successful, still have some limitations in their ability to extract long-term trends in stellar flux. They also leave poorly correlated bias sources, such as drifting moiré pattern, uncorrected. We will illustrate several approaches where applying systematic error correction algorithms to the pixel time series, rather than the co-added raw flux time series, provide significant advantages. Examples include, spatially localized determination of time varying moiré pattern biases, greater sensitivity to radiation-induced pixel sensitivity drops (SPSDs), improved precision of co-trending basis vectors (CBV), and a means of distinguishing the stellar variability from co-trending terms even when they are correlated. For the last item, the approach enables physical interpretation of appropriately scaled coefficients derived in the fit of pixel time series to the CBV as linear combinations of various spatial derivatives of the pixel response function (PRF). We demonstrate that the residuals of a fit of soderived pixel coefficients to various PRF-related components can be deterministically interpreted in terms of physically meaningful quantities, such as the component of the stellar flux time series which is correlated with the CBV, as well as, relative pixel gain, proper motion and parallax. The approach also enables us to parameterize and assess the limiting factors in the uncertainties in these quantities.

Kolodzijczak, Jeffrey J.↗

COS-B observations of localized sources of gamma-ray emission

In October 1975, the high-energy gamma-ray flux from the Vela pulsar measured by COS-B was found to be 1.6 to 2.1 times higher than the flux measured by SAS-2 in 1973. This factor is too large to be accounted for by error in the COS-B calibration or analysis. This is supported by a comparison of the COS-B measurement of the narrow-line component from the galactic center region with the flux derived from the measurements of SAS-2; the COS-B flux comes out about 15 percent lower than the SAS-2 figure. It is interesting to note that a glitch in the pulsar period took place about 1 month prior to the COS-B observation; the previous glitch occurred about 1.5 years before the SAS-2 observation. The increased rotational energy loss after the glitch cannot simply explain the increased gamma-ray luminosity. If the two phenomena are related, the gamma-ray emission, absorption, or beaming process must be extremely sensitive to changes in rotational parameters. The existence is confirmed of a second region of enhanced radiation in the galactic anticenter in addition to that from the Crab pulsar.

Source record↗

HEMCO v1.0: A Versatile, ESMF-Compliant Component for Calculating Emissions in Atmospheric Models

We describe the Harvard-NASA Emission Component version 1.0 (HEMCO), a stand-alone software component for computing emissions in global atmospheric models. HEMCO determines emissions from different sources, regions, and species on a user-defined grid and can combine, overlay, and update a set of data inventories and scale factors, as specified by the user through the HEMCO configuration file. New emission inventories at any spatial and temporal resolution are readily added to HEMCO and can be accessed by the user without any preprocessing of the data files or modification of the source code. Emissions that depend on dynamic source types and local environmental variables such as wind speed or surface temperature are calculated in separate HEMCO extensions. HEMCO is fully compliant with the Earth System Modeling Framework (ESMF) environment. It is highly portable and can be deployed in a new model environment with only few adjustments at the top-level interface. So far, we have implemented HEMCO in the NASA Goddard Earth Observing System (GEOS-5) Earth system model (ESM) and in the GEOS-Chem chemical transport model (CTM). By providing a widely applicable framework for specifying constituent emissions, HEMCO is designed to ease sensitivity studies and model comparisons, as well as inverse modeling in which emissions are adjusted iteratively. The HEMCO code, extensions, and the full set of emissions data files used in GEOS-Chem are available at http: //wiki.geos-chem.org/HEMCO.

GEOS-5↗

A non-canonical fungal peroxisome PTS-1 signal, SYM, and its evolutionary aspects

Abstract Proteins localized to peroxisomes, particularly those expressed under specific conditions or in low abundance, are often undetected by routine proteomics methods due to detection sensitivity limits. In silico identification and experimental validation of peroxisomal targeting signals (PTSs) offer a reliable alternative. We demonstrate that SYM, a non-canonical plant PTS-1 signal, functions similarly inAspergillus nidulans, as GFP tagged with a SYM C-terminal tripeptide localizes to peroxisomes. One of two nativeA. nidulansproteins with C-terminal SYM tripeptide shows weak peroxisomal localization alongside cytoplasmic presence, indicating that only a subset of proteins with non-canonical signals access peroxisomes.In silicoanalysis of 1,010 fungal genomes identified diverse SYM-proteins with variable functions, suggesting that non-canonical PTS-1 signals may evolve spontaneously. Two-thirds of SYM-proteins are predicted to localize to specific intracellular compartments other than the peroxisome. We propose that despite their predicted localization, these proteins possessing SYM as a non-canonical peroxisomal signal might also have peroxisomal presence. Among SYM-proteins, pectinesterases, known plant pathogen virulence factors, were frequent. Notably, 25% of fungal pectinesterases harbor non-canonical PTS-1 signals, suggesting that partial peroxisomal localization of pectinesterases has evolved convergently. This suggests that partial peroxisomal localization may enhance protein functional flexibility, contributing to the organism’s adaptability.

Science & Technology - Other Topics↗

The FIELDS Instrument Suite on MMS: Scientific Objectives, Measurements, and Data Products

The FIELDS instrumentation suite on the Magnetospheric Multiscale (MMS) mission provides comprehensive measurements of the full vector magnetic and electric fields in the reconnection regions investigated by MMS, including the dayside magnetopause and the night-side magnetotail acceleration regions out to 25 Re. Six sensors on each of the four MMS spacecraft provide overlapping measurements of these fields with sensitive cross-calibrations both before and after launch. The FIELDS magnetic sensors consist of redundant flux-gate magnetometers (AFG and DFG) over the frequency range from DC to 64 Hz, a search coil magnetometer (SCM) providing AC measurements over the full whistler mode spectrum expected to be seen on MMS, and an Electron Drift Instrument (EDI) that calibrates offsets for the magnetometers. The FIELDS three-axis electric field measurements are provided by two sets of biased double-probe sensors (SDP and ADP) operating in a highly symmetric spacecraft environment to reduce significantly electrostatic errors. These sensors are complemented with the EDI electric measurements that are free from all local spacecraft perturbations. Cross-calibrated vector electric field measurements are thus produced from DC to 100 kHz, well beyond the upper hybrid resonance whose frequency provides an accurate determination of the local electron density. Due to its very large geometric factor, EDI also provides very high time resolution (∼1 ms) ambient electron flux measurements at a few selected energies near 1 keV. This paper provides an overview of the FIELDS suite, its science objectives and measurement requirements, and its performance as verified in calibration and cross-calibration procedures that result in anticipated errors less than 0.1 nT in B and 0.5 mV/m in E. Summaries of data products that result from FIELDS are also described, as well as algorithms for cross-calibration. Details of the design and performance characteristics of AFG/DFG, SCM, ADP, SDP, and EDI are provided in five companion papers.

Electromagnetic field measurements↗

Epoxide Alcoholysis over M-BEA Zeolites: Effects of Alcohol Chain Length on Rates and Regioselectivities

The structures of nucleophilic reactants affect their coordination behavior among solvent molecules and kinetics of reactions with surface intermediates within the confines of fluidfilled pores of zeolites and other microporous materials. Consequently, rates and regioselectivities of diverse chemistries may depend sensitively on nucleophile identity in manners not observed for classic fluid phase reactions. Here, we examine the impact of varying the primary alcohol (ROH) chain length on the kinetics of 1,2-epoxybutane (C 4 H 8 O) ring-opening within Brønsted (Al-BEA) and Lewis acid (Zr-BEA) zeolites. Turnover rates increase by factors of ∼6 (Al-BEA) and 4-fold (Zr-BEA) between methanol and 1-hexanol, yet the reaction mechanisms remain comparable. Despite modest rate differences, apparent activation enthalpies calculated from rates and activities of solvated reactants decrease linearly by 12 (Al-BEA) to 33 kJ mol −1 (Zr- BEA) with increased proton affinity, which suggests bond formation energies for the nucleophile strongly influence rate increases. The molecular interpretation of these trends demonstrates, however, that the solvation of ring-opening transition states by zeolite pore structures and solvent molecules also governs rates. The impact of local solvating interactions appears most directly as changes in regioselectivities, which tend to enhance terminal alcohol formation with increasing ROH chain length. Regioselectivities largely do not vary with differences in fluid composition for a given ROH. The addition of H 2 O increases the number of hydrogen bonds among reactive species, and trends in regioselectivities imply that the decreased hydrogen bonding ability of longer chain ROH, and not the nucleophile strength or steric bulk, determines the regioselectivities of the resulting products. This work provides direct experimental evidence that nucleophilicity and hydrogen bonding influence reaction barriers and regioselectivities in zeolitecatalyzed epoxide ring-opening, offering pathways to better control reaction kinetics.

acidic zeolites↗

Arctic Shrub Expansion, Plant Functional Trait Variation, and Effects on Belowground Carbon Cycling (Final Technical Report)

Terrestrial ecosystems are undergoing dramatic changes in response to climate warming, and these changes are expected to feedback to the atmosphere, potentially altering the trajectory of future climate change. Feedbacks from Arctic ecosystems are a major concern because the Arctic is projected to warm significantly in the 21 st century and because >50% of global belowground organic carbon is stored in permafrost and overlying soils. Warming-driven release of this carbon could drastically increase atmospheric greenhouse gas concentrations and accelerate climate warming. Plant communities are also responding to warming, as evidenced by the widely documented increase in woody-shrub growth and “greening” across much of the Arctic tundra biome. This vegetation shift may offset or amplify warming by altering carbon cycling. The direction and magnitude of shrub effects remain highly uncertain, however, due to limited understanding of the consequences of shrub expansion for belowground carbon cycling and simplification of these relationships in models. The major shrubs expanding in the Arctic (Betula, Salix, and Alnus) vary widely with respect to aboveground and belowground traits (e.g., tissue production and chemistry, rooting depth, microbial symbionts), and may also exhibit substantial intraspecific variation in these traits in response to environmental conditions. Such variation is likely to have profound implications for soil carbon cycling. The overarching goal of this project was to improve process-based understanding of the influence of shrub expansion on carbon cycling to enable improved representation of carbon dynamics in ecosystem and Earth system models. We investigated how plant functional traits vary among shrub genera, respond to environmental conditions, and affect belowground carbon and nutrient cycling by quantifying relationships among functional traits and biogeochemical cycling along edaphic gradients nested within a climate gradient in the Alaskan tundra. We found consistent differences in leaf and root traits among shrub genera and between shrubs and a widespread sedge species, indicating diverse nutrient acquisition strategies and belowground impacts among different arctic shrubs. We also found striking differences in trait values among individuals within the same species or genera within sites. Soil parameters were more important than climate parameters for predicting size and leaf trait variation, and root trait responses were less dependent on climate overall. For all but one root trait, including parameters representing aboveground traits improved the predictive ability of models. These results demonstrate that tundra shrub traits vary considerably at local scales and soil factors drive this variation, especially belowground. Furthermore, leveraging information about aboveground traits and soil conditions can improve predictions of how belowground traits will respond to climate change. Despite these differences, soil carbon and nitrogen pools in the active layer did not vary among plots dominated by different shrub or sedge genera. Instead, pool sizes generally decreased from warmer to colder sites, consistent with a productivity gradient. Patterns of isotopic N composition indicate that shrubs tighten nitrogen cycling via nitrogen resorption or immobilization of shrub litter. Overall, these results suggest that further identifying the specific shrub genera in the tundra landscape will ultimately provide better predictions of belowground dynamics across the changing arctic. We also performed simulation experiments with the Terrestrial Ecosystem Model (TEM) incorporated in the Predictive Ecosystem Analyzer (PEcAn) framework, treats model parameters as probability distributions, estimates parameters based on a synthesis of available field data, and then quantifies both model sensitivity and uncertainty to a given parameter or suite of parameters. We performed simulations across different types of tundra, including shrub tundra. One key finding was that both model sensitivity and uncertainty to a given parameter could vary within the same type of tundra, but in a different geographical location, such as over the climate gradient of shrub tundra described above. We organized a special session at the annual meeting of the Ecological Society of America in August 2019 to disseminate our results, refine recommendations for model improvement, and initiate collaborations to implement these recommendations in existing models of tundra carbon dynamics at ecosystem to Earth system scales. Our results support DOE near-term priorities by providing mechanistic insights into the role of vegetation change in the terrestrial carbon cycle in a region that is inadequately represented in Earth system models. Current models reduce the complexity of Arctic vegetation to a small number of plant functional types (PFTs). This approach implicitly assumes that each PFT represents the average ecological function of its constituent species, thus ignoring the effects of trait variation on biogeochemical cycling and potentially leading to large uncertainty in the sign and magnitude of ecosystem feedbacks to climate. By quantifying variation of plant functional traits across broad gradients of climatic and edaphic conditions and elucidating the linkages of such variation with carbon and nutrient cycling, our results illustrate the need and create a foundation for further developing trait-based modeling approaches that allow the traits of PFTs to vary as a function of environmental conditions. These approaches should improve the capacity of simulation models to offer insights into ecosystem carbon dynamics associated with novel plant communities in a rapidly changing Arctic.

54 ENVIRONMENTAL SCIENCES↗

NMF-Based Anomaly Detection in CMS 2D Tracking Occupancy Histograms

The CMS experiment relies on Data Quality Monitoring (DQM) to ensure that recorded collision data are suitable for physics analysis. During LHC Run 3, each run contains many lumisections and tracking monitoring elements, making offline inspection challenging, especially for localized detector effects that may appear only for short periods of time. This poster presents an unsupervised machine-learning approach to identify anomalous lumisections in CMS tracking occupancy histograms using Non-Negative Matrix Factorization (NMF). The workflow uses offline CMS DQMIO tracking histograms retrieved with the CMS DIALS API and organized as two-dimensional occupancy maps for each lumisection. After selecting stable lumisections, the occupancy maps are normalized and arranged into a non-negative data matrix. The NMF model learns a compact set of basis patterns describing normal tracking occupancy. Each lumisection is then reconstructed from these learned components, and the reconstruction error is used as an anomaly score. Large residuals indicate occupancy patterns that deviate from normal detector behavior and are flagged for further inspection. This NMF-based approach provides a fast and interpretable way to flag lumisections whose tracking occupancy patterns differ from normal detector behavior. Preliminary studies show sensitivity to known tracking anomalies, and ongoing work is focused on validating the method across additional Run 3 Pixel and Strip detector issues.

Rodríguez Ramos, Iliomar [Puerto Rico U., Mayaguez↗

Low- n stability and plasma response to RMP in various STEP scenarios

The low-n (n is the toroidal mode number) magnetohydrodynamic (MHD) stability and plasma response are numerically investigated for various scenarios designed for STEP, that are relevant for the H-mode pedestal analysis. Control of the edge-localized modes (ELMs) with externally applied resonant magnetic perturbations (RMPs) is considered. Optimization of the ELM control coil current configuration, based on the computed plasma MHD response and well-established figures of merit validated on present-day experiments, finds reasonable robustness of a fixed coil phasing (for a given n-number) to control ELMs in all five STEP plasmas considered. Based on certain semi-empirical criteria, the required coil current to achieve ELM suppression is estimated to be about 10–20 kAt with the n = 1 or 2 RMP configuration and about 100–200 kAt for the n = 4 RMP. Systematic linear stability calculations are used to map out stability windows for the low-n kink-peeling modes, in terms of the ideal-wall location and variation of the edge safety factor q 95 with respect to the target design. The kink-peeling stability boundary is found to be generally sensitive to the q 95 variation, which has implications for achieving the quiescent H-mode regime in STEP. Full toroidal quasilinear initial-value simulations for these STEP plasmas find that generation of the edge-harmonic oscillations (EHOs) depends sensitively on the plasma scenario, the initial linear stability of the kink-peeling modes, the initial plasma toroidal flow and q 95 . In general, it is easier (more robust) to access the EHO-regime for two of the cases considered with smaller plasma volume and higher on-axis safety factor. Finally, quasilinear simulations find robust density pumpout due to applied RMPs in these STEP plasmas, but the effect on the plasma toroidal flow varies among different cases.

EHO↗