Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “solar 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 703 records · Page 39

Enabling the Assimilation of CrIS Shortwave Infrared Observations in Global NWP at NOAA. Part II: OSEs and Results

The assimilation of data from hyperspectral infrared sounders in global data assimilation systems has historically been focused on observations in the longwave infrared (LWIR) region of the spectrum (650-1100 cm−1), despite the often concurrent availability of measurements from the shortwave infrared (SWIR) region of the spectrum (2150-2550 cm−1), because issues (like solar effects) have generally prevented the assimilation of SWIR observations. Recent advances in radiative transfer models have worked to address some of the previous challenges in simulating SWIR observations, and the assimilation of SWIR data (e.g. from potential future small-satellites) is now a feasible prospect. Still, a better understanding of how these observations perform in a data assimilation system and impact resulting analyses and NWP forecasts is necessary. In this study, the value of SWIR observations in global NWP is assessed by assimilating SWIR observations from the Cross-track Infrared Sounder (CrIS) in NOAA’s Global Data Assimilation System (GDAS). The methodologies used to enable the assimilation of these observations, including the implementation of a scene-dependent observation error and the enhancement of quality control procedures, are discussed, as are the results of Observing System Experiments (OSEs) conducted to evaluate the impact of assimilating SWIR observations on forecast skill. The overall results show that SWIR assimilation produces similar forecast impacts to LWIR assimilation. The ability to demonstrate that the assimilation or SWIR observations in NWP is a realistic prospect may help to shape future constellations of small-satellites to serve as a beneficial complement to the current constellation if hyperspectral IR sounders.

Erin Jones↗

High-Resolution WRF-Based Downscaling of Earth System Model Projections for Energy Applications across CONUS

Evaluating energy resources under future scenarios requires meteorological information that adequately resolves regional-scale variability and is suitable for regional energy system studies. Although Earth system model (ESM) outputs provide essential large-scale context, their coarse resolution and inherent biases limit direct use in energy system applications. This work presents a high-resolution dynamical downscaling framework using the Weather Research and Forecasting (WRF) model to generate energy-relevant regional fields for future scenarios over the contiguous United States (CONUS). The framework first identifies an optimal WRF configuration through sensitivity experiments, then evaluates raw and bias-corrected ESM initial and boundary conditions, with soil moisture and soil temperature bias correction implemented as an integral component of the bias-corrected ESM atmospheric forcing prior to WRF dynamical downscaling to improve land-atmosphere interactions. Simulations performed at 4-km resolution show that uncorrected ESM forcing leads to systematically dry and cold soil states, which propagate into elevated near-surface air temperature and solar irradiance biases, particularly during summer for the period 2000-2014. Incorporating bias-corrected atmospheric forcing together with soil state bias correction substantially reduces these errors and improves the representation of surface energy processes in WRF simulations. The results highlight the importance of bias-aware initialization strategies in high-resolution dynamical downscaling for future energy system analysis and planning.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Substorm warnings - An ISEE-3 real time data system

The use of solar wind measurements made by ISEE-3 in its halo orbit around the L1 libration point to predict the onsets of magnetospheric substorms and geomagnetic storms is discussed. Consideration is given to the limitations on the predictive ability of the satellite measurements set by the bulk solar wind velocity, the elliptical orbit of the satellite and the correlation lengths of the magnetic field and the solar wind plasma. The ISEE-3 real-time data system is presented, with attention given to the ground receiving stations, the NASCOM communications system, the Multisatellite Operations Control Center and Information Processing Division at the Goddard Space Flight Center, the link between Goddard and the NOAA Space Environmental Services Center, and the NOAA Space Environment Laboratory data acquisition and display data system, which includes displays allowing storm forecasts. It is noted that the entire system should be operational by March, 1980.

Tsurutani, B.↗

AVHRR-Based Polar Pathfinder Products: Evaluation, Enhancement and Transition to MODIS

