Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “horizontal”

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 1,675 records · Page 93
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Dynamically Downscaled Hourly Future Weather Data with 12-km Resolution Covering Most of North America

This is an hourly future weather dataset for energy modeling applications. The dataset is primarily based on the output of a regional climate model (RCM), i.e., the Weather Research and Forecasting (WRF) model version 3.3.1. The WRF simulations are driven by the output of a general circulation model (GCM), i.e., the Community Climate System Model version 4 (CCSM4). This dataset is in the EPW format, which can be read or translated by more than 25 building energy modeling programs (e.g., EnergyPlus, ESP-r, and IESVE), energy system modeling programs (e.g., System Advisor Model (SAM)), indoor air quality analysis programs (e.g., CONTAM), and hygrothermal analysis programs (e.g., WUFI). It contains 13 weather variables, which are the Dry-Bulb Temperature, Dew Point Temperature, Relative Humidity, Atmospheric Pressure, Horizontal Infrared Radiation Intensity from Sky, Global Horizontal Irradiation, Direct Normal Irradiation, Diffuse Horizontal Irradiation, Wind Speed, Wind Direction, Sky Cover, Albedo, and Liquid Precipitation Depth. The weather data is created for two emissions scenarios: RCP4.5 and RCP8.5 and spans two 10-year time slices in the future: 2045 - 2054 and 2085 - 2094. It offers a spatial resolution of 12 km by 12 km with extensive coverage across most of North America. Due to the enormous size of the entire dataset, in the first stage of its distribution, we provide 20 years of future weather data for the centroid of each Public Use Microdata Area (PUMA), excluding Hawaii. PUMAs are non-overlapping, statistical geographic areas that partition each state or equivalent entity into geographic areas containing no fewer than 100,000 people each. The 2,378 PUMAs as a whole cover the entirety of the U.S. The weather data can be utilized alongside the large-scale energy analysis tools, ResStock and ComStock, developed by National Renewable Energy Laboratory, whose smallest resolution is at the PUMA scale. The data for RCP4.5 is still being processed and will be published soon.

Array↗
DOE OSTI2023

Urbanization may enhance tornado potential: A single case report

Tornadoes pose a risk of catastrophic economic loss and casualty in the United States. Modification of land use by urbanization alters the meteorological conditions that may impact tornado formation and intensification processes. Here we explored the simulated impact of Kansas City urbanization on the tornado potential of a supercell storm. In this studied case, we found that urbanization might enhance tornado potential by a) strengthening the low-level streamwise vorticity in the storm inflow region, thus forming stronger rotating updrafts; and b) intensifying near-surface horizontal vorticity near the boundary of the forward-flank cold pool which increases the ingestion, tilting, and stretching of streamwise horizontal vorticity into vertical vorticity. The former results from the stronger east-to-west pressure perturbation gradient due to the faster, stronger outflow boundary, and the latter is mainly a result of stronger cold fronts and a better alignment of storm-relative inflow with the horizontal vorticity vector. We emphasize that our conclusions only represent one possibility of how urbanization would affect tornado potential, and a more systematic examination is needed to achieve a more general conclusion.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

The O{sub 2} A-Band in the Fluxes and Polarization of Starlight Reflected by Earth-Like Exoplanets

Earth-like, potentially habitable exoplanets are prime targets in the search for extraterrestrial life. Information about their atmospheres and surfaces can be derived by analyzing the light of the parent star reflected by the planet. We investigate the influence of the surface albedo A {sub s}, the optical thickness b {sub cloud}, the altitude of water clouds, and the mixing ratio of biosignature O{sub 2} on the strength of the O{sub 2} A-band (around 760 nm) in the flux and polarization spectra of starlight reflected by Earth-like exoplanets. Our computations for horizontally homogeneous planets show that small mixing ratios ( η < 0.4) will yield moderately deep bands in flux and moderate-to-small band strengths in polarization, and that clouds will usually decrease the band depth in flux and the band strength in polarization. However, cloud influence will be strongly dependent on properties such as optical thickness, top altitude, particle phase, coverage fraction, and horizontal distribution. Depending on the surface albedo and cloud properties, different O{sub 2} mixing ratios η can give similar absorption-band depths in flux and band strengths in polarization, especially if the clouds have moderate-to-high optical thicknesses. Measuring both the flux and the polarization is essential to reduce the degeneracies, although it will not solve them, especially not for horizontally inhomogeneous planets. Observations at a wide range of phase angles and with a high temporal resolution could help to derive cloud properties and, once those are known, the mixing ratio of O{sub 2} or any other absorbing gas.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Challenge Solutions for a Significant Water Cut Reduction History Match on Heavy Oil Polymer EOR

