SEARCH · Engineering Papers

Results for “point clouds”

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.

SourceResearch connectionFrom yearThrough year

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

At least 1,693 records · Page 94

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Performance Evaluation of Intelligent Solar Control Software Through Hardware-in-the-Loop (CRADA Final Report)

Recent research has highlighted the potential for solar to act as a zero-marginal-cost and zero-emission flexibility resource on the bulk power system when operated with advanced control systems. To increase the performance of these systems, leading technologies, including machine learning (ML) and hierarchical inverter set point allocation, have been developed by Latimer Controls, Inc. to estimate the headroom of large PV plants for grid operation and control; however, these technologies lack comprehensive validation under real-world application scenarios. Latimer Controls, Inc. received two voucher awards for research at a national laboratory from the Department of Energy American Made Solar Prize Round 6. The National Renewable Energy Laboratory (NREL) was selected to collaborate with Latimer staff to conduct a performance evaluation of Latimer PV control software. The NREL team will develop a hardware-in-the-loop (HIL) testbed to perform testing and validation of the Latimer PV control technology in a de-risked yet realistic testbed environment. Latimer and NREL worked together to analyze the test data, draw conclusions from the results, and disseminate the resulting scientific findings. In this CRADA work, we propose to test and validate the real-world application of the Latimer Control solution in an HIL environment. We evaluate the performance of different flexible solar technologies in responding to automatic generation control signals in a closed-loop fashion. In particular, a data-driven potential high limit (PHL) estimation is developed for large solar plants to accurately estimate their headroom so that they have fast and short-time regulation and control capability to participate in grid services and respond to grid signals in real time (e.g., AGC). This PHL estimation algorithm is embedded in a hardware power plant controller (PPC) and tested with an IEEE-39 bus system model developed in RTDS. To account for the varying cloud conditions and diverse inverter dispatches, we developed a 135-MW PV plant with detailed modeling of 27 individual PV modules and inverters using RTDS. The real-world communications used in such big plants, such as ModBus TCP/IP for inverter level and DNP3 for plant level, were developed to emulate the real-world applications in big PV plants. The ML-based PHL estimation method is tested under nine separate weather scenarios against the ‘reference-control’ solution, hereafter referred to as the baseline solution. The baseline method reserves a subset of inverters (reference group) to operate at their PHL at all times and dispatches only the remaining inverters (control group) at curtailed levels to fulfill the flexibility need. Despite being successfully piloted by NREL in California in 2017 and Chile in 2020, there exist two gaps in the state of the art to fully unlock the flexibility of PV plants: a. There is a trade-off between the PHL estimation accuracy and the flexibility range. b. There lacks granularity in the PHL estimation to capture the variation across inverters. The Latimer solution seeks to address these gaps by applying machine learning methods to improve PHL estimation accuracy while accounting for variability at every inverter. Performance metrics were taken from the 2023 Georgia Power CARES utility-scale RFP. The results demonstrate that the ML-based approach outperforms the traditional baseline method in PHL estimation accuracy for 7 of 9 scenarios. The average PHL error across the nine scenarios was 7.40% for the ML-based method, 2.06% less than the 9.46% PHL error average across scenarios that was exhibited by the baseline method. Additionally, the PHL error was below 5% for at least 95% of the testing interval for 3 of 9 tested intervals with the ML approach, whereas it did not achieve this metric for any of the baseline tests. Overall, simulation results indicate the superior performance of an ML-based approach compared to the conventional baseline reference-control approach, showcasing its potential to support grid stability and operational efficiency. This laboratory HIL testing using real PPC, representative power system simulation models in real-time with detailed PV plant and inverter models, and real-world communication protocols gives us confidence that this machine learning based PHL estimation algorithm works well in the hardware PPC and therefore de-risks future field commissioning. The end goal of this project is to advance grid technology to address the grid operation challenges brought by solar plant’s variability and uncertainties in power generation.

14 SOLAR ENERGY↗

NASA NTRS2016

Recent Variability Observations of Solar System Giant Planets: Fresh Context for Understanding Exoplanet and Brown Dwarf Weather

Over the past several years a number of high cadence photometric observations of solar system giant planets have been acquired by various platforms. Such observations are of interest as they provide points of comparison to the already expansive set of brown dwarf variability observations and the small, but growing, set of exoplanet variability observations. By measuring how rapidly the integrated light from solar system giant planets can evolve, variability observations of substellar objects that are unlikely to ever be resolved can be placed in a fuller context. Examples of brown dwarf variability observations include extensive work from the ground (e.g., Radigen et al. 2014), Spitzer (e.g., Metchev et al. 2015), Kepler (Gizis et al. 2015), and HST (Yang et al. 2015).Variability has been measured on the planetary mass companion to the brown dwarf 2MASS 1207b (Zhou et al. 2016) and further searches are planned in thermal emission for the known directly imaged planets with ground based telescopes (Apai et al. 2016) and in reflected light with future space based telescopes. Recent solar system variability observations include Kepler monitoring of Neptune (Simon et al. 2016) and Uranus, Spitzer observations of Neptune (Stauffer et al. 2016), and Cassini observations of Jupiter (West et al. in prep). The Cassini observations are of particular interest as they measured the variability of Jupiter at a phase angle of approximately 60 deg, comparable to the viewing geometry expected for space based direct imaging of cool extrasolar Jupiters in reflected light. These solar system analog observations capture many of the characteristics seen in brown dwarf variability, including large amplitudes and rapid light curve evolution on timescales as short as a few rotation periods. Simon et al. (2016) attribute such variations at Neptune to a combination of large scale, stable cloud structures along with smaller, more rapidly varying, cloud patches. The observed brown dwarf and exoplanet variability may well arise from comparable cloud structures. In my presentation I will compare and contrast the nature of the variability observed for the various solar system and other substelar objects and present a wish list for future observations.

giant planets↗

DOE OSTI2024

EPCAPE-PT-LANL Measurements: Single Particle Soot Photometer

