Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spatial constraints”

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 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Multiscale design of nonlinear materials using a Eulerian shape optimization scheme

Motivated by recent advances in manufacturing, the design of materials is the focal point of interest in the material research community. One of the critical challenges in this field is finding optimal material microstructure for a desired macroscopic response. This work presents a computational method for the mesoscale-level design of particulate composites for an optimal macroscale-level response. The method relies on a custom shape optimization scheme to find the extrema of a nonlinear cost function subject to a set of constraints. Three key “modules” constitute the method: multiscale modeling, sensitivity analysis, and optimization. Multiscale modeling relies on a classical homogenization method and a nonlinear NURBS-based generalized finite element scheme to efficiently and accurately compute the structural response of particulate composites using a nonconformal discretization. A three-parameter isotropic damage law is used to model microstructure-level failure. An analytical sensitivity method is developed to compute the derivatives of the cost/constraint functions with respect to the design variables that control the microstructure's geometry. The derivation uncovers subtle but essential new terms contributing to the sensitivity of finite element shape functions and their spatial derivatives. Several structural problems are solved to demonstrate the applicability, performance, and accuracy of the method for the design of particulate composites with a desired macroscopic nonlinear stress-strain response.

42 ENGINEERING↗
NASA NTRS2003

Enhanced Mixing in a Rectangular Duct

An experimental investigation of the mixing of non-reacting opposed rows of jets injected normal to a confined rectangular crossflow has been conducted. Planar Mie-scattering was used to measure the time-average concentration distribution of the jet fluid in planes perpendicular to the duct axis. Particular emphasis was placed on the study of closely spaced orifice configurations applicable to the mixing zone of an RQL combustor. Baseline studies were performed of mixing under "ideal" conditions, i.e., plenum fed jets injecting into a crossflow uniform in velocity and turbulence intensity. In addition, more practical ("non-ideal") issues encountered during hardware design were also studied. As in other studies, mixing effectiveness, determined using a spatial unmixedness parameter based on the variance of mean jet concentration distributions, was found to be optimum when the spacing-to-duct-height ratio was inversely proportional to the square root of the jet-to-mainstream momentum-flux ratio. This relationship is suitable for design under ideal flow conditions. Inlet flow boundary conditions of the jet and approach flow (mainstream) were found to strongly influence mixing performance, but no attempt was made to determine optimum performance under non-ideal conditions. The tests performed do offer some guidance as to expected mixing behavior for several common variables likely to be imposed by hardware constraints. Additionally, in this study it was found that for rows of orifices with opposite centerlines inline, mixing was similar for blockages up to 89 percent (previous crossflow mixing studies concerned with dilution zone configurations, blockages were typically less than 50 percent). Lower levels of unmixedness were obtained as a function of downstream location when axial injection length was minimized. Mixing may be enhanced if orifice centerlines of opposed rows are staggered, but blockage must be =50 percent in this configuration. Round hole and "square" orifice shapes had similar performance. Other variations in orifice shape did not substantially augment overall mixing performance. Furthermore an isothermal mixing data set was generated and used by CFDRC as input to a NO(x) inference code.

Liscinsky, D. S.↗
DOE OSTI2021

Simulation of radon-222 with the GEOS-Chem global model: Emissions, seasonality, and convective transport