The Advanced Very High Resolution Radiometer (AVHRR)-Based Polar Pathfinder (APP) products include calibrated AVHRR channel data, surface temperatures, albedo, satellite scan and solar geometries, and cloud mask, all composited into twice-per-day images, and daily averaged fields of sea ice motion, for regions poleward of 50 latitude. Our general goals under this grant: (1) Quantify the APP accuracy and sources of error by comparing Pathfinder products with field measurements; (2) Determine the consistency of mean fields and trends in comparison with longer time series of available station data and forecast model output; (3) Investigate the consistency of the products between the different AVHRR instruments over the 1982-present period of the NOAA program; and (4) Compare and annual cycle of the APP products with MODIS to establish a baseline for extending Pathfinder-type products into the new ESE period.

Fowler, Charles↗

Application of IRTAM to Support ISS Program Safety

The International Space Station (ISS) orbits near the F-peak of the ionosphere (approximately 400 km altitude). Generally, satellites orbiting at this altitude would have a floating potential (FP) of approximately -1 V due to the electron temperature (Te). However, the ISS has 8 large negatively grounded 160 V solar array wings (SAW) that collect a significant electron current from the ionosphere. This current drives the ISS FP much more negative during insolation and is highly dependent on the electron density (Ne). Also, due to the size of the ISS, magnetic inductance caused by the geomagnetic field produces a delta potential up to 40 V across the truss, possibly producing positive potentials. During Extravehicular Activity (EVA) the negative FP can lead to an arcing hazard when it exceeds -45.5 V, and the positive FP can produce a DC current high enough to stimulate the astronaut's muscles and also cause a hazard. Data collected from the Floating Potential Monitoring Unit (FPMU) have shown that the probability of either of these hazards occurring during times with quiet to moderately disturbed geomagnetic activity is low enough to no longer be considered a risk. However, a study of the ionosphere Ne during severe geomagnetic storm activity has shown that the Ne can be enhanced by a factor of 6 in the ISS orbit. As a result, the ISS Safety Review Panel (SRP) requires that ionospheric conditions be monitored using the FPMU in conjunction with the ISS Plasma Interaction Model (PIM) to determine if a severe geomagnetic storm could result in a plasma environment that could produce a hazard. A 'Real-Time' plasma hazard assessment process was developed to support ISS Program real-time decision making providing constraint relief information for EVAs planning and operations. This process incorporates 'real time' ionospheric conditions, ISS solar arrays' orientation, ISS flight attitude, and where the EVA will be performed on the ISS. This assessment requires real time data that is presently provided by the FPMU including ISS floating potential, along with ionospheric Ne and Te, in order to determine the present environment. Once the present environment conditions are known to be either above, below, or near the current IRI values, the IRI is used to forecast what the environment could become in the event of a severe geomagnetic storm. If the FPMU should fail, the Space Environments team needs another source of data which is utilized to support a short-term forecast for EVAs. The IRI Real-Time Assimilative Mapping (IRTAM) model is an ionospheric model that uses real time measurements from approximately 70 digisondes to produce foF2 and hmF2 global maps in 15 minute cadence. The Boeing Space Environments team has used the IRI coefficients produced in IRTAM to calculate the Ne along the ISS orbital track. The results of the IRTAM model have been compared to FPMU measurements and show excellent agreement (figure 1). IRTAM has been identified as a potential FPMU back-up system will be used as a backup for the FPMU to support the ISS Program following completion of an FPMU/IRTAM validation campaign.

Hartman, William A.↗

SPHINX: A Generalized Tool for Solar Energetic Particle Model Validation

The Integrated Solar Energetic Proton Alert/Warning (ISEP) project was established to aid in the transition of SEP models from the research realm to operational use. This project is a collaboration between the NASA Community Coordinated Model Center (CCMC), the NASA Moon to Mars Space Weather Analysis Office (M2M) and the NASA Space Radiation Analysis Group (SRAG). Through ISEP, the SEP Scoreboards (Intensity, Probability, and All Clear) were developed to visualize real time forecasts from SEP models for use in SRAG operations. To accomplish this, SRAG, CCMC, and M2M work closely with modelers on model development, CCMC works on Scoreboard implementation, and M2M provides human-in-the-loop space weather analysis and ensures that models run in real time.

space weather↗

Statistical theory for estimating sampling errors of regional radiation averages based on satellite measurements

