SEARCH · Engineering Papers

Results for “wells”

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,711 records · Page 95

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Modeling and Assessment of Integrated Hydroelectric Power and Hydrogen Energy Storage

Environmental and regulatory constraints are increasingly posing challenges to the optimal operation of hydroelectric facilities. Hydrogen energy storage (HES) is well-suited to address these challenges.Using excess or low-cost hydroelectric power during off-peak hours to produce hydrogen can help reduce water spills and unlock opportunities for energy arbitrage. Moreover, the oxygen byproduct of electrolysis can be used to improve the dissolved oxygen (DO) levels in tailrace water. %, a key water quality metric for hydroelectric facilities.This paper presents an innovative framework for modeling and assessing the integrated hydroelectric-hydrogen system, considering subsystem constraints and coupling, benefits in spill water reduction and DO mitigation, as well as diverse hydrogen utilization, including re-electrification, blending with natural gas, and direct sales. Comprehensive case studies were conducted to assess the techno-economic performance of the integrated system in different deployment scenarios.

Energy storage, Hydrogen, Hydropower, Optimization↗

DOE OSTI2021

High-Voltage CMOS Active Pixel Sensor

The high-voltage CMOS (HVCMOS) sensors are a novel type of CMOS active pixel sensors for ionizing particles that can be implemented in CMOS processes with deep n-well option. The pixel contains one sensor electrode formed with a deep n-well implanted in a p-type substrate. CMOS pixel electronics, embedded in shallow wells, are placed inside the deep n-well. By biasing the substrate with a high negative voltage and by the use of a lowly doped substrate, a depleted region depth of at least 30 µm can be achieved. The electrons generated by a particle are collected by drift, which induces fast detectable signals. This publication presents a 4.2-cm 2 large HVCMOS pixel sensor implemented in a commercial 180-nm process on a lowly doped substrate and its characterization.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Short-Term Load Forecasting Considering EV Charging Loads with Prediction Interval Evaluation

Short-term load forecasting plays a critical role in power system planning and operation. Along with the electrification of various loads, electricity demands are becoming increasingly hard to predict. Notably, the recent rise in electric vehicles (EVs) has further contributed to this unpredictability. To address this issue, this paper proposes a probabilistic load forecasting strategy utilizing Gaussian process regression, structured in a day-ahead manner. While many works focus on deterministic prediction, probabilistic forecasting offers additional insights into variability and uncertainty, enabling more flexible and reliable operation for power systems. To enhance the accuracy of the load forecasting model, the inputs include features related to EV charging habits as well as commonly used weather information. The load forecasting results are evaluated using various metrics, including conventional ones that assess the accuracy of point forecasts, as well as additional metrics that test the reliability of prediction intervals. The proposed load forecasting method is finally tested on real residential power consumption data and EV charging data sampled from real-world sources. The results prove that the new features can greatly improve the performance of the load forecasting method.

electrical vehicle↗

DOE OSTI2022

The Latent Dirichlet Allocation model applied to airborne LiDAR data: A case study on mapping forest degradation associated with fragmentation and fire in the Amazon region

1. LiDAR data are being increasingly used to provide a detailed characterization of the vertical profile of forests. This characterization enables the generation of new insights on the influence of environmental drivers and anthropogenic disturbances on forest structure as well as on how forest structure influences important ecosystem functions and services. Unfortunately, extracting information from LiDAR data in a way that enables the spatial visualization of forest structure, as well as its temporal changes, is challenging due to the high dimensionality of these data. 2. In this study, we show how the Latent Dirichlet Allocation model applied to LiDAR data (LidarLDA) can be used to identify forest structural types and how the relative abundance of these forest types changes throughout the landscape. The code to fit this model is made available through the open-source r package LidarLDA in github. We illustrate the use of LidarLDA both with simulated data and data from a large-scale fire experiment in the Brazilian Amazon region. 3. Using simulated data, we demonstrate that LidarLDA accurately identifies the number of forest types as well as their spatial distribution and absorptance probabilities. For the empirical data, we found that LidarLDA detects both landscape-level patterns in forest structure as well as the strong interacting effect of fire and forest fragmentation on forest structure based on the experimental fire plots. More specifically, LidarLDA reveals that proximity to forest edge exacerbates the impact of fires, and that burned forests remain structurally different from unburned areas for at least 7 years, even when burned only once. Importantly, LidarLDA generates insights on the 3D structure of forest that cannot be obtained using more standard approaches that just focus on top-of-the-canopy information (e.g. canopy height models based on LiDAR data). 4. By enabling the mapping of forest structure and its temporal changes, we believe that LidarLDA will be of broad utility to the ecological research community.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Greenhouse gas emissions from the global transportation of crude oil: Current status and mitigation potential