Radon-222 (222Rn) is a short-lived radioactive gas naturally emitted from land surface, and has long been used to assess convective transport in atmospheric models. In this study, we simulate 222Rn using the GEOS-Chem chemical transport model with aims to improve our understanding of 222Rn emissions and surface concentration seasonality, and to characterize convective transport associated with two Goddard Earth Observing System (GEOS) meteorological products, MERRA and GEOS-FP. We evaluate four available 222Rn emission scenarios by comparing model results with surface observations at 51 global surface sites. The default emission scenario in GEOS-Chem yields a moderate agreement with global surface observations (< 70% data within a factor of 2) and a large underestimate of wintertime surface 222Rn concentrations at Northern Hemisphere mid- and high-latitudes due to an oversimplified formulation of 222Rn emission fluxes (1 atom cm-2 s-1 over land with a reduction by a factor of 3 under freezing conditions). We compose a new global 222Rn emission scenario based on Zhang et al. (2011) and show its potential to improve simulated surface 222Rn concentrations and seasonality. The regional components of this emission scenario include spatially and temporally varying emission fluxes derived from previous measurements of soil radium content and soil exhalation models, which are key factors to 222Rn emission flux rates. However, large model underestimates of surface 222Rn concentrations still exist in Asia, suggesting unusually high regional 222Rn emissions. We propose a conservative up-scaling factor of 1.2 for 222Rn emission fluxes in China, as also constrained by the observed deposition fluxes of 210Pb (decay daughter of 222Rn). With this modification, the model shows better agreement with the observations in Europe and North America (>80% data within a factor of 2), and reasonable agreement in Asia (close to 70%). Further constraints on 222Rn emissions would require additional observations of surface 222Rn concentrations and emission fluxes in central U.S., Canada, Africa, and Asia. We also compare and assess convective transport in model simulations driven by MERRA and GEOS-FP using observed 222Rn vertical profiles during northern mid-latitude summertime and from two short-term airborne campaigns. While the simulations with both GEOS products are able to capture the observed vertical gradient of 222Rn concentrations in the lower troposphere (0-4 km), neither correctly represents the level of convective detrainment, resulting in biases in the middle and upper troposphere. Compared to GEOS-FP, MERRA leads to stronger convective transport of 222Rn, which is partially compensated by its weaker large-scale vertical advection, resulting in similar global vertical distributions of 222Rn concentrations between the two simulations.

Zhang, Bo↗
DOE OSTI2025

Steady-state mixing state of black carbon aerosols from a particle-resolved model

Abstract. Black carbon (BC) exerts a notable warming effect due to its strong light absorption, largely influenced by its “mixing state”. However, due to computational constraints, the mixing state is challenging to accurately represent in large-scale models. In this study, we employ a particle-resolved model to simulate the evolution of BC mixing state based on field observation. Our result shows that aerosol compositions, coating thickness (CT) distribution, and optical properties of BC aerosols all exhibit a tendency toward a steady state with a characteristic timescale of less than 1 d, considerably shorter than the BC atmospheric lifetime. The rapid attainment of a steady state suggests that it is reasonable to disregard this pre-steady-state period and instead concentrate on the average properties of BC across extensive spatial and temporal scales. The distribution of CT follows an exponential linear distribution and can be characterized by a single slope parameter k. This distribution is independent of the BC core's distribution. In the model simulation, the mean CT, equivalent to the 1/k, is 62 nm, which is consistent with the statistical results indicating a mean CT of 63 nm. Utilizing the slope parameter k, which effectively characterizes the CT distribution under the steady-state simplifying assumption, the BC absorption enhancement closely corresponds to the results obtained via the particle-resolved method. This study simplifies the BC mixing state description and yields a precise evaluation of the BC optical properties, which has the potential utility for modeling efforts in the refinement of the assessment of BC's radiative effects.

Zhang, Zhouyang↗
NASA NTRS2008

Effective Interpolation of Incomplete Satellite-Derived Leaf-Area Index Time Series for the Continental United States

Many earth science modeling applications employ continuous input data fields derived from satellite data. Environmental factors, sensor limitations and algorithmic constraints lead to data products of inherently variable quality. This necessitates interpolation of one form or another in order to produce high quality input fields free of missing data. The present research tests several interpolation techniques as applied to satellite-derived leaf area index, an important quantity in many global climate and ecological models. The study evaluates and applies a variety of interpolation techniques for the Moderate Resolution Imaging Spectroradiometer (MODIS) Leaf-Area Index Product over the time period 2001-2006 for a region containing the conterminous United States. Results indicate that the accuracy of an individual interpolation technique depends upon the underlying land cover. Spatial interpolation provides better results in forested areas, while temporal interpolation performs more effectively over non-forest cover types. Combination of spatial and temporal approaches offers superior interpolative capabilities to any single method, and in fact, generation of continuous data fields requires a hybrid approach such as this.

Jasinski, Michael F.↗
DOE OSTI2020

