Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spatial analysis”

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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Object-oriented analysis as a foundation for building climate storylines of compounding short-term drought and crop heat stress

Introduction: Crops are vulnerable to precipitation and heat extremes during late spring through summer. Methods: We analyzed for a north-central U.S. region short-term drought and agricultural heat stress during April-May-June-July. We used the 4-km Parameter Elevation Regression on Independent Slopes Model (PRISM) for observations, aggregated to a 25-km grid, and two 25-km Regional Climate Model version 4 (RegCM4) simulns used either GFDL- or MPI-GCM boundary conditions. We chose 1981-2000 as our contemporary time period, and 2041- 2060 as our scenario time period, which used the Representative Concentration Pathway 8.5 emissions scenario. We used object-oriented analysis to identify events of interest in observations and simulations by identifying objects in a space-time domain that meet specified criteria, such as exceeding a heat-stress temperature threshold. The event diagnosis allowed analysis of compound events, occurring when temperature and drought objects overlap. Results: Identified objects yielded events that can undermine agricultural productivity and which are thus relevant to decision makers, making them building blocks for possible climate storylines. The observations and simulations showed similar spatial distributions of event frequencies across the analysis region. However, the simulations attained this distribution by having fewer events that tend to cover larger areas compared to observed events, suggesting that the effective resolution of the simulations was coarser than their 25-km grids. Short-term drought frequency increased and heat-stress frequency decreased in transitioning to the scenario climate. When compounding occurred heat-stress events generally preceded the short-term drought events. The overlapping, compound events tended to be more extreme compared to non-overlapping events of either type. Discussion: The information yielded projected changes in these agriculturally motivated events. One prominent conditional behavior emerging from the work was that a heat-stress event should be a warning to watch for potential drought, as both could compound each other to more intense levels.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

High-Resolution South American Wind Resource Data Downscaled with Generative Machine Learning Conditioned on Near-Surface Observations

High-resolution historical wind data was developed for the entirety of South America using the innovative Super-Resolution for Renewable Resource Data (sup3r) machine learning framework. The publicly available Sup3rWind South America dataset represents a significant advancement in wind resource data generation, leveraging generative machine learning conditioned on near-surface observations from the Meteorological Assimilation Data Ingest System (MADIS) to efficiently and accurately downscale coarse reanalysis data from the European Centre for Medium-Range Weather Forecasts (ERA5). This approach produces fine-scale, spatially and temporally coherent wind and meteorological fields hundreds of times more computationally efficient than traditional numerical weather modeling methods, enabling access to high-fidelity wind information across both continental and offshore regions. Sup3rWind South America builds on the earlier Sup3rWind Ukraine dataset through improvements in model architecture and outputs conditioned on near-surface observation inputs. As with the Ukraine data release, this dataset includes wind speed, wind direction, temperature, relative humidity, and pressure at a horizontal resolution of ~2 km, representing a 15x spatial enhancement relative to the 31 km ERA5 grid. Wind speed and direction are provided at 5-minute resolution, a 12x temporal refinement compared to the hourly ERA5 data, while temperature, relative humidity, and pressure remain at hourly resolution. The data covers all years from 2005 to 2024. Before downscaling, ERA5 inputs were bias-corrected using long-term monthly means and a limited number of quality-controlled observations to align large-scale statistics with regional conditions. The resulting dataset is the first publicly available high-resolution timeseries wind record that provides full spatial coverage of South America. Model validation demonstrates strong agreement with observations across several statistical metrics, consistent with other state-of-the-art high-resolution wind resource datasets. The potential applications of Sup3rWind South America span renewable energy resource assessment, energy system modeling, and grid resilience analysis. The 20-year record and high spatial and temporal resolution support accurate estimation of long-term energy yield and the economic feasibility of potential wind development sites. Continuous coverage across both continental and offshore regions enables comprehensive site prospecting within exclusive economic zones. The 2 km, 5-minute resolution data provide the spatial and temporal variability required for power system simulation, operational planning, and regional risk assessments.

17 WIND ENERGY↗
NASA NTRS1983

Comparison of phase and amplitude structures of Pc 5 pulsations in the morning and afternoon sectors