The processes which determine the weather and climate are driven by the radiation received by the earth and the radiation subsequently emitted. A knowledge of the absorbed and emitted components of radiation is thus fundamental for the study of these processes. In connection with the desire to improve the quality of long-range forecasting, NASA is developing the Earth Radiation Budget Experiment (ERBE), consisting of a three-channel scanning radiometer and a package of nonscanning radiometers. A set of these instruments is to be flown on both the NOAA-F and NOAA-G spacecraft, in sun-synchronous orbits, and on an Earth Radiation Budget Satellite. The purpose of the scanning radiometer is to obtain measurements from which the average reflected solar radiant exitance and the average earth-emitted radiant exitance at a reference level can be established. The estimate of regional average exitance obtained will not exactly equal the true value of the regional average exitance, but will differ due to spatial sampling. A method is presented for evaluating this spatial sampling error.

Smith, G. L.↗

Electric Access System Enhancement (EASE): Assessment of a Distributed Energy Resource Management System for Enabling Dynamic Hosting Capacity

The Electric Access System Enhancement (EASE) project demonstrated a scalable, interoperable, and cost-effective means of integrating high penetration of distributed energy resources (DERs). The control architecture developed leveraged a Distributed Energy Resource Management System, Distribution System Operator for transacting energy, and 3rd Party DER Aggregator platform. The project identified ways to enhance the customer interconnection process to the grid and improve access to information from DERs and optimize the usage of DERs to provide energy services in a simulated day-ahead shadow market while maintain grid reliability. By integrating these capabilities into a scalable system of systems, the DSO can effectively balance DER generation and customer demand on the distribution network. This capability allows the grid to host more DER than traditionally possible on wires alone. Hosting more DER has the added benefit of supplying increased demand growth, sometimes beyond the capacity limits of the distribution network itself. This is known as a dynamic hosting capacity, which could help utilities manage the forecasted growth in electricity demand as California switch to electric vehicles and appliances. This could help to establish energy storage and PV generation as a “pseudo firm” generation resource mix if managed appropriately and provide sufficient resource adequacy for distribution capacity upgrade deferrals.

14 SOLAR ENERGY↗

GEWEX Water and Energy Budget Study

Closing the global water and energy budgets has been an elusive Global Energy and Water-cycle Experiment (GEWEX) goal. It has been difficult to gather many of the needed global water and energy variables and processes, although, because of GEWEX, we now have globally gridded observational estimates for precipitation and radiation and many other relevant variables such as clouds and aerosols. Still, constrained models are required to fill in many of the process and variable gaps. At least there are now several atmospheric reanalyses ranging from the early National Centers for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR) and NCEP/Department of Energy (DOE) reanalyses to the more recent ERA40 and JRA-25 reanalyses. Atmospheric constraints include requirements that the models state variables remain close to in situ observations or observed satellite radiances. This is usually done by making short-term forecasts from an analyzed initial state; these short-term forecasts provide the next guess, which is corrected by comparison to available observations. While this analysis procedure is likely to result in useful global descriptions of atmospheric temperature, wind and humidity, there is no guarantee that relevant hydroclimate processes like precipitation, which we can observe and evaluate, and evaporation over land, which we cannot, have similar verisimilitude. Alternatively, the Global Land Data Assimilation System (GLDAS), drives uncoupled land surface models with precipitation, surface solar radiation, and surface meteorology (from bias-corrected reanalyses during the study period) to simulate terrestrial states and surface fluxes. Further constraints are made when a tuned water balance model is used to characterize the global runoff observational estimates. We use this disparate mix of observational estimates, reanalyses, GLDAS and calibrated water balance simulations to try to characterize and close global and terrestrial atmospheric and surface water and energy budgets to within 10-20% for long term (1986-1995), large-scale global to regional annual means.

Roads, J.↗

Higher Throughput, Lower Cost Processing of Flexible Perovskite Solar Cells by Photonic Curing