NSFnets (Navier-Stokes flow nets): Physics-informed neural networks for the incompressible Navier-Stokes equations

In the last 50 years there has been a tremendous progress in solving numerically the Navier-Stokes equations using finite differences, finite elements, spectral, and even meshless methods. Yet, in many real cases, we still cannot incorporate seamlessly (multi-fidelity) data into existing algorithms, and for industrial-complexity applications the mesh generation is time consuming and still an art. Moreover, solving ill-posed problems (e.g., lacking boundary conditions) or inverse problems is often prohibitively expensive and requires different formulations and new computer codes. Here, we employ physics-informed neural networks (PINNs), encoding the governing equations directly into the deep neural network via automatic differentiation, to overcome some of the aforementioned limitations for simulating incompressible laminar and turbulent flows. We develop the Navier-Stokes flow nets (NSFnets) by considering two different mathematical formulations of the Navier-Stokes equations: the velocity-pressure (VP) formulation and the vorticity-velocity (VV) formulation. Since this is a new approach, we first select some standard benchmark problems to assess the accuracy, convergence rate, computational cost and flexibility of NSFnets; analytical solutions and direct numerical simulation (DNS) databases provide proper initial and boundary conditions for the NSFnet simulations. The spatial and temporal coordinates are the inputs of the NSFnets, while the instantaneous velocity and pressure fields are the outputs for the VP-NSFnet, and the instantaneous velocity and vorticity fields are the outputs for the VV-NSFnet. This is unsupervised learning and, hence, no labeled data are required beyond boundary and initial conditions and the fluid properties. The residuals of the VP or VV governing equations, together with the initial and boundary conditions, are embedded into the loss function of the NSFnets. No data is provided for the pressure to the VP-NSFnet, which is a hidden state and is obtained via the incompressibility constraint without extra computational cost. Unlike the traditional numerical methods, NSFnets inherit the properties of neural networks (NNs), hence the total error is composed of the approximation, the optimization, and the generalization errors. Here, we empirically attempt to quantify these errors by varying the sampling (“residual”) points, the iterative solvers, and the size of the NN architecture. For the laminar flow solutions, we show that both the VP and the VV formulations are comparable in accuracy but their best performance corresponds to different NN architectures. The initial convergence rate is fast but the error eventually saturates to a plateau due to the dominance of the optimization error. For the turbulent channel flow, we show that NSFnets can sustain turbulence at , but due to expensive training we only consider part of the channel domain and enforce velocity boundary conditions on the subdomain boundaries provided by the DNS data base. We also perform a systematic study on the weights used in the loss function for balancing the data and physics components, and investigate a new way of computing the weights dynamically to accelerate training and enhance accuracy. In the last part, we demonstrate how NSFnets should be used in practice, namely for ill-posed problems with incomplete or noisy boundary conditions as well as for inverse problems. We obtain reasonably accurate solutions for such cases as well without the need to change the NSFnets and at the same computational cost as in the forward well-posed problems. As a result, we also present a simple example of transfer learning that will aid in accelerating the training of NSFnets for different parameter settings.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Long-baseline quantum sensor network as dark matter haloscope

Ultralight dark photons constitute a well-motivated candidate for dark matter. A coherent electromagnetic wave is expected to be induced by dark photons when coupled with Standard-Model photons through kinetic mixing mechanism, and should be spatially correlated within the de Broglie wavelength of dark photons. Here we report the first search for correlated dark-photon signals using a long-baseline network of 15 atomic magnetometers, which are situated in two separated meter-scale shield rooms with a distance of about 1700 km. Both the network’s multiple sensors and the shields large size significantly enhance the expected dark-photon electromagnetic signals, and long-baseline measurements confidently reduce many local noise sources. Using this network, we constrain the kinetic mixing coefficient of dark photon dark matter over the mass range 4.1 feV-2.1 peV, which represents the most stringent constraints derived from any terrestrial experiments operating over the aforementioned mass range. Our prospect indicates that future data releases may go beyond the astrophysical constraints from the cosmic microwave background and the plasma heating.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Prompt Detection of Fast Optical Bursts with the Vera C. Rubin Observatory

