Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “performance 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 1,711 records · Page 95

A Clustering-Based Scenario Generation Framework for Power Market Simulation with Wind Integration

A critical step in stochastic optimization models of power system analysis is to select a set of appropriate scenarios and significant numbers of scenario generation methods exist in the literature. This paper develops a clustering based scenario generation method, which aims to improve the performance of existing scenario generation techniques by grouping a set of correlated wind sites into clusters according to their cross-correlations. Copula based models are utilized to model spatiotemporal correlations and the Gibbs sampling is then used to generate scenarios for day-ahead markets. Our results show that the generated scenarios based on clustered wind sites outperform existing approaches in terms of reliability and sharpness and can reduce the total computational time for scenario generation and reduction significantly. The clustering-based framework can therefore provide a better support for real-world market simulations with high wind penetration.

data visualization↗

Uncertainty quantification and optimization of precipitating hydrometeor parameters for winter precipitation in a cloud microphysics scheme

The precipitating hydrometeor parameters used in cloud microphysics schemes carry inherent uncertainties. The quantification of these uncertainties, together with parameter optimization, can significantly improve precipitation forecasts. This study investigates the effects of 13 parameters in the Weather Research and Forecasting (WRF) Double-Moment 6-class (WDM6) microphysics scheme, which define the hydrometeor characteristics such as fall velocity–diameter and mass–diameter relationships, as well as the shape parameter of the drop size distribution for precipitating particles such as rain, snow, and graupel on simulated winter precipitation. A comparison between the model's pre-defined parameters and observations from the International Collaborative Experiments for the PyeongChang 2018 Olympic and Paralympic winter games (ICE-POP 2018) field campaign reveals that the fall velocity–diameter relationship for rain, the mass–diameter relationships for snow and graupel, and the shape parameters for all precipitating particles in the WDM6 scheme deviate from the median values observed by the two-dimensional video disdrometer (2DVD). To quantify parameter sensitivities, a perturbed parameter ensemble (PPE) of 256 simulations was conducted within parameter ranges constrained by 2DVD observations for three winter precipitation cases. Bayesian optimization was then applied to identify parameter sets that minimized the root mean square error (RMSE) for each case, achieving reductions of up to 30.2 %. These results demonstrate that ensemble-based uncertainty quantification and parameter optimization can help identify key parameters and provide a pathway to improving precipitation simulation performance. In addition, measurement sites can be strategically selected based on regions that show high sensitivity to variations in hydrometeor characteristic parameters.

Bayesian optimization↗

Divide and conquer: Learning chaotic dynamical systems with multistep penalty neural ordinary differential equations

Forecasting high-dimensional dynamical systems is a fundamental challenge in various fields, such as geosciences and engineering. Neural Ordinary Differential Equations (NODEs), which combine the power of neural networks and numerical solvers, have emerged as a promising algorithm for forecasting complex nonlinear dynamical systems. However, classical techniques used for NODE training are ineffective for learning chaotic dynamical systems. In this work, we propose a novel NODE-training approach that allows for robust learning of chaotic dynamical systems. Here, our method addresses the challenges of non-convexity and exploding gradients associated with underlying chaotic dynamics. Training data trajectories from such systems are split into multiple, non-overlapping time windows. In addition to the deviation from the training data, the optimization loss term further penalizes the discontinuities of the predicted trajectory between the time windows. The window size is selected based on the fastest Lyapunov time scale of the system. Multi-step penalty(MP) method is first demonstrated on Lorenz equation, to illustrate how it improves the loss landscape and thereby accelerates the optimization convergence. MP method can optimize chaotic systems in a manner similar to least-squares shadowing with significantly lower computational costs. Our proposed algorithm, denoted the Multistep Penalty NODE, is applied to chaotic systems such as the Kuramoto-Sivashinsky equation, the two-dimensional Kolmogorov flow, and ERA5 reanalysis data for the atmosphere. It is observed that MP-NODE provide viable performance for such chaotic systems, not only for short-term trajectory predictions but also for invariant statistics that are hallmarks of the chaotic nature of these dynamics.

Chaotic dynamical systems↗

Natural Hazard Forecast Alert Grid Risk System

