Search NASA⌕ Search

SEARCH · Search NASA

Results for “spatial variability”

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

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

At least 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

An Accurate SUPG-stabilized Continuous Galerkin Discretization for Anisotropic Heat Flux in Magnetic Confinement Fusion

We present a novel spatial discretization for the anisotropic heat conduction equation, aimed at improved accuracy at the high levels of anisotropy seen in a magnetized plasma, for example, for magnetic confinement fusion. The new discretization is based on a mixed formulation, introducing a form of the directional derivative along the magnetic field as an auxiliary variable and discretizing both the temperature and auxiliary fields in a continuous Galerkin (CG) space. Both the temperature and auxiliary variable equations are stabilized using the streamline upwind Petrov–Galerkin (SUPG) method, ensuring a better representation of the directional derivatives and therefore an overall more accurate solution. This approach can be seen as the CG-based version of our previous work (Wimmer, Southworth, Gregory, Tang, 2024), where we considered a mixed discontinuous Galerkin (DG) spatial discretization including DG-upwind stabilization. We prove consistency of the novel discretization, and demonstrate its improved accuracy over existing CG-based methods in test cases relevant to magnetic confinement fusion. This includes a long-run tokamak equilibrium sustainment scenario, demonstrating a 35% and 32% spurious heat loss for existing primal and mixed CG-based formulations versus 4% for our novel SUPG-stabilized discretization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Multi-scale integration of satellite remote sensing improves characterization of dry-season green-up in an Amazon tropical evergreen forest

In tropical forests, leaf phenology-particularly the pronounced dry-season green-up-strongly regulates biogeochemical cycles of carbon and water fluxes. However, uncertainties remain in the understanding of tropical forest leaf phenology at different spatial scales. Phenocams accurately characterize leaf phenology at the crown and ecosystem scales but are limited to a few sites and time spans of a few years. Time-series satellite observations might fill this gap, but the commonly used satellites (e.g. MODIS, Landsat and Sentinel-2) have resolutions too coarse to characterize single crowns. To resolve this observational challenge, we used the PlanetScope constellation with a 3m resolution and near daily nadir-view coverage. We first developed a rigorous method to cross-calibrate PlanetScope surface reflectance using daily BRDF-adjusted MODIS as the reference. We then used linear spectral unmixing of calibrated PlanetScope to obtain dry-season change in the fractional cover of green vegetation (GV) and non-photosynthetic vegetation (NPV) at the PlanetScope pixel level. We used the Central Amazon Tapajos National Forest k67 site, as all necessary data (from field to phenocam and satellite observations) was available. For this proof of concept, we chose a set of 22 dates of PlanetScope measurements in 2018 and 16 in 2019, all from the six drier months of the year to provide the highest possible cloud-free temporal resolution. Our results show that MODIS-calibrated dry-season PlanetScope data (1) accurately assessed seasonal changes in ecosystem-scale and crown-scale spectral reflectance; (2) detected an increase in ecosystem-scale GV fraction (and a decrease in NPV fraction) from June to November of both years, consistent with local phenocam observations with R 2 around 0.8; and (3) monitored large seasonal trend variability in crown-scale NPV fraction. Finally, our results highlight the potential of integrating multi-scale satellite observations to extend fine-scale leaf phenology monitoring beyond the spatial limits of phenocams.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

3D optical diagnostics for explosively driven deformation and fragmentation

High-speed, optical imaging diagnostics are presented for three-dimensional (3D) quantification of explosively driven metal fragmentation. At early times after detonation, Digital Image Correlation (DIC) provides non-contact measures of 3D case velocities, strains, and strain rates, while a proposed stereo imaging configuration quantifies in-flight fragment masses and velocities at later times. Experiments are performed using commercially obtained RP-80 detonators from Teledyne RISI, which are shown to create a reproducible fragment field at the benchtop scale. DIC measurements are compared with 3D simulations, which have been ‘leveled’ to match the spatial resolution of DIC. Overall, results demonstrate improved ability to identify predicted quantities-of-interest that fall outside of measurement uncertainty and shot-to-shot variability. Similarly, video measures of fragment trajectories and masses allow rapid experimental repetition and provide correlated fragment size-velocity measurements. Measured and simulated fragment mass distributions are shown to agree within confidence bounds, while some statistically meaningful differences are observed between the measured and predicted conditionally averaged fragment velocities. Together these techniques demonstrate new opportunities to improve future model validation.