Milne Point Field is conducting the first-ever active polymer flood pilot program designed to enhance oil recovery (EOR) of heavy oils on the Alaska North Slope (ANS). The experimental field program consists of two pairs of horizontal injection and production wells deployed in an isolated fault block in the Schrader Bluff reservoir. The EOR method combines polymer flooding and low salinity water flooding to improve recovery. Prior average water cut reduction in polymer EOR fields ordinarily ranges between 10% to 15% for most light oil and heavy oil reservoirs. When comparing other field applications, the long-lasting (24 months) and much lower water cut poses a challenge for numerical history matching when applying typical simulation methods. This abstract discusses two practical new approaches used for history matching of the long-lasting ultralow water cut. These new applications include 1) multiple relative permeability curves with J-function incorporation with low endpoints of relative permeability of water (Krw). This approach integrated the permeability differential adjustment in the geometries of horizontal and vertical directions based on geological review. 2) A specific permeability channel strips along with the horizontal or vertical planar assumption. The current history matching results indicate success in achieving an excellent agreement with the field observed low water-cut.

Namie, S.↗
DOE OSTI2014

Quantitative Methods for Comparing Different Polyline Stream Network Models

Two techniques for exploring relative horizontal accuracy of complex linear spatial features are described and sample source code (pseudo code) is presented for this purpose. The first technique, relative sinuosity, is presented as a measure of the complexity or detail of a polyline network in comparison to a reference network. We term the second technique longitudinal root mean squared error (LRMSE) and present it as a means for quantitatively assessing the horizontal variance between two polyline data sets representing digitized (reference) and derived stream and river networks. Both relative sinuosity and LRMSE are shown to be suitable measures of horizontal stream network accuracy for assessing quality and variation in linear features. Both techniques have been used in two recent investigations involving extracting of hydrographic features from LiDAR elevation data. One confirmed that, with the greatly increased resolution of LiDAR data, smaller cell sizes yielded better stream network delineations, based on sinuosity and LRMSE, when using LiDAR-derived DEMs. The other demonstrated a new method of delineating stream channels directly from LiDAR point clouds, without the intermediate step of deriving a DEM, showing that the direct delineation from LiDAR point clouds yielded an excellent and much better match, as indicated by the LRMSE.

58 GEOSCIENCES↗
DOE OSTI2022

Behavior and mechanisms of Doppler wind lidar error in varying stability regimes

Abstract. Wind lidars are widespread and important tools in atmospheric observations. An intrinsic part of lidar measurement error is due to atmospheric variability in the remote-sensing scan volume. This study describes and quantifies the distribution of measurement error due to turbulence in varying atmospheric stability. While the lidar error model is general, we demonstrate the approach using large ensembles of virtual WindCube V2 lidar performing a profiling Doppler-beam-swinging scan in quasi-stationary large-eddy simulations (LESs) of convective and stable boundary layers. Error trends vary with the stability regime, time averaging of results, and observation height. A systematic analysis of the observation error explains dominant mechanisms and supports the findings of the empirical results. Treating the error under a random variable framework allows for informed predictions about the effect of different configurations or conditions on lidar performance. Convective conditions are most prone to large errors (up to 1.5 m s−1 in 1 Hz wind speed in strong convection), driven by the large vertical velocity variances in convective conditions and the high elevation angle of the scanning beams (62∘). Range-gate weighting induces a negative bias into the horizontal wind speeds near the surface shear layer (−0.2 m s−1 in the stable test case). Errors in the horizontal wind speed and direction computed from the wind components are sensitive to the background wind speed but have negligible dependence on the relative orientation of the instrument. Especially during low winds and in the presence of large errors in the horizontal velocity estimates, the reported wind speed is subject to a systematic positive bias (up to 0.4 m s−1 in 1 Hz measurements in strong convection). Vector time-averaged measurements can improve the behavior of the error distributions (reducing the 10 min wind speed error standard deviation to <0.3 m s−1 and the bias to <0.1 m s−1 in strong convection) with a predictable effectiveness related to the number of decorrelated samples in the time window. Hybrid schemes weighting the 10 min scalar- and vector-averaged lidar measurements are shown to be effective at reducing the wind speed biases compared to cup measurements in most of the simulated conditions, with time averages longer than 10 min recommended for best use in some unstable conditions. The approach in decomposing the error mechanisms with the help of the LES flow field could be extended to more complex measurement scenarios and scans.