Abstract The transient optical sky has remained largely unexplored on very short timescales. While there have been some experiments searching for optical transients from minutes to years, none have had the capability to distinguish millisecond fast optical bursts (FOBs). Such very fast transients could be the optical counterparts of fast radio bursts, the prompt emission from γ -ray bursts, or other previously unknown phenomena. Here, we investigate a novel approach to the serendipitous detection of FOBs, which relies on searching for anomalous spatial images. In particular, due to their short duration, the seeing-distorted images of FOBs should look characteristically different than those of steady sources in a standard optical exposure of finite duration. We apply this idea to simulated observations with the Vera C. Rubin Observatory, produced by tracing individual photons through a turbulent atmosphere, and down through the optics and camera of the Rubin telescope. We compare these simulated images to steady-source star simulations in 15 s integrations, the nominal Rubin exposure time. We report the classification accuracy results of a neural network classifier for distinguishing FOBs from steady sources. From this classifier, we derive constraints in duration–intensity parameter space for unambiguously identifying FOBs in Rubin observations. We conclude with estimates of the total number of detections of FOB counterparts to FRBs expected during the 10 yr Rubin Legacy Survey of Space and Time.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Mesoporous Thin Film Architectures: Addressing Material Demands through Molecular Self-Assembly

Mesoporous thin films spark interest across a wide range of disciplines due to their tunable nanostructures, large internal surface areas, and strong compatibility with planar optical, electronic, and microfluidic devices. While attention in the porous materials community has shifted toward macroporous or disordered nanoporous systems, a resurgence in mesoporous thin film research is underway, driven by new molecular self-assembly methods, advanced materials chemistry, and improved characterization techniques. The integration of high-χN block copolymer design, kinetically persistent micelle templating, and postdeposition processing protocols now allows control over structural parameters such as pore size, wall thickness, porosity, and connectivity. These advances have overcome many of the thermodynamic and processing constraints that previously limited widespread adoption. Rather than serving only as high-surface-area supports, mesoporous thin films are engineered as active interfaces where responsive chemistries and nanoscale confinement act in tandem. Embedding switchable ligands, thermoresponsive polymers, redox mediators, or ion-selective groups directly within the pore walls enables real-time control over transport, optical, and electrochemical properties. These capabilities open up new directions in adaptive coatings, gated membranes, and fast-response biosensors. To further expand their functional scope, mesoporous films are integrated into hierarchical and multicomponent architectures. Techniques such as triblock terpolymer templating, crack-directed assembly, and nanoimprint lithography allow for control over spatial organization on the micron and submicron scale and pore system orientation. This enables programmable anisotropy, enhanced molecular diffusion, and wavelength-selective photonic behavior, essential for next-generation sensing, catalysis, and energy applications. Such structural and functional complexity requires equally sophisticated characterization. Multimodal and in situ techniques can track material dynamics under operational conditions. Recent progress includes extended-range ellipsometric porosimetry (EP) for hierarchical architectures, vacuum EP for interface energetics, time-resolved EP for diffusion kinetics, and correlative AFM-SAXS mapping. The introduction of advanced neutron-based spectroscopies, particularly quasielastic neutron scattering (QENS), promises to provide real-time access to ion transport dynamics and segmental motion under nanoscale confinement, offering a path toward deeper mechanistic understanding of structure-performance correlations in mesoporous systems. This Account reflects the technical advances made and the interdisciplinary collaborations that have shaped our collective vision. The particular dimensions of mesopores enable us to subtly tune interactions at the molecular, interfacial, and mesoscopic levels that permit us to harness nanoconfinement. What emerges is a versatile, modular platform capable of chemical gating, energy transduction, and sensing with a level of tunability unmatched by other porous materials. We highlight critical challenges including the need for more robust large-area processing, a deeper understanding of dynamic behavior under cycling, and better integration with device-level architectures. Our strategies support the transition of mesoporous thin films into active high-performance components in next-generation energy, environmental, and biomedical systems.

oxides↗
DOE OSTI2020

General approaches for shear-correcting coordinate transformations in Bragg coherent diffraction imaging. Part I