42 ENGINEERING↗
DOE OSTI2023

Understanding the hydrological response of a headwater-dominated catchment by analysis of distributed surface–subsurface interactions

Abstract We computationally explore the relationship between surface–subsurface exchange and hydrological response in a headwater-dominated high elevation, mountainous catchment in East River Watershed, Colorado, USA. In order to isolate the effect of surface–subsurface exchange on the hydrological response, we compare three model variations that differ only in soil permeability. Traditional methods of hydrograph analysis that have been developed for headwater catchments may fail to properly characterize catchments, where catchment response is tightly coupled to headwater inflow. Analyzing the spatially distributed hydrological response of such catchments gives additional information on the catchment functioning. Thus, we compute hydrographs, hydrological indices, and spatio-temporal distributions of hydrological variables. The indices and distributions are then linked to the hydrograph at the outlet of the catchment. Our results show that changes in the surface–subsurface exchange fluxes trigger different flow regimes, connectivity dynamics, and runoff generation mechanisms inside the catchment, and hence, affect the distributed hydrological response. Further, changes in surface–subsurface exchange rates lead to a nonlinear change in the degree of connectivity—quantified through the number of disconnected clusters of ponding water—in the catchment. Although the runoff formation in the catchment changes significantly, these changes do not significantly alter the aggregated streamflow hydrograph. This hints at a crucial gap in our ability to infer catchment function from aggregated signatures. We show that while these changes in distributed hydrological response may not always be observable through aggregated hydrological signatures, they can be quantified through the use of indices of connectivity.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

TDCOSMO - XVI. Measurement of the Hubble constant from the lensed quasar WGD 2038–4008

Time-delay cosmography is a powerful technique to constrain cosmological parameters, particularly the Hubble constant (H0). The TDCOSMO Collaboration is performing an ongoing analysis of lensed quasars to constrain cosmology using this method. In this work, we obtain constraints from the lensed quasar WGD 2038−4008 using new time-delay measurements and previous mass models by TDCOSMO. This is the first TDCOSMO lens to incorporate multiple lens modeling codes and the full time-delay covariance matrix into the cosmological inference. The models are fixed before the time delay is measured, and the analysis is performed blinded with respect to the cosmological parameters to prevent unconscious experimenter bias. We obtain DΔ t = 1.68−0.38+0.40 Gpc using two families of mass models, a power-law describing the total mass distribution, and a composite model of baryons and dark matter, although the composite model is disfavored due to kinematics constraints. In a flat ΛCDM cosmology, we constrain the Hubble constant to be H0 = 65−14+23 km s−1 Mpc−1. The dominant source of uncertainty comes from the time delays, due to the low variability of the quasar. Future long-term monitoring, especially in the era of the Vera C. Rubin Observatory’s Legacy Survey of Space and Time, could catch stronger quasar variability and further reduce the uncertainties. This system will be incorporated into an upcoming hierarchical analysis of the entire TDCOSMO sample, and improved time delays and spatially-resolved stellar kinematics could strengthen the constraints from this system in the future.Key words: gravitational lensing: strong / cosmological parameters / distance scale⋆ Corresponding author; kcwong19@gmail.com.⋆⋆ NHFP Einstein fellow.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Amplified bottom water acidification rates on the Bering Sea shelf from 1970–2022

The Bering Sea shelf supports a highly productive marine ecosystem that is vulnerable to ocean acidification (OA) due to the cold, carbon-rich waters. Previous observational evidence suggests that bottom waters on the shelf are already seasonally undersaturated with respect to aragonite (i.e. Ω arag <1) and that OA will continue to increase the spatial extent, duration, and intensity of these conditions. Here, we use a regional ocean biogeochemical model to simulate changes in ocean carbon chemistry for the Bering Sea shelf from 1970–2022. Over this timeframe, model results suggest that surface Ω arag decreases by −0.043 per decade and surface pH by −0.014 per decade, comparable to observed global rates of OA. However, bottom water pH decreases at twice the rate of surface pH, while bottom [H + ] decreases at nearly 3 times the rate of surface [H + ]. This amplified bottom water acidification has emerged over the past 25 years and is likely driven by a combination of anthropogenic carbon accumulation and increasing primary productivity and subsurface respiration and remineralization. Due to this enhanced bottom water acidification, the spatial extent of bottom waters with Ω arag <1 has greatly expanded over the past 2 decades, along with pH conditions harmful to red king crab. Interannual variability in surface and bottom Ω arag , pH, and [H + ] has also increased over the past 2 decades, resulting in part from the increased physical climate variability. We also find that the Bering Sea shelf is a net annual carbon sink of 1.1–7.9 Tg C yr −1 , with the range resulting from the difference in the two different atmospheric forcing reanalysis products used. Seasonally, the shelf is a significant carbon sink from April–October but a somewhat weaker carbon source from November–March.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Design and expected performance of a variable-radii sinusoidal spiral x-ray spectrometer for the National Ignition Facility

