Search NASA⌕ Search

SEARCH · Search NASA

Results for “Seasonal forecasting”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

How Snow Aggregate Ellipsoid Shape and Orientation Variability Affects Fall Speed and Self-Aggregation Rates

Snow aggregate shapes and orientations have long been known to exhibit substantial variability. Despite this observed variability, most weather and climate prediction models use fixed power-law functions that deterministically map particle size to mass and fall speed. As such, integrated quantities like precipitation and self-aggregation rates currently ignore nonlinear effects resulting from variation in shape and orientation for aggregates of the same size. This study therefore develops an analytic framework that couples an empirically based bivariate distribution of ellipsoid shapes to classical hydrodynamic theory so as to capture an appropriate dispersion of masses, projected areas, and fall speeds for an assumed size distribution. For a fixed aggregate size, shape variations produce approximately ±0.13 m s -1 standard deviation of fall speed which increases the mass flux fall speed dispersion by more than 100% over traditional microphysics models. This increased fall speed dispersion results predominantly from shape-induced mass dispersion whereas orientation and drag dispersion play a lesser role. Shape variations can increase mass- and reflectivity-weighted fall speeds by up to 60% of traditional models whereas self-aggregation rates can increase by a factor of 100 for very small slope parameters. This implies that aggregate shape variations effectively forestall the theorized onset of fall speed distribution narrowing and subsequent quenching of the aggregation process. As a result, it is likely that secondary ice formation is necessary to prevent an ever decreasing slope parameter. The mathematical theory presented in this study is used to develop simple correction factors for snow forecast and climate models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Evaluation of the Rapid Refresh Numerical Weather Prediction Model over Arctic Alaska

Abstract Despite a need for accurate weather forecasts for societal and economic interests in the U.S. Arctic, thorough evaluations of operational numerical weather prediction in the region have been limited. In particular, the Rapid Refresh Model (RAP), which plays a key role in short-term forecasting and decision-making, has seen very limited assessment in northern Alaska, with most evaluation efforts focused on lower latitudes. In the present study, we verify forecasts from version 4 of the RAP against radiosonde, surface meteorological, and radiative flux observations from two Arctic sites on the northern Alaskan coastline, with a focus on boundary layer thermodynamic and dynamic biases, model representation of surface inversions, and cloud characteristics. We find persistent seasonal thermodynamic biases near the surface that vary with wind direction, and may be related to the RAP’s handling of sea ice and ocean interactions. These biases seem to have diminished in the latest version of the RAP (version 5), which includes refined handling of sea ice, among other improvements. In addition, we find that despite capturing boundary layer temperature profiles well overall, the RAP struggles to consistently represent strong, shallow surface inversions. Further, while the RAP seems to forecast the presence of clouds accurately in most cases, there are errors in the simulated characteristics of these clouds, which we hypothesize may be related to the RAP’s treatment of mixed-phase clouds.

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

A New Paradigm for Observing and Modeling of Air-Sea Interactions to Advance Earth System Prediction

The protection of people, property, and environmental resources from extreme weather, seasonal patterns, and climate change drives the need for predictions of weather, ocean, and climate patterns that have skill and value at timescales longer than traditional 1-10-day forecasts, including outlooks spanning weeks to decades. Advancing Earth System Prediction (ESP) skill at this range of timescales requires improved observations, understanding, and modeling of the processes in the ocean boundary layer, the atmospheric boundary layer, and their interface. A new way of referring to this coupled system is the Air-Sea Transition Zone (ASTZ). The report that follows is framed by the paradigm that the ASTZ is a single entity that regulates the flow of energy and matter between the ocean and the atmosphere. The ASTZ is thus the medium through which the ocean and atmosphere respond to and influence one another across their often disparate scales of variability. ASTZ modeling, observing, and understanding needs are particularly acute because very few measurements exist over oceans, and even fewer span the entire ASTZ, even though oceans cover 70% of Earth’s surface and are the source of most of the rain and snow that falls on both the land and the oceans.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Harnessing land-atmosphere interactions to enhance subseasonal-to-seasonal predictability