Coastal cities offer a unique environment for studying aerosol-cloud interactions and the effects of urban emissions on cloud properties. As part of the Eastern Pacific Cloud Aerosol Precipitation Experiment (EPCAPE), the Partitioning Thrust by Los Alamos National Laboratory (EPCAPE-PT-LANL) was conducted. Our campaign focused on measuring the optical and chemical properties of aerosols and their interactions within marine stratocumulus clouds in La Jolla, California. EPCAPE-PT-LANL enhances the primary goals of EPCAPE through innovative observations of vapor-phase transitions between aerosols and cloud droplets, the impact of black carbon on aerosol-cloud dynamics, and the effects of cloud processing on aerosol optical properties. Instrument: Single Particle Soot Photometer (Droplet Measurements Technology) Data Notes: Contact us if you want additional data products from this instrument. Reported data is the black carbon (rBC) number and mass concentration with diagnostic flags. The SP2 measures incandescence from particles that is induced with a 1064 Nd-YAG laser. Particles that absorb the laser energy and then emit radiation is assumed to contain black carbon. The amount of intensity of radiation is related to the mass of absorbing material in the particle. Here, we calibrate the incandescence intensity to size selected Regal Black (Cabot) nebulized from solution. Data provided for the EPCAPE campaign used the combined broadband high-gain + low-gain channels (BHBL). Each channel has a lower threshold of detection equivalent to 2-s of the respective channel noise. Thresholding has been applied to select between the high-gain and low-gain channels. The detection limit is 80-540 nm (Deq) or 0.36- 55 fg. It is assumed that rBC is the only aerosol type that is in significant concentration in the sampled atmosphere that absorbs the laser energy. Particle data are integrated over 10 second windows to calculate a rBC number and mass concentration. QC/QA: • Diagnostic flags that impacted measured concentration: - Sample, Sheath, and Purge Flow Rates: Despite observing fluctuations in all flows, the rBC detection and mass quantification are generally observed to be stable, although large changes in sample flow rate did impact detection efficiency. A flag was implemented for sample flow deviations >12 cm3/min from the set point averaged over 30 seconds - Laser Power: Detection efficiency decreases with laser power and deviations in laser power also affect mass quantification. Flag for laser power is set for deviations in laser current from the set value >5 mA. - Primary Detector Threshold: Primary thresholding is the signal value which determines whether a “particle” is recorded. Thresholding errors occur when the threshold value is too HIGH and real particles are ignored. • Several periods without data: - 16-18 Nov: Ultra Zero Air generator failure, flows deviated significantly from set points. rBC # conc recovered but questionable. - 20-21 Nov: Offline for calibrations for several hours each day. - 25-27 Nov: Data is missing. - 30 Nov: Power outage 3 Dec: SP2 hard-drive full. - 3 Dec: Power Outage Header: - BHBL_NumbConc[#/cc]: Refers to the number concentration of refractory black carbon (rBC) measured from combined broadband high-gain and low-gain channels, expressed in particles per cubic centimeter. - BHBL_BCMass_Conc[ug/m3]: Refers to the mass concentration of refractory black carbon (rBC) measured from combined broadband high-gain and low-gain channels, expressed in micrograms per cubic meter. - NoData_Flag[bool]: A boolean flag that indicates whether no data was recorded during a measurement. - NoBC_Flag[bool]: A boolean flag indicating whether no black carbon particles were recorded during the measurement. - Laser_Flag[bool]: A boolean flag indicating deviations in laser current during the measurement. - SampleFlow_Flag[bool]: A boolean flag indicating deviations from the set sample flow rate during the measurement. - PrimThresh_Flag[bool]: A boolean flag indicating deviations from the set primary threshold during the measurement, potentially ignoring real particles. - Manual_Flag[bool]: A boolean flag indicating manual intervention or adjustments during the measurement. - CVI_Flag[bool]: A boolean flag indicating whether the Counterflow Virtual Impactor (CVI) was active or inactive during the measurement.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

High temporal resolution estimates of Arctic snowfall rates emphasizing gauge and radar-based retrievals from the MOSAiC expedition

This article presents the results of snowfall rate and accumulation estimates from a vertically pointing 35-GHz radar and other sensors deployed during the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition. The radar-based retrievals are the most consistent in terms of data availability and are largely immune to blowing snow. The total liquid-equivalent accumulation during the snow accumulation season is around 110 mm, with more abundant precipitation during spring months. About half of the total accumulation came from weak snowfall with rates less than approximately 0.2 mmh–1. The total snowfall estimates from a Vaisala optical sensor aboard the icebreaker are similar to those from radar retrievals, though their daily and monthly accumulations and instantaneous rates varied significantly. Compared to radar retrievals and the icebreaker optical sensor data, measurements from an identical optical sensor at an ice camp are biased high. Blowing snow effects, in part, explain differences. Weighing gauge measurements significantly overestimate snowfall during February–April 2020 as compared to other sensors and are not well suited for estimating instantaneous snowfall rates. The icebreaker optical disdrometer estimates of snowfall rates are, on average, relatively little biased compared to radar retrievals when raw particle counts are available and appropriate snowflake mass-size relations are used. These counts, however, are not available during periods that produced more than a third of the total snowfall. While there are uncertainties in the radar-based retrievals due to the choice of reflectivity-snowfall rate relations, the major error contributor is the uncertainty in the radar absolute calibration. The MOSAiC radar calibration is evaluated using comparisons with other radars and liquid water cloud–drizzle processes observed during summer. Overall, this study describes a consistent, radar-based snowfall rate product for MOSAiC that provides significant insight into Central Arctic snowfall and can be used for many other purposes.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1996

Stratospheric Heterogeneous Chemistry and Microphysics: Model Development, Validation and Applications

The objectives of this project are to: define the chemical and physical processes leading to stratospheric ozone change that involve polar stratospheric clouds (PSCS) and the reactions occurring on the surfaces of PSC particles; study the formation processes, and the physical and chemical properties of PSCS, that are relevant to atmospheric chemistry and to the interpretation of field measurements taken during polar stratosphere missions; develop quantitative models describing PSC microphysics and heterogeneous chemical processes; assimilate laboratory and field data into these models; and calculate the extent of chemical processing on PSCs and the impact of specific microphysical processes on polar composition and ozone depletion. During the course of the project, a new coupled microphysics/physical-chemistry/ photochemistry model for stratospheric sulfate aerosols and nitric acid and ice PSCs was developed and applied to analyze data collected during NASA's Arctic Airborne Stratospheric Expedition-II (AASE-II) and other missions. In this model, detailed treatments of multicomponent sulfate aerosol physical chemistry, sulfate aerosol microphysics, polar stratospheric cloud microphysics, PSC ice surface chemistry, as well as homogeneous gas-phase chemistry were included for the first time. In recent studies focusing on AASE measurements, the PSC model was used to analyze specific measurements from an aircraft deployment of an aerosol impactor, FSSP, and NO(y) detector. The calculated results are in excellent agreement with observations for particle volumes as well as NO(y) concentrations, thus confirming the importance of supercooled sulfate/nitrate droplets in PSC formation. The same model has been applied to perform a statistical study of PSC properties in the Northern Hemisphere using several hundred high-latitude air parcel trajectories obtained from Goddard. The rates of ozone depletion along trajectories with different meteorological histories are presently being systematically evaluated to identify the principal relationships between ozone loss and aerosol state. Under this project, we formulated a detailed quantitative model that predicts the multicomponent composition of sulfate aerosols under stratospheric conditions, including sulfuric, nitric, hydrochloric, hydrofluoric and hydrobromic acids. This work defined for the first time the behavior of liquid ternary-system type-1b PSCS. The model also allows the compositions and reactivities of sulfate aerosols to be calculated over the entire range of environmental conditions encountered in the stratosphere (and has been incorporated into a trajectory/microphysics model-see above). Important conclusions that derived from this work over the last few years include the following: the HNO3 content of liquid-state aerosols dominate PSCs below about 195 K; the freezing of nitric acid ice from sulfate aerosol solutions is likely to occur within a few degrees K of the water vapor frost point; the uptake and reactions of HCl in liquid aerosols is a critical component of PSC heterogeneous chemistry. In a related application of this work, the inefficiency of chlorine injection into the stratosphere during major volcanic eruptions was explained on the basis of nucleation of sulfuric acid aerosols in rising volcanic plumes leading to the formation of supercooled water droplets on these aerosols, which efficiently scavenges HCl via precipitation.

Turco, Richard P.↗

DOE OSTI2025

