Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “constrained regions”

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 163 records · Page 9

Ensemble‐Based Spatially Distributed CLM5 Hydrological Parameter Estimation for the Continental United States

Abstract One of the major challenges in large‐domain hydrological modeling efforts lies in the estimation of spatially distributed hydrological parameters while simultaneously accounting for their associated uncertainties. Addressing this challenge is particularly difficult in ungauged locations. With growing societal demands for large‐scale streamflow projections to inform water resource management and long‐term planning, evaluating and constraining hydrological parameter uncertainty is increasingly vital. This study introduces a hybrid regionalization approach to enhance hydrological predictions of the Community Land Model version 5 (CLM5) across the Continental United States (CONUS), with a total of 50,629 1/8° grid cells. This hybrid method combines the strengths of two existing techniques: parameter regionalization and streamflow signature regionalization. It identifies ensemble behavioral parameters for each 1/8° grid cell across the CONUS domain, tailored to three distinct streamflow signatures focused on low flows, high flows, and annual water balance. Evaluating this hybrid method for 464 CAMELS (Catchment Attributes and Meteorology for Large‐sample Studies) basins demonstrates a significant improvement in CLM5 hydrological predictions, even in challenging arid regions. In CONUS applications, the derived spatially distributed parameter sets capture both spatial continuity and variation of parameters, highlighting their heterogeneous nature within specific regions. Overall, this hybrid regionalization approach offers a promising solution to the complex task of improving hydrological modeling over large domains for important hydrological applications.

CLM5↗

Holography of the photon ring

Space-based next-generation interferometers propose to measure the Lyapunov exponents of the nearly bound geodesics that comprise the photon ring surrounding the black hole M87*. We argue that these classical Lyapunov exponents equal the quantum Ruelle resonances describing the late-time approach to thermal equilibrium of the quantum microstate holographically dual to any Kerr black hole such as M87*. Moreover, we identify 'near-ring regions' in the phase space of fields propagating on Kerr that exhibit critical behavior, including emergent conformal symmetries. These are analogues for sub-extremal Kerr of the much-studied 'near-horizon regions' of (near-)extremal black holes. Here, the emergent conformal symmetries greatly constrain the observational predictions for the fine photon ring substructure around M87* and for quasinormal gravitational-wave ringdowns, as well as any proposal for a quantum holographic dual to the Kerr black hole. More generally, we hope that our identification of several universal features of Kerr spectroscopy provides a useful starting point for a bottom-up approach to holography for astrophysical black holes.

79 ASTRONOMY AND ASTROPHYSICS↗

Evolutionary constraints and climate variability jointly shape starch–sugar balance in woody plants

Nonstructural carbohydrates (NSC) buffer plants against carbon imbalances, yet their partitioning between storage and soluble pools remains elusive at global scales. Here, we compiled a dataset of starch to soluble sugar ratio (St : Su) for 308 woody species across 220 sites world-wide and introduce a dimensionless index that integrates storage and demand while minimizing methodological artifacts. St : Su was strongly associated with growth, identifying it as a key axis of carbon allocation. Foliage consistently exhibited lower St : Su than lignified organs, reflecting a division between transient and conservative pools. Conifers accumulated more starch in foliage but less in stems relative to angiosperms, while leaf habits and mycorrhizal associations further modulated organ-specific strategies. Contrary to expectation, foliar and root St : Su varied little among biomes, but stems exhibited higher ratios in tropical rainforests than in boreal or arid regions, reflecting differences in species composition and adaptive storage under disturbance. Phylogeny constrained stem storage, whereas climatic variability, rather than mean conditions, dominated allocation in leaves and roots. These findings establish St : Su as a robust functional trait linking allocation strategies, growth, and resilience, which can be used to improve vegetation model prediction of forest productivity and mortality under climate variability.