North American IMS magnetometer network data are used in comparisons of the spatial phase and amplitude structures of the morning and afternoon Pc 5 pulsations, using cross spectral analysis of phase and amplitude readings between different stations. The spatial phase structures are found to be very similar in the morning and afternoon sectors for each of the three components. Amplitudes are larger in the magnetic north than in the magnetic east components for both sectors. Latitudinal amplitude and phase difference variations between the magnetic north and east components exhibit dissimilarities in the morning and afternoon sectors. These results suggest that compressional mode magnetic field oscillations are not observable on the ground in this time sector. The pulsations in both time sectors are instead seen as caused by transverse mode magnetospheric field oscillations.

Saka, O.↗
DOE OSTI2020

Optical Characterization of the Keck Array and BICEP3 CMB Polarimeters from 2016 to 2019

The BICEP/Keck experiment (BK) is a series of small-aperture refracting telescopes observing degree-scale cosmic microwave background (CMB) polarization from the South Pole in search of a primordial B-mode signature. This B-mode signal arises from primordial gravitational waves interacting with the CMB and has amplitude parametrized by the tensor-to-scalar ratio r. Since 2016, BICEP3 and the Keck Array have been observing with 4800 total antenna-coupled transition-edge sensor detectors, with frequency bands spanning 95, 150, 220, and 270 GHz. Here we present the optical performance of these receivers from 2016 to 2019, including far-field beams measured in situ with an improved chopped thermal source and instrument spectral response measured with a field-deployable Fourier transform spectrometer. As a pair differencing experiment, an important systematic that must be controlled is the differential beam response between the co-located, orthogonally polarized detectors. We generate per-detector far-field beam maps and the corresponding differential beam mismatch that is used to estimate the temperature-to-polarization leakage in our CMB maps and to give feedback on detector and optics fabrication. The differential beam parameters presented here were estimated using improved low-level beam map analysis techniques, including efficient removal of non-Gaussian noise as well as improved spatial masking. These techniques help minimize systematic uncertainty in the beam analysis, with the goal of constraining the bias on r induced by temperature-to-polarization leakage to be subdominant to the statistical uncertainty. This is essential as we progress to higher detector counts in the next generation of CMB experiments.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2002

Plasmaspheric Density Troughs: Global IMAGE EUV Observations and Analysis via Global Core Plasma Modeling

To date, the IMAGE EUV camera has observed several plasmaspheric density trough features inside the plasmapause under a wide range of geomagnetic activity. From the perspective of EUV, a density trough feature appears as a channel of diminished pixel counts which spans a width of L-shell (DELTA L) and magnetic local time (MLT) inside the plasmapause. Plasmaspheric density troughs are found to be morphologically complex possessing considerable spatial and temporal variability. We present an analysis of the evolution of trough DELTA L and MLT extent as functions of associated D (sub ST) and K (sub p) history. Trough features range in size from 0.16 less than or equal to DELTA L less than or equal to 1.2 with azimuthal extent from 1500 less than or equal to MLT less than or equal to 1200. All cases of plasmaspheric density troughs studied to date appear to have evolved as a result of the inner edge of the afternoon/evening plasma drainage plume being wrapped around through the nightside plasmasphere. The structure of plasmaspheric density trough features is further probed by analyzing simulated EUV images produced by forward modeling artificially introduced regions of depleted density into both static and dynamic global core plasmaspheric models. Forward modeling suggests that (1) L-shell refilling of density troughs during storm recovery can be modeled as filling from the ionosphere toward the equator (i.e., bottom-up refilling), and (2) that an erosion process is operating within flux tubes beyond the outer L-shell wall of the observed density troughs.

Adrian, M. L.↗
DOE OSTI2021

Aquifer Monitoring Using Ambient Seismic Noise Recorded With Distributed Acoustic Sensing (DAS) Deployed on Dark Fiber