Weather events cause most power outages. Often, we even get notifications on our phones to take cover or be prepared for an imminent event. If electric grid utilities had a similar warning that also included probable scenarios and the equipment involved, they could prepare and minimize the effects. Idaho National Laboratory had a project with the U.S. Department of Energy’s Cybersecurity, Energy Security, and Emergency Response program to develop a grid alert application that receives messages from the existing emergency alert system, filters and determines components possibly affected by the emergency event, calculates probable scenarios using MASTERRI (Modeling And Simulation for Targeted Reliability and Resilience Improvement). For high-risk events, the application can then send alert links to subscribed electric distribution utility operations staff to allow them to see and evaluate the scenarios and the impact in a web based interactive map tool. This proof of concept application used data from utilities and organizations, such as the international regulatory body North American Electric Reliability Corporation, which have complied historical failure data of elements that comprise the U.S. electric grid. Nominal failure rates are obtained from this data. To make this tool possible, estimated failure rates were calculated for different component types given the alert type, severity, and location. Historic weather-related grid element failures were correlated with historic weather events from the Integrated Public Alert & Warning System. These correlated events and failures are used along with Bayesian updates from the historical norms to provide a modified failure rate for grid elements in the alert areas and calculate probable scenarios. Working with an industry collaborator, actual grid models and data were used for demonstration cases. This report outlines the work performed for this project.

24 - POWER TRANSMISSION AND DISTRIBUTION↗

Inclusion of Linearized Moist Physics in Nasa's Goddard Earth Observing System Data Assimilation Tools

Inclusion of moist physics in the linearized version of a weather forecast model is beneficial in terms of variational data assimilation. Further, it improves the capability of important tools, such as adjoint-based observation impacts and sensitivity studies. A linearized version of the relaxed Arakawa-Schubert (RAS) convection scheme has been developed and tested in NASA's Goddard Earth Observing System data assimilation tools. A previous study of the RAS scheme showed it to exhibit reasonable linearity and stability. This motivates the development of a linearization of a near-exact version of the RAS scheme. Linearized large-scale condensation is included through simple conversion of supersaturation into precipitation. The linearization of moist physics is validated against the full nonlinear model for 6- and 24-h intervals, relevant to variational data assimilation and observation impacts, respectively. For a small number of profiles, sudden large growth in the perturbation trajectory is encountered. Efficient filtering of these profiles is achieved by diagnosis of steep gradients in a reduced version of the operator of the tangent linear model. With filtering turned on, the inclusion of linearized moist physics increases the correlation between the nonlinear perturbation trajectory and the linear approximation of the perturbation trajectory. A month-long observation impact experiment is performed and the effect of including moist physics on the impacts is discussed. Impacts from moist-sensitive instruments and channels are increased. The effect of including moist physics is examined for adjoint sensitivity studies. A case study examining an intensifying Northern Hemisphere Atlantic storm is presented. The results show a significant sensitivity with respect to moisture.

Convection↗

Multiple Beam Triode Driven RF Sources for Accelerator Applications Phase I Final Report

Calabazas Creek Research, Inc, (CCR), in collaboration with Microwave Power Products, Inc. (MPP), formerly Communications & Power Industries, LLC (CPI,) and JP Accelerator Works, Inc. (JPAW), embarked on a program to develop multiple beam triodes to produce RF power from 350 – 800 MHz with an average power exceeding 200 kW. The effort was motivated by the performance of a triode-based RF source which produces 25 kW of UHF power at 90% efficiency. The CCR effort focused on implementing this technology into a multiple beam device to increase the output power while retaining the low cost, compact size, and high efficiency. The program performed extensive simulations indicating that the goals could be achieved, and a prototype multiple beam triode was built, baked, and tested. Unfortunately, a grid to cathode short terminated the testing before the tube could generate RF power. Nevertheless, the effort demonstrated that a multiple beam triode could be designed, built, baked, and energized to high voltage. The multiple beam triode used oxide cathodes, which are only capable of pulsed operation. The multiple beam triode will be rebuilt using dispenser cathodes, which will allow high duty or continuous operation. The grids were also modified to be more robust to avoid previous issues. The MB triode will provide the beam power for RF generation. The RF is generated by surrounding the triode with input and output cavities to convert beam power to RF power. RF cavities to generate 200 kW CW at 350-450 MHz using the MB triode with dispenser cathodes was assembled during the program. The next Phase of this effort is to assemble the multiple beam triode using the subassemblies built in the Phase I program and test with the RF cavities. The Phase I program also initiated design of a higher frequency, higher power multiple beam triode. That design is forecast to produce approximately 500 kW CW from 350 - 500 MHz.