2025 Advancing Understanding of Land-Atmosphere Interactions and Processes on S2S Predictability Workshop What: 227 registered workshop participants gathered in person (43%) and online (57%) to discuss state-of-the-art scientific understanding and modeling of land-atmosphere interactions and related processes in the context of subseasonal-to-seasonal (S2S) predictability. Topics covered sources of S2S predictability, land model initialization methods, model diagnosis and evaluation metrics, AI/ML analysis and applications, and coordination of future community multi-model S2S forecast focused experiments. To advance the science, this community workshop, organized by NSF NCAR, NOAA, NASA, and DOE, aimed to 1) identify process- and application-oriented metrics for assessing S2S prediction skill and 2) develop experimental protocols for coordinated experiments to isolate, quantify, and understand the role of land-atmosphere interactions in S2S predictability. When: June 16-18, 2025 Where: Boulder, CO, USA, and online.

Land-Atmosphere Interaction↗
DOE OSTI2022

Impacts, processes and projections of the quasi-biennial oscillation

In the tropical stratosphere, deep layers of eastward and westward winds encircle the globe and descend regularly from the upper stratosphere to the tropical tropopause. With a complete cycle typically lasting almost 2.5 years, this quasi-biennial oscillation (QBO) is arguably the most predictable mode of atmospheric variability that is not linked to the changing seasons. Here, the QBO affects climate phenomena outside the tropical stratosphere, including ozone transport, the North Atlantic Oscillation and the Madden–Julian Oscillation, and its high predictability could enable better forecasts of these phenomena if models can accurately represent the coupling processes. Climate and forecasting models are increasingly able to simulate stratospheric oscillations resembling the QBO, but exhibit common systematic errors such as weak amplitude in the lowermost tropical stratosphere. Uncertainties about the waves that force the oscillation, particularly the momentum fluxes from small-scale gravity waves excited by deep convection, make its simulation challenging. Improved representation of the processes governing the QBO is expected to lead to better forecasts of the oscillation and its impacts, increased understanding of unusual events such as the two QBO disruptions observed since 2016, and more reliable future projections of QBO behaviour under climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Muon Flux Variance from Severe Atmospheric Conditions

We analyze data from the European Centre for Medium-Range Weather Forecasts (ECMWF) and the QuarkNet Research Group at Fermilab to study the variance in muon flux due to extreme weather conditions. We considered three high-intensity storms for each season: winter, spring, and summer. The most significant observation was during winter, with a ~125% increase in muon flux corresponding to a temperature drop of ~20K. Additionally, we are developing CosmicWatch muon detectors to be used in tests for the MMBC/MAMBA Bubble Chamber. A pair of these detectors will be employed to trigger on muons, providing extremely precise and inexpensive tracking granularity in the form of a muon telescope.

Aamir, Yusuf↗
DOE OSTI2023

Modeling the Impact of the Organic Aerosol Phase State on Multiphase OH Reactive Uptake Kinetics and the Resultant Heterogeneous Oxidation Timescale of Organic Aerosol in the Amazon Rainforest