A novel high-resolution x-ray spectrometer for point-like emission sources has been developed using a crystal shape having both a variable major and a variable minor radius of curvature. This variable-radii sinusoidal spiral spectrometer (VR-Spiral) allows three common spectrometer design goals to be achieved simultaneously: 1. reduction of aberrations and improved spectral (energy) resolution, 2. reduction of source size broadening, and 3. use of large crystals to improve total throughput. The VR-Spiral concept and its application to practical spectrometer design are described in detail. This concept is then used to design a spectrometer for an extreme extended x-ray absorption fine structure experiment at the National Ignition Facility looking at the Pb L 3 absorption edge at 13.0352 keV. The expected performance of this VR-Spiral spectrometer, both in terms of energy resolution and spatial resolution, is evaluated through the use of a newly developed raytracing tool, xicsrt. Finally, the expected performance of the VR-Spiral concept is compared to that of spectrometers based on conventional toroidal and variable-radii toroidal crystal geometries showing a greatly improved energy resolution.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Spatially explicit terrestrial carbon densities for calibrating the carbon cycle in human-Earth system Models

Soil and vegetation carbon stocks play a critical role in human-Earth system models. These stocks (denominated as densities in MgC/ha) affect variables such as land use change emissions and also influence land use change pathways under climate forcing scenarios where terrestrial carbon is assigned a carbon price. Here we present reharmonized soil and vegetation carbon densities both at the 5-arcmin resolution grid cell level and also aggregated to 235 water sheds for 4 land use types (Cropland, Grazed land, Urban land and unmanaged vegetation) and 15 unmanaged land cover types. Moreover, we use the distribution of carbon within and across pixels to define statistical "states" of carbon, once again differentiated by land type. These statistical states are used to define a range of possible carbon values that can be used for defining initial conditions of soil and vegetation carbon in human-Earth system models. We implement these data in a state-of-the-art multi sector dynamics model, namely the Global Change Analysis Model (GCAM), and show that these new data improve several land use responses, especially when terrestrial carbon is assigned a carbon price.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Mitigating drought mortality by incorporating topography into variable forest thinning strategies

Abstract Drought-induced productivity reductions and tree mortality have been increasing in recent decades in forests around the globe. Developing adaptation strategies hinges on an adequate understanding of the mechanisms governing the drought vulnerability of forest stands. Prescribed reduction in stand density has been used as a management tool to reduce water stress and wildfire risk, but the processes that modulate fine-scale variations in plant water supply and water demand are largely missing in ecosystem models. We used an ecohydrological model that couples plant hydraulics with groundwater hydrology to examine how within-stand variations in tree spatial arrangements and topography might mitigate forest vulnerability to drought at individual-tree and stand scales. Our results demonstrated thinning generally ameliorated plant hydraulic stress and improved carbon and water fluxes of the remaining trees, although the effectiveness varied by climate and topography. Variable thinning that adjusted thinning intensity based on topography-mediated water availability achieved higher stand productivity and lower mortality risk, compared to evenly-spaced thinning at comparable intensities. The results from numerical experiments provided mechanistic evidence that topography mediates the effectiveness of thinning and highlighted the need for an explicit consideration of within-stand heterogeneity in trees and abiotic environments when designing forest thinning to mitigate drought impacts.

Environmental Sciences & Ecology↗
DOE OSTI2024

Science Area 1: Standard Award: Model-Data Fusion to Examine Multiscale Dynamical Controls on Snow Cover and Critical Zone Moisture Inputs (Final Report)