This two-part article series provides a generalized description of the scattering geometry of Bragg coherent diffraction imaging (BCDI) experiments, the shear distortion effects inherent in the 3D image obtained from presently used methods and strategies to mitigate this distortion. Part I starts from fundamental considerations to present the general real-space coordinate transformation required to correct this shear, in a compact operator formulation that easily lends itself to implementation with available software packages. Such a transformation, applied as a final post-processing step following phase retrieval, is crucial for arriving at an undistorted, correctly oriented and physically meaningful image of the 3D crystalline scatterer. As the relevance of BCDI grows in the field of materials characterization, the available sparse literature that addresses the geometric theory of BCDI and the subsequent analysis methods are generalized here. This geometrical aspect, specific to coherent Bragg diffraction and absent in 2D transmission CDI experiments, gains particular importance when it comes to spatially resolved characterization of 3D crystalline materials in a reliable nondestructive manner. This series of articles describes this theory, from the diffraction in Bragg geometry to the corrections needed to obtain a properly rendered digital image of the 3D scatterer. Part I of this series provides the experimental BCDI community with the general form of the 3D real-space distortions in the phase-retrieved object, along with the necessary post-retrieval correction method. Part II builds upon the geometric theory developed in Part I with the formalism to correct the shear distortions directly on an orthogonal grid within the phase-retrieval algorithm itself, allowing more physically realistic constraints to be applied. Taken together, Parts I and II provide the X-ray science community with a set of generalized BCDI shear-correction techniques crucial to the final rendering of a 3D crystalline scatterer and for the development of new BCDI methods and experiments.

36 MATERIALS SCIENCE↗
NASA NTRS1993

Galaxies with spectral energy distributions peaking near 60 micrometers. 1: Optical spectroscopy, infrared photometry, and radio continuum data

We present and discuss broadband infrared photometry in the 1-100 micrometer wavelength range, optical spectroscopy, and radio continuum observations of a sample of IRAS galaxies with unusual spectral energy distributions that peak near 60 micrometers. For inclusion in this sample of '60 PKs', the galaxies must have satisfied the following criteria: IRAS flux ratios f(sub 60)/f(sub 100) greater than 1 and 1 less than f(sub 60)/f(sub 25) less than 4, and galactic latitude absolute value of b greater than 10 deg. In this paper, which is the first of a series, we show that the 60PKs are relatively scarce objects that represent about 2% of the spatial density of 60 micrometer-selected galaxies in the range L(60 micrometers) = 10(exp 9) to 10(exp 12) solar luminosity, but have a far-infrared luminosity function of nearly identical shape. They are detected up to redshifts of 0.2. Besides having the usual high percentage of active galaxies (approximately 30% H II-region like, 50% Seyfert 2, 10% Seyfert 1, 10% unknown) associated with a flat 25 to 60 micrometer spectral index, the sample also includes most of the galaxies that have been found to have a dust-obscured broadline region. We show that the additional f(sub 60)/f(sub 100) greater than 1 constraint selects galaxies with dust that is more centrally concentrated and exposed to a more intense radiation field than in most other IRAS galaxies. In particular, the cirrus component of the far-infrared radiation, which is typically the dominant contributor to the 100 micrometer emission from spiral galaxies, is negligible or missing. This is consistent with the fact that the most distinctive optical signature of 60PKs is the absence of spiral structure: they tend to be peculiar and/or amorphous objects.

Vader, J Patricia↗
NASA NTRS1995

A diagnostic analysis of the VVP single-doppler retrieval technique

A diagnostic analysis of the VVP (volume velocity processing) retrieval method is presented, with emphasis on understanding the technique as a linear, multivariate regression. Similarities and differences to the velocity-azimuth display and extended velocity-azimuth display retrieval techniques are discussed, using this framework. Conventional regression diagnostics are then employed to quantitatively determine situations in which the VVP technique is likely to fail. An algorithm for preparation and analysis of a robust VVP retrieval is developed and applied to synthetic and actual datasets with high temporal and spatial resolution. A fundamental (but quantifiable) limitation to some forms of VVP analysis is inadequate sampling dispersion in the n space of the multivariate regression, manifest as a collinearity between the basis functions of some fitted parameters. Such collinearity may be present either in the definition of these basis functions or in their realization in a given sampling configuration. This nonorthogonality may cause numerical instability, variance inflation (decrease in robustness), and increased sensitivity to bias from neglected wind components. It is shown that these effects prevent the application of VVP to small azimuthal sectors of data. The behavior of the VVP regression is further diagnosed over a wide range of sampling constraints, and reasonable sector limits are established.

