Heating of the Magnetically Closed Corona and Physical Models of Solar and Stellar Spectral Irradiances
Explore the source record for details and available documents.
SEARCH · Engineering Papers
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.
Explore the source record for details and available documents.
Machine learning (ML) models have been applied to forecast solar energy; however, they often lack clarity of interpretability and underlying physics. This work addresses such challenges by developing a hierarchy of ML models that gradually introduce predictors to improve the forecast accuracy based on a physics-based framework. Three ML models (ARIMA, LSTM, and XGBoost) are examined and compared with four physics-informed persistence models reported in Part I and the simple persistence model to assess the improvement of different models. The 7-year measurements at the U.S. Department of Energy's Atmospheric Radiation Measurement's Southern Great Plains Central Facility site are used for forecasts and evaluations. The results reveal that the step-by-step introduction of predictors leads to different improvements for models at different hierarchical levels. Comparison of the ML models with persistence models shows that LSTM and XGBoost outperform all the persistence models, with LSTM having the overall best performance; however, ARIMA underperforms the four physics-informed persistence models. This study demonstrates the importance and utility of incorporating physics into ML models in improving forecast accuracy by introducing a hierarchy of physics-based predictors, distinguishing predictor contributions, and enhancing the ML interpretability. The combined use of Global Horizontal Irradiance (GHI) and Direct Normal Irradiance (DNI) significantly improves the forecast accuracy compared to using individual irradiances alone because the pair contains more information on cloud-radiation interactions.
The National Solar Radiation Database (NSRDB) provides solar resource data across the globe at a high temporal and spatial resolution. This data is primarily used in solar energy modeling. The NSRDB is updated annually for the United States and North, Central and South America and the data is currently available from 1998-2021. In 2022 the NSRDB was updated using the latest version of the underlying Physical Solar Model (PSM). This update includes improved surface albedo and gap-filling of cloud properties. The inclusion of these updates reduced the uncertainty in the data compared to previous versions of the NSRDB. The Himawari and Meteosat Indian Ocean Data Coverage (IODC) satellites were added to the Geostationary Operational Environmental Satellite (GOES) and made our coverage global. While standard data from the GOES continues to be served at an hourly 4km x 4km resolution, full resolution data has also been made available to the user. The NSRDB now contains over 200Tb of data with nearly 40Tb being added annually. We provide significant flexibility for data download depending on the amount of data required by the users. In this paper we provide an update on the current status on the NSRDB.
The National Solar Radiation Database (NSRDB) provides high-resolution spatiotemporal solar irradiance data for the entire globe. The NSRDB uses a two-step Physical Solar Model (PSM) to compute the effects of clouds and other atmospheric variables on the solar radiation reaching the surface of the Earth. Physical and optical cloud properties are fundamental inputs to the PSM and are derived from the National Oceanic and Atmospheric Administration's Geostationary Operational Environmental Satellites. This paper describes recent improvements to the NSRDB driven by physics-guided machine learning methods for cloud property retrieval. The impacts of these new methods on the NSRDB irradiance data are validated using an extensive set of ground measurement sites, showing significant improvement for all sites. We report on average, the mean absolute percentage error for global horizontal irradiance and direct normal irradiance show reductions of 2.16 and 3.95 percentage points respectively for all daylight conditions, 5.92 and 17.39 percentage points respectively for cloudy conditions, and 9.00 and 22.59 percentage points respectively for gap-filled cloudy conditions. These new methods will help improve the quality and accuracy of the irradiance and cloud data in the NSRDB.
The NASA Clouds and the Earth’s Radiant Energy System (CERES) SYN1deg(Ed4.1) provides hourly global horizontal irradiance (GHI), diffuse horizontal irradiance (DHI) and direct horizontal irradiance (DirHI) at 1-degree latitude by 1-degree longitude resolution, and the time span is from March 2000 to near present. While the GHI agrees well with the ground-based Baseline Surface Radiation Network (BSRN) data, the DHI and DirHI are, respectively, positively and negatively biased against the BSRN data, and so is the direct normal irradiance (DNI) derived simply by dividing the DirHI by the cosine of the hourly solar zenith angle (SZA). We found that the biases tend to be well-defined functions of cos(SZA) and the cloud fraction (CLFR). We thus performed a bias-correction on the hourly DHI and DNI in the latitude-cos(SZA)-CLFR phase space, and therefrom, we derived the hourly global tilted irradiance (GTI) at a number of tilt angles using the isotropic diffuse irradiance model, and the results agree well with the GTI derived from the original BSRN near-instantaneous records. Meanwhile, the National Renewable Energy Laboratory (NREL) has produced the National Solar Radiation Database (NSRDB) with finer spatiotemporal resolution using its Physical Solar Model (PSM) and surface-based measured data, model results and satellite-based data as inputs. Now the NSRDB data covers not only the United States and its neighboring regions, but Europe, Africa, Asia and Oceania. The available temporal resolutions are 5, 10, 15, 30 and 60 minutes and the spatial resolutions are 2, 4 and 10 km, depending on the selected region. In this presentation, we compare both the CERES-based and NSRDB GHI, DHI and DNI with their BSRN counterparts and show how finer spatial resolution may give us an advantage in data accuracy and usability.
In an effort to develop computational tools for predicting radiation hazards from solar energetic particles (SEPs), we have created a data-driven physics-based particle transport model to calculate the injection, acceleration, and propagation of SEPs from coronal mass ejection (CME) shocks traversing through the solar corona and interplanetary magnetic fields. The model runs on an input of corona and heliospheric plasma and magnetic field configuration from a magnetohydrodynamic model driven by solar photospheric magnetic field measurements superposed with observed CME shocks determined from coronagraph images. SEP source particles are injected at the shock using the result of diffusive shock acceleration formulation from a characteristic obliquity-dependent injection from a heated solar wind thermal tail population. With several advanced computation techniques involving stochastic simulation and integration, the model obtains the particle intensity at any location in interplanetary space through the rigorous solution to the time-dependent 5D focus transport equation in the phase space that includes perpendicular diffusion. We apply the model to the 2011 November 3 CME event. The calculation results reproduce multispacecraft SEP observations at Earth and STEREO-B reasonably well without normalization of particle flux. The observations at STEREO-A can be reproduced by rescaling particle energy or modified energy dependence of particle diffusion coefficients. This circumsolar SEP event seen by spacecraft at Earth, STEREO-A, and STEREO-B at widely separated longitudes can be explained by diffusive shock acceleration by a single CME shock with a moderate speed.
Abstract In an effort to develop computational tools for predicting radiation hazards from solar energetic particles (SEPs), we have created a data-driven physics-based particle transport model to calculate the injection, acceleration, and propagation of SEPs from coronal mass ejection (CME) shocks traversing through the solar corona and interplanetary magnetic fields. The model runs on an input of corona and heliospheric plasma and magnetic field configuration from a magnetohydrodynamic model driven by solar photospheric magnetic field measurements superposed with observed CME shocks determined from coronagraph images. SEP source particles are injected at the shock using the result of diffusive shock acceleration formulation from a characteristic obliquity-dependent injection from a heated solar wind thermal tail population. With several advanced computation techniques involving stochastic simulation and integration, the model obtains the particle intensity at any location in interplanetary space through the rigorous solution to the time-dependent 5D focus transport equation in the phase space that includes perpendicular diffusion. We apply the model to the 2011 November 3 CME event. The calculation results reproduce multispacecraft SEP observations at Earth and STEREO-B reasonably well without normalization of particle flux. The observations at STEREO-A can be reproduced by rescaling particle energy or modified energy dependence of particle diffusion coefficients. This circumsolar SEP event seen by spacecraft at Earth, STEREO-A, and STEREO-B at widely separated longitudes can be explained by diffusive shock acceleration by a single CME shock with a moderate speed.
The physical modeling of active regions (ARs) and of the global coronal is receiving increasing interest lately. Recent attempts to model ARs using static equilibrium models were quite successful in reproducing AR images of hot soft X-ray (SXR) loops. They however failed to predict the bright EUV warm loops permeating ARs: the synthetic images were dominated by intense footpoint emission. We demonstrate that this failure is due to the very weak dependence of loop temperature on loop length which cannot simultaneously account for both hot and warm loops in the same AR. We then consider time-dependent AR models based on nanoflare heating. We demonstrate that such models can simultaneously reproduce EUV and SXR loops in ARs. Moreover, they predict radial intensity variations consistent with the localized core and extended emissions in SXR and EUV AR observations respectively. We finally show how the AR morphology can be used as a gauge of the properties (duration, energy, spatial dependence, repetition time) of the impulsive heating.
Direct Normal Irradiance (DNI) is one of the most used quantities to quantify the magnitude of solar energy resource. The concept of DNI is often interpreted differently for ground measurements and solar forecasting by numerical weather prediction (NWP) models, leading to substantial bias during evaluation of DNI forecasts especially under cloudy-sky conditions. To eliminate the bias, we use the Fast All-sky Radiation Model for Solar applications with DNI (FARMS-DNI) to provide a physics-based solution of solar radiation in the circumsolar region. The FARMS-DNI is implemented in the Weather Research and Forecasting model with solar extensions (WRF-Solar) to forecast day-ahead DNI in the north America. By comparing with conventional predictions from WRF-Solar and satellite observations from the National Solar Radiation Data Base (NSRDB), we found significant improvements in our prediction of DNI.
The National Solar Radiation Database (NSRDB) has significantly evolved since the first release of the point source database in 1993. The NSRDB has been widely used by the solar energy industry to provide long-term time-series solar resource data for various project phases. The NSRDB represents the state of the art in the satellite-based estimation of solar resource information and uses a unique physics-based modeling approach that allows improvements in accuracy with the deployment of the next-generation geostationary satellites. This poster provides an overview of (1) the improved spatiotemporal resolution; (2) on-demand services and their applications; (3) future Improvements, such as a new direct normal irradiance model and new methods to gap-fill missing data using physics-guided machine learning; (4) data quality; and (5) data dissemination.
A theory is proposed to explain the generation mechanism of type II solar radio bursts. It is suggested that the shock wave formed at the leading edge of a coronal transient can accelerate electrons. Because of the nature of the acceleration process, the energized electrons can possess a 'hollow-beam' type distribution function. When the electron beam propagates along the ambient magnetic field to lower altitudes and attains larger pitch angles, a synchrotron-maser instability can set in. This instability leads to the amplification of unpolarized or weakly polarized radiation. The present discussion incorporates a model which describes the ambient magnetic field and background plasma by means of MHD simulation. The potential emission regions may be located approximately, according to the time-dependent MHD simulation. Since the average local plasma frequency in the source region can be evaluated from the MHD model, the frequent drift associated with the radiation may be estimated. The result seems to be in good agreement with that derived from observations.
The discrepancy between the observed and computed flux of solar neutrinos can be eliminated if the temperature gradient in the core of the model, and hence its central temperature, is reduced. This could in principle be accomplished by reducing the core opacity; alternatively, it has been proposed that a fraction of the energy transport in the solar core results from the motion of hypothetical particles (the so-called WIMPs). The resulting changes in solar structure have measurable effects on the frequencies of the solar 5-minute oscillations. Here, models are considered where the opacity in the solar core has been decreased in a manner which roughly simulates the effects of the WIMPs. The analysis of the models and their frequencies provides insight into the consequences of modifications to the physics of solar models. In particular, it is found that models with the observed neutrino flux are inconsistent with observations of low-degree solar oscillations.
This presentation provides and overview of the current status of the National Solar Radiation Database (NSRDB).
EUV spectra of three active region loops observed above the solar limb with the SO55 spectrometer on Skylab are analyzed. It is noted that the lengths, peak temperatures, and pressures of the loops are typical of the X-ray coronal loops to which static models have been applied. It is found that the physical parameters of the coronal loop plasma derived from EUV spectra and raster pictures are not well represented by the static models. Although the loops also contain a significant quantity of cool plasma, no physical reason is found to differentiate them from other active region loops of similar length, pressure, and temperature. Several line ratios in the loop spectrum suggest departures from ionization equilibrium caused by rapid cooling. The source of this cooling material is discussed with reference to several models of loop dynamics.
The paper presents the results from a multi-year effort to develop and validate image processing methods forselecting the best physical models based on solar image observations. The approach consists of selecting the physical models based on their agreement with coronal holes extracted from the images. Ultimately, the goal is to use physical models to predict geomagnetic storms. We decompose the problem into three subproblems: (i) coronal hole segmentation based on physical constraints, (ii) matching clusters of coronal holes between different maps, and (iii) physical map classification. For segmenting coronal holes, we develop a multi-modal method that uses segmentation maps from three different methods to initialize a level-set method that evolves the initial coronal hole segmentation to the magnetic boundary. Then, we introducea new method based on Linear Programming for matchingclusters of coronal holes. The final matching is then performedusing Random Forests. The methods were carefully validatedusing consensus maps derived from multiple readers, manualclustering, manual map classification, and method validation for50 maps. The proposed multi-modal segmentation method significantly outperformed SegNet, U-net, Henney-Harvey, and FCNby providing accurate boundary detection. Overall, the methodgave a 95.5% map classification accuracy.
Primitive solar nebula physical model considered under angular momentum conservation in collapsing fragment
The goals of this research are to provide scientific and technical capabilities in the areas of solar and plasma physics contained in research programs and instrumentation development relative to current rocket and balloon projects; to develop flight instrumentation design, flight hardware, and flight program objectives and participate in peer reviews as appropriate; and to participate in solar-terrestrial physics modeling studies and analysis of flight data and provide theoretical investigations as required by these studies.
The paper presents the g-mode oscillation for a set of modern solar models. Each solar model is based on a single modification or improvement to the physics of a reference solar model. Improvements were made to the nuclear reaction rates, the equation of state, the opacities, and the treatment of the atmosphere. The error in the predicted g-mode periods associated with the uncertainties in the model physics is predicted and the specific sensitivities of the g-mode periods and their period spacings to the different model structures are described. In addition, these models are compared to a sample of published observations. A remarkably good agreement is found between the 'best' solar model and the observations of Hill and Gu (1990).