The objectives of this SETO project are to explore whether photonic curing can be used to perform annealing in perovskite solar cell (PSC) fabrication, and if so, whether it can enable high-throughput PSC manufacturing. The successful outcome of the project will reduce manufacturing cost towards realizing the SETO goal of the Levelized cost of electricity (LCOE) $0.03/kWh by 2030. The project is a collaboration between Hsu’s group at the University of Texas at Dallas (UTD) and NovaCentrix, an Austin TX based company that pioneered the PulseForge® photonic curing tools and has expertise in applying this technology to printed electronics. Thin-film synthesis typically involves a thermal annealing step to convert from precursors to final material phases, to obtain the desired crystalline phase, or to improve materials structural or electrical properties. However, thermal annealing typically requires high temperatures that are incompatible with inexpensive plastic substrates and annealing times lasting tens-to-hundreds of minutes, making it incompatible with roll-to-roll (R2R) manufacturing at a desirable web speed. Photonic curing uses a flash Xe lamp to deliver short (20 µs to 100 ms) but high intensity (up to 50 kW/cm 2 ) pulses of broadband light (200 – 1500 nm) to the sample. It has been successfully applied to sinter printed metal nanoparticle inks into conductive patterns. This project aimed to extend the applications of photonic curing to convert halide perovskite films and transparent metal oxide transport layers. Despite the interruption by COVID-19, we have achieved significant accomplishments to establish photonic curing as a viable technology for the high-throughput manufacturing of PSCs. (1) We demonstrated that photonic curing can convert both metal oxide transport layer and halide perovskite absorber in the PSCs and successfully made flexible PSCs on Corning Willow glass® (WG)/ indium tin oxide (ITO) substrates without any thermal annealing steps. These photonically cured devices achieve performance comparable to conventional thermally annealed devices but were fabricated with a total processing time reduced by six orders of magnitude. We also show that photonic curing can achieve desirable web speed (up to 26 m/min) and large-area uniformity. (2) By establishing photonic curing processing phase spaces for the perovskite, we show that the longer pulses are more forgiving in small variation in the photonic curing conditions, making it more suitable for R2R manufacturing. (3) Our results clearly illustrated the importance of improving mechanical and optical properties of transparent electrodes on flexible substrates. (4) We also performed fundamental studies that elucidate the interactions between perovskite absorber and metal oxide transport layer that can undermine the PSC stability. Five refereed journal papers in top journals (ACS Energy Letters, npj Flexible Electronics, ACS Applied Energy Materials, Materials Advances, and Frontiers in Energy Research), one conference proceeding (PVSC-47), six conference presentations (PVSC and Materials Research Society meetings), and a Ph.D. dissertation (Trey B. Daunis) were produced under this project. Energy Materials Corporation (EMC) has forecast that with a 1.5 m wide web moving at a web speed of 30 m/min, R2R processing of PSC could provide up to 4 GW/year for a single R2R assembly line. Replacing thermal annealing with photonic curing, as demonstrated by the accomplishments of this project, is necessary to realize such an energy production goal.

14 SOLAR ENERGY↗

Carbon‐Neutral Pathways for the United States

Abstract The Intergovernmental Panel on Climate Change (IPCC) Special Report on Global Warming of 1.5°C points to the need for carbon neutrality by mid‐century. Achieving this in the United States in only 30 years will be challenging, and practical pathways detailing the technologies, infrastructure, costs, and tradeoffs involved are needed. Modeling the entire U.S. energy and industrial system with new analysis tools that capture synergies not represented in sector‐specific or integrated assessment models, we created multiple pathways to net zero and net negative CO 2 emissions by 2050. They met all forecast U.S. energy needs at a net cost of 0.2–1.2% of GDP in 2050, using only commercial or near‐commercial technologies, and requiring no early retirement of existing infrastructure. Pathways with constraints on consumer behavior, land use, biomass use, and technology choices (e.g., no nuclear) met the target but at higher cost. All pathways employed four basic strategies: energy efficiency, decarbonized electricity, electrification, and carbon capture. Least‐cost pathways were based on >80% wind and solar electricity plus thermal generation for reliability. A 100% renewable primary energy system was feasible but had higher cost and land use. We found multiple feasible options for supplying low‐carbon fuels for non‐electrifiable end uses in industry, freight, and aviation, which were not required in bulk until after 2035. In the next decade, the actions required in all pathways were similar: expand renewable capacity 3.5 fold, retire coal, maintain existing gas generating capacity, and increase electric vehicle and heat pump sales to >50% of market share. This study provides a playbook for carbon neutrality policy with concrete near‐term priorities.