Boccippio, Dennis J.↗
NASA NTRS1995

The stellar content of 30 doradus derived from spatially integrated ultraviolet spectra: A test of spectral synthesis models

Using the IUE satellite, we have obtained spatially integrated ultraviolet spectra of three areas within the giant H II region 30 Dor in the Large Magellanic Cloud. The spectra correspond to spatial reginswith sizes of 20 sec x 20 sec, 1 min x 1 min, and 3 min x 3 min, all of which are approximately centered on R136. We have performed a spectral synthesis analysis of the spectra of the two larger regions and compared the results with the known stellar content in these regions. The spectral synthesis models are sensitive to the ultraviolet continuum level, the P Cygni profile of the C Iv wavelength 1550 line, the absorption strength of the Si IV wavelength 1400 line, and the emission strength of the He II wavelength 1640 line. The intrinsic continuum levels and the profiles of these stellar wind lines provide constraints on the age and duration of the starburst episode within a region, as well as on the upper curoff mass of the initial mass function. From our analysis we find that the present-day value of the upper cutoff mass in the 1 min x 1 min and 3 min x 3 min regions has a lower limit of approximately 50 solar mass, a result which is in good agreement with several other recent determinations. The age of the starburst episode must be less than approximately 3 Myr, also in agreement with other estimates. Comparison of the observed total numbers of O and W-R stars with those predicted from the various models favors an instantaneous burst of star formation in the regions. However, the differences between the two burst scenarios we investigated (instantaneous and continuous) are small at such a young age, and distinguishing between the two is difficult. We are now confident that these spectral synthesis models can be used to determine the stellar content of more distant star-forming regions.

Vacca, William D.↗
DOE OSTI2020

Three-dimensional traction microscopy accounting for cell-induced matrix degradation

Tractions exerted by cells on the extracellular matrix (ECM) are critical in many important physiological and pathological processes such as embryonic morphogenesis, wound healing, and cancer metastasis. Three-dimensional Traction Microscopy (3DTM) is a tool to quantify cellular tractions by first measuring the displacement field in the ECM in response to these tractions, and then using this measurement to infer tractions. Most applications of 3DTM have assumed that the ECM has spatially-uniform mechanical properties, but cells secrete enzymes that can locally degrade the ECM. In this work, a novel computational method is developed to quantify both cellular tractions and ECM degradation. In particular, the ECM is modeled as a hyperelastic, Neo-Hookean solid, whose material parameters are corrupted by a single degradation parameter. The feasibility of determining both the traction and the degradation parameter is first demonstrated by showing the existence and uniqueness of the solution. An inverse problem is then formulated to determine the nodal values of the traction vector and the degradation parameter, with the objective of minimizing the difference between a predicted and measured displacement field, under the constraint that the predicted displacement field satisfies the equation of equilibrium. The inverse problem is solved by means of a gradient-based optimization approach, and the gradient is computed efficiently using appropriately derived adjoint fields. The computational method is validated in-silico using a geometrically realistic neuronal cell model and synthetic traction and degradation fields. It is found that the method accurately recovers both the traction and degradation fields. Moreover, it is found that neglecting ECM degradation can yield significant errors in traction measurements. Finally, our method can extend the range of context where tractions can be appropriately measured.

42 ENGINEERING↗
NASA NTRS2021

Simulation of Radon-222 with the GEOS-Chem Global Model: Emissions, Seasonality, and Convective Transport