Abstract Groundwater is a critical resource for human activities worldwide, and a vital component of many natural ecosystems. However, the state and dynamics of water‐bearing aquifers remain uncertain, mostly due to the paucity of subsurface data at high spatial and temporal resolution. Here, we show that analysis of infrastructure‐generated ambient seismic noise acquired on distributed acoustic sensing (DAS) arrays has potential as a tool to track variations in seismic velocities (d v / v ) caused by groundwater level fluctuations. We analyze 5 months of ambient noise acquired along an unused, 23 km‐long telecommunication fiber‐optic cable in the Sacramento Valley, CA, a so‐called “dark fiber." Three array subsections, ∼6 km apart, are processed and the stretching technique is applied to retrieve daily d v / v beneath each location. Near the Sacramento river, d v / v variations in the order of 2%–3% correlate with precipitation events and fluctuations in river stage of ∼1.5 m. In contrast, regions away (2.5 km) from the river do not experience large d v / v variations. These observations reveal short‐scale spatial variability in aquifer dynamics captured by this approach. Dispersion analysis and surface wave inversion of noise gathers reveal that seismic velocity perturbations occur at depths of 10–30 m. Rock physics modeling confirms that observed d v / v are linked to pore pressure changes at these depths, caused by groundwater table fluctuations. Our results suggest that DAS combined with ambient noise interferometry provides a means of tracking aquifer dynamics at high spatial and temporal resolutions at local to regional scales, relevant for effective groundwater resource management.

58 GEOSCIENCES↗
NASA NTRS1985

Residual stresses of thin, short rectangular plates

The analysis of the residual stresses in thin, short rectangular plates is presented. The analysis is used in conjunction with a shadow moire interferometry technique by which residual stresses are obtained over a large spatial area from a strain measurement. The technique and analysis are applied to a residual stress measurement of polycrystalline silicon sheet grown by the edge-defined film growth technique.

Andonian, A. T.↗
DOE OSTI2025

Disentangling the Impacts of Microtopography and Shrub Distribution on Snow Depth in a Subarctic Watershed: Toward a Predictive Understanding of Snow Spatial Variability: Supporting Data and Code

This repository contains R code and associated datasets for reproducing the analysis described in the manuscript titled “Disentangling the Impacts of Microtopography and Shrub Distribution on Snow Depth in a Subarctic Watershed: Toward a Predictive Understanding of Snow Spatial Variability” (DOI: 10.1029/2024JG008604). The provided scripts facilitate a comprehensive analysis of snow depth variability influenced by microtopography and vegetation distribution in a subarctic watershed. Included datasets are high-resolution spatial maps of snow depth, terrain elevation, vegetation height, and distance from shrubs taller than 1 meter, all formatted as text files (.txt). These data are fully describe in doi:10.15485/2316038. Users can adapt the provided R scripts to accommodate different data formats or larger spatial domains, noting that some output files may require modification due to their size.The code includes implementations for boosted regression tree analysis adapted from methods outlined in Elith et al. (2008). Users interested in understanding or modeling landscape-scale snow distribution patterns, particularly in Arctic or subarctic ecosystems, will find this package useful. The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research.The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska.Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Annihilation-gamma-based diagnostic techniques for magnetically confined electron–positron pair plasma

Efforts are underway to magnetically confine electron–positron pair plasmas to study their unique behaviour, which is characterized by significant changes in plasma time and length scales, supported waves and unstable modes. However, use of conventional plasma diagnostics presents challenges with these low-density and annihilating matter–antimatter plasmas. To address this problem, we propose to develop techniques based on the distinct emission provided by annihilation. This emission exhibits two spatial correlations: the distance attenuation of isotropic sources and the back-to-back propagation of momentum-preserving 2 $\gamma$ annihilation. We present the results of our analysis of the $\gamma$ emission rate and the spatial profile of the annihilation in a magnetized pair plasma from direct pair collisions, from the formation and decay of positronium as well as from transport processes. In order to demonstrate the effectiveness of annihilation-based techniques, we tested them on annular $\gamma$ emission profiles produced by a $\beta ^+$ radioisotope on a rotating turntable. Direct and positronium-mediated annihilation result in overlapping volumetric $\gamma$ sources, and the 2 $\gamma$ emission from these volumetric sources can be tomographically reconstructed from coincident counts in multiple detectors. Transport processes result in localized annihilation where field lines intersect walls, limiters or internal magnets. These localized sources can be identified by the fractional $\gamma$ counts on spatially distributed detectors.

Physics↗
DOE OSTI2021

A Spatially Explicit Individual-Based Modeling Approach to Evaluate the Cumulative Effects of Wind Energy Development on the Greater Sage-Grouse: Sensitivity Analysis and Validation of Key Parameters