Toward Drilling the Perfect Geothermal Well: An International Research Coordination Network for Geothermal Drilling Optimization Supported by Deep Machine Learning and Cloud Based Data Aggregation

The EDGE project, supported by the U.S. Department of Energy Geothermal Technologies Office under award DE-EE0008793, established a data-driven framework for improving the efficiency, cost-effectiveness, and reliability of geothermal well drilling. The project focused on developing scalable data infrastructure, advanced machine learning and probabilistic models, and integrated analytics tools to support continuous drilling optimization. A central objective was to reduce geothermal drilling costs by up to seventy percent while minimizing the risk of well failure through predictive diagnostics and adaptive planning. Over the project period, a comprehensive data repository was designed and deployed, incorporating records from over one hundred geothermal wells across varied geological settings. This repository supported both structured and unstructured data and adhered to FAIR data principles, enabling provenance tracking, quality control, and standardized metadata. The project introduced automated ingestion pipelines and a cloud-hosted platform that facilitated access to raw, processed, and derived datasets. This infrastructure served as the foundation for model development and analysis. Machine learning workflows were developed to predict key drilling metrics including rate of penetration, non-productive time, and total drilling costs. Self-organizing maps and dimensionality reduction methods were used to uncover operational patterns and outliers, while supervised learning algorithms such as random forests and deep neural networks were applied to forecast performance outcomes. The models were validated on heterogeneous datasets from both U.S. and Icelandic fields, demonstrating variable but significant predictive accuracy. The results indicated that finer temporal resolution, inclusion of lithological data, and consistency in operational annotations could substantially improve model performance. The project also implemented process mining techniques to reconstruct state-transition models from drilling event logs. These models enabled the identification of deviations from optimal workflows and provided insights into recurring failure modes. Analysis of non-productive time highlighted the impact of equipment failures, geological challenges, and human factors, offering opportunities for targeted mitigation strategies. The EDGE Dashboard was developed as a web-based expert system integrating data visualization, model outputs, and user-driven queries. It provided an accessible interface for operators to explore historical data, evaluate predicted outcomes, and compare drilling scenarios. Initial feedback from project partners suggested that the dashboard could serve as a foundation for more advanced advisory and optimization tools. Overall, the EDGE project demonstrated the feasibility and value of applying modern data science techniques to geothermal drilling. It delivered a set of interoperable tools and models that can support more efficient, lower-risk well development. The findings point toward a viable path for transitioning from advisory analytics to semi-autonomous drilling systems, contingent on continued collaboration, expanded datasets, and field validation. The project results have immediate relevance for drilling operations, data management practices, and future geothermal R&D efforts aimed at achieving reliable, cost-competitive geothermal energy at scale.

15 GEOTHERMAL ENERGY↗

DOE OSTI2026

Latency Analysis of the Nexus Digital Twin Framework

Real-time digital catalogs are increasingly relied upon to track metadata and connect disparate data sources for cloud-based data integration efforts. One such tool, Deeplynx Nexus is supporting real-time digital twin efforts through event-driven data integration and time-series queries. Nexus’s usefulness for these applications depends critically on how quickly individual records can be uploaded and downloaded, since delays directly affect the responsiveness of any system built on top of it. However, the actual latency a user should expect from Nexus has not been systematically measured before, particularly for the small, frequent transactions typical of live sensor feeds. Here we show that single-record round-trip latency is 61.1 ms on a local Nexus instance and 391.7 ms on the hosted production infrastructure, a roughly 6.4x difference driven primarily by fixed per-request overhead rather than data volume. This overhead dominates at small scale: comparing single-record and ten-record trials suggests approximately 56 ms of each single-record request is fixed connection and authentication cost rather than data-transfer time, meaning batching even a handful of records is substantially more efficient than transmitting them individually. At large batch sizes, this pattern reverses for uploads, which converge to near parity between local and hosted environments by 25,000-50,000 records, while download latency remains persistently 5.7-6.4x slower on hosted infrastructure even at scale. These results suggest that Nexus deployments intended for real-time digital twin applications should prioritize record batching over single-record transactions, and that download-path optimization on hosted infrastructure offers the largest remaining opportunity to reduce latency at scale. We anticipate these baseline measurements will serve as a reference point for future digital twin projects evaluating whether Nexus’s latency profile meets their real-time requirements, and as a benchmark for tracking the effect of future infrastructure or API changes.

99 - GENERAL AND MISCELLANEOUS↗

NASA NTRS1994

The Origins of Polycyclic Aromatic Hydrocarbons: Are They Everywhere?

During the past 15 years considerable progress in observational techniques has been achieved in the middle-infrared region (5000-500 per centimeter, 2-20 micron), the region where most diagnostic molecular vibrations occur. Spectra of many different astronomical infrared sources, some deeply embedded in dark molecular clouds and others at their edges, are now available. These spectra provide a powerful probe, not only for the identification of interstellar molecules in both the gas and solid phases, but also of the physical and chemical conditions which prevail in these two very different domains. The two lectures will focus on the evidence that polycyclic aromatic hydrocarbons (PAHs) are an important, ubiquitous and abundant interstellar species. PAHs are. extremely stable species which can range in size from a few angstroms across to several hundred angstroms (PAHs are also the building blocks of amorphous carbon particles). This identification rests on the suggestive agreement between the laboratory spectra of PAHs with a set of IR emission bands which emanate from many different sources where ultraviolet starlight impinges on a "dusty" region. The picture is that individual PAHs are first pumped into highly vibrationally excited states and relax by fluorescence at their fundamental vibrational frequencies. That PAHs are a ubiquitous interstellar component has serious ramifications in other spectral regions as well, including the strong extinction in the ultraviolet, and the classic visible diffuse interstellar bands discovered more than 50 years ago (but unexplained to this day) The first part of the course will focus on the interpretation of astronomical spectra. The second lecture will start by showing how recent laboratory data on PAHs taken under realistic interstellar conditions has con borated the PAH hypothesis and led to great insight into the conditions in the PAH containing regions. This lecture will end by reviewing the ever-increasing evidence for interstellar PAHs in meteorites and interplanetary dust particles. This in conjunction with the recent suggestion that PAHs are abundant in Jupiter's atmosphere will make the point that "PAHs are Everywhere".

Allamandola, L. J.↗

DOE OSTI2022

Three-Dimensional Convective–Stratiform Echo-Type Classification and Convectivity Retrieval from Radar Reflectivity

The Echo Classification from COnvectivity (ECCO) algorithm identifies convective and stratiform types of radar echo in three dimensions. It is based on the calculation of reflectivity texture—a combination of the intensity and the heterogeneity of the radar echoes on each horizontal plane in a 3D Cartesian volume. Reflectivity texture is translated into convectivity, which is designed to be a quantitative measure of the convective nature of each 3D radar grid point. It ranges from 0 (100% stratiform) to 1 (100% convective). By thresholding convectivity, a more traditional qualitative categorization is obtained, which classifies radar echoes as convective, mixed, or stratiform. In contrast to previous algorithms, these echo-type classifications are provided on the full 3D grid of the reflectivity field. The vertically resolved classifications, in combination with temperature data, allow for subclassifications into shallow, mid-, deep, and elevated convective features, and low, mid-, and high stratiform regions—again in three dimensions. The algorithm was validated using datasets collected over the U.S. Great Plains during the PECAN field campaign. An analysis of lightning counts shows ~90% of lightning occurring in regions classified as convective by ECCO. A statistical comparison of ECCO echo types with the well-established GPM radar precipitation-type categories show 84% (88%) of GPM stratiform (convective) echo being classified as stratiform (convective) or mixed by ECCO. ECCO was applied to radar grids for the continental United States, the United Arab Emirates, Australia, and Europe, illustrating its robustness and adaptability to different radar grid characteristics and climatic regions.