Li, Weibin [Lanzhou Univ. (China)] (ORCID:00000001↗

Global Methane Budget 2000–2020

Abstract. Understanding and quantifying the global methane (CH4) budget is important for assessing realistic pathways to mitigate climate change. CH4 is the second most important human-influenced greenhouse gas in terms of climate forcing after carbon dioxide (CO2), and both emissions and atmospheric concentrations of CH4 have continued to increase since 2007 after a temporary pause. The relative importance of CH4 emissions compared to those of CO2 for temperature change is related to its shorter atmospheric lifetime, stronger radiative effect, and acceleration in atmospheric growth rate over the past decade, the causes of which are still debated. Two major challenges in quantifying the factors responsible for the observed atmospheric growth rate arise from diverse, geographically overlapping CH4 sources and from the uncertain magnitude and temporal change in the destruction of CH4 by short-lived and highly variable hydroxyl radicals (OH). To address these challenges, we have established a consortium of multidisciplinary scientists under the umbrella of the Global Carbon Project to improve, synthesise, and update the global CH4 budget regularly and to stimulate new research on the methane cycle. Following Saunois et al. (2016, 2020), we present here the third version of the living review paper dedicated to the decadal CH4 budget, integrating results of top-down CH4 emission estimates (based on in situ and Greenhouse Gases Observing SATellite (GOSAT) atmospheric observations and an ensemble of atmospheric inverse-model results) and bottom-up estimates (based on process-based models for estimating land surface emissions and atmospheric chemistry, inventories of anthropogenic emissions, and data-driven extrapolations). We present a budget for the most recent 2010–2019 calendar decade (the latest period for which full data sets are available), for the previous decade of 2000–2009 and for the year 2020. The revision of the bottom-up budget in this 2025 edition benefits from important progress in estimating inland freshwater emissions, with better counting of emissions from lakes and ponds, reservoirs, and streams and rivers. This budget also reduces double counting across freshwater and wetland emissions and, for the first time, includes an estimate of the potential double counting that may exist (average of 23 Tg CH4 yr−1). Bottom-up approaches show that the combined wetland and inland freshwater emissions average 248 [159–369] Tg CH4 yr−1 for the 2010–2019 decade. Natural fluxes are perturbed by human activities through climate, eutrophication, and land use. In this budget, we also estimate, for the first time, this anthropogenic component contributing to wetland and inland freshwater emissions. Newly available gridded products also allowed us to derive an almost complete latitudinal and regional budget based on bottom-up approaches. For the 2010–2019 decade, global CH4 emissions are estimated by atmospheric inversions (top-down) to be 575 Tg CH4 yr−1 (range 553–586, corresponding to the minimum and maximum estimates of the model ensemble). Of this amount, 369 Tg CH4 yr−1 or ∼ 65 % is attributed to direct anthropogenic sources in the fossil, agriculture, and waste and anthropogenic biomass burning (range 350–391 Tg CH4 yr−1 or 63 %–68 %). For the 2000–2009 period, the atmospheric inversions give a slightly lower total emission than for 2010–2019, by 32 Tg CH4 yr−1 (range 9–40). The 2020 emission rate is the highest of the period and reaches 608 Tg CH4 yr−1 (range 581–627), which is 12 % higher than the average emissions in the 2000s. Since 2012, global direct anthropogenic CH4 emission trends have been tracking scenarios that assume no or minimal climate mitigation policies proposed by the Intergovernmental Panel on Climate Change (shared socio-economic pathways SSP5 and SSP3). Bottom-up methods suggest 16 % (94 Tg CH4 yr−1) larger global emissions (669 Tg CH4 yr−1, range 512–849) than top-down inversion methods for the 2010–2019 period. The discrepancy between the bottom-up and the top-down budgets has been greatly reduced compared to the previous differences (167 and 156 Tg CH4 yr−1 in Saunois et al. (2016, 2020) respectively), and for the first time uncertainties in bottom-up and top-down budgets overlap. Although differences have been reduced between inversions and bottom-up, the most important source of uncertainty in the global CH4 budget is still attributable to natural emissions, especially those from wetlands and inland freshwaters. The tropospheric loss of methane, as the main contributor to methane lifetime, has been estimated at 563 [510–663] Tg CH4 yr−1 based on chemistry–climate models. These values are slightly larger than for 2000–2009 due to the impact of the rise in atmospheric methane and remaining large uncertainty (∼ 25 %). The total sink of CH4 is estimated at 633 [507–796] Tg CH4 yr−1 by the bottom-up approaches and at 554 [550–567] Tg CH4 yr−1 by top-down approaches. However, most of the top-down models use the same OH distribution, which introduces less uncertainty to the global budget than is likely justified. For 2010–2019, agriculture and waste contributed an estimated 228 [213–242] Tg CH4 yr−1 in the top-down budget and 211 [195–231] Tg CH4 yr−1 in the bottom-up budget. Fossil fuel emissions contributed 115 [100–124] Tg CH4 yr−1 in the top-down budget and 120 [117–125] Tg CH4 yr−1 in the bottom-up budget. Biomass and biofuel burning contributed 27 [26–27] Tg CH4 yr−1 in the top-down budget and 28 [21–39] Tg CH4 yr−1 in the bottom-up budget. We identify five major priorities for improving the CH4 budget: (i) producing a global, high-resolution map of water-saturated soils and inundated areas emitting CH4 based on a robust classification of different types of emitting ecosystems; (ii) further development of process-based models for inland-water emissions; (iii) intensification of CH4 observations at local (e.g. FLUXNET-CH4 measurements, urban-scale monitoring, satellite imagery with pointing capabilities) to regional scales (surface networks and global remote sensing measurements from satellites) to constrain both bottom-up models and atmospheric inversions; (iv) improvements of transport models and the representation of photochemical sinks in top-down inversions; and (v) integration of 3D variational inversion systems using isotopic and/or co-emitted species such as ethane as well as information in the bottom-up inventories on anthropogenic super-emitters detected by remote sensing (mainly oil and gas sector but also coal, agriculture, and landfills) to improve source partitioning. The data presented here can be downloaded from https://doi.org/10.18160/GKQ9-2RHT (Martinez et al., 2024).

54 ENVIRONMENTAL SCIENCES↗

Measurement of Charged-Pion Production in Deep-Inelastic Scattering off Nuclei with the CLAS Detector

Background: Energetic quarks in nuclear deep-inelastic scattering propagate through the nuclear medium. Processes that are believed to occur inside nuclei include quark energy loss through medium-stimulated gluon bremsstrahlung and intranuclear interactions of forming hadrons. More data are required to gain a more complete understanding of these effects. Purpose: To test the theoretical models of parton transport and hadron formation, we compared their predictions for the nuclear and kinematic dependence of pion production in nuclei. Methods: We have measured charged-pion production in semi-inclusive deep-inelastic scattering off D, C, Fe, and Pb using the CLAS detector and the CEBAF 5.014-GeV electron beam. We report results on the nuclear-to-deuterium multiplicity ratio for π + and π - as a function of energy transfer, four-momentum transfer, and pion energy fraction or transverse momentum—the first three-dimensional study of its kind. Results: The π + multiplicity ratio is found to depend strongly on the pion fractional energy z and reaches minimum values of 0.67 ± 0.03, 0.43 ± 0.02, and 0.27 ± 0.01 for the C, Fe, and Pb targets, respectively. The z dependencies of the multiplicity ratios for π + and π - are equal within uncertainties for C and Fe targets but show differences at the level of 10% for the Pb-target data. The results are qualitatively described by the GIBUU transport model, as well as with a model based on hadron absorption, but are in tension with calculations based on nuclear fragmentation functions. Conclusions: These precise results will strongly constrain the kinematic and flavor dependence of nuclear effects in hadron production, probing an unexplored kinematic region. They will help to reveal how the nucleus reacts to a fast quark, thereby shedding light on its color structure and transport properties and on the mechanisms of the hadronization process.

Moran, S.↗

A system of shared autonomous vehicles for Chicago: Understanding the effects of geofencing the service

With autonomous vehicles (AVs) still in the testing phase, researchers and planners must resort to simulation techniques to explore possible futures regarding shared and automated mobility. An agent-based discrete-event transport simulator, POLARIS, is used in this study to simulate travel in the 20-county Chicago region with a shared AV (SAV) mobility option. Using this framework, the effect of an SAV fleet on system performance when constrained to serve within geofences is studied under four distinct scenarios: service restricted to the city, to the city plus suburban core, to the core plus exurban areas, and to the entire region — along with the choice of dynamic ridesharing (DRS) versus solo travel in an SAV. Results indicate that service areas need a balanced mix of trip generators and attractors, and an SAV fleet’s empty VMT (eVMT) can be noticeably reduced through suitable geofencing and DRS. Geofences can also help lower response times, reduce systemwide VMT across all modes, and ensure uniform access to SAVs. DRS is most useful in lowering VMT and %eVMT that arises from sprawled land development, but with insufficient demand to share rides, savings from the use of geofences is higher. Geofences targeting neighborhoods with high trip density bring about low response times and %eVMT, but fleet sizes in these regions need to be designed for uniformly low response times throughout a large region, as opposed to maximizing vehicle use in a 24-hour day.

33 ADVANCED PROPULSION SYSTEMS↗

Embedded Sensors in Additively Manufactured Silicon Carbide

Silicon carbide (SiC) components are being considered for a wide range of nuclear applications due to their high-temperature strength retention, low neutron absorption, chemical inertness, and dimensional stability under neutron irradiation. However, machining and joining of SiC components have traditionally limited its application to relatively simple geometries. Recent work has demonstrated additive manufacturing of complex, high-purity, crystalline SiC components using a combination of binder jet printing and densification via chemical vapor infiltration (CVI). The process lends itself to embedding of fuel, absorbers, moderators, and sensors at strategic locations within a component. The latter could allow for enhanced in situ performance monitoring of limiting fuel temperatures, self-shielded neutron flux, and potentially spatially distributed strain within complex SiC components if sensors can be successfully embedded during CVI. This work describes (1) methods for embedding sensors; (2) thermodynamic analyses and material compatibility testing for identifying sensors capable of surviving high temperatures and exposure to hydrogen and hydrogen chloride during CVI; and (3) nuclear applications for embedded sensors, including potential failure modes during fabrication and during reactor operation. Molybdenum-sheathed thermocouples were successfully embedded in a complex SiC component, whereas niobium-sheathed high-temperature irradiation-resistant thermocouples started to drift as soon as the reactant gases were introduced and ultimately failed during CVI due to severe constrained expansion, potentially resulting from niobium hydride formation in the low-temperature region of the CVI system. Optical fibers were successfully embedded in SiC, but further work is needed to protect the fragile fiber leads after their protective coatings are removed during CVI.

thermocouples↗

Wind‐Associated Melt Trends and Contrasts Between the Greenland and Antarctic Ice Sheets

Föhn and katabatic winds (downslope winds) can increase ice sheet surface melt, run-off, and ice-shelf vulnerability to hydrofracture and are poorly constrained on the Greenland and Antarctic ice sheets (GIS and AIS). We use regional climate model simulations of the GIS and AIS to quantify and intercompare trends in downslope winds and associated melt since 1960. Results reveal surface melt associated with downslope wind is significant on both the GIS and AIS representing 27.5 ± 4.5% and 19.7 ± 3.8% of total surface melt respectively. Wind-associated melt has decreased 31.8 ± 5.3% on the AIS while total melt decreased 15.4 ± 2.4% due to decreased föhn-induced melt on the Antarctic Peninsula and increasing stratospheric ozone. Wind-associated melt has increased 10.3 ± 2.5% on the GIS, combining with a more positive North Atlantic Oscillation and warmer surface to increase total melt 34 ± 5.8%.

58 GEOSCIENCES↗

ERF: Energy Research and Forecasting Model

High performance computing (HPC) architectures have undergone rapid development in recent years. As a result, established software suites face an ever increasing challenge to remain performant on and portable across modern systems. Many of the widely adopted atmospheric modeling codes cannot fully (or in some cases, at all) leverage the acceleration provided by General-Purpose Graphics Processing Units, leaving users of those codes constrained to increasingly limited HPC resources. Energy Research and Forecasting (ERF) is a regional atmospheric modeling code that leverages the latest HPC architectures, whether composed of only Central Processing Units (CPUs) or incorporating GPUs. ERF contains many of the standard discretizations and basic features needed to model general atmospheric dynamics. The modular design of ERF provides a flexible platform for exploring different physics parameterizations and numerical strategies. ERF is built on a state-of-the-art, well-supported, software framework (AMReX) that provides a performance portable interface and ensures ERF's long-term sustainability on next generation computing systems. This paper details the numerical methodology of ERF, presents results for a series of verification/validation cases, and documents ERF's performance on current HPC systems. The roughly 5× speed up of ERF (using GPUs) over Weather Research and Forecasting (CPUs only) for a 3D squall line test case highlights the significance of leveraging GPU acceleration.

17 WIND ENERGY↗

COMBLE ARM Mobile Facility (AMF) Measurements of Ice Nucleating Particles Field Campaign Report

During the Cold-Air Outbreaks in the Marine Boundary Layer Experiment (COMBLE), a project with the overarching goal to shed light on the dynamics and microphysical properties of clouds and precipitation in the high-latitude marine boundary layer during cold-air outbreaks (CAOs), we were tasked with providing and assisting the collection of aerosol filter samples for measuring ice nucleating particle (INP) concentrations. Little is known about the properties of CAO cloud systems, how they vary with surface, environmental, and aerosol conditions, the role of cold-air outbreaks in the global atmospheric and ocean circulation, and the accuracy of the treatment of this atmospheric regime in climate models. Uncertainties in climate projections for the arctic and sub-arctic region are large, in part due to the lack of observational guidance to constrain the treatment of aerosol-cloud-precipitation linkages in climate models.

54 ENVIRONMENTAL SCIENCES↗

Dark Photons from Perturbative Decay of a Misaligned Higgs Field

We reconsider the production of dark photons $A'$ as dark matter, from the perturbative decay of a dark Higgs field $h$, that is stochastically misaligned from the minimum of its potential during inflation. This is a simple and predictive framework for generating the $A'$ relic abundance. It is constrained by structure formation, since the $A'$ are initially boosted, and inflationary isocurvature fluctuations, which require small quartic couplings $λh^4$. We identify $A'$ masses between 100 eV and 1 GeV and gauge couplings $g\sim 10^{-15}-10^{-10}$ that are consistent in this scenario, and which become more tightly constrained if a generic level of kinetic mixing is present. The favored parameter region could be tested through future CMB or Lyman-$α$ observations, and, in the presence of kinetic mixing, by direct detection experiments or diffuse soft gamma-ray searches.

Cline, James M. [McGill U., Montreal (main)]↗

Polymorphism of gold under laser-based ramp compression to 690 GPa

The high-pressure behavior and polymorphism of gold has attracted strong recent interest from theory as well as from both static and dynamic experiments. In this study, gold was dynamically ramp-compressed to 690 GPa via laser ablation using the Omega Laser Facility (U. Rochester) in order to investigate the stability regions of this noble metal's polymorphs. Here, stresses were determined from measurements of surface or interface velocities by laser interferometry. In situ X-ray diffraction was carried out to constrain the structure and density. The ambient face-centered-cubic (fcc) phase was observed at 162 GPa, followed by a mixed phase region consisting of the fcc and body-centeredcubic (bcc). The presence of only the bcc phase was observed between 377 and 690 GPa. Our results are consistent with recent shock compression experiments and demonstrate that gold transforms to the bcc phase under conditions of both ramp and shock loading.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Evaluating seasonal and regional distribution of snowfall in regional climate model simulations in the Arctic

In this study, we investigate how the regional climate model HIRHAM5 reproduces the spatial and temporal distribution of Arctic snowfall when compared to CloudSat satellite observations during the examined period of 2007–2010. For this purpose, both approaches, i.e., the assessments of the surface snowfall rate (observation-to-model) and the radar reflectivity factor profiles (model-to-observation), are carried out considering spatial and temporal sampling differences. The HIRHAM5 model, which is constrained in its synoptic representation by nudging to ERA-Interim, represents the snowfall in the Arctic region well in comparison to CloudSat products. The spatial distribution of the snowfall patterns is similar in both identifying the southeastern coast of Greenland and the North Atlantic corridor as regions gaining more than twice as much snowfall as the Arctic average, defined here for latitudes between 66 and 81°N. Excellent agreement (difference less than 1%) in the Arctic-averaged annual snowfall rate between HIRHAM5 and CloudSat is found, whereas ERA-Interim reanalysis shows an underestimation of 45% and significant deficits in the representation of the snowfall rate distribution. From the spatial analysis, it can be seen that the largest differences in the mean annual snowfall rates are an overestimation near the coastlines of Greenland and other regions with large orographic variations as well as an underestimation in the northern North Atlantic Ocean. To a large extent, the differences can be explained by clutter contamination, blind zone or higher resolution of CloudSat measurements, but clearly HIRHAM5 overestimates the orographic-driven precipitation. The underestimation of HIRHAM5 within the North Atlantic corridor south of Svalbard is likely connected to a poor description of the marine cold air outbreaks which could be identified by separating snowfall into different circulation weather type regimes. By simulating the radar reflectivity factor profiles from HIRHAM5 utilizing the Passive and Active Microwave TRAnsfer (PAMTRA) forward-modeling operator, the contribution of individual hydrometeor types can be assessed. Looking at a latitude band at 72–73°N, snow can be identified as the hydrometeor type dominating radar reflectivity factor values across all seasons. The largest differences between the observed and simulated reflectivity factor values are related to the contribution of cloud ice particles, which is underestimated in the model, most likely due to the small sizes of the particles. The model-to-observation approach offers a promising diagnostic when improving cloud schemes, as illustrated by comparison of different schemes available for HIRHAM5.

54 ENVIRONMENTAL SCIENCES↗

Simulation of Physics-Based 0-10Hz Strong Motion Using High Performance Computing Supporting Refinements to Regional Ground Motion Models for the Central Eastern US

In collaboration with the U.S. Nuclear Regulatory Commission (NRC) the LLNL has developed a computationally efficient simulation platform designed to perform physics-based ground motion simulations for crustal earthquakes in the Stable Continental Regions of Central and Eastern US (CEUS), using high-performance computing. The main objective of the earthquake simulations was to use synthetic ground motion to provide constrains to refinements of existing ergodic Ground Motion Models (GMMs), for large magnitude earthquakes and near-fault distances, for which these models are less reliable. Physics-based broadband (0-10Hz) ground motion simulations were used to estimate the near-fault ground motion amplitudes and within event and between-event variabilities associated with fault rupture characteristics. In our simulations we used a 3D regional velocity model that was based on Saikia’s 1D velocity model (1994). In simulations performed during the first stage of this project the Saikia’s velocity model demonstrated better performance in modelling high frequency regional wave propagation for the CEUS region recorded during the Mw5.0 November 7, 2016, Cushing Oklahoma (Taylor et al., 2017), and Mw5.8 September 3, 2016, Pawnee Oklahoma earthquakes. The proposed regional 3D model includes random perturbations to the 1D background model using the stochastic scheme of Pitarka and Mellors (2021). In addition, validation analysis of the rupture generator and regional wave propagation models, using comparisons with different GMMs for Mw6.5 and Mw7.0 scenario earthquakes in the CEUS region resulted in a very good match between the simulated and empirical ground motion models. For the purposes of seismic hazard assessment at the existing and planned nuclear power plants, NRC is interested in studies aimed at improving the current ground motion models (GMM) for both Stable Continental Regions (SCR) in the Central and Eastern US and Active Crustal Regions (ACR) in the Western US. Due to lack of recorded data, these improvements require synthetic data for short fault distances and large magnitude earthquakes for which the existing recorded data is not enough to uniquely constrain the GMMs. The need for simulations and strong motion data is especially critical for the CEUS region where we do not have recorded data from potentially large damaging earthquakes with moment magnitudes 6.0 and higher. In this the project, we focused on 10Hz simulations of Mw7.0 scenario earthquakes with strike slip and thrust faulting mechanisms. We used more than 50 Mw7.0 earthquake rupture scenarios to investigate the ground motion uncertainty due to unknown earthquake rupture parameters, in particular, the slip distribution, rupture velocity, and faulting mechanism, and their implication on ground motion amplification due to forward rupture directivity effects.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Hard X-Ray Emission from the Eastern Jet of SS 433 Powering the W50 “Manatee” Nebula: Evidence for Particle Reacceleration

We present a broadband X-ray study of W50 (the “Manatee” nebula), the complex region powered by the microquasar SS 433, that provides a test bed for several important astrophysical processes. The W50 nebula, a Galactic PeVatron candidate, is classified as a supernova remnant but has an unusual double-lobed morphology likely associated with the jets from SS 433. Using NuSTAR, XMM-Newton, and Chandra observations of the inner eastern lobe of W50, we have detected hard nonthermal X-ray emission up to ~30 keV, originating from a few-arcminute-sized knotty region (“Head”) located ≲18' (29 pc for a distance of 5.5 kpc) east of SS 433, and constrained its photon index to 1.58 ± 0.05 (0.5–30 keV band). The index gradually steepens eastward out to the radio “ear” where thermal soft X-ray emission with a temperature kT ~ 0.2 keV dominates. The hard X-ray knots mark the location of acceleration sites within the jet and require an equipartition magnetic field of the order of ≳12 μG. The unusually hard spectral index from the “Head” region challenges classical particle acceleration processes and points to particle injection and reacceleration in the subrelativistic SS 433 jet, as seen in blazars and pulsar wind nebulae.

79 ASTRONOMY AND ASTROPHYSICS↗

Comparison of Conventional and Constrained Variational Methods for Computing Large-Scale Budgets and Forcing Fields

Analyses of atmospheric heat and moisture budgets serve as an effective tool to study convective characteristics over a region and to provide large-scale forcing fields for various modeling applications. This paper examines two popular methods for computing large-scale atmospheric budgets: the conventional budget method (CBM) using objectively gridded analyses based primarily on radiosonde data and the constrained variational analysis (CVA) approach which supplements vertical profiles of atmospheric fields with measurements at the top of the atmosphere and at the surface to conserve mass, water, energy, and momentum. Successful budget computations are dependent on accurate sampling and analyses of the thermodynamic state of the atmosphere and the divergence field associated with convection and the large-scale circulation that influences it. Utilizing analyses generated from data taken during Dynamics of the Madden-Julian Oscillation (DYNAMO) field campaign conducted over the central Indian Ocean from October to December 2011, we evaluate the merits of these budget approaches and examine their limitations. While many of the shortcomings of the CBM, in particular effects of sampling errors in sounding data, are effectively minimized with CVA, accurate large-scale diagnostics in CVA are dependent on reliable background fields and rainfall constraints. Finally, for the DYNAMO analyses examined, the operational model fields used as the CVA background state provided wind fields that accurately resolved the vertical structure of convection in the vicinity of Gan Island. However, biases in the model thermodynamic fields were somewhat amplified in CVA resulting in a convective environment much weaker than observed.

54 ENVIRONMENTAL SCIENCES↗

The eROSITA Final Equatorial-Depth Survey (eFEDS)

Here, we investigate the physical properties – such as the stellar mass (M * ), star-formation rate, infrared (IR) luminosity (L IR ), X-ray luminosity (L X ), and hydrogen column density (N H ) – of mid-IR (MIR) galaxies and active galactic nuclei (AGN) at z < 4 in the 140 deg 2 field observed by eROSITA on SRG using the Performance-and-Verification-Phase program named the eROSITA Final Equatorial Depth Survey (eFEDS). By cross-matching the WISE 22 μm (W4)-detected sample and the eFEDS X-ray point-source catalog, we find that 692 extragalactic objects are detected by eROSITA. We have compiled a multiwavelength dataset extending from X-ray to far-IR wavelengths. We have also performed (i) an X-ray spectral analysis, (ii) spectral-energy-distribution fitting using X-CIGALE, (iii) 2D image-decomposition analysis using Subaru Hyper Suprime-Cam images, and (iv) optical spectral fitting with QSFit to investigate the AGN and host-galaxy properties. For 7088 WISE 22 μm objects that are undetected by eROSITA, we have performed an X-ray stacking analysis to examine the typical physical properties of these X-ray faint and probably obscured objects. We find that (i) 82% of the eFEDS–W4 sources are classified as X-ray AGN with log L X > 42 erg s –1 ; (ii) 67 and 24% of the objects have log(L IR /L ⊙ ) > 12 and 13, respectively; (iii) the relationship between L X and the 6 μm luminosity is consistent with that reported in previous works; and (iv) the relationship between the Eddington ratio and N H for the eFEDS–W4 sample and a comparison with a model prediction from a galaxy-merger simulation indicates that approximately 5.0% of the eFEDS–W4 sources in our sample are likely to be in an AGN-feedback phase, in which strong radiation pressure from the AGN blows out the surrounding material from the nuclear region. Thanks to the wide area coverage of eFEDS, we have been able to constrain the ranges of the physical properties of the WISE 22 μm-selected sample of AGNs at z < 4, providing a benchmark for forthcoming studies on a complete census of MIR galaxies selected from the full-depth eROSITA all-sky survey.

79 ASTRONOMY AND ASTROPHYSICS↗

Variance-Reduced Accelerated First-Order Methods: Central Limit Theorems and Confidence Statements

In this paper, we consider a strongly convex stochastic optimization problem and propose three classes of variable sample-size stochastic first-order methods: (i) the standard stochastic gradient descent method, (ii) its accelerated variant, and (iii) the stochastic heavy-ball method. In each scheme, the exact gradients are approximated by averaging across an increasing batch size of sampled gradients. We prove that when the sample size increases at a geometric rate, the generated estimates converge in mean to the optimal solution at an analogous geometric rate for schemes (i)–(iii). Based on this result, we provide central limit statements, whereby it is shown that the rescaled estimation errors converge in distribution to a normal distribution with the associated covariance matrix dependent on the Hessian matrix, the covariance of the gradient noise, and the step length. If the sample size increases at a polynomial rate, we show that the estimation errors decay at a corresponding polynomial rate and establish the associated central limit theorems (CLTs). Under certain conditions, we discuss how both the algorithms and the associated limit theorems may be extended to constrained and nonsmooth regimes. As a result, we provide an avenue to construct confidence regions for the optimal solution based on the established CLTs and test the theoretical findings on a stochastic parameter estimation problem.

Lei, Jinlong↗