17 WIND ENERGY↗
DOE OSTI2022

Sizing ice hydrometeor populations using the dual-wavelength radar ratio

Dual-wavelength (3.2 and 0.32 cm, i.e., X- and W-radar bands) radar ratio (DWR) measurements in ice clouds and precipitation using Canada's National Research Council Institute for Aerospace Research airborne radar are compared to closely collocated particle microphysical in situ sampling data in order to develop relations between DWR and characteristic hydrometeor size. This study uses the radar and in situ data sets collected during the In-Cloud ICing and Large-drop Experiment (ICICLE) campaign in midlatitude frontal clouds. Since atmospheric particle scattering at X band is predominantly in the Rayleigh regime and the W-band frequency is the highest frequency usually used for hydrometeor remote sensing, the X–W-band combination provides a relatively strong dual-wavelength reflectivity difference. This study considers radar and in situ measurements conducted in relatively homogeneous cloud and precipitation conditions. Measurements show that under these conditions, the difference between the X-band radar reflectivities observed with vertical and horizontal pointing of the radar beam are generally small and often negligible. However, W-band reflectivities at vertical beam pointing are, on average, larger than those for horizontal beam pointing by about 4 dB, which is a non-Rayleigh scattering effect from preferentially oriented non-spherical particles. A horizontal radar beam DWR–mean volume particle size relation, D v , provides robust estimates of this characteristic size for populations of particles with different habits. Uncertainties of D v retrievals using DWR are around 0.6 mm when D v is greater than approximately 1 mm. Size estimates using vertical radar beam DWRs have larger uncertainties due to smaller dual-wavelength signals and stronger influences of hydrometeor habits and orientations at this geometry of beam pointing. Mean relations among different characteristic sizes, which describe the entire particle size distribution (PSD), such as D v , and other sizes used in various applications (e.g., the mean, effective, and median sizes) are derived, so the results of this study can be used for estimating different PSD characteristic sizes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

The CU 2-D-MAX-DOAS instrument – Part 1: Retrieval of 3-D distributions of NO 2 and azimuth-dependent OVOC ratios

Here, we present an innovative instrument telescope and describe a retrieval method to probe three-dimensional (3-D) distributions of atmospheric trace gases that are relevant to air pollution and tropospheric chemistry. The University of Colorado (CU) two-dimensional (2-D) multi-axis differential optical absorption spectroscopy (CU 2-D-MAX-DOAS) instrument measures nitrogen dioxide (NO 2 ), formaldehyde (HCHO), glyoxal (CHOCHO), oxygen dimer (O 2 –O 2 , or O 4 ), and water vapor (H 2 O); nitrous acid (HONO), bromine monoxide (BrO), and iodine monoxide (IO) are among other gases that can in principle be measured. Information about aerosols is derived through coupling with a radiative transfer model (RTM). The 2-D telescope has three modes of operation: mode 1 measures solar scattered photons from any pair of elevation angle (–20° < EA < +90° or zenith; zero is to the horizon) and azimuth angle (–180° < AA < +180°; zero being north); mode 2 measures any set of azimuth angles (AAs) at constant elevation angle (EA) (almucantar scans); and mode 3 tracks the direct solar beam via a separate view port. Vertical profiles of trace gases are measured and used to estimate mixing layer height (MLH). Horizontal distributions are then derived using MLH and parameterization of RTM (Sinreich et al., 2013). NO 2 is evaluated at different wavelengths (350, 450, and 560 nm), exploiting the fact that the effective path length varies systematically with wavelength. The area probed is constrained by O 4 observations at nearby wavelengths and has a diurnal mean effective radius of 7.0 to 25 km around the instrument location; i.e., up to 1960 km 2 can be sampled with high time resolution. The instrument was deployed as part of the Multi-Axis DOAS Comparison campaign for Aerosols and Trace gases (MAD-CAT) in Mainz, Germany, from 7 June to 6 July 2013. We present first measurements (modes 1 and 2 only) and describe a four-step retrieval to derive (a) boundary layer vertical profiles and MLH of NO 2 ; (b) near-surface horizontal distributions of NO 2 ; (c) range-resolved NO 2 horizontal distribution measurements using an "onion-peeling" approach; and (d) the ratios HCHO to NO 2 ( R FN ), CHOCHO to NO 2 ( R GN ), and CHOCHO to HCHO ( R GF ) at 14 pre-set azimuth angles distributed over a 360° view. Three-dimensional distribution measurements with 2-D-MAX-DOAS provide an innovative, regional perspective of trace gases as well as their spatial and temporal concentration gradients, and they maximize information to compare near-surface observations with atmospheric models and satellites.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