Radars/Radar observations↗

DOE OSTI2023

Characterization of Orography-Influenced Riming and Secondary Ice Production and Their Effects on Precipitation Rates Using Radar Polarimetry and Doppler Spectra (CORSIPP-SAIL)

The Characterization of Orography-Influenced Riming and Secondary Ice Production and Their Effects on Precipitation Rates Using Radar Polarimetry and Doppler Spectra (CORSIPP) project was conducted to help improve the understanding of precipitation formation in orographically influenced terrain. Special focus is put on the two processes of riming and secondary ice production and their external drivers. Two instruments, a polarimetric W-band simultaneous transmission simultaneous reception (STSR) Doppler cloud radar manufactured by Radiometer Physics GmbH (RPG, instrument type RPG-FMCW-94-DP), from now on named LIMRAD94, and the video in situ snowfall sensor (VISSS), were deployed at the U. S. Department of Energy Atmospheric Radiation Measurement (ARM) user facility’s Surface Atmosphere Integrated Field Laboratory (SAIL) site in Gothic, Colorado between November 2022 and June 2023 during the second SAIL winter. Note that the exact dates of data availability differ between the instruments. Both instruments arrived at Gothic on November 2, 2022, 09:20 local time. VISSS, described by Maahn et al., is equipped with two camera systems with telecentric lenses. The two cameras are at a 90° angle to each other. This configuration allows for size-independent measurements by capturing images of hydrometeors from two sides at a high frame rate of 250 Hz. With a minimum detection size of 200 μm, VISSS provides valuable insights into particle size, number, shape, complexity, and fall velocity. The VISSS was deployed on the grassland next to the ARM facility with the amazing help of the ARM employees on site. The setup started on November 2, 2022, and was finished on November 5, 2022, without major problems. VISSS measurements were started on November 6, 2022. LIMRAD94 was installed on a scaffolding platform near Orehouse (Gothic) on November 9, 2022, with the great help of RMBL staff. LIMRAD94 was mounted on a cold temperature scanner prototype. After a short test of the setup on November 9, 2022, the digital control of the scanner elevation stopped working for (at that time) unknown reasons. All attempts to resolve the problem failed. This malfunction made it impossible to operate LIMRAD94 in scanning mode. The scanner was then manually moved to zenith pointing mode and between November 10 and November 15, 2022, vertical observations for polarimetric calibration were performed. On November 15, 2022, after the polarimetric calibration was applied, the scanner was manually moved to 40° elevation with azimuthal view towards the Ka-band ARM Zenith Radar (KAZR) and measurements were continued at constant elevation. Investigation of the scanner malfunction on February 6, 2023, by Benn Schmatz revealed a disengagement between the cogwheel of the elevation motor and the cogwheel moving the scanner in elevation. This mechanical issue was temporarily solved by re-engaging the cogwheels. This made the scanner operational again for about four weeks, until mechanical force disengaged the cogwheels again on March 15, 2023. This repeated scanner failure remained undetected for about three weeks until April 8, 2023; during this time the scanner was stuck at 72° elevation. However, the radar software continued to produce data files falsely indicating that the scanner was still operational. After the scanner failure was noticed, scanning was stopped again, and we returned to constant elevation measurements. On May 15, 2023, the cogwheels of the elevation motor were secured with additional screws sent by the manufacturer. At some point in May, the cogwheels of the azimuth motor were also disengaged by mechanical force, which still allowed for range height indicator (RHI) but no plan position indicator (PPI) scans in the last weeks of the campaign. The azimuth motor was repaired in Germany after the end of the campaign. Throughout the campaign, RMBL and ARM staff kept the radar and the VISSS free of snow.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2008

The EXIST Mission Concept Study

EXIST is a mission designed to find and study black holes (BHs) over a wide range of environments and masses, including: 1) BHs accreting from binary companions or dense molecular clouds throughout our Galaxy and the Local Group, 2) supermassive black holes (SMBHs) lying dormant in galaxies that reveal their existence by disrupting passing stars, and 3) SMBHs that are hidden from our view at lower energies due to obscuration by the gas that they accrete. 4) the birth of stellar mass BHs which is accompanied by long cosmic gamma-ray bursts (GRBs) which are seen several times a day and may be associated with the earliest stars to form in the Universe. EXIST will provide an order of magnitude increase in sensitivity and angular resolution as well as greater spectral resolution and bandwidth compared with earlier hard X-ray survey telescopes. With an onboard optical-infra red (IR) telescope, EXIST will measure the spectra and redshifts of GRBs and their utility as cosmological probes of the highest z universe and epoch of reionization. The mission would retain its primary goal of being the Black Hole Finder Probe in the Beyond Einstein Program. However, the new design for EXIST proposed to be studied here represents a significant advance from its previous incarnation as presented to BEPAC. The mission is now less than half the total mass, would be launched on the smallest EELV available (Atlas V-401) for a Medium Class mission, and most importantly includes a two-telescope complement that is ideally suited for the study of both obscured and very distant BHs. EXIST retains its very wide field hard X-ray imaging High Energy Telescope (HET) as the primary instrument, now with improved angular and spectral resolution, and in a more compact payload that allows occasional rapid slews for immediate optical/IR imaging and spectra of GRBs and AGN as well as enhanced hard X-ray spectra and timing with pointed observations. The mission would conduct a 2 year full sky survey in scanning mode, interrupted for several orbits per day by GRB follow-ups, followed by a combined pointing-scanning mission phase for optical/IR spectroscopy and redshifts for the large AGN sample found in the survey as well as GRBs and LSST transients. A Team of university, NASA, and industry investigators will conduct the study to determine the full sensitivity and capabilities of this new configuration for EXIST. It will build on the extensive studies of the prior design for the mission and the HET and will incorporate the optical/IR telescope (hereafter IRT) now fully developed by our ITT partner for the NextView Commercial Remote Sensing mission (early 2008 launch) with a focal plane to be developed at GSFC based in part on JWST/NIRSPEC designs. No new technology is needed for either the IRT or HET instruments. The study will pay close attention to full mission cost and present a design for the Decadal Survey Workshop to ensure this even more capable EXIST mission is once again part of the next Decadal Survey.

Fishman, Gerald J.↗

NASA NTRS1998

Multiwavelength Observations of a Dramatic High-Energy Flare in the Blazar 3C 279