In many mountain watersheds of the world, seasonal snowpacks play an important role as natural reservoirs of water. Seasonal snowpacks accumulate water during cold, wet winter months that subsequently melts. Downstream communities depend on water from melting seasonal snowpacks to support agricultural, industrial, and municipal water needs. Rapidly melting snowpacks can also present a flooding hazard, particularly if snowpacks melt at rates faster than anticipated and where adequate reservoir capacity is unavailable to buffer river flows associated with melt. The spatial and temporal dynamics of snow accumulation and melt also play an important role in supporting upland ecosystems in mountain landscapes. Snowmelt provides soil moisture that enable terrestrial ecosystem productivity and exert control on soil microorganisms that play important roles in global carbon cycles. Climate warming is gradually decreasing the amount of precipitation in mountain watersheds arriving as snow, presenting potentially profound disruptions to mountain ecosystems, as well as downstream delivery of water. The overarching goal of this project was to understand how interactions between the near-surface atmosphere and surface topography control the input, accumulation, retention, and release of water from mountain snowpacks. Over a 5-year period, this project pursued an approach combining high-resolution regional climate modeling, satellite and airborne remote sensing data, and ground-based observations to develop and analyze benchmark datasets to address overarching science questions and hypotheses. Key products include a continuous, long-term, high spatiotemporal resolution (1 km/1 hr) dataset characterizing key climate variables in the Upper Colorado River Basin. The dataset included historical estimates of precipitation, temperature, humidity, solar and longwave radiation, and wind speeds. Additionally, the project developed a 20+ year long, 30 m spatial, daily temporal multi-sensor dataset characterizing snow presence/absence in the East/Taylor River watersheds in the Central Rocky Mountains of Colorado. The project supported training of 1 postdoctoral scholar, 1 Ph.D. student, and 1 M.S. student.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Spectral Fidelity of Earth's Terrestrial and Aquatic Ecosystems

The Surface Biology and Geology (SBG) investigation will create global maps of spectral surface reflectance and emissivity at a cadence of 16 days or better, with coverage to address global questions about Earth's geology, cryosphere and ecosystems. The revolutionary potential poses a commensurate challenge: creating contiguous maps free from regional biases induced by atmosphere, observation geometry, or inversion error. This will require an accurate calibration with precise knowledge of each channel's spectral response. In this work, we quantify the impact of spectral calibration on SBG's aquatic and terrestrial ecosystem objectives. We find that contemporary algorithms for ecosystem trait retrieval demand more accurate spectral calibration than historical missions. Errors due to drift or spatial nonuniformity in the wavelength calibration that have previously been considered acceptable can cause systematic errors larger than the instrument noise and of the same order as the variability SBG aims to measure. Moreover, their impact on atmospheric correction can induce climate-dependent systematic errors that thwart comparisons between ecosystems. These results underscore the importance of spectral response accuracy in SBG mission design.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Assessment of turbulent combustion models for simulating pre-chamber ignition in a natural gas engine

Lean combustion in an internal combustion engine is a promising strategy to increase thermal efficiency by leveraging a more favorable specific heat ratio of the fresh mixture and simultaneously suppressing the heat losses to the cylinder wall. However, unstable ignition events and slow flame propagation at fuel-lean condition lead to high cycle-to-cycle variability and hence limit the high-efficiency engine operating range. Pre-chamber ignition is considered an effective concept to extend the lean operating limit, by providing spatially distributed ignition with multiple turbulent flame-jets and enabling faster combustion rate compared to the conventional spark ignition approach. From a numerical modeling perspective, to date, still the science base and available simulation tools are inadequate for understanding and predicting the combustion processes in pre-chamber ignited engines. In this paper, conceptually different RANS combustion models widely adopted in the engine modeling community were used to simulate the ignition and combustion processes in a medium-duty natural gas engine with a pre-chamber spark-ignition system. A flamelet-based turbulent combustion model, i.e., G-equation, and a multi-zone well-stirred reactor model were employed for the multi-dimensional study. Simulation results were compared with experimental data in terms of in-cylinder pressure and heat release rate. Finally, the analysis of the performance of the two models is carried out to highlight the strengths and limitations of the two formulations respectively.

combustion model evaluation↗
DOE OSTI2020

Effect of Computational Schemes on Coupled Flow and Geo-Mechanical Modeling of CO 2 Leakage through a Compromised Well