43 PARTICLE ACCELERATORS↗

Karhunen–Loève deep learning method for surrogate modeling and approximate Bayesian parameter estimation

We evaluate the performance of the Karhunen-Loève Deep Neural Network (KL-DNN) framework for surrogate modeling and approximate Bayesian parameter estimation in partial differential equation models. In the surrogate model, the Karhunen-Loève (KL) expansions are used for the dimensionality reduction of the number of unknown parameters and variables, and a deep neural network is employed to relate the reduced space of parameters to that of the state variables. The KL-DNN surrogate model is used to formulate a maximum-a-posteriori-like least-squares problem, which is randomized to draw samples of the posterior distribution of the parameters. We test the proposed framework for a hypothetical unconfined aquifer via comparison with the forward MODFLOW and inverse PEST++ iterative ensemble smoother (IES) solutions as well as the state-of-the-art Fourier neural operator (FNO) and deep operator networks (DeepONets) operator learning surrogate models. Our results show that the KL-DNN surrogate model outperforms FNO and DeepONet for forward predictions. For solving inverse problems, the randomized algorithm provides the same or more accurate Bayesian predictions of the parameters than IES as evidenced by the higher log-predictive probability of both the estimated parameter field and the forecast hydraulic head. The posterior mean obtained from the randomized algorithm is closer to the reference parameter field than that obtained with FNO as the maximum a posteriori estimate.

Approximate Bayesian inference↗

Porting the WAVEWATCH III (v6.07) wave action source terms to GPU

Abstract. Surface gravity waves play a critical role in several processes, including mixing, coastal inundation, and surface fluxes. Despite the growing literature on the importance of ocean surface waves, wind–wave processes have traditionally been excluded from Earth system models (ESMs) due to the high computational costs of running spectral wave models. The development of the Next Generation Ocean Model for the DOE’s (Department of Energy) E3SM (Energy Exascale Earth System Model) Project partly focuses on the inclusion of a wave model, WAVEWATCH III (WW3), into E3SM. WW3, which was originally developed for operational wave forecasting, needs to be computationally less expensive before it can be integrated into ESMs. To accomplish this, we take advantage of heterogeneous architectures at DOE leadership computing facilities and the increasing computing power of general-purpose graphics processing units (GPUs). This paper identifies the wave action source terms, W3SRCEMD, as the most computationally intensive module in WW3 and then accelerates them via GPU. Our experiments on two computing platforms, Kodiak (P100 GPU and Intel(R) Xeon(R) central processing unit, CPU, E5-2695 v4) and Summit (V100 GPU and IBM POWER9 CPU) show respective average speedups of 2× and 4× when mapping one Message Passing Interface (MPI) per GPU. An average speedup of 1.4× was achieved using all 42 CPU cores and 6 GPUs on a Summit node (with 7 MPI ranks per GPU). However, the GPU speedup over the 42 CPU cores remains relatively unchanged (∼ 1.3×) even when using 4 MPI ranks per GPU (24 ranks in total) and 3 MPI ranks per GPU (18 ranks in total). This corresponds to a 35 %–40 % decrease in both simulation time and usage of resources. Due to too many local scalars and arrays in the W3SRCEMD subroutine and the huge WW3 memory requirement, GPU performance is currently limited by the data transfer bandwidth between the CPU and the GPU. Ideally, OpenACC routine directives could be used to further improve performance. However, W3SRCEMD would require significant code refactoring to make this possible. We also discuss how the trade-off between the occupancy, register, and latency affects the GPU performance of WW3.

58 GEOSCIENCES↗

Ice storage model-predictive control in an office building with PV: scenario, error and sensitivity analysis