Radon-222 (Rn-222) is a short-lived radioactive gas naturally emitted from land surfaces and has long been used to assess convective transport in atmospheric models. In this study, we simulate Rn-222 using the GEOS-Chem chemical transport model to improve our understanding of Rn-222 emissions and surface concentration seasonality and characterize convective transport associated with two Goddard Earth Observing System (GEOS) meteorological products, the Modern-Era Retrospective analysis for Research and Applications (MERRA) and GEOS Forward Processing (GEOS-FP). We evaluate four global Rn-222 emission scenarios by comparing model results with observations at 51 surface sites. The default emission scenario in GEOS-Chem yields a moderate agreement with surface observations globally (68.9 % of data within a factor of 2) and a large underestimate of winter surface Rn-222 concentrations at Northern Hemisphere midlatitudes and high latitudes due to an oversimplified formulation of Rn-222 emission fluxes (1 atom cm−2 s−1 over land with a reduction by a factor of 3 under freezing conditions). We compose a new global Rn-222 emission scenario based on Zhang et al. (2011) and demonstrate its potential to improve simulated surface Rn-222 concentrations and seasonality. The regional components of this scenario include spatially and temporally varying emission fluxes derived from previous measurements of soil radium content and soil exhalation models, which are key factors in determining Rn-222 emission flux rates. However, large model underestimates of surface Rn-222 concentrations still exist in Asia, suggesting unusually high regional Rn-222 emissions. We therefore propose a conservative upscaling factor of 1.2 for Rn-222 emission fluxes in China, which was also constrained by observed deposition fluxes of 210Pb (a progeny of Rn-222). With this modification, the model shows better agreement with observations in Europe and North America (> 80 % of data within a factor of 2) and reasonable agreement in Asia (close to 70 %). Further constraints on Rn-222 emissions would require additional concentration and emission flux observations in the central United States, Canada, Africa, and Asia. We also compare and assess convective transport in model simulations driven by MERRA and GEOS-FP using observed Rn-222 vertical profiles in northern midlatitude summer and from three short-term airborne campaigns. While simulations with both GEOS products are able to capture the observed vertical gradient of Rn-222 concentrations in the lower troposphere (0–4 km), neither correctly represents the level of convective detrainment, resulting in biases in the middle and upper troposphere. Compared with GEOS-FP, MERRA leads to stronger convective transport of Rn-222, which is partially compensated for by its weaker large-scale vertical advection, resulting in similar global vertical distributions of Rn-222 concentrations between the two simulations. This has important implications for using chemical transport models to interpret the transport of other trace species when these GEOS products are used as driving meteorology.

Bo Zhang↗
NASA NTRS2016

On the Distance to the North Polar Spur and the Local CO-H2 Factor

Aims. Most models identify the X-ray bright North Polar Spur (NPS) with a hot interstellar (IS) bubble in the Sco-Cen star-forming region at approx. = to 130 pc. An opposite view considers the NPS as a distant structure associated with Galactic nuclear outflows. Constraints on the NPS distance can be obtained by comparing the foreground IS gas column inferred from X-ray absorption to the distribution of gas and dust along the line of sight. Absorbing columns towards shadowing molecular clouds simultaneously constrain the CO-H2 conversion factor. Methods. We derived the columns of X-ray absorbing matter N(sub Habs) from spectral fitting of dedicated XMM-Newton observations towards the NPS southern terminus (l(sup II) approx. = to 29 degrees, b(sup II) approx. = to +5 to +11 degrees). The distribution of the IS matter was obtained from absorption lines in stellar spectra, 3D dust maps and emission data, including high spatial resolution CO measurements recorded for this purpose. Results. N(sub Habs) varies from approx. = to 4.3 to approx. = to 1.3 x 10(exp 21) per sq cm along the 19 fields. Relationships between X-ray brightness, absorbing column and hardness ratio demonstrate a brightness decrease with latitude governed by increasing absorption. The comparison with absorption data, local and large-scale dust maps rules out a NPS near side closer than 300 pc. The correlation between N(sub Habs) and the reddening increases with the sightline length from 300 pc to 4 kpc and is the tightest with Planck tau(sub 353GHz) -based reddening, suggesting a much larger distance. N(H)/E(B-V)(sub tao) approx. = to 4.1 x 10(exp 21) per sq cm mag, close to Fermi-Planck determinations. N(sub Habs) absolute values are compatible with HI-CO clouds at -5 less than or equal to V(sub LSR) less than or equal to +25 to +45 km/s and a NPS potentially far beyond the Local Arm. A shadow cast by a b=+9 degree molecular cloud constrains X(sub CO) in that direction to less than or equal to 1.0 x 10(exp 20) sec/((sq cm)K km). The average X(sub CO) over the fields is less than or equal to 0.75 x 10(exp 20) sec/((sq cm)K km).