Williams, James H.↗

Agent-Based Coordination Scheme for PV Integration (ABC4PV)

Renewables and especially photovoltaics (PV) have benefitted significantly from a host of incentives and policies targeted toward enhanced integration and adoption of specific energy technologies. However, with the push to move forward into a subsidy-free market framework, behind-the-meter residential PV applications have generally struggled to retain their value (unlike utility scale and commercial projects) [1]. This project focused on developing control-theoretic solutions aimed at improving the integration and interaction of behind-the-meter residential PV with other distribution system assets (controllable and non-controllable) to enhance the integrated value of residential PV. To this end, a suite of decentralized control methodologies have been developed to enable effective coordination and control of behind-the-meter residential load customers’ PV, battery storage systems (BSS), controllable loads and other similar assets within a distribution feeder. This interaction aims at procuring energy savings and, thus, energy bill savings. The main source of savings is drawn from reducing the effect of demand charge pricing and is realized at the feeder level, assuming community level interaction and management among the aforementioned assets. Optimal control of the assets is implemented with a distributed optimization methodology, leveraging consensus-based algorithms. The results gathered from the optimal control simulations demonstrates that the savings can be duly achieved and the algorithm decision times (to dynamically control asset set points, for example) are fast. As for the overall efficiency of PV+BSS systems, to procure energy savings from curtailment of the demand charge pricing effects, the optimal control is set up so as to minimize the variance of the load for all customers, throughout a feeder and throughout time in a rolling horizon scheduling with model predictive control. The control takes into account inter-temporal electrochemical storage (battery) degradation costs: specifically, we have developed a long-term lifetime model for the BSS that weighs in the effect of the degradation factor in the dispatch formulations, thus, a considerable operating cost that affects energy decision making. The levelized cost of energy (LCOE – redefined for the purpose of quantifying asset integration effectiveness through the customers’ energy cost) is shown to be below the threshold set for the combined PV+BSS topology of $ 0.14/kWh for multiple cases of PV penetration all the way up to 50%, provided that a policy of shared ownership of and savings is in place. Further, the LCOE calculated for the case before the deployment PV+BSS systems is also achievable, i.e. the deployment of PV+BSS, if planned and scheduled optimally. will have no effect on customers’ energy costs. From the control methodology viewpoint, the developed consensus-based algorithms are shown to converge for a wide range of problem cases (spanning normal operating scenarios and contingencies), guaranteeing dispatch solutions under forecasting errors, communication break-downs and cyber-security attacks. The proposed control solutions are scalable and real-time implementable, with dispatch computations and device set-point updates converging in less than 2s in most practical instances of the above events.

14 SOLAR ENERGY↗

Yarkovsky Drift Detections for 247 Near-Earth Asteroids

The Yarkovsky effect is a thermal process acting upon the orbits of small celestial bodies, which can cause these orbits to slowly expand or contract with time. The effect is subtle (〈da/dt〉∼10{sup −4} au My{sup −1} for a 1 km diameter object) and is thus generally difficult to measure. We analyzed both optical and radar astrometry for 600 Near-Earth Asteroids (NEAs) for the purpose of detecting and quantifying the Yarkovsky effect. We present 247 NEAs with measured drift rates, which is the largest published set of Yarkovsky detections. This large sample size provides an opportunity to examine the Yarkovsky effect in a statistical manner. In particular, we describe two independent population-based tests that verify the measurement of Yarkovsky orbital drift. First, we provide observational confirmation for the Yarkovsky effect’s theoretical size dependence of 1/D, where D is diameter. Second, we find that the observed ratio of negative to positive drift rates in our sample is 2.34, which, accounting for bias and sampling uncertainty, implies an actual ratio of 2.7{sub −0.7}{sup +0.3}. This ratio has a vanishingly small probability of occurring due to chance or statistical noise. The observed ratio of retrograde to prograde rotators is two times lower than the ratio expected from numerical predictions from NEA population studies and traditional assumptions about the sense of rotation of NEAs originating from various main belt escape routes. We also examine the efficiency with which solar energy is converted into orbital energy and find a median efficiency in our sample of 12%. We interpret this efficiency in terms of NEA spin and thermal properties.

