Search NASA⌕ Search

SEARCH · Search NASA

Results for “spatial data visualization”

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.

427 records · Page 24

Geometric validation of self-gating k-space-sorted 4D-MRI vs 4D-CT using a respiratory motion phantom

Purpose: MRI is increasingly being used for radiotherapy planning, simulation, and in-treatment-room motion monitoring. To provide more detailed temporal and spatial MR data for these tasks, we have recently developed a novel self-gated (SG) MRI technique with advantage of k-space phase sorting, high isotropic spatial resolution, and high temporal resolution. The current work describes the validation of this 4D-MRI technique using a MRI- and CT-compatible respiratory motion phantom and comparison to 4D-CT. Methods: The 4D-MRI sequence is based on a spoiled gradient echo-based 3D projection reconstruction sequence with self-gating for 4D-MRI at 3 T. Respiratory phase is resolved by using SG k-space lines as the motion surrogate. 4D-MRI images are reconstructed into ten temporal bins with spatial resolution 1.56 × 1.56 × 1.56 mm{sup 3}. A MRI-CT compatible phantom was designed to validate the performance of the 4D-MRI sequence and 4D-CT imaging. A spherical target (diameter 23 mm, volume 6.37 ml) filled with high-concentration gadolinium (Gd) gel is embedded into a plastic box (35 × 40 × 63 mm{sup 3}) and stabilized with low-concentration Gd gel. The phantom, driven by an air pump, is able to produce human-type breathing patterns between 4 and 30 respiratory cycles/min. 4D-CT of the phantom has been acquired in cine mode, and reconstructed into ten phases with slice thickness 1.25 mm. The 4D images sets were imported into a treatment planning software for target contouring. The geometrical accuracy of the 4D MRI and CT images has been quantified using target volume, flattening, and eccentricity. The target motion was measured by tracking the centroids of the spheres in each individual phase. Motion ground-truth was obtained from input signals and real-time video recordings. Results: The dynamic phantom has been operated in four respiratory rate (RR) settings, 6, 10, 15, and 20/min, and was scanned with 4D-MRI and 4D-CT. 4D-CT images have target-stretching, partial-missing, and other motion artifacts in various phases, whereas the 4D-MRI images are visually free of those artifacts. Volume percentage difference for the 6.37 ml target ranged from 5.3% ± 4.3% to 10.3% ± 5.9% for 4D-CT, and 1.47 ± 0.52 to 2.12 ± 1.60 for 4D-MRI. With an increase of respiratory rate, the target volumetric and geometric deviations increase for 4D-CT images while remaining stable for the 4D-MRI images. Target motion amplitude errors at different RRs were measured with a range of 0.66–1.25 mm for 4D-CT and 0.2–0.42 mm for 4D-MRI. The results of Mann–Whitney tests indicated that 4D-MRI significantly outperforms 4D-CT in phase-based target volumetric (p = 0.027) and geometric (p < 0.001) measures. Both modalities achieve equivalent accuracy in measuring motion amplitude (p = 0.828). Conclusions: The k-space self-gated 4D-MRI technique provides a robust method for accurately imaging phase-based target motion and geometry. Compared to 4D-CT, the current 4D-MRI technique demonstrates superior spatiotemporal resolution, and robust resistance to motion artifacts caused by fast target motion and irregular breathing patterns. The technique can be used extensively in abdominal targeting, motion gating, and toward implementing MRI-based adaptive radiotherapy.

60 APPLIED LIFE SCIENCES↗

High spatial resolution and contrast radiography of hydrodynamic instabilities at the National Ignition Facility

In this work, we are developing techniques for studying the Rayleigh–Taylor (RT) and Richtmyer–Meshkov (RM) instabilities in a planar geometry at high-energy-densities at the National Ignition Facility (NIF). In particular, through the improvement of experimental imaging quality, we are progressing toward the study of the turbulent regime of the mixing regions in capsule implosion experiments for inertial confinement fusion, which requires few micrometers resolution. Using 60 NIF beams, a solid shock tube is driven launching a shock wave that crosses the interface between a dense and a light material pre-machined in the target to obtain sinusoidal ripples, which results in RM and RT instabilities that are imaged using the NIF Crystal Backlighter Imager. High-quality images were obtained with a mean resolution of 7 μm and improved contrast. While the obtained resolution does not allow the observation of the smallest scale of the “turbulent” energy spectrum, the generated image encompasses 63% of the total flow energy, a 50% improvement over previous studies, which is observed for the first time a roll-up feature in a high energy density-type RT experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Differences in soil biological activity by terrain types at the sub-field scale in central Iowa US

Soil microbial communities are structured by biogeochemical processes that occur at many different spatial scales, which makes soil sampling difficult. Because soil microbial communities are important in nutrient cycling and soil fertility, it is important to understand how microbial communities function within the heterogeneous soil landscape. In this study, a self-organizing map was used to determine whether landscape data can be used to characterize the distribution of microbial biomass and activity in order to provide an improved understanding of soil microbial community function. Points within a row crop field in south-central Iowa were clustered via a self-organizing map using six landscape properties into three separate landscape clusters. Twelve sampling locations per cluster were chosen for a total of 36 locations. After the soil samples were collected, the samples were then analysed for various metabolic indicators, such as nitrogen and carbon mineralization, extractable organic carbon, microbial biomass, etc. It was found that sampling locations located in the potholes and toe slope positions had significantly greater microbial biomass nitrogen and carbon, total carbon, total nitrogen and extractable organic carbon than the other two landscape position clusters, while locations located on the upslope did not differ significantly from the other landscape clusters. However, factors such as nitrate, ammonia, and nitrogen and carbon mineralization did not differ significantly across the landscape. Altogether, this research demonstrates the effectiveness of a terrain-based clustering method for guiding soil sampling of microbial communities.

54 ENVIRONMENTAL SCIENCES↗

SOUTHERN MASSIVE STARS AT HIGH ANGULAR RESOLUTION: OBSERVATIONAL CAMPAIGN AND COMPANION DETECTION