The blazar 3C 279, one of the brightest identified extragalactic objects in the gamma-ray sky, underwent a large (factor of approximately 10 in amplitude) flare in gamma-rays toward the end of a 3 week pointing by Compton Gamma Ray Observatory (CGRO), in 1996 January-February. The flare peak represents the highest gamma-ray intensity ever recorded for this object. During the high state, extremely rapid gamma-ray variability was seen, including an increase of a factor of 2.6 in approximately 8 hr, which strengthens the case for relativistic beaming. Coordinated multifrequency observations were carried out with Rossi X-Ray Timing Explorer (RXTE), Advanced Satellite for Cosmology and Astrophysics (ASCA; or, Astro-D), Roentgen Satellite (ROSAT), and International Ultraviolet Explorer (IUE) and from many ground-based observatories, covering most accessible wavelengths. The well-sampled, simultaneous RXTE light curve shows an outburst of lower amplitude (factor of approximately = 3) well correlated with the gamma-ray flare without any lag larger than the temporal resolution of approximately 1 day. The optical-UV light curves, which are not well sampled during the high-energy flare, exhibit more modest variations (factor of approximately 2) and a lower degree of correlation. The flux at millimetric wavelengths was near a historical maximum during the gamma-ray flare peak, and there is a suggestion of a correlated decay. We present simultaneous spectral energy distributions of 3C 279 prior to and near to the flare peak. The gamma-rays vary by more than the square of the observed IR-optical flux change, which poses some problems for specific blazar emission models. The synchrotron self-Compton (SSC) model would require that the largest synchrotron variability occurred in the mostly unobserved submillimeter/far-infrared region. Alternatively, a large variation in the external photon field could occur over a time-scale of a few days. This occurs naturally in the "mirror" model wherein the flaring region in the jet photoionizes nearby broad emission line clouds, which, in turn, provide soft external photons that are Comptonized to gamma-ray energies.

Wehrle, A. E.↗

NASA NTRSDate not supplied

Imaging Science Requirements for a Uranus Flagship Mission

Introduction: Our presentation will summarize the requirements for UV/Visible/Near-IR remote-sensing imaging science instruments (including mapping spectrometers) envisioned for the Uranus Orbiter and Probe (UOP) mission recommended by the recent planetary decadal survey. The Uranus system offers many targets for imaging investigations. Imaging science investigations are crucial to answering fundamental questions about various physical processes that shape the Uranian system and drive its evolution. Each science objective addressed through imaging has requirements for observation and illumination geometries as well as coverage and resolution in the spatial, spectral and temporal dimensions. By combining the requirements for the diverse science targets in the Uranian system, we present a set of overarching imaging science requirements for UOP. Atmospheric Science: The study of Uranus’ atmosphere is crucial in understanding the planet's energy balance and interior heat flux. Imaging observations map how the atmosphere scatters incident solar radiation. Imaging is key to characterizing planetary-scale energy/material transport though global atmospheric dynamics and local meteorology. Hyperspectral mapping combined with radiative transfer analysis reveals the vertical distribution of clouds and hazes. Chronicling temporal variabilities of atmospheric features is key to understanding the various processes that operate within the dynamic Uranian atmosphere. A comprehensive understanding of Uranus' atmosphere is essential in deciphering the planet's origins and evolution as well as providing insight into the broader field of planetary and exoplanetary science. Magnetospheric Interactions: Imaging science is critical to understand the elusive Uranian aurorae. Mapping the aurorae and their temporal dynamics will reveal the interaction of Uranus’ complex magnetosphere with the solar wind. The apparent decades-long cooling of Uranus’ upper atmosphere is thought to be due to changes in energy input from the overlying magnetosphere, and so mapping and quantifying the auroral evidence for this energy flux is particularly relevant for upper atmospheric science In addition, Uranus’ large moons may generate auroral footprints, which can be used to constrain the planet’s unusual magnetic field structure. Satellite Interiors: Imaging science can address the potential habitability of the larger moons of Uranus by measuring their shapes and rotational states (including librations), which combined with gravity data are crucial to understanding their internal density structure. Satellite Origins: Origin markers such as organics and volatiles may be preserved on the satellite surfaces, which may be sensed by UV/IR mapping. Imaging combined with the moons’ densities may reveal whether the moons formed in Uranus’ circumplanetary disk, rings, or from ejecta of a giant impact on Uranus. Satellite Geology: Imaging investigations, including stereo imaging to obtain topography, will enable constraining each satellite’s geologic history by investigating the distributions and morphologies of impact craters, tectonic and potential cryovolcanic features, and regolith and mass wasting deposits. These records will help unravel the past complex history of orbital resonances that occurred between Uranian satellites, which may have driven geologic activity. Analyses of these features will also point to resurfacing styles and modification processes important for understanding each moon’s thermal evolution, and any present-day geological activity, including potential subsurface deep oceans or possible current activity. Additionally, hyperspectral imaging will enable mapping the surface composition and identifying the potential presence of material that may point to a subsurface ocean, such as ammonia-bearing species. Rings and Small Moons: Imaging observations elucidate the physical processes that contribute to the complex ring-moon system’s inferred unstable state. Imaging surveys will seek smaller moons and other debris that could provide evidence of past collisions, while precise astrometric data of moons will determine the current dynamical state of the system. Many of Uranus' rings are exceptionally narrow, and images of the rings covering a broad range of longitudes, and high-resolution images of features within the rings will reveal how these rings are confined and maintained, as well as illustrating the rings' internal structure and dynamics. Images of the rings and moons obtained at a variety of wavelengths and lighting conditions are also needed to quantify the composition and the size distribution of the small satellites and ring particles.

K M Sayanagi↗

NASA NTRS1998

Multiwavelength Observations of a Dramatic High-Energy Flare in the Blazar 3C 279

The blazar 3C 279, one of the brightest identified extragalactic objects in the gamma-ray sky, underwent a large (factor of approx. 10 in amplitude) flare in gamma-rays toward the end of a 3 week pointing by Compton Gamma Ray Observatory (CGRO), in 1996 January-February. The flare peak represents the highest gamma-ray intensity ever recorded for this object. During the high state, extremely rapid gamma-ray variability was seen, including an increase of a factor of 2.6 in approx. 8 hrs., which strengthens the case for relativistic beaming. Coordinated multifrequency observations were carried out with Rossi X-Ray Timing Explorer (RXTE), Advanced Satellite for Cosmology and Astrophysics (ASCA; or, Astro-D), Roentgen Satellite (ROSAT), and International Ultraviolet Explorer (IUE) and from many ground-based observatories, covering most accessible wavelengths. The well-sampled, simultaneous RXTE light curve shows an outburst of lower amplitude (factor of approx. or equal to 3) well correlated with the gamma-ray flare without any lag larger than the temporal resolution of approx. 1 day. The optical-UV light curves, which are not well sampled during the high-energy flare, exhibit more modest variations (factor of approx. 2) and a lower degree of correlation. The flux at millimetric wavelengths was near a historical maximum during the gamma-ray flare peak, and there is a suggestion of a correlated decay. We present simultaneous spectral energy distributions of 3C 279 prior to and near to the flare peak. The gamma-rays vary by more than the square of the observed IR-optical flux change, which poses some problems for specific blazar emission models. The synchrotron self-Compton (SSC) model would require that the largest synchrotron variability occurred in the mostly unobserved submillimeter/far-infrared region. Alternatively, a large variation in the external photon field could occur over a time scale of a few days. This occurs naturally in the "mirror" model wherein the flaring region in the jet photoionizes nearby broad emission line clouds, which, in turn, provide soft external photons that are Comptonized to gamma-ray energies.

Wehrle, A. E.↗

DOE OSTI2021

Combined Mesonet and Tracker