Implicit–explicit (IMEX) Runge–Kutta methods for non-hydrostatic atmospheric models

Abstract. The efficient simulation of non-hydrostatic atmospheric dynamics requires time integration methods capable of overcoming the explicit stability constraints on time step size arising from acoustic waves. In this work, we investigate various implicit–explicit (IMEX) additive Runge–Kutta (ARK) methods for evolving acoustic waves implicitly to enable larger time step sizes in a global non-hydrostatic atmospheric model. The IMEX formulations considered include horizontally explicit – vertically implicit (HEVI) approaches as well as splittings that treat some horizontal dynamics implicitly. In each case, the impact of solving nonlinear systems in each implicit ARK stage in a linearly implicit fashion is also explored. The accuracy and efficiency of the IMEX splittings, ARK methods, and solver options are evaluated on a gravity wave and baroclinic wave test case. HEVI splittings that treat some vertical dynamics explicitly do not show a benefit in solution quality or run time over the most implicit HEVI formulation. While splittings that implicitly evolve some horizontal dynamics increase the maximum stable step size of a method, the gains are insufficient to overcome the additional cost of solving a globally coupled system. Solving implicit stage systems in a linearly implicit manner limits the solver cost but this is offset by a reduction in step size to achieve the desired accuracy for some methods. Overall, the third-order ARS343 and ARK324 methods performed the best, followed by the second-order ARS232 and ARK232 methods.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Impact of numerical choices on water conservation in the E3SM Atmosphere Model version 1 (EAMv1)

Abstract. The conservation of total water is an important numerical feature for global Earth system models. Even small conservation problems in the water budget can lead to systematic errors in century-long simulations. This study quantifies and reduces various sources of water conservation error in the atmosphere component of the Energy Exascale Earth System Model. Several sources of water conservation error have been identified during the development of the version 1 (V1) model. The largest errors result from the numerical coupling between the resolved dynamics and the parameterized sub-grid physics. A hybrid coupling using different methods for fluid dynamics and tracer transport provides a reduction of water conservation error by a factor of 50 at 1∘ horizontal resolution as well as consistent improvements at other resolutions. The second largest error source is the use of an overly simplified relationship between the surface moisture flux and latent heat flux at the interface between the host model and the turbulence parameterization. This error can be prevented by applying the same (correct) relationship throughout the entire model. Two additional types of conservation error that result from correcting the surface moisture flux and clipping negative water concentrations can be avoided by using mass-conserving fixers. With all four error sources addressed, the water conservation error in the V1 model becomes negligible and insensitive to the horizontal resolution. The associated changes in the long-term statistics of the main atmospheric features are small. A sensitivity analysis is carried out to show that the magnitudes of the conservation errors in early V1 versions decrease strongly with temporal resolution but increase with horizontal resolution. The increased vertical resolution in V1 results in a very thin model layer at the Earth's surface, which amplifies the conservation error associated with the surface moisture flux correction. We note that for some of the identified error sources, the proposed fixers are remedies rather than solutions to the problems at their roots. Future improvements in time integration would be beneficial for V1.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