To address the need for an approach to evaluate cumulative ecological impacts of wind energy development, particularly those on critical wildlife habitats and species, the authors developed a proof-of-concept, landscape-based, spatially explicit individual-based modeling (IBM) framework for populations of the greater sage-grouse (Centrocercus urophasianus) in Albany County, Wyoming, based on published and other available information on the life history of the species. This sage-grouse IBM estimates the spatial and temporal movements of sage-grouse based on modeled effects of habitat characteristics and developments on sage-grouse condition, survivorship, and reproduction. After completing the initial version of the model, we conducted a series of sensitivity analyses of key model parameters and validation of lek occurrences and spatial distributions to better understand performance of the model relative to observed values of important life history characteristics. The parameter sensitivity analysis was conducted as a series of one-factor-at-a-time tests in which the values of selected parameters were individually varied to evaluate the corresponding effects on model output. We examined the model’s response to changes in parameter values to gain insight into the overall robustness of the model, identify key driver parameters, and initiate qualitative validation. Five functions selected for sensitivity analyses were (1) the rate of change in body condition in relationship to habitat suitability index (HSI) of occupied habitat (condition change rate), (2) the effect of competition between individuals with overlapping home ranges expressed as a percentage reduction in HSI (competition factor), (3) value of condition at which survivorship probability is equal to zero (zero survivorship condition), (4) the relationship between female body condition and clutch size, and (5) the relationship between female body condition and nest success. We varied the parameter values to be tested for the sensitivity analysis around the original value used in the model (i.e., tested values were higher and lower than the model’s original value) to examine the effect of this variation on the predicted output of the model. Model outputs evaluated in the sensitivity analysis included (1) predicted age-class distribution (as measured by the percentage of immature individuals within the population); (2) predicted sex ratio; (3) predicted age-sex class distribution (percent of population in each of six age-sex classes, i.e., male and female juveniles, yearlings, and adults); (4) lifespan; (5) population size; and (6) spatial distribution. Of the five parameters analyzed, condition change rate, survivorship, clutch size, and nest success did not substantially affect any of the evaluated model outputs. For the predicted age-class distribution, the percentage of immature individuals predicted by the model for these four parameters ranged from 61.9% to 65.0%, which was reasonable with respect to the reference values of 51.4 to 57.8% in published studies. The ratio of females to males ranged from 1.09 to 1.14, which was comparable to the reference values of 1.2 to 3.0 presented in an existing study. The predicted mean lifespan ranged from 1.62 to 1.74 years, which appears to be reasonable with respect to the reference values of 0.9 to 1.1 years for sharp-tail grouse and greater prairie chickens presented by published studies. The predicted population size ranged from 1,877 to 3,140. This prediction would be reasonable with respect to our estimate of 5,000 yearlings and adults, which was based on the USFWS estimate for the state of Wyoming scaled to the number of leks in the county. The predicted spatial distribution was not substantially affected by the parameter values examined. In our sensitivity analyses, the only tested parameter that had a noticeable effect on model output was competition factor. Over the range of values tested for this parameter, predicted population size exhibited the largest range (from 1,833 to 4,737) of all parameters tested. As the competition factor increased, the predicted population size decreased. This reflects the effect competition had on the number of individuals that could co occur in high-quality habitat patches. Despite this apparently significant effect, the range of values appears to be reasonable with respect to our estimate of 5,000 yearlings and adults, which was based on the USFWS estimate for the state of Wyoming scaled to the number of leks in the county.

17 WIND ENERGY↗
NASA NTRS2007

Theorizing Land Cover and Land Use Change: The Peasant Economy of Colonization in the Amazon Basin