Title: Combined Mesonet and Trackers (UNL Mobile Mesonets) Authors University of Nebraska PI: Adam Houston, UNL Professor (ahouston2@unl.edu) Mailing Address: 126 Bessey Hall P.O. Box 880340 Lincoln, NE 68588-0340 CoMeT Overview The University of Nebraska-Lincoln operates three Combined Mesonet and Tracker (CoMeTs). CoMeTs are Ford Explorers (model years 2015, 2017, and 2019) with forward-mounted suites of meteorological sensors and dual moonroofs, combining the capability of a mobile mesonet to collect near-surface observations with the capability of an unmanned aircraft systems (UAS) tracker vehicle, which enables an observer in the second row of seats to see the aircraft and maintain compliance with Federal Aviation Administration policies on UAS operation. The CoMeTs collect observations of slow temperature and humidity at ~2 m above ground level (AGL) using a Vaisala HMP155A, fast temperature at ~2 m AGL using a Campbell Scientific 109SS-L thermistor, pressure at ~2.5 m AGL using a Vaisala PTB210, wind speed and direction at ~3.25 m AGL using an R.M. Young 05103 propeller anemometer, and vehicle heading using a KVH Industries C-100 fluxgate compass (Barbieri et al. 2019). The HMP155A and 109SS-L thermistor are shielded and aspirated within a U-tube (Waugh and Frederickson 2010; Houston et al. 2016). This list of sensors is also included in the CoMeT data file metadata. Manufacturer specifications for these instruments are given in Table 1 of Hanft and Houston (2018). The reported measured quantities are summarized below.CoMeT-3 was funded through an equipment allocation included in the NSF TORUS award (AGS-1824649). Instrument Description The specific sensors included on each CoMeT are summarized in the table at the end of this section. In general each CoMeT collects observations of slow temperature and humidity at ~2 m above ground level (AGL) using a Vaisala HMP155, fast temperature at ~2 m AGL using a Campbell Scientific 109SS-L thermistor, pressure at ~2.5 m AGL using a Vaisala PTB210 barometer with a Gill pressure port, wind speed and direction at ~3.25 m AGL using an R.M. Young 05103 propeller anemometer, position using a Garmin 19x HVS receiver, and vehicle heading using a KVH Industries C-100 fluxgate compass. The HMP155 and 109SS are shielded and aspirated within a U-tube (Waugh and Frederickson 2010; Houston et al. 2016). Fast temperature and corrected RH measurements (using sensors housed within the U-tube) have a time constant of 10-12 s based on data collected across a temperature and RH shock during the CLOUD-MAP 2017 calibration/validation tests on June 26, 2017. Vehicle speed was < 10 kts for this test. CoMeT-1 CoMeT-2 CoMeT-3 Slow Temperature Slow RH Vaisala HMP155A-L20-PT Part #: 22280-7 Vaisala HMP155E Part #: E1AA11A0B1A1A0A Vaisala HMP155E Part #: E1AA11A0B1A1A0A Fast temperature Campbell Scientific 109SS-L20-PT Part #: 21448-3 Campbell Scientific 109SS-L12-PW Part #: 21448-109 Campbell Scientific 109SS-L12-PW Part #: 21448-150 Pressure Vaisala PTB-210 Part #: A1A1B Gill Pressure Port Part #: 61002 Vaisala PTB-210 Part #: A1A1B Gill Pressure Port Part #: 61002 Vaisala PTB-210 Part #: A1A1B Gill Pressure Port Part #: 61002 Wind RM Young 05103-L20-PT Part #: 18435-310 RM Young 05103-L20-PW Part #: 18435-244 RM Young 05103-L20-PW Part #: 18435-244 GPS Garmin GPS 19x HVS (NMEA 0183) Part #: 010-01010-00 Garmin GPS 19x HVS (NMEA 0183) Part #: 010-01010-00 Garmin GPS 19x HVS (NMEA 0183) Part #: 010-01010-00 Compass KVH C-100 Part #: 01-0177-15 KVH C-100 Part #: 01-0177-15 KVH C-100 Part #: 01-0177-15 Logger Campbell Scientific CR6-NA-XT-SW Part #: 28385-9 Campbell Scientific CR6-WIFI-XT-SW Part #: 28385-6 Campbell Scientific CR6-WIFI-XT-SW Part #: 28385-6 Data Collection and Real-Time Processing The reported measured quantities are summarized in the table below. Quantity Units Source Epoch time Seconds GPS Latitude and longitude Degrees GPS Altitude m GPS Pressure hPa PTB210 Temperature (fast) deg C 109SS-L Temperature (slow) deg C HMP155 RH (slow) % HMP155 Vehicle speed m/s GPS Vehicle heading deg C-100 and GPS In addition to the measured variables, several derived variables are calculated. Corrected/fast relative humidity (%) Relative humidity is adjusted to the fast temperature following Richardson et al. (1998) and Houston et al. (2016). Water vapor mixing ratio (g/kg) Dew point temperature (°C) Potential temperature (Kelvin) Virtual potential temperature (Kelvin) Equivalent potential temperature (Kelvin) Regular intercomparisons between all three CoMeTs were performed during TORUS 2019. Comparisons were also conducted between CoMeT-1 and CoMeT-2 during LAPSE-RATE (2018) on 14 July. In these intercomparisons, the vehicles were parked adjacent to each other aligned perpendicular to (and facing into) the wind. To minimize engine heating effects, intercomparisons were only conducted when the wind speed was >10 kts. Data Format Original data files for each deployment are saved as text files and then converted to NetCDF. NetCDF versions have units that are CF compliant and may not match the original units in the txt files. The naming convention for the NetCDF files is as follows: UNL.CoMeT3.{deployment date YYYYMMDD}.{start time of observation collection in UTC HHMM}.L2.{post-processing codes}.cdf example: UNL.CoMeT3.20190627.1931.L2.g1.f1.cdf Post-processing codes are included to track modifications to the raw data. These codes are closely connected to error flags associated with each record. Each letter corresponds to a particular instrument: g: GPS p: Barometer tf: Fast temperature ts: Slow temperature rh: Relative humidity f: Compass w: Wind monitor a: All instruments Each number corresponds to a particular post-processing action described more below. Measured and derived variables are included in the following table. Variable Heading Standard Name Units time Time seconds since 00:00:00, 01-01-1970 Alt Altitude meters lat Latitude degrees north lon Longitude degrees east fast_temp Air Temperature kelvin slow_temp Air Temperature kelvin pressure Air Pressure pascals logger_RH Relative Humidity percent calc_corr_RH Relative Humidity percent wind_speed Wind Speed meters per second wind_dir Wind From Direction degrees vehicle_dir Vehicle Direction degrees dewpoint Dew Point Temperature kelvin mixing_ratio Humidity Mixing Ratio g/g theta Air Potential Temperature kelvin theta_v Virtual Potential Temperature kelvin theta_e Equivalent Potential Temperature Kelvin error_flag The error_flag variable is a string that matches the post-processing codes listed above. All instruments will have an associated code, but will have a “0” if the datum is unchanged from the initial processed value. Error Codes The following table summarizes the error codes for data collected before 2020: Error Code Relevant CoMeT Description g1 1,2,3 Exact correction. GPS position and time reprocessed from raw data g2 1 As far as we can tell this is an exact correction to an error in the GPS time. During the correct time periods the time suddenly went backwards ~250s and stayed at this offset for 750s when it corrected itself. The offset was applied to the “time warp” period. p1 2 Approximate correction. Hole in the pressure tube connecting the pressure port to the barometer. Resulted in erroneously low air pressure measurements when the vehicle was in motion. Derived variables recalculated (dew point temperature [e depends on qv and p], water vapor mixing ratio, potential temperature, virtual potential temperature, equivalent potential temperature) a1 3 Exact correction. Missing data reprocessed from raw data a2 1 Bug fix to bias correction for ts1, ts2, and rh1: water vapor mixing ratio was off by a factor of 10 and virtual potential temperature was wrong because of this. f1 3 No correction, missing data. Fluxgate compass inoperable. Wind speed and direction calculated using GPS-derived vehicle heading instead. rh1 1 Approximate correction. Constant bias of +1.7% removed from relative humidity. Derived variables recalculated (corrected/fast relative humidity, dew point temperature, water vapor mixing ratio, virtual potential temperature, equivalent potential temperature) ts1 1 Approximate correction. Constant bias of +0.6 K removed from slow temperature. Derived variables recalculated (corrected/fast relative humidity, dew point temperature, water vapor mixing ratio, virtual potential temperature, equivalent potential temperature) ts2 1 Approximate correction. Constant bias of +1.0 K removed. Derived variables recalculated (corrected/fast relative humidity, dew point temperature, water vapor mixing ratio, virtual potential temperature, equivalent potential temperature) References Bolton, D., 1980: The Computation of Equivalent Potential Temperature. Mon. Wea. Rev., 108, 1046–1053, https://doi.org/10.1175/1520-0493(1980)108<1046:TCOEPT>2.0.CO;2. Hanft, W., and A. L. Houston, 2018: An Observational and Modeling Study of Mesoscale Air Masses with High Theta-E. Mon. Wea. Rev., 146, 2503–2524, https://doi.org/10.1175/MWR-D-17-0389.1.Wexler Houston, A. L., R. J. Laurence III, T. W. Nichols, S. Waugh, B. Argrow, and C. L. Ziegler, 2016: Intercomparison of unmanned aircraft-borne and mobile mesonet atmospheric sensors. Journal of Atmospheric and Oceanic Technology. 33, 1569-1582, doi: 10.1175/JTECH-D-15-0178.1. Lowe, P. R., 1977: An Approximating Polynomial for the Computation of Saturation Vapor Pressure. J. Applied Meteorology, 16, 100–103. Richardson, S. J., S. E. Frederickson, F. V. Brock, and J. A. Brotzge, 1998: Combination temperature and relative humidity probes: Avoiding large air temperature errors and associated relative humidity errors. Preprints, 10th Symp. On Meteorological Observations and Instrumentation, Phoenix, AZ, Amer. Meteor. Soc., 278–283. Waugh, S., and S. E. Frederickson, 2010: An improved aspirated temperature system for mobile meteorological observations, especially in severe weather. 25th Conf. on Severe Local Storms, Denver, CO, Amer. Meteor. Soc., P5.2. [Available online at https://ams.confex.com/ams/25SLS/techprogram/paper_176205.htm.]

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2018