X-rays: ISM / radio lines: ISM / local insterstell↗
DOE OSTI2023

Temperature Uncertainty Modeling with Proxy Structural Data as Geostatistical Constraints for Well Siting: An Example Applied to Granite Springs Valley, NV, USA

Utilizing existing temperature and structural information around Granite Springs Valley, Nevada, we build 3D stochastic temperature models with the aim of evaluating the 3D uncertainty of temperature and choosing between candidate exploration well locations . The data used to support the modeling are measured temperatures and structural proxies from 3D geologic modeling, the latter considered "secondary" data. Two stochastic geostatistical techniques are explored for incorporating the structural proxies: cosimulation and local varying mean. With both the cosimulation and local varying mean methods, many equally likely temperature models (i.e., realizations) are produced, from which temperature probability profiles are calculated at candidate well locations. To aid in choosing between the candidate locations, two quantities summarize the temperature probabilities: Vprior and entropy. Vprior quantifies the likelihood for economic temperatures at each candidate location, whereas entropy identifies where new information has the most potential to reduce uncertainty. In general, the cosimulation realizations have smoother spatial structure, and extrapolate high temperatures at candidate locations that are located along the direction of the longest spatial correlation, which are down dip from existing temperature logs. The smooth realizations result in tight temperature probability profiles that are easier to interpret, but they have unrealistic temperature reversals in some locations because the cosimulation technique does not enforce a conductive geothermal gradient as a baseline (i.e., linearly increasing temperature with depth). The local varying mean results produce realizations with more realistic geothermal gradients, with temperatures increasing downward since a depth-temperature relationship is included. However, because they have much noisier spatial nature compared to cosimulation, it is harder to interpret the temperature probability profiles. The different local varying mean results allow the geologist to determine which proxy (e.g., dilation versus distance to fault termination) should be used given the specific geothermal system. In general, Vprior from local varying mean results identify locations that are close to high values for the structural proxies: areas with highe r probabilities for higher temperatures. The entropy results identify where uncertainty is greatest and therefore new drilling information could be most useful. Though these techniques provide useful information, even when applied to areas of sparse data, our comp arison of these two techniques demonstrates the need for new geothermal geostatistics techniques that combine the advantages of these two methods and that are tailored to the spatial uncertainty issues inherent in geothermal exploration.

3D temperature modeling↗
DOE OSTI2022

Rising CO 2 and warming reduce global canopy demand for nitrogen

Nitrogen (N) limitation has been considered as a constraint on terrestrial carbon uptake in response to rising CO 2 and climate change. By extension, it has been suggested that declining carboxylation capacity (V cmax ) and leaf N content in enhanced-CO 2 experiments and satellite records signify increasing N limitation of primary production. We predicted Vcmax using the coordination hypothesis and estimated changes in leaf-level photosynthetic N for 1982–2016 assuming proportionality with leaf-level V cmax at 25°C. The whole-canopy photosynthetic N was derived using satellite-based leaf area index (LAI) data and an empirical extinction coefficient for V cmax , and converted to annual N demand using estimated leaf turnover times. The predicted spatial pattern of V cmax shares key features with an independent reconstruction from remotely sensed leaf chlorophyll content. Predicted leaf photosynthetic N declined by 0.27% yr -1 , while observed leaf (total) N declined by 0.2–0.25% yr -1 . Predicted global canopy N (and N demand) declined from 1996 onwards, despite increasing LAI. Leaf-level responses to rising CO 2 , and to a lesser extent temperature, may have reduced the canopy requirement for N by more than rising LAI has increased it. This finding provides an alternative explanation for declining leaf N that does not depend on increasing N limitation.

54 ENVIRONMENTAL SCIENCES↗