Global crude-oil transportation contributes a significant portion of greenhouse gas (GHG) emissions in the marine transportation sector. In this work, we first compile a detailed country-level global crude-oil transportation network in 2018 and estimate that the direct and well-to-hull GHG emissions related to crude transportation were 97 and 109 million metric tons, respectively. Combining with the country-specific crude recovery GHG intensities, the consumption-based well-to-country-gate crude-oil GHG intensities are derived for individual countries, ranging from 2.99 to 27.32 g CO 2 eq/MJ, with a global crude-volume-weighted average of 8.67 g CO 2 eq/MJ. We then project the global crude transportation emissions at the regional level in 2050 under a static (no change) scenario (based on current ship energy efficiency) and a sustainable-development (SD) scenario (based on the International Energy Agency's projections of ship energy efficiency and penetration of alternative marine fuels). Results show that the global well-to-hull GHG emissions related to crude transportation would be 82 and 59 million metric tons in 2050 in the static and SD scenarios, respectively. To further evaluate the impact of potential fuel-switching on decarbonizing the crude oil transportation sector, we estimate the GHG emissions for 20 fuel/production options in 2050. We find that, in comparison to the static scenario, ~50% reduction in global well-to-hull GHG emissions from crude transportation could be achieved under the SD scenario if green ammonia further replaces conventional ammonia. The methodology developed here can be applied to other commodities to estimate the emissions associated with their global marine transportation and to evaluate the potential emission mitigation options.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Evaluation of Energy Storage Potential of Unconventional Shale Reservoirs Using Numerical Simulation of Cyclic Gas Injection

Compressed air energy storage (CAES) stores energy as compressed air in underground formations, typically salt dome caverns. When electricity demand grows, the compressed air is released through a turbine to produce electricity. CAES in the US is limited to one plant built in 1991, due in part to the inherent risk and uncertainty of developing subsurface storage reservoirs. As an alternative to CAES, we propose using some of the hundreds of thousands of hydraulically fractured horizontal wells to store energy as compressed natural gas in unconventional shale reservoirs. To store energy, produced or “sales” natural gas is injected back into the formation using excess electricity and is later produced through an expander to generate electricity. To evaluate this concept, we performed numerical simulations of cyclic natural gas injection into unconventional shale reservoirs using cmg-gem commercial reservoir modeling software. We tested short-term (diurnal) and long-term (seasonal) energy storage potential by modeling well injection and production gas flowrates as a function of bottom-hole pressure. First, we developed a conceptual model of a single fracture stage in an unconventional shale reservoir to characterize reservoir behavior during cyclic injection and production. Next, we modeled cyclic injection in the Marcellus shale gas play using published data. Results indicate that Marcellus unconventional shale reservoirs could support both short- and long-term energy storage at capacities of 100–1000 kWe per well. The results indicate that energy storage in unconventional shale gas wells may be feasible and warrants further investigation.

25 ENERGY STORAGE↗

DOE OSTI2021

Large Eddy Simulation of Lean Mixed-Mode Combustion Assisted by Partial Fuel Stratification in a Spark-Ignition Engine

Partial fuel stratification (PFS) is a promising fuel injection strategy to improve the stability of lean combustion by applying a small amount of pilot injection right before spark timing. Mixed-mode combustion, which makes use of end-gas autoignition following conventional deflagration-based combustion, can be further utilized to speed up the overall combustion. In this study, PFS-assisted mixed-mode combustion in a lean-burn direct injection spark-ignition (DISI) engine is numerically investigated using multi-cycle large eddy simulation (LES). A previously developed hybrid G-equation/well-stirred reactor combustion model for the well-mixed operation is extended to the PFS-assisted operation. The experimental spray morphology is employed to derive spray model parameters for the pilot injection. The LES-based model is validated against experimental data and is further compared with the Reynolds-averaged Navier–Stokes (RANS)-based model. Overall, both RANS and LES predict the mean pressure and heat release rate traces well, while LES outperforms RANS in capturing the cycle-to-cycle variation (CCV) and the combustion phasing in the mass burned space. Liquid and vapor penetrations obtained from the simulations agree reasonably well with the experiment. Detailed flame structures predicted from the simulations reveal the transition from a sooting diffusion flame to a lean premixed flame, which is consistent with experimental findings. LES captures more wrinkled and stretched flames than RANS. Finally, the LES model is employed to investigate the impacts of fuel properties, including heat of vaporization (HoV) and laminar burning speed (SL). Combustion phasing is found more sensitive to SL than to HoV, with a larger fuel property sensitivity of the heat release rate from autoignition than that from deflagration. Moreover, the combustion phasing in the PFS-assisted operation is shown to be less sensitive to SL compared with the well-mixed operation.

fuel property sensitivity↗

DOE OSTI2025

Multidimensional Modeling of Mixture Formation in a Hydrogen-Fueled Heavy-Duty Optical Engine With Direct Injection

Hydrogen (H 2 ), as a carbon-free fuel, is considered as one of the most promising solutions to reduce the carbon footprint of hard-to-decarbonize energy and transportation sectors. As such, hydrogen-fueled internal combustion engines (H 2 ICEs) have recently been receiving increasing attention, particularly in applications such as on-road/off-road heavy-duty transport and combined heat and power. The direct injection (DI) of gaseous hydrogen into the combustion chamber offers great potential for achieving high power density and high engine efficiency, while mitigating the risk of backfire and reducing pre-ignition. However, the numerical simulation of H 2 DI system remains a formidable challenge associated with the high computational cost of reproducing compressible supersonic flow and shocks in narrow injector passages and in near-nozzle regions. In general, there is a lack of well-established and validated practices for the modeling of high-pressure H 2 DI in large-bore engines. Here, to this end, this study focuses on computational fluid dynamics (CFD) modeling of the mixture formation process in a heavy-duty optical engine employing a medium-pressure H 2 DI system. Both large eddy simulations (LES) and Reynolds Averaged Navier–Stokes (RANS) simulations are performed and evaluated against optical data. Gaseous hydrogen is injected into the combustion chamber via a centrally located outward opening hollow-cone injector at a pressure of 40 bar. Simulations are carried out for two injection timings, namely, −120 and −60 °CA. The numerical predictions for H 2 distribution in different horizontal and vertical planes during the compression stroke are systematically compared against optical data obtained through planar laser-induced fluorescence (PLIF) measurements. Overall, the LES approach using the Dynamic Structure model is found to have good predictive capabilities for the early jet penetration in terms of length and shape, as well as the later H 2 distributions. However, the unsteady RANS approach with the renormalization group $k - ϵ$ model, which is widely used by industry to model heavy-duty ICEs, significantly underpredicts the H 2 mixing, even at similar mesh resolution to that used in LES. These results indicate that there is a need for the improvement of mixing submodels within the RANS approach when applied to H 2 DI simulations.

LES↗

DOE OSTI2024

Response of Graphite to Dynamic Loading and Hypervelocity Jet Impacts

The compressive strengths of three varieties of high purity graphite, PCEA, NBG-18, and NBG-25, as well as the depth of penetration of small-scale charges into these materials was experimentally determined. These grades are similar in density, ranging from 1.80 – 1.85 g/cc, and nominal apparent porosity, ranging from 18% to 20%, but provide a wide range in maximum grain or particle size from 10s to 1000s of µm. Two very different manufacturing methods are also represented; PCEA is extruded while NBG-18 and NBG-25 are iso-molded. The quasistatic and dynamic strengths of each grade were determined on a load frame and split-Hopkinson pressure bar, respectively. The depth of penetration (DOP) of two small-scale shaped charges, the Teledyne RP-1 and RP-4, was determined against graphite. The global response of the RP-4 impacts was markedly different as the PCEA samples remained intact while all the NBG-25 samples split into 2 or 3 pieces after the jet penetration had completed. However, for all tests, the trusted DOPs fell within 2 cm. Preliminary hydrocode modeling of the penetration events used existing models that were not designed for graphite. The results can be tuned to reasonably reproduce the DOP, but the wound channel geometry is not reproduced well. A model designed for graphite would need to represent graphite’s non-linear and energy dissipation characteristics.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Multichannel spectral estimation in acoustics: A state-space approach

Spectral estimation is a necessary methodology to analyze the frequency content of noisy data sets especially in acoustic applications. Many spectral techniques have evolved starting with the classical Fourier transform methods based on the well-known Wiener-Khintchine relationship relating the covariance-to-spectral density as a transform pair culminating with more elegant model-based parametric techniques that apply prior knowledge of the data to produce a high-resolution spectral estimate. Multichannel spectral representations are a class of both nonparametric, as well as parametric, estimators that provide improved spectral estimates. In any case, classical nonparametric multichannel techniques can provide reasonable estimates when coupled with peak-peaking methods as long as the signal levels are reasonably high. Furthermore, parametric multichannel methods can perform quite well in low signal level environments even when applying simple peak-picking techniques. In this paper, the performance of both nonparametric (periodogram) and parametric (state-space) multichannel spectral estimation methods are investigated when applied to both synthesized noisy structural vibration data as well as data obtained from a sounding rocket flight. It is demonstrated that for the multichannel problem, state-space techniques provide improved performance, offering a parametric alternative compared to classical methods.

42 ENGINEERING↗

DOE OSTI2020

Absence of evidence for chiral Majorana modes in quantum anomalous Hall-superconductor devices

A quantum anomalous Hall (QAH) insulator coupled to an s-wave superconductor is predicted to harbor chiral Majorana modes. A recent experiment interprets the half-quantized two-terminal conductance plateau as evidence for these modes in a millimeter-size QAH-niobium hybrid device. However, non-Majorana mechanisms can also generate similar signatures, especially in disordered samples. Here, we studied similar hybrid devices with a well-controlled and transparent interface between the superconductor and the QAH insulator. When the devices are in the QAH state with well-aligned magnetization, the two-terminal conductance is always half-quantized. Our experiment provides a comprehensive understanding of the superconducting proximity effect observed in QAH-superconductor hybrid devices and shows that the half-quantized conductance plateau is unlikely to be induced by chiral Majorana fermions in samples with a highly transparent interface.

Science & Technology - Other Topics↗

DOE OSTI2026

Beyond $\rho ^{2/3}$ scaling: microscopic origins and multimessengers of high-density nuclear symmetry energy

The nature and Equation of State (EOS) of dense neutron-rich matter are still very poorly known, while they have broad impacts on many interesting issues in both astrophysics and nuclear physics. In particular, the nuclear symmetry energy $E_{\textrm{sym}}(\rho )$ encoding the cost to make nuclear matter more neutron-rich has been the most uncertain component of the EOS of dense neutron-rich nucleonic matter. It significantly affects the radii, tidal deformations, cooling rates, and frequencies of various oscillation modes of isolated neutron stars as well as the strain amplitude and frequencies of gravitational waves from their mergers, besides its many effects on the structures of nuclei as well as the dynamics and observables of their collisions. Siemens (1970s) observed that $E_{\textrm{sym}}(\rho )$ scales as $(\rho /\rho _0)^{2/3}$ near the saturation density $\rho _0$ of nuclear matter, since both the kinetic part and the potential contribution (quadratic in momentum) exhibit this dependence. The scaling holds if: (1) the nucleon isoscalar potential is quadratic in momentum, and (2) the isovector interaction is weakly density-dependent. After examining many empirical evidences and understanding theoretical findings in the literature, we conclude that: (1) Siemens’ $\rho ^{2/3}$ scaling is robust and serves as a valuable benchmark for both nuclear theories and experiments up to $2\rho _0$ but breaks down at higher densities, (2) Experimental and theoretical findings about $E_{\textrm{sym}}(\rho )$ up to $2\rho _0$ are broadly consistent, but uncertainties remain large for its curvature $K_{\textrm{sym}}(\rho )$ and higher-order parameters, (3) Above $2\rho _0$ , uncertainties grow due to poorly constrained spin-isospin-dependent tensor and three-body forces as well as the resulting nucleon short-range correlations. Looking forward, combining signals from both observations of neutron stars and terrestrial heavy-ion reaction experiments is the most promising path to finally constraining the high-density $E_{\textrm{sym}}(\rho )$ and the EOS of supradense neutron-rich matter. Multiple examples of community efforts to further constrain the high-density $E_{\textrm{sym}}(\rho )$ using both real and mocked data of present and future high-precision observations of neutron stars, as well as heavy-ion collisions involving high-energy rare isotopes, are briefly reviewed.

Li, Bao-An [East Texas A & M Univ., Commerce, TX (↗

DOE OSTI2021

Influence of the proton relay spacer on hydrogen electrocatalysis by cobalt hangman porphyrins

A cobalt hangman porphyrin system with a phenyl spacer between the porphyrin ring and an internal carboxylic acid group as well as its non-hangman analogue were synthesized and utilized for the study of the proton-coupled electron transfer (PCET) kinetics attendant to electrocatalytic hydrogen evolution. Cyclic voltammetry (CV) together with simulations show that a short distance between the proton relay and the redox active cobalt center as well as the increased proton donating strength results in superior catalytic activity. Here, the mechanism of hydrogen generation is at the nexus of proton transfer–electron transfer (PTET) and concerted proton–electron transfer (CPET), as opposed to an ETPT mechanism that is characteristic of hangman systems with longer proton relay networks.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Potential Dependent Mn Oxidation and Its Role in Passivation of Ni 38 Fe 20 Cr 22 Mn 10 Co 10 Multi-Principal Element Alloy Using Multi-Element Resolved Atomic Emission Spectroelectrochemistry

The identity of passivating oxides on multi-principal element alloys is of great interest as their optimization offers the potential for exceptional corrosion resistance in aqueous solutions over a broad range of potential and pH. This study focuses on a non-equiatomic Ni 38 Fe 20 Cr 22 Mn 10 Co 10 solid solution alloy and tracks the fate of each alloying element during linear sweep voltammetry, low and intermediate potential holds in the passive potential domain as well as during open circuit relaxation after anodic polarization in slightly acidified Cl − solution. Ni dissolves at all potentials investigated in this work, Fe and Co are incorporated into oxides or hydroxides in low concentrations whilst Cr and Mn are enriched at passive potentials. At low passivating potentials, Mn(II) dissolves and is incorporated in minor amounts in oxides containing large concentrations of Cr(III). Considerable enrichment in Mn(II)-species occurs relative to Cr(III) in the oxide at 0.1 V vs SCE. Electrochemical impedance spectroscopy suggests the presence of layered oxides with marginal passivation at high Mn(II) levels. The formation of these oxides depends on a combination of thermodynamic and kinetic factors as well as the sequence of passivation.

Han, Junsoo (ORCID:0000000192952954)↗

DOE OSTI2023

Catalyst Layer Resistance and Utilization in PEM Electrolysis

Improving utilization, performance, and stability of low iridium (Ir)-loaded anodes is a key goal to enable widespread adoption of polymer electrolyte membrane water electrolysis (PEMWE) for clean hydrogen production. A potential limitation is high ionic or electronic resistance of the anode catalyst layer, which leads to poor catalyst utilization, increased voltage losses, and high local overpotentials that can accelerate degradation. While catalyst layer resistance is relatively well-understood in fuel cells and other porous electrode systems, characterization of these effects is not as well established in PEMWE research. Here we present in-situ methods for measuring catalyst layer resistance in electrolysis cells using a non-faradaic H 2 /H 2 O condition as well as methods for calculating the associated voltage losses. These methods are applied to anode catalyst layers based on IrO 2 nanoparticles as well as dispersed nano-structured thin film (NSTF) Ir catalysts. Trends with anode catalyst loading and interactions between the porous transport layer and catalyst layer are investigated for IrO 2 anodes. Post-mortem microscopic analysis of durability-tested anodes is also presented, showing uneven degradation of the catalyst layer caused by catalyst layer resistance.

08 HYDROGEN↗

DOE OSTI2020

Analysis of Electric Vehicle Charging Behavior Patterns with Function Principal Component Analysis Approach

This manuscript focused on analyzing electric vehicles’ (EV) charging behavior patterns with a functional data analysis (FDA) approach, with the goal of providing theoretical support to the EV infrastructure planning and regulation, as well as the power grid load management. 5-year real-world charging log data from a total of 455 charging stations in Kansas City, Missouri, was used. The focuses were placed on analyzing the daily usage occupancy variability, daily energy consumption variability, and station-level usage variability. Compared with the traditional discrete-based analysis models, the proposed FDA modeling approach had unique advantages in preserving the smooth function behavior of the data, bringing more flexibility in the modeling process with little required assumptions or background knowledge on independent variables, as well as the capability of handling time series data with different lengths or sizes. In addition to the patterns revealed in the EV charging station’s occupancy and energy consumption, the differences between EV driver’s charging time and parking time were analyzed and called for the needs for parking regulation and enforcement. The different usage patterns observed at charging stations located on different land-use types were also analyzed.

Engineering↗

DOE OSTI2026

fp-tools

The output of a molecular dynamics simulation is a data file containing properties of the motion of the simulated particles as a function of time (e.g., positions, velocities). From this output, scientists can extract various properties which give insight into the physics of the system; for example, the radial distribution function can be computed from the particle positions, which gives insight into how close the system is to melting or freezing. Many publicly available codes to perform molecular dynamics simulations exist and have been well adopted by the scientific community (e.g., LAMMPS, VASP). However, it remains standard practice for scientists to write their own post-processing scripts to extract various properties of interest from the molecular dynamics simulation data. Since in many cases scientists are interested in the same set of “textbook” properties, there is a large duplication of effort in writing these post-processing scripts. Our proposed code, FP-Tools, is a C++ toolkit designed to help scientists extract and analyze quantities of interest from the output of a molecular dynamics simulation, eliminating the need for them to write these codes themselves and thus reducing this duplication of effort. The quantities which are computed by our code are well documented in the literature (either in publications or textbooks), and the algorithms we implement to compute these properties are also well known in the field. We are not introducing new science or methods here; rather, our primary goal is to provide a useful tool to the community in the form of a well-documented, easy to use package.

Hartman, Leah↗

DOE OSTI2024

Improved Diagnosis of Precipitation Type with LightGBM Machine Learning

Abstract Existing precipitation-type algorithms have difficulty discerning the occurrence of freezing rain and ice pellets. These inherent biases are not only problematic in operational forecasting but also complicate the development of model-based precipitation-type climatologies. To address these issues, this paper introduces a novel light gradient-boosting machine (LightGBM)-based machine learning precipitation-type algorithm that utilizes reanalysis and surface observations. By comparing it with the Bourgouin precipitation-type algorithm as a baseline, we demonstrate that our algorithm improves the critical success index (CSI) for all examined precipitation types. Moreover, when compared with the precipitation-type diagnosis in reanalysis, our algorithm exhibits increased F1 scores for snow, freezing rain, and ice pellets. Subsequently, we utilize the algorithm to compute a freezing-rain climatology over the eastern United States. The resulting climatology pattern aligns well with observations; however, a significant mean bias is observed. We interpret this bias to be influenced by both the algorithm itself and assumptions regarding precipitation processes, which include biases associated with freezing drizzle, precipitation occurrence, and regional synoptic weather patterns. To mitigate the overall bias, we propose increasing the precipitation cutoff from 0.04 to 0.25 mm h −1 , as it better reflects the precision of precipitation observations. This adjustment yields a substantial reduction in the overall bias. Finally, given the strong performance of LightGBM in predicting mixed precipitation episodes, we anticipate that the algorithm can be effectively utilized in operational settings and for diagnosing precipitation types in climate model outputs. Significance Statement Freezing rain can have significant impacts on transportation and infrastructure, making accurate prediction of precipitation types crucial. In this study, we use a machine learning method known as LightGBM to predict precipitation types. We show that the new algorithm performs better than the existing methods for all precipitation types examined. Additionally, we compute a freezing-rain climatology over the eastern United States. Although the resulting climatology pattern corresponds well to observations, the algorithm overpredicts freezing-rain occurrence. We argue that this bias can be substantially reduced by increasing the precipitation cutoff from 0.04 to 0.25 mm h −1 . Overall, this work highlights the potential of the LightGBM algorithm for both weather forecasting and diagnosing precipitation types in climate models.

Meteorology & Atmospheric Sciences↗