Carbon capture, utilization, and storage (CCUS) describes a set of technically viable processes to separate carbon dioxide (CO 2 ) from industrial byproduct streams and inject it into deep geologic formations for long-term storage. Legacy wells located within the spatial domain of new injection and production activities represent potential pathways for fluids (i.e., CO 2 and aqueous phase) to leak through compromised components (e.g., through fractures or micro-annulus pathways). The finite element (FE) method is a well-established numerical approach to simulate the coupling between multi-phase fluid flow and solid phase deformation interactions that occur in a compromised well system. We assumed the spatial domain consists of a three-phases system: a solid, liquid, and gas phase. For flow in the two fluids phases, we considered two sets of primary variables: the first considering capillary pressure and gas pressure (PP) scheme, and the second considering liquid pressure and gas saturation (PS) scheme. Fluid phases were coupled with the solid phase using the full coupling (i.e., monolithic coupling) and iterative coupling (i.e., sequential coupling) approaches. The challenge of achieving numerical stability in the coupled formulation in heterogeneous media was addressed using the mass lumping and the upwinding techniques. Numerical results were compared with three benchmark problems to assess the performance of coupled FE solutions: 1D Terzaghi’s consolidation, Liakopoulos experiments, and the Kueper and Frind experiments. We found good agreement between our results and the three benchmark problems. For the Kueper and Frind test, the PP scheme successfully captured the observed experimental response of the non-aqueous phase infiltration, in contrast to the PS scheme. These exercises demonstrate the importance of fluid phase primary variable selection for heterogeneous porous media. We then applied the developed model to the hypothetical case of leakage along a compromised well representing a heterogeneous media. Considering the mass lumping and the upwinding techniques, both the monotonic and the sequential coupling provided identical results, but mass lumping was needed to avoid numerical instabilities in the sequential coupling. Additionally, in the monolithic coupling, the magnitude of primary variables in the coupled solution without mass lumping and the upwinding is higher, which is essential for the risk-based analyses.

deformation flow↗
DOE OSTI2020

A new spatially distributed added value index for regional climate models: the EURO-CORDEX and the CORDEX-CORE highest resolution ensembles

The added value of using regional climate models (RCMs) to downscale data from general circulation models (GCMs) has often been questioned and researched. Although several studies have used different methods to identify (and in some cases quantify) the added value, there is still a need to find a general metric that quantifies the added value of any variable. This paper builds on past studies to propose a new metric of added value in the simulation of present-day climate which measures the difference in the probability density functions (PDFs) at each grid-cell between a model and an observation source, and then compares the results of the RCM and GCM in order to spatially compute the added value index. The same method is also adapted to quantify the climate change downscaling signal in a way that is consistent with the present-day metric. These new metrics are tested on the daily precipitation output from the EURO-CORDEX and CORDEX-CORE projection ensembles and reveal an overall positive added value of RCMs, especially at the tail-end of the distribution. Higher added value is obtained in areas of complex topography and coast-lines, as well as in tropical regions. Areas with large added value in present-day climate are consistent with areas of significant climate change downscaling signal in the RCP 8.5 far future simulations, and when the analysis is repeated at a low-resolution. The use of different resolution observations shows that the added value tends to decrease when models are compared to low-resolution observation datasets.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Gridded Sub-daily Climate Forcings for North America Based on Daymet and ERA5 (Daymet-ERA5)

To support high spatial and temporal resolution land surface modeling, this dataset provides hourly time step historic weather forcing at 1-km spatial resolution for the entire North America. The latest Daymet V4 data provides gridded historic daily weather observations at 1-km spatial resolution from 1980 to 2024. Using sub-daily temporal information from the ECMWF ERA5 reanalysis, Daymet was further temporally downscaled to hourly time steps and provided in the format required for land surface model simulations. The process of temporal downscaling preserves the relative magnitude in each hourly time step from ERA5 while maintaining the total and average values from Daymet for each day. This results in a blended 1980-2024 Daymet-ERA5 dataset. Available variables include surface air temperature, precipitation, specific humidity, shortwave and longwave radiation, wind speed, and pressure. These data can be used as a high-resolution meteorological forcing dataset to support high-resolution land surface modeling where accurate meteorological forcing datasets built from historic observations and/or reanalysis datasets are desirable.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Atmospheric teleconnection associated with the Atlantic multidecadal variability in summer: assessment of the CESM1 model