DCMIP2016: the tropical cyclone test case

This paper describes and analyzes the Reed–Jablonowski (RJ) tropical cyclone (TC) test case used in the 2016 Dynamical Core Model Intercomparison Project (DCMIP2016). This intermediate-complexity test case analyzes the evolution of a weak vortex into a TC in an idealized tropical environment. Reference solutions from nine general circulation models (GCMs) with identical simplified physics parameterization packages that participated in DCMIP2016 are analyzed in this study at 50 km horizontal grid spacing, with five of these models also providing solutions at 25 km grid spacing. Evolution of minimum surface pressure (MSP) and maximum 1 km azimuthally averaged wind speed (MWS), the wind–pressure relationship, radial profiles of wind speed and surface pressure, and wind composites are presented for all participating GCMs at both horizontal grid spacings. While all TCs undergo a similar evolution process, some reach significantly higher intensities than others, ultimately impacting their horizontal and vertical structures. TCs simulated at 25 km grid spacings retain these differences but reach higher intensities and are more compact than their 50 km counterparts. These results indicate that dynamical core choice is an essential factor in GCM development, and future work should be conducted to explore how specific differences within the dynamical core affect TC behavior in GCMs.

Willson, Justin L.↗
DOE OSTI2017

Impact of numerical choices on water conservation in the E3SM Atmosphere Model Version 1 (EAM V1)

The conservation of total water is an important numerical feature for global Earth system models. Even small conservation problems in the water budget can lead to systematic errors in century-long simulations for sea level rise projection. This study quantifies and reduces various sources of water conservation error in the atmosphere component of the Energy Exascale Earth System Model. Several sources of water conservation error have been identified during the development of the version 1 (V1) model. The largest errors result from the numerical coupling between the resolved dynamics and the parameterized sub-grid physics. A hybrid coupling using different methods for fluid dynamics and tracer transport provides a reduction of water conservation error by a factor of 50 at 1&deg; horizontal resolution as well as consistent improvements at other resolutions. The second largest error source is the use of an overly simplified relationship between the surface moisture flux and latent heat flux at the interface between the host model and the turbulence parameterization. This error can be prevented by applying the same (correct) relationship throughout the entire model. Two additional types of conservation error that result from correcting the surface moisture flux and clipping negative water concentrations can be avoided by using mass-conserving fixers. With all four error sources addressed, the water conservation error in the V1 model is negligible and insensitive to the horizontal resolution. The associated changes in the long-term statistics of the main atmospheric features are small. A sensitivity analysis is carried out to show that the magnitudes of the conservation errors decrease strongly with temporal resolution but increase with horizontal resolution. The increased vertical resolution in the new model results in a very thin model layer at the Earth’s surface, which amplifies the conservation error associated with the surface moisture flux correction. We note that for some of the identified error sources, the proposed fixers are remedies rather than solutions to the problems at their roots. Future improvements in time integration would be beneficial for this model.

Zhang, Kai↗
DOE OSTI2016

High Resolution Model Intercomparison Project (HighResMIP v1.0) for CMIP6