79 ASTRONOMY AND ASTROPHYSICS↗

AVHRR-Based Polar Pathfinder Products: Evaluation, Enhancement, and Transition to MODIS

The AVHRR-Based Polar Pathfinder (APP) products include calibrated AVHRR channel data, surface temperatures, albedo, satellite scan and solar geometries, and a cloud mask composited into twice- per-day images, and daily averaged fields of sea ice motion, for regions poleward of 50 deg. latitude. Our goals under this grant, in general, are four-fold: 1. To quantify the APP accuracy and sources of error by comparing Pathfinder products with field measurements. 2. To determine the consistency of mean fields and trends in comparison with longer time series of available station data and forecast model output. 3. To investigate the consistency of the products between the different AVHRR instruments over the 1982-present period of the NOAA program. 4. To compare an annual cycle of the AVHRR Pathfinder products with MODIS to establish a baseline for extending Pathfinder-type products into the new ESE period. Year One tasks include intercomparisons of the Pathfinder products with field measurements, testing of algorithm assumptions, collection of field data, and further validation and possible improvement of the multi-sensor ice motion fields. Achievements for these tasks are summarized below.

Fowler, Charles↗

Bias Correcting NOAA's High-Resolution Rapid Refresh (HRRR) Wind Resource Data for Grid Integration Applications [Slides]

Many weather years of high-quality wind data are widely accepted in the grid integration community to be important for studying wind energy technical potential, energy system operations, and grid resilience. NREL makes high-quality wind and solar resource data available. NREL's Grid-Atmosphere workshop (March 2024) identified NREL National Solar Radiation Database as widely used in grid integration modeling, but there is less agreement on commonly used wind datasets. One important factor identified by ESIG's 2023 report 'Weather Dataset Needs for Planning and Analyzing Modern Power Systems' for gold standard wind data is regular updates. To address the need for regular updates, NREL's team can now process all currently available and regularly updated High-Resolution Rapid Refresh (HRRR) outputs. HRRR is an hourly-updated operational forecast product produced by the National Oceanic and Atmospheric Administration (NOAA) (Dowell et al., 2022). One barrier to NREL using HRRR is systematic bias and consistency with NREL's existing wind datasets (e.g. WIND Toolkit, 'WTK') across weather years. To address this barrier, we show that the HRRR can be interpolated and bias-corrected to be consistent with NRE's existing datasets. We call the new dataset BC-HRRR (bias-corrected HRRR). As with historical datasets like the WTK, BC-HRRR is intended for use in grid integration modeling (e.g., capacity expansion, production cost, and resource adequacy modeling). BC-HRRR's (2015-present) consistency with WTK (2007-2013) allows NREL to extend internal grid integration tooling with 15+ weather years of wind data with low-overhead extensibility to future years as they are made available by NOAA. The rest of this slide deck documents the BC-HRRR processing methods, validation, and its implications for intended use.

17 WIND ENERGY↗

Space Environment (Natural and Induced)

Considerable effort and improvement have been made in the study of ionizing radiation exposure occurring in various regions of space. Satellites and spacecrafts equipped with innovative instruments are continually refining particle data and providing more accurate information on the ionizing radiation environment. The major problem in accurate spectral definition of ionizing radiation appears to be the detailed energy spectra, especially at high energies, which is important parameter for accurate radiation risk assessment. Magnitude of risks posed by exposure to radiation in future space missions is subject to the accuracies of predictive forecast of event size of SPE, GCR environment, geomagnetic fields, and atmospheric radiation environment. Although heavy ion fragmentations and interactions are adequately resolved through laboratory study and model development, improvements in fragmentation cross sections for the light nuclei produced from HZE nuclei and their laboratory validation are still required to achieve the principal goal of planetary GCR simulation at a critical exposure site. More accurate prediction procedure for ionizing radiation environment can be made with a better understanding of the solar and space physics, fulfillment of required measurements for nuclear/atomic processes, and their validation and verification with spaceflights and heavy ion accelerators experiments. It is certainly true that the continued advancements in solar and space physics combining with physical measurements will strengthen the confidence of future manned exploration of solar system. Advancements in radiobiology will surely give the meaningful radiation hazard assessments for short and long term effects, by which appropriate and effective mitigation measures can be placed to ensure that humans safely live and work in the space, anywhere, anytime.