Multiplicity is one of the most fundamental observable properties of massive O-type stars and offers a promising way to discriminate between massive star formation theories. Nevertheless, companions at separations between 1 and 100 milliarcsec (mas) remain mostly unknown due to intrinsic observational limitations. At a typical distance of 2 kpc, this corresponds to projected physical separations of 2-200 AU. The Southern MAssive Stars at High angular resolution survey (SMaSH+) was designed to fill this gap by providing the first systematic interferometric survey of Galactic massive stars. We observed 117 O-type stars with VLTI/PIONIER and 162 O-type stars with NACO/Sparse Aperture Masking (SAM), probing the separation ranges 1-45 and 30-250 mas and brightness contrasts of ΔH < 4 and ΔH < 5, respectively. Taking advantage of NACO's field of view, we further uniformly searched for visual companions in an 8'' radius down to ΔH = 8. This paper describes observations and data analysis, reports the discovery of almost 200 new companions in the separation range from 1 mas to 8'' and presents a catalog of detections, including the first resolved measurements of over a dozen known long-period spectroscopic binaries. Excluding known runaway stars for which no companions are detected, 96 objects in our main sample (δ < 0°; H < 7.5) were observed both with PIONIER and NACO/SAM. The fraction of these stars with at least one resolved companion within 200 mas is 0.53. Accounting for known but unresolved spectroscopic or eclipsing companions, the multiplicity fraction at separation ρ < 8'' increases to f {sub m} = 0.91 ± 0.03. The fraction of luminosity class V stars that have a bound companion reaches 100% at 30 mas while their average number of physically connected companions within 8'' is f {sub c} = 2.2 ± 0.3. This demonstrates that massive stars form nearly exclusively in multiple systems. The nine non-thermal radio emitters observed by SMaSH+ are all resolved, including the newly discovered pairs HD 168112 and CPD–47°2963. This lends strong support to the universality of the wind-wind collision scenario to explain the non-thermal emission from O-type stars.

79 ASTRONOMY AND ASTROPHYSICS↗

Accounting for hydrological and microbial processes on greenhouse gas budgets from river systems (Final Report)

River systems have a disproportionately large impact on carbon and nitrogen cycling on this planet. For instance, despite their much smaller size in area, it is estimated that global carbon dioxide emission from river systems exceed the emission from lakes and reservoirs, and are on par with the estimated oceanic and land carbon sinks. Closing carbon and nitrogen budgets in river systems requires knowledge of end member greenhouse gas (GHG) products (e.g. carbon dioxide, nitrous oxide, methane), which typically have not been measured in river systems. Moreover, while the production of these gasses is thought to be microbially derived, little is known about microbial activities in sediments of large rivers, despite their potentially high influence on biogeochemical budgets. Typically, when we think of rivers we visualize the channels that carve through the earth, but this only represents a small part of the hydrologic exchange that a river participates in. Rivers are more than conduits; they are dynamically exchanging water and materials with surrounding terrestrial and subsurface environments. In fact, 96% of the carbon dioxide produced within river systems occurs in hyporheic zones, or sediment regions beneath and alongside a river where groundwater and river water mix. In fact, it is widely recognized that microbial processes in hyporheic zones, also known as the river’s liver, strongly influence the fate of carbon, nutrients, and contaminants transported through watersheds. Despite this general knowledge, observational studies linking nitrogen and carbon gas production to microbial activity, and studies of overall gas emissions, are limited across river systems, but especially absent from larger rivers. To identify microbial processes responsible for GHG emissions from large rivers, and to quantify their environmental drivers and dependencies, we selected the Columbia River as a representative field site. This is the largest river west of the North American Continental Divide, and produces a drainage and discharge roughly double of the Nile in Egypt. This river is also an important source for agricultural irrigation and generates more hydroelectric power than any other North American river. Dams manage the release of water in the river, with these changes in river discharge patterns strongly impacting microbial activities in hyporheic zones. Specifically, we target a region down-river from the Priest Rapids Dam (WA) associated with the PNNL SFA, where water discharge from the dam is highly controlled, resulting in 1.5 m fluctuations in river stage daily. This dynamic river stage creates an ideal research opportunity for investigating hydrological perturbations in GHG production, uptake, and emission from large river systems. In this exploratory 1-year project, we conducted data-driven field research to provide observations needed to support predictions about GHG flux from river ecosystems. Using porewater dialysis and surface flux measurements, we observed where and when GHGs were produced and is the emission impacted by river stage. Using meta-omics tools, we determined what microbial processes control production and consumption of GHG in these systems. Finally, we provided key observations to parameterize models that could predict the effect of environmental conditions on GHG emission in the Columbia River near shore environment, at the scale that the modeled processes operate. Closing knowledge gaps in the linkages between river stage, microbial processes, and GHG production and emission, as well as the spatial distribution and ecosystem dependencies of these processes, is necessary to better parameterize land surface models. We provide some of the first, detailed GHG measurements paired microbial activity from large river systems. Our findings have direct relevance to many other large rivers with flows altered by hydroelectric energy, reservoirs, snowmelt, tides, or agriculture.

54 ENVIRONMENTAL SCIENCES↗

Carbon Management Projects (CONNECT) Database and Explorer