Robust projections and predictions of climate variability and change, particularly at regional scales, rely on the driving processes being represented with fidelity in model simulations. The role of enhanced horizontal resolution in improved process representation in all components of the climate system is of growing interest, particularly as some recent simulations suggest both the possibility of significant changes in large-scale aspects of circulation as well as improvements in small-scale processes and extremes. However, such high-resolution global simulations at climate timescales, with resolutions of at least 50 km in the atmosphere and 0.25° in the ocean, have been performed at relatively few research centres and generally without overall coordination, primarily due to their computational cost. Assessing the robustness of the response of simulated climate to model resolution requires a large multi-model ensemble using a coordinated set of experiments. The Coupled Model Intercomparison Project 6 (CMIP6) is the ideal framework within which to conduct such a study, due to the strong link to models being developed for the CMIP DECK experiments and other model intercomparison projects (MIPs). Increases in high-performance computing (HPC) resources, as well as the revised experimental design for CMIP6, now enable a detailed investigation of the impact of increased resolution up to synoptic weather scales on the simulated mean climate and its variability. The High Resolution Model Intercomparison Project (HighResMIP) presented in this paper applies, for the first time, a multi-model approach to the systematic investigation of the impact of horizontal resolution. A coordinated set of experiments has been designed to assess both a standard and an enhanced horizontal-resolution simulation in the atmosphere and ocean. The set of HighResMIP experiments is divided into three tiers consisting of atmosphere-only and coupled runs and spanning the period 1950–2050, with the possibility of extending to 2100, together with some additional targeted experiments. This paper describes the experimental set-up of HighResMIP, the analysis plan, the connection with the other CMIP6 endorsed MIPs, as well as the DECK and CMIP6 historical simulations. Lastly, HighResMIP thereby focuses on one of the CMIP6 broad questions, “what are the origins and consequences of systematic model biases?”, but we also discuss how it addresses the World Climate Research Program (WCRP) grand challenges.

58 GEOSCIENCES↗
DOE OSTI2023

Investigating the physical mechanisms that modify wind plant blockage in stable boundary layers

Abstract. Wind plants slow down the approaching wind, a phenomenon known as blockage. Wind plant blockage undermines turbine performance for front-row turbines and potentially for turbines deeper into the array. We use large-eddy simulations to characterize blockage upstream of a finite-size wind plant in flat terrain for different atmospheric stability conditions and investigate the physical mechanisms modifying the flow upstream of the turbines. To examine the influence of atmospheric stability, we compare simulations of two stably stratified boundary layers using the Weather Research and Forecasting model in large-eddy simulation mode, representing wind turbines using the generalized actuator disk approach. For a wind plant, a faster cooling rate at the surface, which produces stronger stably stratified flow in the boundary layer, amplifies blockage. As a novelty, we investigate the physical mechanisms amplifying blockage by evaluating the different terms in the momentum conservation equation within the turbine rotor layer. The velocity deceleration upstream of a wind plant is caused by an adverse pressure gradient and momentum advection out of the turbine rotor layer. The cumulative deceleration of the flow upstream of the front-row turbines instigates vertical motions. The horizontal flow is diverted vertically, reducing momentum availability in the turbine rotor layer. Although the adverse pressure gradient upstream of the wind plant remains unchanged with atmospheric stability, vertical advection of horizontal momentum is amplified in the more strongly stable boundary layer, mainly by larger shear of the horizontal velocity, thus increasing the blockage effect.

17 WIND ENERGY↗
DOE OSTI2012

SACRADV-3D3C

SACR-ADV-3D3C remaps the outputs of SACRCORR for cross-wind range-height indicator (CW-RHI) scans to a Cartesian grid and reports reflectivity CFAD and best estimate domain averaged cloud fraction. CW-RHI scans consist of multiple consecutive horizon-to-horizon scans (elevations ~10-170°) (Kollias et al. 2014). The CW-RHI scan strategy is typically continuously operated for 25-55 min at the ARM sites. The major inputs of this VAP are generated by the SACRCORR VAP. They include: Reflectivity corrected for gas attenuation, Unfolded mean Doppler velocity, spectral width and linear depolarization ratio, all of which are free of noise, insects and second trip echo and all of which are in native radar polar coordinates. An additional input is the horizontal wind speed from the nearest (in time) sounding. These 2-D grid slices are used to estimate a reflectivity CFAD. This CFAD contains information for each height about the distribution of reflectivity observed for the duration of the scan strategy (e.g. 19 minutes). A profile of maximum observable domain is reported for user guidance. For completeness, the time coordinate is converted to distance using the sounding horizontal wind speed to provide an estimate of x-y-z cloud fraction. The 2-D grids are also used as an input to a 3-D gridding technique. This technique interpolates or extrapolates the 2-D grids in time. The final output is a single NetCDF file containing all aforementioned corrected radar moments remapped on a 3-D Cartesian grid, the SACR reflectivity CFAD, a profile of best estimate cloud fraction, a profile of maximum observable x-domain size (xmax), a profile time to horizontal distance estimate and a profile of minimum observable reflectivity (dBZmin).