Thermal energy storage (TES) can enable more building-sited renewable electricity generation and lower utility bill costs for buildings owners and occupants, especially when there are high demand and variable time-of-use (TOU) charges. A model predictive control (MPC) strategy can offer additional savings over a schedule-based control with added complexity and reliance on forecasts. Here, this study examines savings for medium office buildings with chiller plants in three locations with building-installed solar photovoltaics (PV) to understand the impact of MPC. Control setpoints are fixed by a schedule-based control or optimized by nonlinear MPC. These control setpoints are actuated within EnergyPlus building models to simulate the utility cost of the chiller plant. NLP solutions can be unstable or unrealistic, but our results show that by regularizing the NLP, the solutions can be reasonably followed by the building model. MPC models make simplifications that lead to errors once the controller is participating in and changing the operation of the building. These errors average 9 % across the cases, showing that the most important parts of the system are represented. The no-thermal load costs are computed to show that the optimization can in some cases achieve both the minimum TOU and minimum monthly demand costs by demand management while reducing TOU energy costs by energy arbitrage. The MPC saves 35–66 % in the annual chiller plant operating costs, which is an additional savings above the schedule by 1–33 %. PV and TES are complementary and mostly independent, but a load with PV often results in better performance for the schedule. Our case study and sensitivity analysis show the importance of modeling and optimization for complex rates, but also the circumstances wherein a simpler strategy achieves the same performance with less potential for error.

14 SOLAR ENERGY↗

Controller Strategies for Managing Air Traffic in High Altitude Arrival Sectors

Substantial increases in the volume of air traffic in the National Airspace System (NAS) are forecast for the next decade, with the number of passengers travelling on U.S. airlines expected to increase by as much as 60%. This increased demand on system capacity will be accompanied by increases in traffic complexity as air traffic service providers routinely accommodate user preferred routing requests. Changes to the NAS to meet these new demands are currently underway, including development of new decision support tools to aid controllers in monitoring and managing air traffic, and increased air-to-air and air-to-ground information exchange. Changes in roles and responsibilities of pilots and controllers in flight path management will accompany these changes in traffic patterns and information technology, however the ultimate responsibility for maintaining aircraft separation will remain with the air traffic controller. A thorough understanding of the methods controllers use to manage air traffic will help ensure that changes to the NAS are implemented in a way that maintains the controller's ability to separate aircraft as the system evolves. This presentation describes the strategies controllers use today to manage arrival traffic in its descent from cruise altitude to the Terminal Radar Approach Control (TRACON) boundary. Factors that increase the complexity of this task include the presence of overflight traffic, varying aircraft performance characteristics, winds aloft, ground speed variations with altitude, the need to merge arrival traffic into a single stream, and, when arrival traffic exceeds airport runway capacity, the added task of metering flow into the TRACON. Because of the limited information available to controllers to manage arrival traffic, their strategies are often driven by the need to reduce the task's complexity, which can result in de-optimized flight paths for individual aircraft (e.g., sub-optimal descent or speed profiles). Understanding these strategies and the cognitive demands that drive them will support a safe transition to a NAS that relies on enhanced technologies. In addition, it could enable system developers to identify opportunities for new automation-based procedures or information displays that could reduce the controller's workload and increase operational efficiency.

Smith, Nancy↗

Utilization of Unmanned Aircraft Systems for Environmental Purposes at the Savannah River Site – 26578

Born in the 1950s, the Savannah River Plant was constructed as a plutonium and tritium production plant. As the Cold War era came and went, a change of direction was signaled when the name of the facility changed from Savannah River Plant to Savannah River Site (1989) – the main mission at Savannah River shifted from “production” to “cleanup.” The Department of Energy is responsible for managing the 310 square miles of land that is the Savannah River Site and the cleanup/remediation activities that occur. The Savannah River National Laboratory utilizes unmanned aircraft systems to assist with some of those environmental monitoring and remediation activities. One function of unmanned aircraft systems at the Savannah River National Laboratory is conducting aerial photography and videography. Various infrastructure evaluations have been performed with unmanned aircraft – facility rooftop inspections, water tower lock verifications, earthen dam surveys, etc. An unmanned aircraft system has been used for progress footage of remediation projects – Dunbarton Bay remediation, Saltstone Disposal Unit construction, etc. Thermal cameras on an unmanned aircraft system are used to help identify cracks/leaks in structures from vantage points not readily accessible by personnel. Virtual tours of different waste units are conducted with unmanned aircraft systems for Site groups and federal/state regulators to save costs and travel time. Additionally, advanced remote sensing technologies are used on unmanned aircraft systems at the Savannah River Site. Light Detection and Ranging uses laser pulses to measure distances to the Earth's surface or other targets to create highly detailed topographic maps that accurately portray terrain (e.g., elevation changes, slopes, and contours). Data collection with a Light Detection and Ranging unit mounted on an unmanned aircraft system is quick and efficient – large-area surveys are conducted in shorter time frames. Because Light Detection and Ranging can penetrate through foliage and vegetation to ground level, it is being used in conjunction with a watershed model to study the General Separations Area of the Site to determine stormwater flow volume and direction, monitor stream stability, and forecast river flow behavior. A hyperspectral camera captures a wide spectrum of electromagnetic wavelengths across numerous narrow bands, beyond what the human eye can see. It provides detailed spectral information about the objects or surfaces in an image, which can be used to identify and analyze materials based on their spectral signatures. A hyperspectral camera mounted on an unmanned aircraft system has been used at the Savannah River Site for non-destructive evaluation of concrete and concrete structures. Unmanned aircraft systems at the Savannah River Site are also used for the application of herbicide on identified rooftops. P-Reactor and R-Reactor are entombed nuclear reactors at the Savannah River Site. Unwanted vegetation is often present on these rooftops – over time roots can penetrate seams, cracks, and joints of the general roof and concrete caps, leading to water ingress and deterioration of the structural components. For time and cost efficiencies as well as in the interest of personnel safety, an unmanned aircraft is used to dispense herbicide on rooftop areas where vegetation is noticed to help alleviate the issues/hazards.