This paper addresses deforestation processes in the Amazon basin. It deploys a methodology combining remote sensing and survey-based fieldwork to examine, with regression analysis, the impact household structure and economic circumstances on deforestation decisions made by colonist farmers in the forest frontiers of Brazil. Unlike most previous regression-based studies, the methodology implemented analyzes behavior at the level of the individual property. The regressions correct for endogenous relationships between key variables, and spatial autocorrelation, as necessary. Variables used in the analysis are specified, in part, by a theoretical development integrating the Chayanovian concept of the peasant household with spatial considerations stemming from von Thuenen. The results from the empirical model indicate that demographic characteristics of households, as well as market factors, affect deforestation in the Amazon. Thus, statistical results from studies that do not include household-scale information may be subject to error. From a policy perspective, the results suggest that environmental policies in the Amazon based on market incentives to small farmers may not be as effective as hoped, given the importance of household factors in catalyzing the demand for land. The paper concludes by noting that household decisions regarding land use and deforestation are not independent of broader social circumstances, and that a full understanding of Amazonian deforestation will require insight into why poor families find it necessary to settle the frontier in the first place.

Caldas, Marcellus↗
NASA NTRS1980

Spatial frequency response of an optical heterodyne receiver

The analysis of a passive heterodyne receiver with respect to its imaging performance (transfer function) and its heterodyne efficiency shows departures from the results which are obtained in strictly coherent or incoherent imaging systems. The cascading property of modulation transfer function analysis must be carefully applied, since the coherent transfer function of the optical receiver and that due to the local oscillator-detector combination are not separable but are related by the convolution of their products. Application of these results to a spacelab-type optical heterodyne receiver shows that resolutions of the order of 1.5 to 2.0 km are possible for worst-case type orbital scenarios. In addition to an analysis of the derogatory effects of sampling (aliasing) a comparison of obscured-type receivers (e.g., Cassegrains) with unobscured receivers shows that both resolution and efficiency are severly degraded in an obscured-type receiver and, consequently, should not be used for passive heterodyne detection schemes.

Robinson, D. M.↗
NASA NTRS1980

Evaluation of SLAR and thematic mapper MSS data for forest cover mapping using computer-aided analysis techniques

To facilitate comparison between the four different spatial resolution of the NS-001 MSS data sets, a supervised approach was taken in defining training blocks for each of the different cover types. The training fields representing each cover type category were grouped and this group was clustered to determine the individual spectral classes within each cover type category which would effectively characterize the entire test site. Graphs show the variation in spectral response level with respect to distance in the across track dimension for four sampling intervals. Radar digitization procedures were developd. Flight characteristics and parameters for digitization of radar imagery are tabulated. The statement of work for phase 3 was reviewed and modifications were suggested to meet funding reduction.

Hoffer, R. M.↗
NASA NTRS2013

Global Analysis of Aerosol Properties Above Clouds

The seasonal and spatial varability of Aerosol Above Cloud (AAC) properties are derived from passive satellite data for the year 2008. A significant amount of aerosols are transported above liquid water clouds on the global scale. For particles in the fine mode (i.e., radius smaller than 0.3 m), including both clear sky and AAC retrievals increases the global mean aerosol optical thickness by 25(+/- 6%). The two main regions with man-made AAC are the tropical Southeast Atlantic, for biomass burning aerosols, and the North Pacific, mainly for pollutants. Man-made AAC are also detected over the Arctic during the spring. Mineral dust particles are detected above clouds within the so-called dust belt region (5-40 N). AAC may cause a warming effect and bias the retrieval of the cloud properties. This study will then help to better quantify the impacts of aerosols on clouds and climate.

Seasonal spatial varability of AAC↗
NASA NTRS2019

Alternative Vegetation States in Tropical Forests and Savannas: The Search for Consistent Signals in Diverse Remote Sensing Data