Water Across Synthetic Aperture Radar Data (WASARD): SAR Water Body Classification for the Open Data Cube

The detection of inland water bodies from Synthetic Aperture Radar (SAR) data provides a great advantage over water detection with optical data, since SAR imaging is not impeded by cloud cover. Traditional methods of detecting water from SAR data involves using thresholding methods that can be labor intensive and imprecise. This paper describes Water Across Synthetic Aperture Radar Data (WASARD): a method of water detection from SAR data which automates and simplifies the thresholding process using machine learning on training data created from Geoscience Australia’s WOFS algorithm. Of the machine learning models tested, the Linear Support Vector Machine was determined to be optimal, with the option of training using solely the VH polarization or a combination of the VH and VV polarizations. WASARD was able to identify water in the target area with a correlation of 97% with WOFS. Sentinel-1, Open Data Cube, Earth Observations, Machine Learning, Water Detection 1. INTRODUCTION Water classification is an important function of Earth imaging satellites, as accurate remote classification of land and water can assist in land use analysis, flood prediction, climate change research, as well as a variety of agricultural applications [2]. The ability to identify bodies of water remotely via satellite is immensely cheaper than contracting surveys of the areas in question, meaning that an application that can accurately use satellite data towards this function can make valuable information available to nations which would not be able to afford it otherwise. Highly reliable applications for the remote detection of water currently exist for use with optical satellite data such as that provided by LANDSAT. One such application, Geoscience Australia’s Water Observations from Space (WOFS) has already been ported for use with the Open Data Cube [6]. However, water detection using optical data from Landsat is constrained by its relatively long revisit cycle of 16 days [5], and water detection using any optical data is constrained in that it lacks the ability to make accurate classifications through cloud cover [2]. The alternative solution which solves these problems is water detection using SAR data, which images the Earth using cloud-penetrating microwaves. Because of its advantages over optical data, much research has been done into water detection using SAR data. Traditionally, this has been done using the thresholding method, which involves picking a polarization band and labeling all pixels for which this band’s value is below a certain threshold as containing water. The thresholding method works since water tends to return a much lower backscatter value to the satellite than land [1]. However, this method can be flawed since estimating the proper threshold is often imprecise, complicated, and labor intensive for the end user. Thresholding also tends to use data from only one SAR polarization, when a combination of polarizations can provide insight into whether water is present. [2] In order to alleviate these problems, this paper presents an application for the Open Data Cube to detect water from SAR data using support vector machine (SVM) classification. 2. PLATFORM WASARD is an application for the Open Data Cube, a mechanism which provides a simple yet efficient means of ingesting, storing, and retrieving remote sensing data. Data can be ingested and made analysis ready according to whatever specifications the researcher chooses, and easily resampled to artificially alter a scene’s resolution. Currently WASARD supports water detection on scenes from ESA’s Sentinel-1 and JAXA’s ALOS. When testing WASARD, Sentinel-1 was most commonly used due to its relatively high spatial resolution and its rapid 6 day revisit cycle [5]. With minor alterations to the application's code, however, it could support data from other satellites. 3. METHODOLOGY Using supervised classification, WASARD compares SAR data to a dataset pre-classified by WOFS in order to train an SVM classifier. This classifier is then used to detect water in other SAR scenes outside the training set. Accuracy was measured according to the following metrics:  Precision: a measure of what percentage of the points WASARD labels as water are truly water  Recall: a measure of what percentage of the total water cover WASARD was able to identify.  F1 Score: a harmonic average of the precision and recall scores Both precision and recall are calculated at the end of the training phase, when the trained classifier is compared to a testing dataset. Because the WOFS algorithm’s classifications are used as the truth values when training a WASARD classifier, when precision and recall are mentioned in this paper, they are always with respect to the values produced by WOFS on a similar scene of Landsat data, which themselves have a classification accuracy of 97% [6]. Visual representations of water identified by WASARD in this paper were produced using the function wasard_plot(), which is included in WASARD. 3.1 Algorithm Selection The machine learning model used by WASARD is the Linear Support Vector Machine (SVM). This model uses a supervised learning algorithm to develop a classifier, meaning it creates a vector which can be multiplied by the vector formed by the relevant data bands to determine whether a pixel in a SAR scene contains water. This classifier is trained by comparing data points from selected bands in a SAR scene to their respective labels, which in this case are “water” or “not water” as given by the WOFS algorithm. The SVM was selected over the Random Forest model, which outperformed the SVM in training speed, but had a greater classification time and lower accuracy, and the Multilayer Perceptron Artificial Neural Network, which had a slightly higher average accuracy than the SVM, but much greater training and classification times. Figure 1: Visual representation of the SVM Classifier. Each white point represents a pixel in a SAR scene. In Figure 1, the diagonal line separating pixels determined to be water from those determined not to be water represents the actual classification vector produced by the SVM. It is worth noting that once the model has been trained, classification of pixels is done in a similar manner as in the thresholding method. This is especially true if only one band was used to train the model. 3.1 Feature Selection Sentinel-1 collects data from two bands: the Vertical/Vertical polarization (VV) and the Vertical/Horizontal polarization (VH). When 100 SVM classifiers were created for each polarization individually, and for the combination of the two, the following results were achieved: Figure 2: Accuracy of classifiers trained using different polarization bands. Precision and Recall were measured with respect to the values produced by WOFS. Figure 2 demonstrates that using both the VV and VH bands trades slightly lower recall for significantly greater precision when compared with the VH band alone, and that using the VV band alone is inferior in both metrics. WASARD therefore defaults to using both the VV and VH bands, and includes the option to use solely the VH band. The VV polarization’s lower precision compared to the VH polarization is in contrast to results from previous research and may merit further analysis [4]. 3.2 Training a Classifier The steps in training a classifier with WASARD are 1. Selecting two scenes (one SAR, one optical) with the same spatial extents, and acquired close to each other in time, with a preference that the scenes are taken on the same day. 2. Using the WOFS algorithm to produce an array of the detected water in the scene of optical data, to be used as the labels during supervised learning 3. Data points from the selected bands from the SAR acquisition are bundled together into an array with the corresponding labels gathered from WOFS. A random sample with an equal number of points labeled “Water” and “Not Water” is selected to be partitioned into a training and a testing dataset 4. Using Scikit-Learn’s LinearSVC object, the training dataset is used to produce a classifier, which is then tested against the testing dataset to determine its precision and recall The result is a wasard_classifier object, which has the following attributes: 1. f1, recall, and precision: 3 metrics used to determine the classifier’s accuracy 2. Coefficient: Vector which the SVM uses to make its predictions. The classifier detects water when the dot product of the coefficient and the vector formed by the SAR bands is positive 3. Save(): allows a user to save a classifier to the disk in order to use it without retraining 4. wasard_classify(): Classifies an entire xarray of SAR data using the SVM classifier All of the above steps are performed automatically when the user creates a wasard_classifier object. 3.3 Classifying a Dataset Once the classifier has been created, it can be used to detect water in an xarray of SAR data using wasard_classify(). By taking the dot product of the classifier’s coefficients and the vector formed by the selected bands of SAR data, an array of predictions is constructed. A classifier can effectively be used on the same spatial extents as the ones where it was trained, or on any area with a similar landscape. While