Lorier, Troy [Savannah River National Laboratory (↗

Intelligent Observation Strategies for Geosynchronous Remote Sensing for Natural Hazards

Geosynchronous satellites offer a unique perspective for monitoring environmental factors important to understanding natural hazards and supporting the disasters management life cycle, namely forecast, detection, response, recovery and mitigation. In the NASA decadal survey for Earth science, the GEO-CAPE mission was proposed to address coastal and air pollution events in geosynchronous orbit, complementing similar initiatives in Asia by the South Koreans and by ESA in Europe, thereby covering the northern hemisphere. In addition to analyzing the challenges of identifying instrument capabilities to meet the science requirements, and the implications of hosting the instrument payloads on commercial geosynchronous satellites, the GEO-CAPE mission design team conducted a short study to explore strategies to optimize the science return for the coastal imaging instrument. The study focused on intelligent scheduling strategies that took into account cloud avoidance techniques as well as onboard processing methods to reduce the data storage and transmission loads. This paper expands the findings of that study to address the use of intelligent scheduling techniques and near-real time data product acquisition of both the coastal water and air pollution events. The topics include the use of onboard processing to refine and execute schedules, to detect cloud contamination in observations, and to reduce data handling operations. Analysis of state of the art flight computing capabilities will be presented, along with an assessment of cloud detection algorithms and their performance characteristics. Tools developed to illustrate operational concepts will be described, including their applicability to environmental monitoring domains with an eye to the future. In the geostationary configuration, the payload becomes a networked thing with enough connectivity to exchange data seamlessly with users. This allows the full field of view to be sensed at very high rate under the control of ground infrastructure, resulting in improved efficiencies, accuracy and science benefits. Hence a remote sensing payload and its data may become one of millions of connected objects in the emerging Internet of Things (IoT), and be as easily accessible by a users smart phone as any other smart appliance.

Satellite Soil Moisture Data Assimilation Impacts on Modeling Weather Variables and Ozone in the Southeastern US – Part 2: Sensitivity to Dry-Deposition Parameterizations

Ozone (O3) dry deposition is a major O3 sink. As a follow-up study of Huang et al. (2021), we quantify the impact of satellite soil moisture (SM) on model representations of this process when different dry-deposition parameterizations are implemented, based on which the implications for interpreting O3 air pollution levels and assessing the O3 impacts on human and ecosystem health are provided. The SM data from NASA's Soil Moisture Active Passive mission are assimilated into the Noah-Multiparameterization (Noah-MP) land surface model within the NASA Land Information System framework, semicoupled with Weather Research and Forecasting model with online Chemistry (WRF-Chem) regional-scale simulations covering the southeastern US. Major changes in the modeling system used include enabling the dynamic vegetation option, adding the irrigation process, and updating the scheme for the surface exchange coefficient. Two dry-deposition schemes are implemented, i.e., the Wesely scheme and a “dynamic” scheme, in the latter of which dry-deposition parameterization is coupled with photosynthesis and vegetation dynamics. It is demonstrated that, when the dynamic scheme is applied, the simulated O3 dry-deposition velocities vd and their stomatal and cuticular portions, as well as the total O3 fluxes Ft, are larger overall; vd and Ft are 2–3 times more sensitive to the SM changes due to the data assimilation (DA). Further, through case studies at two forested sites with different soil types and hydrological regimes, we highlight that, applying the Community Land Model type of SM factor controlling stomatal resistance (i.e., β factor) scheme in replacement of the Noah-type β factor scheme reduced the v(d) sensitivity to SM changes by ∼75 % at one site, while it doubled this sensitivity at the other site. Referring to multiple evaluation datasets, which may be associated with variable extents of uncertainty, the model performance of vegetation, surface fluxes, weather, and surface O3 concentrations shows mixed responses to the DA, some of which display land cover dependency. Finally, using model-derived concentration- and flux-based policy-relevant O3 metrics as well as their matching exposure–response functions, the relative biomass/crop yield losses for several types of vegetation/crops are estimated to be within a wide range of 1 %–17 %. Their sensitivities to the model's dry-deposition scheme and the implementation of SM DA are discussed.

Min Huang↗

Model predictive control of heating, ventilation, and air conditioning (HVAC) systems: A state-of-the-art review

Due to the fast advancement of communication and information technology, intelligent buildings have garnered great interest. These buildings can forecast weather, ambient temperature, and sun irradiation and can modify heating, ventilation, and air conditioning (HVAC) operations appropriately, based on current and previous data. This change is intended to reduce HVAC system energy usage while maintaining an appropriate degree of thermal comfort and indoor air quality. Since its inception, model predictive control (MPC) has been one of the prospective solutions for HVAC management systems to reduce both costs and energy usage. Additionally, MPC is becoming increasingly practical as the processing capacity of building automation systems increases and a large quantity of monitored building data becomes available. MPC also provides the potential to improve the energy efficiency of HVAC systems via its capacity to consider limitations, to predict disruptions, and to factor in multiple competing goals such as interior thermal comfort and building energy consumption. Although substantial research has been conducted on MPC in building HVAC systems, there is a shortage of critical reviews and a lack of a comprehensive framework that formulates and defines the applications. Here, this article provides a comprehensive state-of-the-art overview of MPC in HVAC systems. Detailed discussions of modeling approaches and optimization algorithms are included. Numerous design aspects such as prediction horizon, occupancy behavior, building type, and cost function, that impact MPC performance are discussed in detail. The technical characteristics, advantages, and disadvantages of various types of modeling software are discussed. The primary objective of this work is to highlight critical design characteristics for the MPC control scheme and to give improved suggestions for future research. Moreover, numerous prospective scenarios have been suggested that might provide future research direction.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Extracting and Generating PV Soiling Profiles for Analysis, Forecasting, and Cleaning Optimization

The identification and prediction of the daily soiling profiles of a photovoltaic site is essential to plan the optimal cleaning schedule. In this article, we analyze and propose various methods to extract and generate photovoltaic soiling profiles, in order to improve the analysis and the forecast of the losses. New soiling rate extraction methods are proposed to reflect the seasonal variability of the soiling rates and, for this reason, are found to identify the most convenient cleaning day with the highest accuracy for the investigated sites. Also, we present an approach that could be used to predict future soiling losses through the implementation of stochastic weather generation algorithms whose ability to identify in advance the best cleaning schedule is also successfully tested. The methods presented in this article can optimize the operation and maintenance schedule and could make it possible, in the future, to predict soiling losses through analysis based only on environmental parameters, such as rainfall and particulate matter, without the need of long-term soiling data.

14 SOLAR ENERGY↗

Evaluation of Mixed-Phase Microphysics Within Winter Storms Using Field Data and In Situ Observations

Snow prediction within models is sensitive to the snow densities, habits, and degree of riming within the BMPs. Improving these BMPs is a crucial step toward improving both weather forecasting and climate predictions. Several microphysical schemes in the Weather Research and Forecasting (WRF) model down to 1.33-km grid spacing are evaluated using aircraft, radar, and ground in situ data from the Global Precipitation Mission Cold-season Precipitation Experiment (GCPEx) experiment over southern Ontario, as well as a few years (12 winter storms) of surface measurements of riming, crystal habit, snow density, and radar measurements at Stony Brook, NY (SBNY on north shore of Long Island) during the 2009-2012 winter seasons. Surface microphysical measurements at SBNY were taken every 15 to 30 minutes using a stereo microscope and camera, and snow depth and snow density were also recorded. During these storms, a vertically-pointing Ku band radar was used to observe the vertical evolution of reflectivity and Doppler vertical velocities. The GCPex presentation will focus on verification using aircraft spirals through warm frontal snow band event on 18 February 2012. All the BMPs realistically simulated the structure of the band and the vertical distribution of snow/ice aloft, except the SBU-YLIN overpredicted slightly and Thompson (THOM) underpredicted somewhat. The Morrison (MORR) scheme produced the best slope size distribution for snow, while the Stony Brook (SBU) underpredicted and the THOM slightly overpredicted. Those schemes that have the slope intercept a function of temperature (SBU and WSM6) tended to perform better for that parameter than others, especially the fixed intercept in Goddard. Overall, the spread among BMPs was smaller than in other studies, likely because there was limited riming with the band. For the 15 cases at SBNY, which include moderate and heavy riming events, the non-spherical snow assumption (THOM and SBU-YLIN) simulated a more realistic distribution of reflectivity than spherical snow assumptions in the WSM6 and MORR schemes. The MORR, WSM6, and SBU schemes are comparable to the observed velocity distribution in light and moderate riming periods. The THOM is approx. 0.25 m/s too slow with its velocity distribution in these periods. In heavier riming, the vertical Doppler velocities in the WSM6, THOM, and MORR schemes were approx. 0.25 m/s too slow, while the SBU was 0.25 to 0.5 m/s too fast because of some excessive cloud water issues.

Colle, Brian A.↗

Parameter Dependencies of Early-Stage Tangential Discontinuity-Driven Foreshock Bubbles in Local Hybrid Simulations

Foreshock bubbles (FBs) are significant foreshock transients that can accelerate particles and disturb the magnetosphere-ionosphere system. In the kinetic formation model, foreshock ions interact with the discontinuity by performing partial gyrations to generate currents that change the magnetic field topology around the discontinuity. However, how different foreshock ion properties affect the growth of the field variations is not well understood. Therefore, we use 2-D local hybrid simulations to study the effects of different foreshock ion distributions and properties on the growth of tangential discontinuity (TD)-driven FBs. We discover that for a gyrophase-bunched distribution with an initial phase where the guiding center is on the other side of the TD, the foreshock ions gyrate together across the TD, causing more foreshock ions to cross the TD and leading to a faster expansion of the structure than for a Maxwellian distribution. A ring distribution also yields higher expansion speeds because of the higher projected velocity into the new perpendicular direction. For Maxwellian distributions, there are positive and linear correlations of the FB expansion speeds with the initial foreshock ion densities, thermal speeds, parallel speeds, and sine of the TD magnetic shear angles. These parameter dependencies grow in strength as the structures evolve with time. The foreshock ion distributions and properties that lead to stronger currents produce more significant magnetic field variations and higher expansion speeds. Our study helps quantify the formation and expansion of FBs to forecast their space weather effects and contribution to shock acceleration.

Andrew Vu↗

Validation for Global Solar Wind Prediction Using Ulysses Comparison: Multiple Coronal and Heliospheric Models Installed at the Community Coordinated Modeling Center

The prediction of the background global solar wind is a necessary part of space weather forecasting. Several coronal and heliospheric models have been installed and/or recently upgraded at the Community Coordinated Modeling Center (CCMC), including the Wang-Sheely-Arge (WSA)-Enlil model, MHD-Around-a-Sphere (MAS)-Enlil model, Space Weather Modeling Framework (SWMF), and Heliospheric tomography using interplanetary scintillation data. Ulysses recorded the last fast latitudinal scan from southern to northern poles in 2007. By comparing the modeling results with Ulysses observations over seven Carrington rotations, we have extended our third-party validation from the previous near-Earth solar wind to middle to high latitudes, in the same late declining phase of solar cycle 23. Besides visual comparison, wehave quantitatively assessed the models capabilities in reproducing the time series, statistics, and latitudinal variations of solar wind parameters for a specific range of model parameter settings, inputs, and grid configurations available at CCMC. The WSA-Enlil model results vary with three different magnetogram inputs.The MAS-Enlil model captures the solar wind parameters well, despite its underestimation of the speed at middle to high latitudes. The new version of SWMF misses many solar wind variations probably because it uses lower grid resolution than other models. The interplanetary scintillation-tomography cannot capture the latitudinal variations of solar wind well yet. Because the model performance varies with parameter settings which are optimized for different epochs or flow states, the performance metric study provided here can serve as a template that researchers can use to validate the models for the time periods and conditions of interest to them.

Jian, L. K.↗