Globally, the spatial distribution of vegetation is governed primarily by climatological factors (rainfall and temperature, seasonality, and inter-annual variability). The local distribution of vegetation, however, depends on local edaphic conditions (soils and topography) and disturbances (fire, herbivory, and anthropogenic activities). Abrupt spatial or temporal changes in vegetation distribution can occur if there are positive (i.e., amplifying) feedbacks favoring certain vegetation states under otherwise similar climatic and edaphic conditions. Previous studies in the tropical savannas of Africa and other continents using the MODerate Resolution Imaging Spectroradiometer (MODIS) vegetation continuous fields (VCF) satellite data product have focused on discontinuities in the distribution of tree cover at different rainfall levels, with bimodal distributions (e.g., concentrations of high and low tree cover) interpreted as alternative vegetation states. Such observed bimodalities over large spatial extents may not be evidence for alternate states, as they may include regions that have different edaphic conditions and disturbance histories. In this study, we conduct a systematic multi-scale analysis of diverse MODIS data streams to quantify the presence and spatial consistency of alternative vegetation states in Sub-Saharan Africa. The analysis is based on the premise that major discontinuities in vegetation structure should also manifest as consistent spatial patterns in a range of remote sensing data streams, including, for example, albedo and land surface temperature (LST). Our results confirm previous observations of bimodal and multimodal distributions of estimated tree cover in the MODIS VCF. However, strong disagreements in the location of multimodality between VCF and other data streams were observed at 1 km scale. Results suggest that the observed distribution of VCF over vast spatial extents are multimodal, not because of local-scale feedbacks and emergent bifurcations (the definition of alternative states), but likely because of other factors including regional scale differences in woody dynamics associated with edaphic, disturbance, and/or anthropogenic processes. These results suggest the need for more in-depth consideration of bifurcation mechanisms and thus the likely spatial and temporal scales at which alternative states driven by different positive feedback processes should manifest.

Savanna↗
NASA NTRS1986

Hybrid analog/digital, large format, photon counting detectors for astronomy

The development of a new microchannel plate photon-counting detector with an analog readout method based on a resistive anode is reported. This detector exhibits extremely high, stable electron gains of 10 to the 8th. At this gain, the spatial resolution is no longer primarily limited by the noise of the resistive anode, so that digital methods of readout, such as discrete conductors, lose their advantage. These detectors can be readily scaled to 40 mm and 70 mm formats to match plate scales of 2-m (and larger) telescopes. New, high speed digital electronics fully exploit the high spatial and time resolution made possible by gains of this level. Analysis of the theoretical performance of this detector shows that the major limitation to the spatial resolution is the proximity focus of the photocathode and the first microchannel plate. The detector has been mated to an echelle spectrograph developed.

Crocker, J.↗
NASA NTRS1995

Assimilation of TOPEX Sea Level Measurements with a Reduced-Gravity, Shallow Water Model of the Tropical Pacific Ocean

Sea surface height variability measured by TOPEX is analyzed in the tropical Pacific Ocean by way of assimilation into a wind-driven, reduced-gravity, shallow water model using an approximate Kalman filter and smoother. The analysis results in an optimal fit of the dynamic model to the observations, providing it dynamically consistent interpolation of sea level and estimation of the circulation. Nearly 80% of the expected signal variance is accounted for by the model within 20 deg of the equator, and estimation uncertainty is substantially reduced by the voluminous observation. Notable features resolved by the analysis include seasonal changes associated with the North Equatorial Countercurrent and equatorial Kelvin and Rossby waves. Significant discrepancies are also found between the estimate and TOPEX measurements, especially near the eastern boundary. Improvements in the estimate made by the assimilation are validated by comparisons with independent tide gauge and current meter observations. The employed filter and smoother are based on approximately computed estimation error covariance matrices, utilizing a spatial transformation and an symptotic approximation. The analysis demonstrates the practical utility of a quasi-optimal filter and smoother.

Fukumori, Ichiro↗
DOE OSTI2024

Source apportionment of airborne volatile organic compounds near unconventional oil and gas development

Oil and natural gas (ONG) extraction emits volatile organic compounds (VOCs). Certain VOCs are identified as hazardous air pollutants (HAPS) while others contribute to ozone formation. This study examines the impact of ONG operations on VOC levels during the development of multi-well ONG pads in suburban Broomfield, Colorado. From October 2018 to December 2020, weekly VOC measurements were taken at 18 sites across the area. These included spots near well pads, in adjacent neighborhoods, and at a background site, covering various stages of well pad development including drilling, hydraulic fracturing, flowback, and production. Analysis using Positive Matrix Factorization (PMF) identified six factors, including combustion, background/biogenic sources, light and complex alkanes, drilling activities, and ONG acetylene. Factors linked to local ONG activities exhibited clear temporal and spatial correlations with Broomfield well development. Benzene source analysis revealed distinct contribution gradients, with ONG-related sources notably influencing areas near the well pads, particularly in pre-production. ONG-related weekly benzene contributions varied from 9% to 63% at a community background site and 18% to 89% in a neighborhood close to a well pad.

54 ENVIRONMENTAL SCIENCES↗