Accurate predictions of chemical lifetime, i.e., oxidation timescale and change in the mass of organic aerosols (OAs) in atmospheric models, are critical to quantify the impacts of OA on aerosol–cloud interactions, radiative forcing, and air quality. The heterogeneous oxidation of OA by hydroxyl (OH) radicals is a key process governing OA mass changes during their oxidation. Recently, laboratory observations of the OH uptake coefficient (γ OH ) and heterogeneous reaction rate (k het,OH ) for different OA systems with varying phase states have been combined to develop a new parameterization of γ OH and k het,OH as a function of the OA phase state, i.e., viscosity (η OA ). Here, in this work, we use a recently published viscosity prediction framework to analyze the new γ OH and k het,OH parameterization with a box model. Subsequently, we implement this box model to predict γ OH and k het,OH over the Amazon rainforest in the dry-to-wet transition season with both significant biomass burning and biogenic influences. Relevant parameters within our box model are specified based on detailed regional model simulations over the Amazon using the Weather Research and Forecasting model coupled with Chemistry (WRF-Chem). η OA is predicted as a function of species volatility, OA composition (including water uptake) along with ambient relative humidity (RH), and temperature. Based on ambient conditions simulated by WRF-Chem over the entire atmospheric column, we use the box model to predict the upper bounds of the heterogeneous oxidation timescale of OA. Based on previous laboratory measurements, we assume that this heterogeneous OH oxidation causes fragmentation (carbon loss) of OA, followed by the evaporation of the fragmented molecules that results in the exponential decay of OA. We predict that the oxidation timescale of OA is ~1 month near the Earth’s surface in the pristine Amazonian background because of low OH concentrations. But OA is oxidized more rapidly within urban and wildfire plumes near the Earth’s surface with an order of magnitude higher OH concentrations compared to the pristine background, causing the simulated oxidation timescale of OA to be much shorter, ~3–4 days. At 3–5 km altitudes where biomass burning OA is predicted to be semisolid, the heterogeneous oxidation timescale is estimated as ~6 days to 3 weeks and decreases with increasing OH concentrations within plumes. We show that the simulated mass loss of OA is strongly size-dependent, where smaller particles oxidize more rapidly compared to larger particles due to their greater per-particle surface area-to-volume ratio. Increasing urbanization and deforestation in the Amazon in the future might increase OH concentrations in the background Amazon, causing faster oxidation of OA. At higher altitudes above liquid clouds, especially in the upper troposphere where temperatures approach ~250 K, future measurements are needed to reduce the uncertainties related to the mass loss of OA at colder temperatures and low relative humidity conditions due to its heterogeneous OH oxidation.

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

Improved Subseasonal Forecasting of Extreme Polar Vortices Using Machine Learning

Our research was focused on forecasting the position and shape of the winter stratospheric polar vortex at a subseasonal timescale of 15 days in advance. To achieve this, we employed both statistical and neural network machine learning techniques. The analysis was performed on 42 winter seasons of reanalysis data provided by NASA giving us a total of 6,342 days of data. The state of the polar vortex for determined by using geometric moments to calculate the centroid latitude and the aspect ratio of an ellipse fit onto the vortex. Timeseries for thirty additional precursors were calculated to help improve the predictive capabilities of the algorithm. Feature importance of these precursors was performed using random forest to measure the predictive importance and the ideal number of precursors. Then, using the precursors identified as important, various statistical methods were tested for predictive accuracy with random forest and nearest neighbor performing the best. An echo state network, a type of recurrent neural network that features sparsely connected hidden layer and a reduced number of trainable parameters that allows for rapid training and testing, was also implemented for the forecasting problem. Hyperparameter tuning was performed for each methods using a subset of the training data. The algorithms were trained and tuned on the first 41 years of data, then tested for accuracy on the final year. In general, the centroid latitude of the polar vortex proved easier to predict than the aspect ratio across all algorithms. Random forest outperformed other statistical forecasting algorithms overall but struggled to predict extreme values. Forecasting from echo state network suggested a strong predictive capability past 15 days, but further work is required to fully realize the potential of recurrent neural network approaches.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Climate Dynamics Preceding Summer Forest Fires in California and the Extreme Case of 2018

Recent record-breaking wildfire seasons in California prompt an investigation into the climate patterns that typically precede anomalous summer burned forest area. Using burned-area data from the U.S. Forest Service’s Monitoring Trends in Burn Severity (MTBS) product and climate data from the fifth major global reanalysis produced by the European Centre for Medium-Range Weather Forecasts (ERA5) over 1984–2018, relationships between the interannual variability of antecedent climate anomalies and July California burned area are spatially and temporally characterized. Lag correlations show that antecedent high vapor pressure deficit (VPD), high temperatures, frequent extreme high temperature days, low precipitation, high subsidence, high geopotential height, low soil moisture, and low snowpack and snowmelt anomalies all correlate significantly with July California burned area as far back as the January before the fire season. Seasonal regression maps indicate that a global midlatitude atmospheric wave train in late winter is associated with anomalous July California burned area. July 2018, a year with especially high burned area, was to some extent consistent with the general patterns revealed by the regressions: low winter precipitation and high spring VPD preceded the extreme burned area. However, geopotential height anomaly patterns were distinct from those in the regressions. Extreme July heat likely contributed to the extent of the fires ignited that month, even though extreme July temperatures do not historically significantly correlate with July burned area. While the 2018 antecedent climate conditions were typical of a high-burned-area year, they were not extreme, demonstrating the likely limits of statistical prediction of extreme fire seasons and the need for individual case studies of extreme years.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