Kim, Myung-Hee Y.↗

Adaptive Numerical Algorithms in Space Weather Modeling

Space weather describes the various processes in the Sun-Earth system that present danger to human health and technology. The goal of space weather forecasting is to provide an opportunity to mitigate these negative effects. Physics-based space weather modeling is characterized by disparate temporal and spatial scales as well as by different physics in different domains. A multi-physics system can be modeled by a software framework comprising of several components. Each component corresponds to a physics domain, and each component is represented by one or more numerical models. The publicly available Space Weather Modeling Framework (SWMF) can execute and couple together several components distributed over a parallel machine in a flexible and efficient manner. The framework also allows resolving disparate spatial and temporal scales with independent spatial and temporal discretizations in the various models. Several of the computationally most expensive domains of the framework are modeled by the Block-Adaptive Tree Solar wind Roe Upwind Scheme (BATS-R-US) code that can solve various forms of the magnetohydrodynamics (MHD) equations, including Hall, semi-relativistic, multi-species and multi-fluid MHD, anisotropic pressure, radiative transport and heat conduction. Modeling disparate scales within BATS-R-US is achieved by a block-adaptive mesh both in Cartesian and generalized coordinates. Most recently we have created a new core for BATS-R-US: the Block-Adaptive Tree Library (BATL) that provides a general toolkit for creating, load balancing and message passing in a 1, 2 or 3 dimensional block-adaptive grid. We describe the algorithms of BATL and demonstrate its efficiency and scaling properties for various problems. BATS-R-US uses several time-integration schemes to address multiple time-scales: explicit time stepping with fixed or local time steps, partially steady-state evolution, point-implicit, semi-implicit, explicit/implicit, and fully implicit numerical schemes. Depending on the application, we find that different time stepping methods are optimal. Several of the time integration schemes exploit the block-based granularity of the grid structure. The framework and the adaptive algorithms enable physics based space weather modeling and even forecasting.

Toth, Gabor↗

Forecasting Coronal Mass Ejections from Vector Magnetograms

In a 17 vector magnetogram study of 12 bipolar active regions (Falconer, Moore, & Gary, 2002, ApJ in press), we correlated four quantitative global magnetic measures with the Coronal Mass Ejections (CME) productivity of the active region. The global measures included a measure of active region size, the total magnetic flux phi and three measures of an active region global nonpotentiality 1) the net current (I (sub N)), 2) the length of the strong-shear, strong-field main neutral line (L(sub SS)) and 3) and the normalized twist (alpha = muIN/PHI). The CME productivity was determined from YOHKOH/SXT observations, Geostationary Operational Environmental Satellite (GOES), and when possible Solar and Heliospheric Observatory/Large Angle and Spectrometric Coronagraph Experiment (SOHO/LASCO) observations within 12 days of the day of the magnetogram. We found that the three measures of global nonpotentiality (I(sub N), L(sub SS), alpha) were all well correlated (greater than 99% confidence level) with an active region's CME productivity. The sample size was to small to confirm if there was a statistical significant correlation of the globally nonscientist measures with future CME activity (i.e. from the date of the magnetogram forward). We are doubling our sample, and will report on the statistical significance of global nonpotentiality as a predictor of future CME productivity. The new active regions are all from the first year of the upgraded MSFC vector magnetograms. This work, is funded by NSF through the Space Weather Program, by NASA through the Living with the Star, Targeted Research and Technology, and by NASA Solar Physics Supporting Research and Technology Program. The upgrade to the MSFC vector magnetograph was supported by the High Energy Solar Spectroscopic Imager (HESSI) mission.

Falconer, D. A.↗