3-D Cartesian grid, reflectivity CFAD↗
DOE OSTI2023

Produced Fluid Induced Mineralogy and Elemental Alterations of Caney Shale, Southern Oklahoma

ABSTRACT This study involves batch reactor experiments and subsequent analyses of samples from Caney Shale in the Ardmore Basin of South-Central Oklahoma. Samples include mainly rock cores and cuttings recovered from two wells respectively drilled vertically through and horizontally across the Caney Shale. Mineralogical compositions are obtained by X-Ray Diffraction (XRD) measurements whilst microstructure and elemental distribution are acquired by Scanning Electron Microscopy/Energy Dispersive Spectroscopy (SEM/EDS) respectively. Batch experiments are then conducted using selected rock samples and produced fluid from the Caney Formation. Deionized water is also reacted with some samples to serve as standard. Experiments are conducted at 95°C and ambient pressure for 7 and 30 days to assess the geochemical rock-fluid interactions. Results show rock mineralogical compositions are predominantly quartz, feldspar, carbonates, and clay with minor pyrite. Post-experimental mineralogical changes observed in samples include increased amorphous entities especially within the clay portions of XRD plots and dissolution of feldspar and carbonate minerals and formation of new mineral phases, mostly clays and salts. These are corroborated by EDS elemental analyses which show decreased elemental compositions. The implications of reactions mentioned above include but not limited to, scale formation, clay fines migration and shale softening all of which pose significant permeability impairment on formation over time. INTRODUCTION Shale reservoirs account for a large share of unconventional reservoirs in the world (Lyu et al., 2015). However, ultra-low permeability and high clay compositions pose significant challenges when producing from these reservoirs (Dawuda and Srinivasan, 2022, 2023). Producing from these reservoirs therefore requires horizontal drilling and hydraulic fracturing technologies which have proven their efficacy in generating substantial permeability in reservoirs to ensure production (Fujian et al., 2019; Liu et al., 2018). Even after expensive horizontal drilling and hydraulic fracturing, geochemical reactions between engineered fluids and formation leads to fracture constriction and adversely impact petrophysical properties (permeability and porosity) of the reservoir. These technologies are therefore under constant development and improvement in various aspects to ensure fine tuning for specific reservoirs. Under present conditions, much of the hydrocarbon reserves in unconventional shale reservoirs are left unproduced due to rapid decline in permeability following resumption to production after hydraulic fracturing. It is therefore essential to understand the range of geochemical reactions that cause rapid depletion of permeability after hydraulic fracturing and apply these to each shale reservoir to ensure substantial recovery rates.

Awejori, G. A.↗
DOE OSTI2025

Baltimore Social-Environmental Collaborative (BSEC) Doppler Lidar & Derived Products

This repository contains all processed Doppler‐lidar outputs from the PSU lidar deployed for the Baltimore Social‐Environmental Collaborative (BSEC) project. Vertical Stare Scans (fixed‐beam, vertical profiling): 1 Hz backscatter intensity (m⁻¹ sr⁻¹), signal‐to‐noise ratio (unitless), and Doppler vertical‐velocity (m s⁻¹) on ~30 m range gates, stored as CF-compliant NetCDF. Wind Profiles (horizontal‐wind retrieval): daily NetCDF outputs of retrieved horizontal wind speed (m s⁻¹) and direction (degrees), computed from the angled‐scan returns. Profile Statistics (summary statistics on the vertical velocity): 15 min windows (default) of mean, variance, skewness, kurtosis, high-frequency variance, etc., as a function of height; saved as CF-compliant NetCDF files. Boundary Layer Height (BLH) (fuzzy-logic output): 15 min BLH estimates (m), with lower/upper fuzzy bounds (m) and a quality flag (0–4) indicating data status (e.g., no data, good, below range, ran out of signal, cloud-topped). Cloud Base Height (Haar-gradient detection): 15 min estimates of cloud-base height (m) with a cloud-detection quality flag (0–3: none, low, moderate, high). All five product streams are organized by year and date under their own top-level folders (01_Vertical_Stare_Scans/ through 05_Cloud_Height/). Each folder contains a data_ /YYYY/ subdirectory with daily CF-compliant NetCDF outputs (96 windows per day at 15 min intervals). Global attributes in each file include creation history, version (2.0.0), institution, and source. Instrument & MeasurementsThe PSU Doppler Lidar samples aerosol backscatter (m⁻¹ sr⁻¹), signal-to-noise ratio, and radial velocity at ~1 Hz. Vertical stare scans point the beam straight up; after collecting angled scans through multiple elevation angles, the "Wind Profiles" product contains the fully retrieved horizontal wind speed and direction. Data were collected continuously at ~30 m range resolution, with a typical height ceiling of ~12 km. How to Use Open any NetCDF with Python's xarray, MATLAB, or similar CF-compliant tools. Stare scans and angled-scan retrievals (Wind Profiles) are CF-compliant daily NetCDF files. Profile-Statistics, BLH, and Cloud Height files are daily 15 min summaries (96 time steps per file). Inspect the included variables (e.g., vertical_velocity_variance, wind_speed, BLH, cloud_base_height) for your analyses. Use the quality flags (BLH_flag, cloud_flag) to filter out poor-quality retrievals. For more information or questions about processing methods, please contact:Nicholas E. Prince ⟨nec5299@psu.edu⟩Penn State Department of Meteorology & Atmospheric Science