An iterative bidirectional gradient boosting approach for CVR baseline estimation

Here this paper presents a novel Iterative Bidirectional Gradient Boosting Model (IBi-GBM) for estimating the baseline of Conservation Voltage Reduction (CVR) programs. In contrast to many existing methods, we treat CVR baseline estimation as a missing data retrieval problem. The approach involves dividing the load and its corresponding temperature profiles into three periods: pre-CVR, CVR, and post-CVR. To restore the missing load profile during the CVR period, the method employs a three-step process. First, a forward-pass GBM is executed using data from the pre-CVR period as inputs. Subsequently, a backward-pass GBM is applied using data from the post-CVR period. The two restored load profiles are reconciled, considering pre-calculated weights derived from forecasting accuracy, and only the leftmost and rightmost points are retained. The newly restored points are then included as inputs for the subsequent iteration. This iterative procedure continues until the original load data in the CVR period is fully restored. We develop IBi-GBM using actual smart meter and Supervisory Control and Data Acquisition (SCADA) data. Our results demonstrate that IBi-GBM exhibits robust performance across various data resolutions and in different seasons and outperforms existing methods by achieving a 1-2% reduction in normalized Root Mean Square Error (nRMSE).

42 ENGINEERING↗
DOE OSTI2021

Sources of Subseasonal-To-Seasonal Predictability of Atmospheric Rivers and Precipitation in the Western United States

Atmospheric rivers (ARs) account for a large portion of winter precipitation in the western US. To evaluate the sources of AR and precipitation predictability at subseasonal-to-seasonal timescales, we examine the relationships between two climate modes, the El Niño/Southern Oscillation (ENSO) and Madden-Julian oscillation (MJO), and winter hydroclimate in the western US. Our analysis uses a large ensemble of the Weather Research and Forecast (WRF) model simulations from 1981 to 2017, facilitating the assessment of uncertainty in climate mode-AR relationship due to climate variability and short length of observations. Over the North Pacific basin, we find ENSO-related latitudinal shifts of ARs, but ENSO has little effect on basin-wide averaged AR frequency. Over the western US, there is some uncertainty in the ENSO-AR connection due to the impacts of ENSO indices and datasets. However, extreme El Niño events defined by the ENSO longitude index are consistently linked to increased landfalling AR activity. The MJO can lead to significantly enhanced or suppressed landfalling AR activity in the large-ensembles, depending on the phase of MJO and time lag, while observations are too short in length to robustly show this signal. ENSO substantially modulates the MJO-AR relationship, triggering variegated responses of landfalling ARs and AR precipitation in La Niña and El Niño years. Our findings highlight the need to evaluate concurrent effects of different climate modes on ARs and precipitation, and may shed light on a path toward more accurate subseasonal-to-seasonal prediction of ARs and precipitation over the western US.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Do Machine Learning Approaches Offer Skill Improvement for Short-Term Forecasting of Wind Gust Occurrence and Magnitude?