Kreiser, Zachary↗

DOE OSTI2021

Cloud droplet diffusional growth in homogeneous isotropic turbulence: bin microphysics versus Lagrangian super-droplet simulations

The increase in the spectral width of an initially monodisperse population of cloud droplets in homogeneous isotropic turbulence is investigated by applying a finite-difference fluid flow model combined with either Eulerian bin microphysics or a Lagrangian particle-based scheme. The turbulence is forced applying a variant of the so-called linear forcing method that maintains the mean turbulent kinetic energy (TKE) and the TKE partitioning between velocity components. The latter is important for maintaining the quasi-steady forcing of the supersaturation fluctuations that drive the increase in the spectral width. We apply a large computational domain (64 3 m 3 ), one of the domains considered in Thomas et al. (2020). The simulations apply 1 m grid length and are in the spirit of the implicit large eddy simulation (ILES), that is, with small-scale dissipation provided by the model numerics. This is in contrast to the scaled-up direct numerical simulation (DNS) applied in Thomas et al. (2020). Two TKE intensities and three different droplet concentrations are considered. Analytic solutions derived in Sardina et al. (2015), valid for the case when the turbulence integral timescale is much larger than the droplet phase relaxation timescale, are used to guide the comparison between the two microphysics simulation techniques. The Lagrangian approach reproduces the scalings relatively well. Representing the spectral width increase in time is more challenging for the bin microphysics because appropriately high resolution in the bin space is needed. The bin width of 0.5 µm is only sufficient for the lowest droplet concentration (26 cm -3 ). For the highest droplet concentration (650 cm -3 ), an order of magnitude smaller bin size is barely sufficient. The scalings are not expected to be valid for the lowest droplet concentration and the high-TKE case, and the two microphysics schemes represent similar departures. Finally, because the fluid flow is the same for all simulations featuring either low or high TKE, one can compare point-by-point simulation results. Such a comparison shows very close temperature and water vapor point-by-point values across the computational domain and larger differences between simulated mean droplet radii and spectral width. The latter are explained by fundamental differences in the two simulation methodologies, numerical diffusion in the Eulerian bin approach and a relatively small number of Lagrangian particles that are used in the particle-based microphysics.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Influence of Lake Ice Biases in Reanalysis Data on Downscaled Climate Simulations over the Great Lakes Region

This data package contains observation-based and model-simulated datasets (all provided in NetCDF format) for evaluating how wintertime lake-ice representation affects regional weather and climate over the Laurentian Great Lakes (freshwater lake ecosystem) during the high–ice-cover winter of 2009. The observational component includes: (1) Stage IV gridded precipitation at 4 km, hourly resolution for January–February 2009 over the Great Lakes region (radar–gauge multisensor precipitation analyses); (2) Great Lakes Surface Environmental Analysis (GLSEA) satellite-derived lake-ice coverage at 1.3 km, daily resolution for the 2009 winter months, providing ice coverage over Lakes Superior, Michigan, Huron, Erie, and Ontario; and (3) in situ measurements at the Standard Rock site on Lake Superior from the Great Lakes Evaporation Network (GLEN) at hourly resolution, including near-surface atmospheric variables and sensible and latent heat fluxes (air–lake exchange) at a fixed point location. The modeling component provides corresponding fields from two simulations, both archived at 4 km, hourly resolution: a standalone Weather Research Forecasting model (WRF) run driven by the European Centre for Medium-Range Weather Forecasts (ECMWF) Reanalysis 5 (ERA5), and a two-way coupled model using WRF and the Finite Volume Community Ocean Model (WRF-FVCOM, a 3-D hydrodynamic lake model). These outputs include variables relevant to air–lake interaction and lake-effect processes (e.g., near-surface temperature, humidity, wind, precipitation, and surface turbulent fluxes), enabling direct comparison with the observational datasets. Users can analyze and visualize these NetCDF files with common tools such as Python (e.g., xarray, netCDF4, numpy, pandas), NCO/CDO, Panoply, or ncview; NetCDF variables can also be converted to other formats (e.g., CSV, GeoTIFF) using these utilities.

EARTH SCIENCE > ATMOSPHERE > ATMOSPHERIC TEMPERATU↗