Observational analysis shows that the Atlantic multidecadal variability (AMV) is associated with climate variability in the Northern Hemisphere through a zonal atmospheric teleconnection extending from the North Atlantic Ocean and propagating eastward around the Northern Hemisphere. Here, we studied the fidelity of model simulations in reproducing the observed summer AMV and the associated impacts on the mid-latitude climate by analysing simulations using the National Centre for Atmospheric Research Community Earth System Model Version 1 (CESM1), including CESM1 North Atlantic idealized and pacemaker simulations, CESM1 large ensemble twentieth century uninitialized simulations and large ensemble initialized CESM1 decadal predictions. To further compare the fidelity of CESM1, we also analysed large ensemble simulations from three other models. Our results suggest that the uninitialized large ensemble simulations from all models can produce an AMV time evolution and its regional climate impacts similar to the observations to certain degree. By initializing the observed oceanic condition in decadal prediction simulations, the simulated AMV and its regional impacts are closer to the observed ones than those in uninitialized ensemble simulations. In addition, the pacemaker simulations that nudged the time-evolving observed North Atlantic sea surface temperature anomalies produce spatiotemporal characteristics of the AMV and AMV climate impacts closer to the observed ones than the uninitialized simulations. We conclude that although coupled models can produce AMV and its regional impacts similar to observed, proper initialization and bias correction of the sea surface temperature spatial and temporal structure can improve this capability.

Initialization↗
DOE OSTI2023

Verification and validation of the open-source plasma fluid code: Zapdos

Zapdos is an open-source finite element plasma fluid solver based on the MOOSE multiphysics framework. This paper outlines Zapdos verification, benchmarking, and validation efforts for 1D and 2D RF capacitively coupled plasma discharge models for mid-range pressures (0.1 - 1 Torr). The verification process involved using the method of manufactured solutions to assess Zapdos spatial and temporal error convergence. L2 errors ranged from 10 -2 to 10 -4 , while the convergence’s slope were in agreement with the predicted slopes for the tested variable and time integration orders. The benchmarking process involved comparisons to previously results from the validated finite element code, LSODI. These works included 1D and 2D simulations for a range of plasma parameters (densities, temperatures, voltage, etc.). For the 1D cases, Zapdos and LSODI results were in very good agreement. In the 2D cases, variable behaviors matched, with slight discrepancies in peak values. The validation process involved comparisons to experimental works including electron density measurements by microwave interferometry and metastable density measurements by planar laser-induced fluorescence imaging. Results shown reasonable agreement at higher pressure, with results starting to diverge at low pressures. Furthermore, probable causes for this diverges are the limitation of the fluid assumption for plasmas at low pressure, or the need for more robust boundary conditions. Overall, Zapdos shown reasonable results for the verification, benchmarking, and validation efforts, and Zapdos can be downloaded at https://github.com/shannon-lab/zapdos.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Integrated Modeling of U.S. Agricultural Soil Emissions of Reactive Nitrogen and Associated Impacts on Air Pollution, Health, and Climate

Agricultural soils are leading sources of reactive nitrogen (Nr) species including nitrogen oxides (NO x ), ammonia (NH 3 ), and nitrous oxide (N 2 O). Here, the propensity of NO x and NH 3 to generate ozone and fine particulate matter (PM) and associated impacts on health are highly variable, whereas the climate impacts of long-lived N 2 O are independent of emission timing and location. However, these impacts have rarely been compared on a spatially-resolved, monetized basis. In this study, we update the nitrogen scheme in an agroecosystem model to simulate the agricultural emissions of Nr from across the United States. We then apply a reduced-form air pollution health effects model to assess air quality impacts from NO x and NH 3 and a social cost of N 2 O to assess the climate impacts. Assuming an $\$$8.2 million Value of a Statistical Life (VSL) and a $\$$13,100/ton social cost of N 2 O, the air quality impacts are a factor of ~7 to 15 times as large as the climate impacts in heavily populated coastal regions, whereas the ratios are closer to 2.5 in sparsely populated regions. Our results show that air pollution, health, and climate should be considered jointly in future assessments of how farming practices affect Nr emissions.

54 ENVIRONMENTAL SCIENCES↗