Abstract Wind gusts, and in particular intense gusts, are societally relevant but extremely challenging to forecast. This study systematically assesses the skill enhancement that can be achieved using artificial neural networks (ANNs) for forecasting of wind gust occurrence and magnitude. Geophysical predictors from the ERA5 reanalysis are used in conjunction with an autoregressive term in regression and ANN models with different predictors, and varying model complexity. Models are derived and assessed for the warm (April–September) and cold (October–March) seasons for three high passenger volume airports in the United States. Model uncertainty is assessed by deriving models for 1000 different randomly selected training (70%) and testing (30%) subsets. Gust prediction fidelity in independent test samples is critically dependent on inclusion of an autoregressive term. Gust occurrence probabilities derived using five-layer ANNs exhibit consistently higher fidelity than those from regression models and shallower ANNs. Inclusion of the autoregressive term and increasing the number of hidden layers in ANNs from 1 to 5 also improve the model performance for gust magnitudes (lower RMSE, increased correlation, and model standard deviations that more closely approximate observed values). Deeper ANNs (e.g., 20 hidden layers) exhibit higher skill in forecasting strong (17–25.7 m s −1 ) and damaging (≥25.7 m s −1 ) wind gusts. However, such deep networks exhibit evidence of overfitting and still substantially underestimate (by 50%) the frequency of strong and damaging wind gusts at the three airports considered herein. Significance Statement Improved short-term forecasting of wind gusts will enhance aviation safety and logistics and may offer other societal benefits. Here we present a rigorous investigation of the relative skill of models of wind gust occurrence and magnitude that employ different statistical methods. It is shown that artificial neural networks (ANNs) offer considerable skill enhancement over regression methods, particularly for strong and damaging wind gusts. For wind gust magnitudes in particular, application of deeper learning networks (e.g., five or more hidden layers) offers tangible improvements in forecast accuracy. However, deeper networks are vulnerable to overfitting and exhibit substantial variability with the specific training and testing data subset used. Also, even deep ANNs reproduce only half of strong and damaging wind gusts. These results indicate the need for future work to elucidate the dynamical mechanisms of intense wind gusts and advance solutions to their prediction.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

A Comparison of Machine Learning Methods to Forecast Tropospheric Ozone Levels in Delhi

Ground-level ozone is a pollutant that is harmful to urban populations, particularly in developing countries where it is present in significant quantities. It greatly increases the risk of heart and lung diseases and harms agricultural crops. This study hypothesized that, as a secondary pollutant, ground-level ozone is amenable to 24 h forecasting based on measurements of weather conditions and primary pollutants such as nitrogen oxides and volatile organic compounds. We developed software to analyze hourly records of 12 air pollutants and 5 weather variables over the course of one year in Delhi, India. To determine the best predictive model, eight machine learning algorithms were tuned, trained, tested, and compared using cross-validation with hourly data for a full year. The algorithms, ranked by R2 values, were XGBoost (0.61), Random Forest (0.61), K-Nearest Neighbor Regression (0.55), Support Vector Regression (0.48), Decision Trees (0.43), AdaBoost (0.39), and linear regression (0.39). When trained by separate seasons across five years, the predictive capabilities of all models increased, with a maximum R 2 of 0.75 during winter. Bidirectional Long Short-Term Memory was the least accurate model for annual training, but had some of the best predictions for seasonal training. Out of five air quality index categories, the XGBoost model was able to predict the correct category 24 h in advance 90% of the time when trained with full-year data. Separated by season, winter is considerably more predictable (97.3%), followed by post-monsoon (92.8%), monsoon (90.3%), and summer (88.9%). These results show the importance of training machine learning methods with season-specific data sets and comparing a large number of methods for specific applications.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Representations of Precipitation Diurnal Cycle in the Amazon as Simulated by Observationally Constrained Cloud‐System Resolving and Global Climate Models

Abstract The ability of an observationally‐constrained cloud‐system resolving model (Weather Research and Forecasting; WRF, 4‐km grid spacing) and a global climate model (Energy Exascale Earth System Model; E3SM, 1‐degree grid spacing) to represent the precipitation diurnal cycle over the Amazon basin during the 2014 wet season is assessed. The WRF model coupled with a 3‐D variational data assimilation scheme reproduces the spatial variability of the precipitation diurnal cycle over the basin and the lifecycle of westward propagating MCSs initiated by the coastal sea‐breeze front. In contrast, a single morning peak in rainfall is produced by E3SM for simulations despite the nudging of large‐scale winds toward global reanalysis, indicating precipitation in E3SM is largely controlled by local convection associated with diurnal heating. The role of propagating MCS on the environment are discussed by using a multivariate perturbation analysis. We also find that the advection of moisture perturbations from ocean to inland regions have a higher correlation with the occurrence of MCSs in the Amazon than the intensity of colder air intrusion associated with sea breezes along the coast. Moreover, the presence of large cold pools over the central Amazon basin are responsible for the maintenance of propagating deep convection.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Interpretable Net Load Forecasting Using Smooth Multiperiodic Features