Overview The Carbon Management Projects (CONNECT) Toolkit is an online exploratory visualization tool developed by the U.S. Department of Energy's (DOE) Office of Fossil Energy and Carbon Management (FECM) with support from other federal agencies such as the U.S. Environmental Protection Agency (EPA) and the U.S. Department of Transportation (DOT). It provides a single point of access to authoritative information on federal agency investment in a portfolio of research, development, and demonstration (RD&D) projects that have been publicly announced to advance technologies for point source carbon capture, carbon dioxide removal, transport, storage, and conversion, collectively referred to as carbon management. The RD&D programs covered in this tool are authorized by annual congressional appropriations ("Base Program") and the 2021 Infrastructure Investment and Jobs Act (IIJA). The tool also incorporates public information on other federal initiatives, such as the Regional Clean Hydrogen Hubs, and public information released by other government agencies, such as the Environmental Protection Agency's (EPA) and Primacy States’ Underground Injection Control Class VI permits and EPA’s facility level greenhouse gas (GHG) emissions. Developed in a geographic information system, the tool organizes carbon management projects into five groups based on the primary technology that a project aims to advance, each visually represented as a digital layer ("carbon management project layer"). Only federally funded projects are included, which can be awarded projects that are completed or ongoing, or projects that have been selected but are currently under negotiation. Project information can be viewed in the map or in the attribute table below it when turned on. In the map view, each project is displayed at either its host site (for field work), where available, or its performer site (project lead's location, further explained in the table below). Host sites and performer sites are represented in distinct icons. Several reference layers offer additional public information on infrastructural and natural resource environment for carbon management. These reference layers, combined with multiple geographical basemaps, enable users to visualize the carbon management project layers in context. Carbon management project information will be updated monthly based on feedback and information availability. Carbon management project layers Point Source Carbon Capture (PSC) This layer contains DOE-funded projects focused on capturing carbon dioxide (CO2) from power plants or industrial facilities. Carbon Dioxide Removal (CDR) This layer contains DOE-funded projects focused on capturing CO2 from the atmosphere, including direct air capture (DAC) and DAC hubs, direct ocean capture, enhanced mineralization, and biomass carbon removal and storage. For projects with multiple host sites, each of the sites are displayed individually with the project cost and cost sharing information representing the total for the entire project. Carbon Transport This layer contains DOE- and DOT-funded projects focused on CO2 transport. The Transport Research and Development sublayer contains projects that do not involve physical infrastructure; the Proposed Transport Corridor sublayer contains projects for which either a route for the transport infrastructure has been proposed or a general area for the transport infrastructure has been identified. Carbon Storage This layer contains DOE-funded key projects focused on CO2 storage. For projects with multiple field-work sites, each of the sites are displayed individually on the map with the project cost and cost sharing information representing the overall total for the entire project. Carbon Conversion This layer contains DOE-funded projects focused on converting CO2 into economically valuable products. Reference layers The following layers provide additional information in the geographic proximity of carbon management projects. Users should reference the original sources for more details (weblinks provided below and in pop-up windows on the map). Regional Clean Hydrogen Hub and Facility These layers illustrate the approximate areas of the Regional Clean Hydrogen Hubs announced by DOE's Office of Clean Energy Demonstrations (OCED) and the approximate locations of individual facilities that constitute the hubs (see "Where are the H2Hubs located?" on the webpage linked above). EPA Facility Level GHG Emissions (direct emitter) This layer shows direct CO2 emissions from stationary sources in 2022, using data extracted from EPA's Facility Level Information on GreenHouse gases Tool (FLIGHT). Captured and injected CO2 are not deducted from direct emitters’ total emissions. Contact EPA for additional details. Underground Injection Control Class VI permit/permit application This layer shows the locations of CO2 injection wells that are granted or in the process of applying for an Underground Injection Control Class VI permit by EPA or a Primacy State (currently Louisiana, North Dakota, and Wyoming). The URLs for the permits or permit applications are provided in the pop-up windows associated with the well locations. Contact EPA for additional details. Carbon Storage Resource This layer contains information on prospective CO2 storage resources in saline formations and oil and gas reservoirs provided by the National Carbon Sequestration Database and Geographic Information System (NATCARB) spatial database. Contact NETL for additional details. Existing CO2 pipeline This layer shows active CO2 pipelines based on information digitized from the map issued by the Pipeline and Hazardous Materials Safety Administration (PHMSA). Contact PHMSA for additional details.

Carbon Conversion↗

How does C-VIEW image quality compare with conventional 2D FFDM?

Purpose: The FDA approved the use of digital breast tomosynthesis (DBT) in 2011 as an adjunct to 2D full field digital mammography (FFDM) with the constraint that all DBT acquisitions must be paired with a 2D image to assure adequate interpretative information is provided. Recently manufacturers have developed methods to provide a synthesized 2D image generated from the DBT data with the hope of sparing patients the radiation exposure from the FFDM acquisition. While this much needed alternative effectively reduces the total radiation burden, differences in image quality must also be considered. The goal of this study was to compare the intrinsic image quality of synthesized 2D C-VIEW and 2D FFDM images in terms of resolution, contrast, and noise. Methods: Two phantoms were utilized in this study: the American College of Radiology mammography accreditation phantom (ACR phantom) and a novel 3D printed anthropomorphic breast phantom. Both phantoms were imaged using a Hologic Selenia Dimensions 3D system. Analysis of the ACR phantom includes both visual inspection and objective automated analysis using in-house software. Analysis of the 3D anthropomorphic phantom includes visual assessment of resolution and Fourier analysis of the noise. Results: Using ACR-defined scoring criteria for the ACR phantom, the FFDM images scored statistically higher than C-VIEW according to both the average observer and automated scores. In addition, between 50% and 70% of C-VIEW images failed to meet the nominal minimum ACR accreditation requirements—primarily due to fiber breaks. Software analysis demonstrated that C-VIEW provided enhanced visualization of medium and large microcalcification objects; however, the benefits diminished for smaller high contrast objects and all low contrast objects. Visual analysis of the anthropomorphic phantom showed a measureable loss of resolution in the C-VIEW image (11 lp/mm FFDM, 5 lp/mm C-VIEW) and loss in detection of small microcalcification objects. Spectral analysis of the anthropomorphic phantom showed higher total noise magnitude in the FFDM image compared with C-VIEW. Whereas the FFDM image contained approximately white noise texture, the C-VIEW image exhibited marked noise reduction at midfrequency and high frequency with far less noise suppression at low frequencies resulting in a mottled noise appearance. Conclusions: Their analysis demonstrates many instances where the C-VIEW image quality differs from FFDM. Compared to FFDM, C-VIEW offers a better depiction of objects of certain size and contrast, but provides poorer overall resolution and noise properties. Based on these findings, the utilization of C-VIEW images in the clinical setting requires careful consideration, especially if considering the discontinuation of FFDM imaging. Not explicitly explored in this study is how the combination of DBT + C-VIEW performs relative to DBT + FFDM or FFDM alone.

60 APPLIED LIFE SCIENCES↗

EXCHANGE Campaign 1: A Community-Driven Baseline Characterization of Soils, Sediments, and Water Across Coastal Gradients

The EXploration of Coastal Hydrobiogeochemistry Across a Network of Gradients and Experiments (EXCHANGE) program is a consortium of scientists working together to improve our understanding of how the two-way exchange of water between estuaries or large lake lacustuaries and the terrestrial landscape influence the state and function of ecosystems across the coastal interface. EXCHANGE Campaign 1 (EC1) focuses on the spatial variation in biogeochemical structure and function at the coastal terrestrial-aquatic interface (TAI). In the Fall of 2021, the EXCHANGE Consortium gathered samples from 52 TAIs. Samples collected from EC1 were analyzed for bulk geochemical parameters, bulk physicochemical parameters, organic matter characteristics, and redox-sensitive elements.Please download ec1_README.pdf for a complete list of available data in each .zip folder, package version history, and detailed information about the project. This README will serve as the central place for EC1 Data Package updates. Experimental setup and v1 methods are documented in Myers-Pigg and Pennington et al., 2023 (https://doi.org/10.1038/s41597-023-02548-7).EC1 Data Package Structure:ec1_README.pdfec1_methods.pdfec1_metadata_v3.zip...ec1_dd.csv...ec1_flmd.csv...ec1_sample_catalog.csv...ec1_metadata_kitlevel.csv...ec1_metadata_collectionlevel.csv...ec1_data_collectionlevel.csv...ec1_igsn_metadata.csvec1_soil_v3.zipec1_sediment_v3.zipec1_water_v3.zipec1_processingscripts_v3.zipThis data package is on v3 and was originally published May 2023 (v1). Subsequent updates will be published here with new version numbers. Please see the Change History section in ec1_README.pdf for detailed changes.---Acknowledging EXCHANGE: General Support and Data Product UseWe ask that users of EXCHANGE data add the following acknowledgement when publishing data in scholarly articles and data repositories:"This research is based on work supported by COMPASS-FME, a multi-institutional project supported by the U.S. Department of Energy, Office of Science, Biological and Environmental Research as part of the Environmental System Science Program."

54 ENVIRONMENTAL SCIENCES↗

Coastal Ohio Wind Project

The Coastal Ohio Wind Project intends to address problems that impede deployment of wind turbines in the coastal and offshore regions of Northern Ohio. The project evaluates different wind turbine designs and the potential impact of offshore turbines on migratory and resident birds by developing multidisciplinary research, which involves wildlife biology, electrical and mechanical engineering, and geospatial science. Firstly, the project conducts cost and performance studies of two- and three-blade wind turbines using a turbine design suited for the Great Lakes. The numerical studies comprised an analysis and evaluation of the annual energy production of two- and three-blade wind turbines to determine the levelized cost of energy. This task also involved wind tunnel studies of model wind turbines to quantify the wake flow field of upwind and downwind wind turbine-tower arrangements. The experimental work included a study of a scaled model of an offshore wind turbine platform in a water tunnel. The levelized cost of energy work consisted of the development and application of a cost model to predict the cost of energy produced by a wind turbine system placed offshore. The analysis found that a floating two-blade wind turbine presents the most cost effective alternative for the Great Lakes. The load effects studies showed that the two-blade wind turbine model experiences less torque under all IEC Standard design load cases considered. Other load effects did not show this trend and depending on the design load cases, the two-bladed wind turbine showed higher or lower load effects. The experimental studies of the wake were conducted using smoke flow visualization and hot wire anemometry. Flow visualization studies showed that in the downwind turbine configuration the wake flow was insensitive to the presence of the blade and was very similar to that of the tower alone. On the other hand, in the upwind turbine configuration, increasing the rotor blade angle of attack reduced the wake size and enhanced the vortices in the flow downstream of the turbine-tower compared with the tower alone case. Mean and rms velocity distributions from hot wire anemometer data confirmed that in a downwind configuration, the wake of the tower dominates the flow, thus the flow fields of a tower alone and tower-turbine combinations are nearly the same. For the upwind configuration, the mean velocity shows a narrowing of the wake compared with the tower alone case. The downwind configuration wake persisted longer than that of an upwind configuration; however, it was not possible to quantify this difference because of the size limitation of the wind tunnel downstream of the test section. The water tunnel studies demonstrated that the scale model studies could be used to adequately produce accurate motions to model the motions of a wind turbine platform subject to large waves. It was found that the important factors that affect the platform is whether the platform is submerged or surface piercing. In the former, the loads on the platform will be relatively reduced whereas in the latter case, the structure pierces the wave free surface and gains stiffness and stability. The other important element that affects the movement of the platform is depth of the sea in which the wind turbine will be installed. Furthermore, the wildlife biology component evaluated migratory patterns by different monitoring systems consisting of marine radar, thermal IR camera and acoustic recorders. The types of radar used in the project are weather surveillance radar and marine radar. The weather surveillance radar (1988 Doppler), also known as Next Generation Radar (NEXRAD), provides a network of weather stations in the US. Data generated from this network were used to understand general migratory patterns, migratory stopover habitats, and other patterns caused by the effects of weather conditions. At a local scale our marine radar was used to complement the datasets from NEXRAD and to collect additional monitoring parameters such as passage rates, flight paths, flight directions, and flight altitudes of nocturnal migrating species. Our work focused on the design and development of custom built marine radar that used t-bar and parabolic dish antennas. The marine radar used in the project was Furuno (XANK250) which was coupled with a XIR3000B digitizing card from Russell Technologies for collection of the radar data. The radar data was processed by open source radR processing software using different computational techniques and methods. Additional data from thermal IR imaging cameras were collected to detect heat emitted from objects and provide information on movements of birds and bats, data which we used for different animal flight behavior analysis. Lastly, the data from the acoustic recorders were used to provide the number of bird calls for assessing patterns and peak passage rates during migration. The development of the geospatial database included collection of different data sources that are used to support offshore wind turbine development. Many different data sets were collected and organized using initial version of web-based repository software tools that can accommodate distribution of rectified pertinent data sets such as the lake depth, lake bottom engineering parameters, extent of ice, navigation pathways, wind speed, important bird habitats, fish efforts and other layers that are relevant for supporting robust offshore wind turbine developments. Additional geospatial products developed during the project included few different prototypes for offshore wind farm suitability which can involve different stakeholders and participants for solving complex planning problems and building consensus. Some of the prototypes include spatial decision support system (SDSS) for collaborative decision making, a web-based Participatory Geographic Information System (PGIS) framework for evaluating importance of different decision alternatives using different evaluation criteria, and an Android application for collection of field data using mobile and tablet devices . In summary, the simulations of two- and three-blade wind turbines suggested that two-bladed machines could produce comparable annual energy as the three-blade wind turbines but have a lighter tower top weight, which leads to lower cost of energy. In addition, the two-blade rotor configuration potentially costs 20% less than a three blade configuration that produces the same power at the same site. The cost model analysis predicted a potential cost savings of approximately 15% for offshore two-blade wind turbines. The foundation design for a wind turbine in Lake Erie is likely to be driven by ice loads based on the currently available ice data and ice mechanics models. Hence, for Lake Eire, the cost savings will be somewhat smaller than the other lakes in the Great Lakes. Considering the size of cranes and vessels currently available in the Great Lakes, the cost optimal wind turbine size should be 3 MW, not larger. The surveillance data from different monitoring systems suggested that bird and bat passage rates per hour were comparable during heavy migrations in both spring and fall seasons while passage rates were significantly correlated to wind directions and wind speeds. The altitude of migration was higher during heavy migrations and higher over water relative to over land. Notable portions of migration on some spring nights occurred parallel the shoreline, often moving perpendicular to southern winds. The birds approaching the Western basin have a higher propensity to cross than birds approaching the Central basin of Lake Erie and as such offshore turbine development might be a better option further east towards Cleveland than in the Western basin. The high stopover density was more strongly associated with migration volume the following night rather than the preceding night. The processed mean scalar wind speeds with temporal resolutions as fine as 10-minute intervals near turbine height showed that August is the month with the weakest winds while December is the month, which typically has the strongest winds. The ice data suggests that shallow western basin of Lake Erie has higher ice cover duration many times exceeding 90 days during some winters.

17 WIND ENERGY↗

Particle-image Velocimetry Sensitivity through Automatic Differentiation

A commonly encountered difficulty in modern experimental fluid dynamics experiments is the need for full-flow visualization and measurement; traditional techniques such as hot-wire probes and laser Doppler velocimetry cannot provide such information without difficult scanning procedures. Particle-image velocimetry (PIV) is a measurement technique which can overcome this spatial limitation to provide real-time three-dimensional fluid velocity vectors for an entire plane within the flow, as opposed to a single point that is provided by other measurement techniques. To accomplish this, a PIV system seeds the fluid flow with tracer particles; these must have nearly the same density as that of the fluid and a high drag coefficient to ensure that they follow the fluid flow without deviation. A laser is sent through optics to expand its beam into a sheet which illuminates the plane of interest in the fluid. The particles reflect more light than the fluid, and their image is recorded using a high speed camera. Due to limitations in bandwidth, video cannot be taken fast enough to be useful, and instead pairs of images are captured with a known separation in time. The illuminated particles will move between the two images taken by the camera, and their relative motion can be used to compute the fluid velocity at that point. When a single camera is used, only two-dimensional velocity vectors can be computed, but with the addition of another camera this can be improved to full three-dimensional vectors. The flow's time evolution can be measured by taking series of image pairs and processing each pair to get the velocity field at the time between each image pair. Only the two-dimensional version will be considered in this work. To compute the velocity vectors of the flow, the images are broken up into interrogation regions, usually in a uniform grid. In the most basic PIV algorithm, corresponding regions in the image pairs are cross-correlated to compute the most likely displacement from one image to the next. This is used since the cross correlation will achieve a maximum value when the two interrogation regions line up the best. The location of the peak in the correlation map and scaling information for the optical setup is used to compute the particle displacement from one image to the next. This process is similar to holding two semi-transparent images and moving one until they match as well as possible. Sub-pixel displacements between images can be obtained by fitting a Gaussian function to the correlation peak. Knowing the displacement vector for each set of interrogation regions in the image pair, it is a simple matter to divide by the time gap between the two images to compute the fluid velocity at the point in time and space. Since there are many interrogation regions in each image pair, the velocity field for a sheet is sampled according to the arrangement of the interrogation regions. The sensitivity of the final computed vector to the input images is difficult to compute, and has led most considerations of uncertainty in the PIV field to apply Monte Carlo methods and correlations to quantify the systematic and random errors which result from application of various PIV algorithms. In this work, automatic differentiation (AD) has been used to trace the dependence of velocity vectors computed using the correlation-peak PIV algorithm back to the input values given to an artificial image generator. These initial results form the foundation for continued work in the application of AD to algorithms in PIV. They demonstrate that AD methods can be successfully applied to PIV algorithms without undue complexity. The derivatives of each individual vector with regard to the coefficients of the velocity field used to create the artificial images can be computed and the statistics of them considered collectively, or correlated between each other to quantify the effectiveness of a PIV algorithm. With the initial implementation complete, adding extra algorithms is relatively simple and considerable data can be generated for analysis. It is hoped that application of these methods can lead to improved PIV algorithms which demonstrate better sensitivity to small variations in particle displacement. While these methods can be used as an analysis tool for refinement of new algorithms, computation of actual resultant vector sensitivities to input images' individual pixels is the next logical step of this work. Such information will be analyzed with experimental data to produce a better understanding of how output vectors' uncertainty relates to the actual image pairs, and lead to automatic bad vector removal. (authors)

42 ENGINEERING↗

HarDWR - Raw Water Rights Records

For a detailed description of the database of which this record is only one part, please see the HarDWR meta-record. In order to hold a water right in the western United States, an entity, (e.g., an individual, corporation, municipality, sovereign government, or non-profit) must register a physical document with the state's water regulatory agency. State water agencies each maintain their own database containing all registered water right documents within the state, along with relevant metadata such as the point of diversion and place of use of the water. All western U.S. states have digitized their individual water rights databases, along with the geospatial data describing the spatial units where water rights are managed. Each state maintains and provides their own water rights data in accordance with individual state regulations and standards. We collected water rights databases from 11 western United States states either by downloading them from publicly accessible web portals, or by contacting state water management representatives; detailed descriptions of where and when the data was collected is provided in the README.txt, as well as Lisk et al.(in review). This collection of data are those raw water rights. Each state formats their data differently, meaning that file types, field availability, and names vary from state to state. Note, the data provided here reflects the state of the water rights databases at the time we collected the data; updates have likely occurred in many states. Some pieces of information are common among all states. These are: priority date, volume or flow of water allowed by the right, stated water use of the right, and some means of identifying the geography and source of the water pertaining to the right - typically the coordinates of the Point of Diversion (PoD) of a waterbody or well. Arizona regulates water in a different way than the other 10 states. Outside of some relatively small critical agricultural areas called Active Management Areas (AMAs), Arizona does not maintain any water rights. However, the state does require registration of surface and groundwater pumping devices, which includes disclosing the mechanical specifics of the devices. We used these records as a proxy for water rights. Each state, and their respective water right authorities, have made their water right records available for non-commercial reference uses. In addition, the states make no guarantees as to the completeness, accuracy, or timeliness of their respective databases, let alone the modifications which we, the authors of this paper, have made to the collected records. None of the states should be held liable for using this data outside of its intended use. In addition, the following states have requested specifically worded disclaimers to be included with their data. Colorado: "The data made available here has been modified for use from its original source, which is the State of Colorado. THE STATE OF COLORADO MAKES NO REPRESENTATIONS OR WARRANTY AS TO THE COMPLETENESS, ACCURACY, TIMELINESS, OR CONTENT OF ANY DATA MADE AVAILABLE THROUGH THIS SITE. THE STATE OF COLORADO EXPRESSLY DISCLAIMS ALL WARRANTIES, WHETHER EXPRESS OR IMPLIED, INCLUDING ANY IMPLIED WARRANTIES OF MERCHANTABILITY, OR FITNESS FOR A PARTICULAR PURPOSE. The data is subject to change as modifications and updates are complete. It is understood that the information contained in the Web feed is being used at one's own risk." Montana: "The Montana State Library provides this product/service for informational purposes only. The Library did not produce it for, nor is it suitable for legal, engineering, or surveying purposes. Consumers of this information should review or consult the primary data and information sources to ascertain the viability of the information for their purposes. The Library provides these data in good faith but does not represent or warrant its accuracy, adequacy, or completeness. In no event shall the Library be liable for any incorrect results or analysis; any direct, indirect, special, or consequential damages to any party; or any lost profits arising out of or in connection with the use or the inability to use the data or the services provided. The Library makes these data and services available as a convenience to the public, and for no other purpose. The Library reserves the right to change or revise published data and/or services at any time." Oregon: "This product is for informational purposes and may not have been prepared for, or be suitable for legal, engineering, or surveying purposes. Users of this information should review or consult the primary data and information sources to ascertain the usability of the information." The available data is provided as a series of compressed files, which each containing the full data collected from each state. Some of the files have been renamed, to more easily know which state the data belongs to. The file renaming was also required as some files from different states had the same name. In other cases, the data for a state has been placed in a folder indicating which state it belongs to - as the state organized its data by selected subregions. Below is a brief description of the format of the collected data from each state. ArizonaRights_StatementOfClaimants: A folder containing a database of interconnected CSV files. The soc_erd.pdf file contains a visual flowchart of how the various files are connected, beginning with SOC_MAIN.csv in the center of the page. ArizonaRights_SurfaceWaterRightsData: A folder containing a database of a single Shapefile and 10 associated CSVs. SurfaceWater.pdf contains a visual flowchart of how the various files are connected, beginning with ADWR_SW_APPL_REGRY.csv. ArizonaRights_Well55Registry: A folder containing a database of a single Shapefile and 59 associated CSVs. Wells55.pdf contains a visual flowchart of how the various files are connected, beginning with WellRegistry.shp. CaliforniaRights_eWRIMS_directDatabase: A folder containing a collection of four "series" Microsoft Excel files, as either XLS or XLSX. The four "series": byCounty, byEntity (what type of legal entity holds the right), byUse (stated water use), and byWatershed, are various methods by which the California water rights are organized within the state's database. However, it was observed that by only collecting a single series, not all water rights were being provided. So, essentially, the majority of records within each "series" are copies of each other, with each "series" containing some unique records. ColoradoRights_NetAmounts: A folder containing 78 CSV files, with one file per Colorado Water District. IdahoRights_PointOfDiversion: A Shapefile containing the Points of Diversion for the entire state of Idaho. IdahoRights_PlaceOfUse: A Shapefile containing the Place of Use polygons for the entire state of Idaho. MontanaRights_WaterRights: A Geodatabase file containing the Points of Diversion and Places of Use for the entire state of Montana. The name of the Points of Diversion Feature Layer within the Geodatabase is "WRDIV", and the name of the Places of Use Feature Layer is "WRPOU". NevadaRights_POD_Sites: A Shapefile containing the Points of Diversion for the entire state of Nevada. NewMexicoRights_Points_of_Diversion: A Shapefile containing the Points of Diversion for the entire state of New Mexico. OregonRights_state_shp: A folder containing 36 Shapefiles and are split between "pod" (Point of Diversion) and "pou" (Place of Use) for each water management basin within Oregon. In other words, each basin has one "pod" file and one "pou" file. The "pod" files are point shapes, and the "pou" files are polygons. UtahRights_Points_of_Diversion: A Shapefile containing the Points of Diversion for the entire state of Utah. WashingtonRights_WaterDiversions_ECY_NHD: A Geodatabase file containing both the Points of Diversion for the entire state of Washington. The name of the Feature Layer within the Geodatabase is "WaterDiversions_ECY_NHD". WyomingRights: A folder containing four subdirectories, one for each Wyoming Water Division. Each Division directory includes a varying number of subdirectories for each Wyoming Water District. Each District folder contains two copies of the Point of Diversion records for that area, with one copying being in CSV and one copy in Microsoft Excel XLS format.

Lisk, Matthew↗

Evaluating the Effectiveness of an Ultrasonic Acoustic Deterrent in Reducing Bat Fatalities at Wind Energy Facilities

This project was designed to use thermal video cameras and fatality monitoring to evaluate the effectiveness of an ultrasonic acoustic deterrent (UAD) on bat activity and mortality, respectively. Our goals were to redesign the UAD device and installation infrastructure, determine the placement on wind turbines to optimize safety, compatibility and functionality, and to compare the mortality among the following conditions: Control (deterrents off and turbines feathered up to the manufacturer’s cut-in speed of 3.5 m/s), Deterrent (deterrents on and turbines feathered up to the manufacturer’s cut-in speed of 3.5 m/s), Curtailment (deterrents off and turbines feathered up to 5/ m/s), and combination (deterrents on and turbines feathered up to 5 m/s). The project was divided into a Feasibility Study and Comparative Study. The objectives for the Feasibility Study were to develop an installation strategy, redesign a previous iteration of a UAD to improve performance and weatherization, and test the effectiveness of the deterrents on bat activity. The Feasibility Study was intended to work out potential issues using a relatively small number of devices prior to manufacturing and installing numerous devices for a larger-scale comparative study. The division of the project into these separate studies was based on previous experience and the challenges of assessing the capabilities of an untested UAD. During the initial development and manufacturing of the UAD, NRG Systems decided to use a piezoelectric transducer rather than an electrostatic transducer, which was used by a previous device (i.e. Deaton UAD). NRG Systems conducted lab testing (i.e., IP67 or Ingress Protection) to ensure no water or dust ingress. In addition, shock/drop trials and variations in temperature exposure were conducted as part of the reliability testing. NRG Systems also developed a communications system to allow for continuous performance monitoring of the UADs. For the Feasibility Study, we installed 6 UADs on each of 2 Gamesa G90 2-MW wind turbines (Turbine 14 and 8) at the South Chestnut Wind Energy Facility, Pennsylvania and monitored activity under control and treatment (i.e. Deterrent) conditions using thermal video monitoring. There are two major sources of variation in bat activity (beyond the anticipated treatment effect): 1) environment around the turbines might inherently favor more activity at one than the other; 2) weather conditions on any given night or within season difference (e.g., migration later during the study period) might favor more activity on some nights than on other nights. Because we could only monitor two turbines on any night, we sought to control the potential influences of these two sources by alternating the turbine on which deterrents were activated each night. If there were no loss of data due to technical failures, this design would result in an equal number of deterrent and control nights at each turbine through the monitoring period, balancing the effects of both sources of variation. We compared the time bats spent and the number of events that occurred in overlapping cameras FOV as an indicator of risk, since 80% of the overlapping FOV of the cameras was in the RSA. Equipment failures, majority due to lightning, resulted in only 17 nights with useable data, with unbalanced treatment assignment within turbines and uneven distribution of treatment assignment throughout the observational period. This resulted in a confounding of treatment assignment and seasonal change. Deterrent treatment was measured at Turbine 14 on only 2 of the first 8 usable nights (spanning the period from 8/18-9/17), whereas 6 times on Turbine 8. From 9/18-927, deterrent was on at Turbine 14 on 6 of the remaining 10 nights, and 4 on Turbine 8. We recorded a total of 1,057 bats and observed a reduction in number of events and duration of events at the UAD-activated turbine when it was Turbine 14. When Turbine 8 had the UAD activated, we observed no difference in number or duration of events. Variation between turbines is not unusual and can cause issues when study designs have no true replication, i.e., multiple turbines per treatment. Within-turbine differences suggested a trend for reduced activity when UADs were activated, particularly for Turbine 14. These results may be caused by the overall higher bat activity at Turbine 8 and confounding of treatment assignment and seasonal trends. We mapped 58 bat events in 3-dimensional space (3D), 30 and 28 events during control and treatment conditions, respectively. We observed bats crossing the rotor plane under both control and treatment conditions and observed a total of 40 confirmed or near-collisions (i.e., target close to blade but no visual confirmation of a strike) out of a total of 1,491 medium and high confidence bat observations (880 control, 611 at treatment). Twice as many collisions/possible collisions were observed during control conditions. Given the challenges with the equipment and potential confounding of the data (i.e., different activity levels at the two wind turbines), we were unable to determine whether this initial turbine placement and orientation was optimal. Given no new information on how best to install the devices, we elected to use the same placement and orientation for the comparative study. For the Comparative Study, the objectives were to investigate the relative mortality rates among 4 treatments. We searched the area within 90 m of each turbine daily to recover the highest number of fresh fatalities possible. We were unable to detect a clear reduction in mortality from deterrents alone for any individual species. Surprisingly, mortality rate of the eastern red bat (Lasiurus borealis) was estimated to be 1.3–4.2 times as much when turbines were operating normally and UADs were on than when UADs were off. Reduction in mortality of all bat species combined due to curtailment of turbines was estimated to be between 0%–38%. This effect was nullified when, in addition to curtailment, UADs were on, with 95% confidence interval ranging from a 45% reduction to a 36% increase in mortality. This was likely due to the large proportion of eastern red bats in the total carcass population. Mortality of all low-frequency echolocating bats combined (i.e. hoary bat [L. cinereus], big brown bat [Eptesicus fuscus], silver-haired bat [Lasionycteris noctivagans]) relative to control was lower when curtailed (95% CI: 0%–74%), but the addition of UADs had no detectable effect (95%CI: 13%–79%). The combined treatment reduced mortality in silver-haired bats relative to control by 11%–99%, compared to curtailment (81% reduction–67% increase) or deterrent (82% reduction–67% increase) alone. Because silver-haired bats comprised a large proportion of low-frequency calling bats found during this study, a similar effect was seen for that group. The higher mortality observed for eastern red bats at UAD compared to control could have been caused by several factors, such as the effective range of the UAD, particularly at higher frequencies, behavior, positioning of the devices on the nacelle, or a combination of these. We used 3D thermal videography to compare control and UAD bat behavior from two turbines using a total of 203 3D bat-tracks across 34 nights. We recorded a similar number of bat-tracks between treatment groups, with 51% and 49% for control and UAD, respectively. Due to potential differences in bat behavior around spinning vs stationary turbine blades, we examined UAD effectiveness separately for non-operating turbines (feathered below cut-in speed of 3.5 m/s) and operating (normal operation above wind speed of 3.5 m/s). We found a higher proportion of bat-tracks at operating turbines (82%) compared to non-operating turbines (18%), although this does not account for overall time turbines were operating versus not. At non-operating turbines the UAD appears to be effective at reducing the amount of time, flight length, and number of passes through the rotor plane, compared to control. In addition, we found bats approached turbines similarly between control and treatment turbines, with 61% of control and 63% of UAD bat-tracks originating leeward of the hub. In contrast, at operating turbines, we saw little change in bat behavior in response to UADs. For example, we found an increase in the average duration of bat-tracks between non-operating and operating turbines for UAD but at control turbines we found average duration decreased once turbines became operational. Both control and UAD had a high proportion of bat-tracks that crossed the rotor plane (i.e., collision risk) originate from the windward side when turbines were operational 65% to 92%, respectively. Given that the UAD devices closest to the blades were orientated parallel to the blades, its possible bats were not exposed to the signal until they were close to the turbine blades, as suggested by the slightly higher mean duration within 5 meters of the blades for UAD turbines. Future research should consider concentrating UAD intensity on the areas of risk (i.e. blades) with enough buffer to allow bats to react to the sound before entering the rotor-swept area (RSA). In addition, investigating the potential of installing UAD units windward of the turbine blades (e.g. a hub-mounted UAD), particularly since even under control conditions, a relatively high proportion (65%) of crosses through the blade plane originated windward. 3D thermal videography provided valuable information on future testing strategies (e.g. device placement, UAD orientation) to improve UAD effectiveness when bats are at risk (i.e., operating wind turbines). Across the entire project we experienced issues with the operation and communication with the UADs. Most of the issues occurred during the Feasibility Study and were resolved prior to the Comparability Study. Additional challenges surfaced during the Comparability Study but were remedied immediately and are thought to have little impact on the results. We had logistical constraints at the project that limited our ability use traditional methods in our camera calibration. Several calibrations showed inaccurate scales, which may have been related to inadequate spatial coverage of “points” in the camera calibration volume. Because we were using actual video recordings of bats at a wind turbine, the behavior of bats could have concentrated “points” in specific areas of the turbine (i.e. leeward of nacelle), and limited “points” in other areas, resulting in camera calibration issues. We have plans to address these inconsistencies and improving the software and related methodologies by early 2020.