Air Quality↗
DOE OSTI2026

Baltimore Social-Environmental Collaborative (BSEC) Doppler Lidar & Derived Products

This repository contains all processed Doppler‐lidar outputs from the PSU lidar deployed for the Baltimore Social‐Environmental Collaborative (BSEC) project. Vertical Stare Scans (fixed‐beam, vertical profiling): 1 Hz backscatter intensity (m⁻¹ sr⁻¹), signal‐to‐noise ratio (unitless), and Doppler vertical‐velocity (m s⁻¹) on 30 m range gates (and 3 m range gates), stored as CF-compliant NetCDF. Wind Profiles (horizontal‐wind retrieval): daily NetCDF outputs of retrieved horizontal wind speed (m s⁻¹) and direction (degrees), computed from the angled‐scan returns. Profile Statistics (summary statistics on the vertical velocity): 15 min windows (default) of mean, variance, skewness, kurtosis, high-frequency variance, etc., as a function of height; saved as CF-compliant NetCDF files. Boundary Layer Height (BLH) (fuzzy-logic output): 15 min BLH estimates (m), with lower/upper fuzzy bounds (m) and a quality flag (0–4) indicating data status (e.g., no data, good, ran out of signal, below range, cloud-topped). Cloud Base Height (Haar-gradient detection): 10 min estimates of cloud-base height (m). All five product streams are organized by year and date under their own top-level folders (01_Vertical_Stare_Scans/ through 05_Cloud_Height/). Each folder contains a data_ /YYYY/ subdirectory with daily CF-compliant NetCDF outputs (96 windows per day at 15 min intervals). Global attributes in each file include creation history, version (3.0.0), institution, and source. Instrument & MeasurementsThe PSU Doppler Lidar samples aerosol backscatter (m⁻¹ sr⁻¹), signal-to-noise ratio, and radial velocity at ~1 Hz. Vertical stare scans point the beam straight up; after collecting angled scans through multiple elevation angles, the "Wind Profiles" product contains the fully retrieved horizontal wind speed and direction. Data were collected continuously at ~30 m range resolution (and 3 m for the year of 2025), with a typical height ceiling of ~12 km. How to Use Open any NetCDF with Python's xarray, MATLAB, or similar CF-compliant tools. Stare scans and angled-scan retrievals (Wind Profiles) are CF-compliant daily NetCDF files. Profile-Statistics, BLH, and Cloud Height files are daily 15 min (10 min for Cloud Heights) summaries (96 time steps per file). Inspect the included variables (e.g., vertical_velocity_variance, wind_speed, BLH, cloud_base_height) for your analyses. Use the quality flags (BLH_flag, cloud_flag) to filter out poor-quality retrievals. For more information or questions about processing methods, please contact:Nicholas E. Prince ⟨nec5299@psu.edu⟩Penn State Department of Meteorology & Atmospheric Science

Air Quality↗