We consider the problem of forecasting net load over a horizon such as one day, using a trailing window of past net load values as well as date and time. We focus on three variations on this problem: point forecasts, marginal quantile forecasts, and generating conditional samples of the future value. We propose a method that relies on linear regression using some custom engineered time-based features to capture multiple periodicities, such as daily, weekly, and seasonal, and their interactions. Our proposed models are readily interpretable, and rely on efficient and reliable convex optimization [1] to fit. We illustrate our method on four years worth of hourly net load data, comparing predictions made with various subsets of the features.

Ogut, Mehmet G↗
DOE OSTI2024

Seasonal variability of wake impacts on US mid-Atlantic offshore wind plant power production

Abstract. The mid-Atlantic will experience rapid wind plant development due to its promising wind resource located near large population centers. Wind turbines and wind plants create wakes, or regions of reduced wind speed, that may negatively affect downwind turbines and plants. We evaluate wake variability and annual energy production with the first yearlong modeling assessment using the Weather Research and Forecasting model, deploying 12 MW turbines across the domain at a density of 3.14 MW km−2, matching the planned density of 3 MW km−2. Using a series of simulations with no wind plants, one wind plant, and complete build-out of lease areas, we calculate wake effects and distinguish the effect of wakes generated internally within one plant from those generated externally between plants. We also provide a first step towards uncertainty quantification by testing the amount of added turbulence kinetic energy (TKE) by 0 % and 100 %. We provide a sensitivity analysis by additionally comparing 25 % and 50 % for a short case study period. The strongest wakes, propagating 55 km, occur in summertime stable stratification, just when New England's grid demand peaks in summer. The seasonal variability of wakes in this offshore region is much stronger than the diurnal variability of wakes. Overall, yearlong simulated wake impacts reduce power output by a range between 38.2 % and 34.1 % (for 0 %–100 % added TKE). Internal wakes cause greater yearlong power losses, from 29.2 % to 25.7 %, compared to external wakes, from 14.7 % to 13.4 %. The overall impact is different from the linear sum of internal wakes and external wakes due to non-linear processes. Additional simulations quantify wake uncertainty by modifying the added amount of turbulent kinetic energy from wind turbines, introducing power output variability of 3.8 %. Finally, we compare annual energy production to New England grid demand and find that the lease areas can supply 58.8 % to 61.2 % of annual load. We note that the results of this assessment are not intended to make nor are they suitable to make commercial judgments about specific wind projects.

17 WIND ENERGY↗
DOE OSTI2022

Detecting and characterizing simulated sea breezes over the US northeastern coast with implications for offshore wind energy

Abstract. With the planned construction of vast offshore wind farms along the US East Coast, identifying and understanding key coastal processes, such as sea breezes, has become a critical need for the sustainability and development of US offshore wind energy. In this study, a new two-step identification method is proposed to detect and characterize three types of sea breezes (pure, corkscrew and backdoor) over the US northeastern coast from a year-long WRF (Weather Research and Forecasting) simulation. The results suggest that the proposed detection method can identify the three different types of sea breezes in the model simulation. Key sea breeze features, such as the calm zone associated with pure sea breezes and coastal jets associated with corkscrew sea breezes, are evident in the sea breeze composite imagery. In addition, the simulated sea breeze events indicate a seasonal transition from pure to corkscrew sea breeze between March and August as the land–sea thermal contrast increases. Furthermore, the location and extension of the sea breeze front are different for each type of sea breeze, suggesting that the coastal impact of sea breeze varies with sea breeze type. From the wind energy perspective, the power production associated with a 10 MW offshore wind turbine would be approximately 3 to 4 times larger during a corkscrew sea breeze event than the other two types of sea breezes. This highlights the importance of identifying the correct type of sea breeze in numerical weather/wind energy forecasting.

17 WIND ENERGY↗