17 WIND ENERGY↗

Small Scale Field Test Demonstrating CO 2 Sequestration In Arbuckle Saline Aquifer And By CO 2 -Eor At Wellington Field, Sumner County, Kansas

The objectives of this project are to understand the processes that occur when a maximum of 70,000 metric tonnes of CO 2 are injected into two different formations to evaluate the response in different lithofacies and depositional environments. The evaluation will be accomplished through the use of both in situ and indirect MVA (monitoring, verification, and accounting) technologies. The project will optimize for carbon storage accounting for 99% of the CO 2 using lab and field testing and comprehensive characterization and modeling techniques. Site characterization and CO 2 injection should demonstrate state-of-the-art MVA tools and techniques to monitor and visualize the injected CO 2 plume and to refine geomodels developed using nearly continuous core, exhaustive wireline logs, and well tests and a multi-component 3-D seismic survey. Reservoir simulation studies will map the injected CO 2 plume and estimate tonnage of CO 2 stored in solution, as residual gas, and by mineralization and integrate MVA results and reservoir models shall be used to evaluate CO 2 leakage. A rapid-response mitigation plan was developed to minimize CO 2 leakage and provide a comprehensive risk management strategy. The CO 2 was intended to be supplied from a reliable facility and have an adequate delivery and quality of CO 2 . However, several unforeseen circumstances complicated this plan: (1) the initially negotiated CO 2 supply facility went offline and contracts associated with CO 2 supply had to be renegotiated, (2) a UIC Class VI permit proved to be difficult to obtain due to the experimental nature of the project. Both subjects are detailed in separate deliverables attached to this report. The CO 2 enhanced oil recovery (EOR) and geologic storage in Mississippian carbonate reservoir was sucessully deployed. Approximately 20,000 metric tons of CO 2 was injected in the upper part of the Mississippian reservoir to verify CO 2 EOR viability in carbonate reservoirs and evaluate a potential of transitioning to geologic CO 2 storage through EOR. A total of 1,101 truckloads, 19,803 metric tons—an average of 120 tonnes per day—were delivered over the course of injection that lasted from January 9 to June 21, 2016. After cessation of CO 2 injection, the KGS 2-32 well was converted to water injector and continues to operate. CO 2 EOR progression in the field was monitored weekly with fluid level, temperature, and production recording and formation fluid composition sampling. It is important to note that normally, CO 2 EOR pilots are less efficient than commercial operations due to lack of directional and precise well control, lack of surface facilities for CO 2 recycling, and other factors. As a result of this pilot CO 2 injection, the observed incremental average oil production increase was ~68% with only ~18% of injected CO 2 produced back. Decline curve analysis forecasts of additional cumulative oil produced were 32.44M STB to the end of 2027. Wellington Mississippian pilot efficiency by the end of forecast calculations is 11 MCF per barrel of produced oil. Using 32M STB oil production and $1,964,063 cost of CO 2 , CO 2 EOR cost per barrel of oil production is ~$60. Simple but robust monitoring technologies proved to be very efficient in detecting and locating CO 2 . High CO 2 reservoir retentions with low yields within an actively producing field could help to estimate real-world risks of CO 2 geological storage for future projects. The Wellington Field CO 2 EOR was executed in a controlled environment with high efficiency. This case study proves that CO 2 EOR could be successfully applied in Kansas carbonate reservoirs if CO 2 sources and associated infrastructure are available. Recent developments in unconventional resources development in Mid-Continent USA and associated large volume disposal of backflow water and the resulting seismic activity have brought more focus and attention to the Arbuckle Group in southern Kansas. Despite the commercial interest, limited essential information about reservoir properties and structural elements has impeded the management and regulation of disposal, an issue brought to the forefront by recent seismicity in and near areas of large volumes and rates of brine disposal. The Kansas Geological Survey (KGS) collected, compiled, and analyzed available data, including well logs, core data, step rate tests, drill stem tests, 2-D and 3-D seismic data, water level measurements, and others types of data. Several exploratory wells were drilled and core was collected and modern suites of logs were analyzed. Reservoir properties were populated into several site-specific geological models. The geological models illustrate the highly heterogeneous nature of the Arbuckle Group. Vertical and horizontal variability results in several distinct hydro-stratigraphic units that are the result of both depositional and diagenetic processes. During the course of this project, it has been demonstrated that advanced seismic interpretation methods can be used successfully for characterization of the Mississippian reservoir and Arbuckle saline aquifer. Analysis of post-stack 3-D seismic data at the Mississippian reservoir showed the response of a gradational velocity transition. Pre-stack gather analysis showed that porosity zones of the Mississippian and Arbuckle reservoirs exhibit characteristic amplitude versus offset (AVO) response. Simultaneous AVO inversion estimated P- and S-impedances. The 3-D survey gather azimuthal anisotropy analysis (AVAZ) provided information about the fault and fracture network and showed good agreement to the regional stress field and well data. Mississippian reservoir porosity and fracture predictions agreed well with the observed mobility of injected CO 2 in KGS well 2-32. Fluid substitution modeling predicted acoustic impedance reduction in the Mississippian carbonate reservoir introduced by the presence of CO 2 . Seismicity in the United States midcontinent has increased by orders of magnitude over the past decade. Spatiotemporal correlations of seismicity to wastewater injection operations have suggested that injection-related pore fluid pressure increases are inducing the earthquakes. In this investigation, we examine earthquake occurrence in southern Kansas and northern Oklahoma and its relation to the change in pore pressure. The main source of data comes from the Wellington Array in the Wellington oil field, in Sumner County, Kansas, which has monitored for earthquakes in central Sumner County, Kansas, since early 2015. The seismometer array was established to monitor CO 2 injection operations at Wellington Field. Although no seismicity was detected in association with the spring 2016 Mississippian CO 2 injection, the array has recorded more than 2,500 earthquakes in the region and is providing valuable understanding to induced seismicity. A catalog of earthquakes was built from this data and was analyzed for spatial and temporal changes, stress information, and anisotropy information. The region of seismic concern has been shown to be expanding through use of the Wellington earthquake catalog, which has revealed a northward progression of earthquake activity reaching the metropolitan area of Wichita. The stress orientation was also calculated from this earthquake catalog through focal mechanism inversion. The calculated stress orientation was confirmed through comparison to other stress measurements from well data and previous earthquake studies in the region. With this knowledge of the stress orientation, the anisotropy in the basement could be understood. This allowed for the anisotropy measurements to be correlated to pore pressure increases. The increase in pore pressure was monitored through time-lapse shear-wave anisotropy analysis. Since the onset of the observation period in 2010, the orientation of the fast shear wave has rotated 90°, indicating a change associated with critical pore pressure build up. The time delay between fast and slow shear wave arrivals has increased, indicating a corresponding increase in anisotropy induced by pore pressure rise. In-situ near-basement fluid pressure measurements corroborate the continuous pore pressure increase revealed by the shear-wave anisotropy analysis over the earthquake monitoring period. This research is the first to identify a change in pore fluid pressure in the basement using seismological data and it was recently published in the AAAS journal Science Advances (Nolte et al., 2017). The shear-wave splitting analysis is a novel application of the technique, which can be used in other regions to identify an increase in pore pressure. This increasing pore fluid pressure has become more regionally extensive as earthquakes are occurring in southern Kansas, where they previously were absent. These monitoring techniques and analyses provide new insight into mitigating induced seismicity’s impact on society.

02 PETROLEUM↗