Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical Weather Prediction Model”

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 595 records · Page 33

A New Ensemble of Perturbed-Input-Parameter Simulations by the Community Atmosphere Model

Uncertainty quantification (UQ) is a fundamental challenge in the numerical simulation of Earth's weather and climate, and other complex systems. It entails much more than attaching defensible error bars to predictions: in particular it includes assessing low-probability but high-consequence events. To achieve these goals with models containing a large number of uncertain input parameters, structural uncertainties, etc., raw computational power is needed. An automated, self-adapting search of the possible model configurations is also useful. Our UQ initiative at the Lawrence Livermore National Laboratory has produced the most extensive set to date of simulations from the US Community Atmosphere Model. We are examining output from about 3,000 twelve-year climate simulations generated with a specialized UQ software framework, and assessing the model's accuracy as a function of 21 to 28 uncertain input parameter values. Most of the input parameters we vary are related to the boundary layer, clouds, and other sub-grid scale processes. Our simulations prescribe surface boundary conditions (sea surface temperatures and sea ice amounts) to match recent observations. Fully searching this 21+ dimensional space is impossible, but sensitivity and ranking algorithms can identify input parameters having relatively little effect on a variety of output fields, either individually or in nonlinear combination. Bayesian statistical constraints, employing a variety of climate observations as metrics, also seem promising. Observational constraints will be important in the next step of our project, which will compute sea surface temperatures and sea ice interactively, and will study climate change due to increasing atmospheric carbon dioxide.

54 ENVIRONMENTAL SCIENCES↗

Comparison of Mesospheric Winds From a High-Altitude Meteorological Analysis System and Meteor Radar Observations During the Boreal Winters of 2009-2010 and 2012-2013

We present a study of horizontal winds in the mesosphere and lower thermosphere (MLT) during the boreal winters of 2009-2010 and 2012-2013 produced with a new high-altitude numerical weather prediction (NWP) system. This system is based on a modified version of the Navy Global Environmental Model (NAVGEM) with an extended vertical domain up to approximately 116 km altitude coupled with a hybrid four-dimensional variational (4DVAR) data assimilation system that assimilates both standard operational meteorological observations in the troposphere and satellite-based observations of temperature, ozone and water vapor in the stratosphere and mesosphere. NAVGEM-based MLT analyzed winds are validated using independent meteor radar wind observations from nine different sites ranging from 69 deg N-67 deg S latitude. Time-averaged NAVGEM zonal and meridional wind profiles between 75 and 95 km altitude show good qualitative and quantitative agreement with corresponding meteor radar wind profiles. Wavelet analysis finds that the 3-hourly NAVGEM and 1-hourly radar winds both exhibit semi-diurnal, diurnal, and quasi-diurnal variations whose vertical profiles of amplitude and phase are also in good agreement. Wavelet analysis also reveals common time-frequency behavior in both NAVGEM and radar winds throughout the Northern extra tropics around the times of major stratospheric sudden warmings (SSWs) in January 2010 and January 2013, with a reduction in semi-diurnal amplitudes beginning around the time of a mesospheric wind reversal at 60 deg N that precedes the SSW, followed by an amplification of semi-diurnal amplitudes that peaks 10-14 days following the onset of the mesospheric wind reversal. The initial results presented in this study demonstrate that the wind analyses produced by the high altitude NAVGEM system accurately capture key features in the observed MLT winds during these two boreal winter periods.

McCormack, J.↗

The potential impact of scatterometry on oceanography - A wave forecasting case

A series of observing system simulation experiments have been performed in order to assess the potential impact of marine surface wind data on numerical weather prediction. In addition to conventional data, the experiments simulated the time-continuous assimilation of remotely sensed marine surface wind or temperature sounding data. The wind data were fabricated directly for model grid points intercepted by a Seasat-1 scatterometer swath and were assimilated into the lowest active level (945 mb) of the model using a localized successive correction method. It is shown that Seasat wind data can greatly improve numerical weather forecasts due to better definition of specific features. The case of the QE II storm is examined.

Cane, M. A.↗

Assimilating microwave cloudy observations into NASA GEOS model using a novel Bayesian Monte Carlo technique

Despite the importance of clouds and their influence on atmospheric water and energy balance, Numerical Weather Prediction (NWP) centers systematically exclude cloud information from the assimilation process and only assimilate clear-sky radiances (Janiskov´a et al. 2012). In order to ensure that only clear sky radiances are assimilated, strict cloud detection thresholds are applied before radiances are fed into data assimilation (DA) systems. This process not only excludes a large portion of satellite radiances, but causes loss of information in the regions that are of high interest to meteorologists and are most challenging for weather forecasts (Errico et al. 2007; Haddad et al. 2015). Although, in recent years there has been great advances in the operational weather forecasting, the prediction of tropical cyclones (TC), especially the intensity of TCs, remains challenging. According to Aksoy et al. (2013), in addition to the model deficiencies, another important factor that contributes to this challenge includes lack of observations in the peripheral environment (rain- bands) of TCs mainly because of the selective assimilation of existing observations. Satellite observations provide more than 90 % of the input data for the initialization of NWP models but more than 75 % of satellite observations are discarded due to the cloud contamination as well as land, snow, and ice emissivity issues (Bauer et al. 2010).

Rainband↗

Numerical Simulation of Aircraft Trailing Vortices

The increase in air traffic is currently outpacing the development of new airport runways. This is leading to greater air traffic congestion, resulting in costly delays and cancellations. The National Aeronautics and Space Administration (NASA) under its Terminal Area Productivity (TAP) program is investigating new technologies that will allow increased airport capacity while maintaining the present standards for safety. As an element of this program, the Aircraft Vortex Spacing System (AVOSS) is being demonstrated in July 2000, at Dallas Ft-Worth Airport. This system allows reduced aircraft separations, thus increasing the arrival and departure rates, while insuring that wake vortices from a leading aircraft do not endanger trailing aircraft. The system uses predictions or wake vortex position and strength based on input from the current weather state. This prediction is accomplished by a semi-empirical model developed from theory, field observations, and relationships derived from numerical wake vortex simulations. Numerical experiments with a Large Eddy Simulation (LES) model are being conducted in order to provide guidance for the enhancement of these prediction algorithms. The LES Simulations of wake vortices are carried out with NASA's Terminal Area Simulation System (TASS). Previous wake vortex investigations with TASS are described. The primary objective of these numerical studies has been to quantify vortex transport and decay in relation to atmospheric variables. This paper summarizes many of the previous investigations with the TASS model and presents some new results regarding the onset of wake vortex decay.

Proctor, Fred H.↗

Graphics Processing Unit (GPU) Acceleration of the Goddard Earth Observing System Atmospheric Model

The Goddard Earth Observing System 5 (GEOS-5) is the atmospheric model used by the Global Modeling and Assimilation Office (GMAO) for a variety of applications, from long-term climate prediction at relatively coarse resolution, to data assimilation and numerical weather prediction, to very high-resolution cloud-resolving simulations. GEOS-5 is being ported to a graphics processing unit (GPU) cluster at the NASA Center for Climate Simulation (NCCS). By utilizing GPU co-processor technology, we expect to increase the throughput of GEOS-5 by at least an order of magnitude, and accelerate the process of scientific exploration across all scales of global modeling, including: The large-scale, high-end application of non-hydrostatic, global, cloud-resolving modeling at 10- to I-kilometer (km) global resolutions Intermediate-resolution seasonal climate and weather prediction at 50- to 25-km on small clusters of GPUs Long-range, coarse-resolution climate modeling, enabled on a small box of GPUs for the individual researcher After being ported to the GPU cluster, the primary physics components and the dynamical core of GEOS-5 have demonstrated a potential speedup of 15-40 times over conventional processor cores. Performance improvements of this magnitude reduce the required scalability of 1-km, global, cloud-resolving models from an unfathomable 6 million cores to an attainable 200,000 GPU-enabled cores.

Putnam, Williama↗

Accelerating Radiation Computations for Dynamical Models With Targeted Machine Learning and Code Optimization

Abstract Atmospheric radiation is the main driver of weather and climate, yet due to a complicated absorption spectrum, the precise treatment of radiative transfer in numerical weather and climate models is computationally unfeasible. Radiation parameterizations need to maximize computational efficiency as well as accuracy, and for predicting the future climate many greenhouse gases need to be included. In this work, neural networks (NNs) were developed to replace the gas optics computations in a modern radiation scheme (RTE+RRTMGP) by using carefully constructed models and training data. The NNs, implemented in Fortran and utilizing BLAS for batched inference, are faster by a factor of 1–6, depending on the software and hardware platforms. We combined the accelerated gas optics with a refactored radiative transfer solver, resulting in clear‐sky longwave (shortwave) fluxes being 3.5 (1.8) faster to compute on an Intel platform. The accuracy, evaluated with benchmark line‐by‐line computations across a large range of atmospheric conditions, is very similar to the original scheme with errors in heating rates and top‐of‐atmosphere radiative forcings typically below 0.1 K day −1 and 0.5 W m −2 , respectively. These results show that targeted machine learning, code restructuring techniques, and the use of numerical libraries can yield material gains in efficiency while retaining accuracy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Tropical Atlantic wind field variations during Sequal Preliminary results

The Sequal (Seasonal Response of the Equatorial Atlantic) program has the objective to model correctly the dominant physical processes producing the seasonal redistribution of heat, mass, and salt in the upper ocean. The efforts of Sequal are greatly enhanced by the Focal (Francais Ocean et Climat dans l'Atlantique Equatorial) program. A knowledge of fields of wind stress over the tropical Atlantic is needed for the achievement of the Sequal objectives. The Sequal program makes use of the Fleet Numerical Oceanography Central (FNOC) 'Nogaps' wind fields as a first-guess field which will subsequently be modified. Nogaps is a global weather prediction model consisting of independent global wind and mass analysis submodels. The present investigation is concerned with some preliminary comparisons between Nogaps predictions and surface observations, taking into account monthly mean Nogaps fields and monthly means from a century-long ship report-based climatology.

Harrison, D. E.↗

Preventive Power Outage Estimation Based on a Novel Scenario Clustering Strategy

The increasing occurrence of extreme weather events is challenging power grid operation. For extreme weather events, the system operator is responsible for estimating the power outages and scheduling the restoration resources. This paper proposes an outage evaluation framework to identify the possible unserved load profiles, vulnerable areas, and mobile energy adequacy. The outputs of an outage prediction model tool are used to generate numerous faulted line scenarios. Next, each scenario's nodal unserved load profile is obtained by solving a three-phase restoration model that considers repair crews and mobile energy resources (MERs). Then, a novel scenario clustering strategy is developed to cluster the unserved load profiles into multiple representative profiles which the system operator can focus on. Finally, case studies on a distribution system evaluate the damage caused by an extreme weather event and verify the effectiveness of the proposed scenario clustering strategy.

MATHEMATICS AND COMPUTING,POWER TRANSMISSION AND D↗

GEFSv12 High- and Low-Skill Day-10 Tornado Forecasts

Abstract On average, modern numerical weather prediction forecasts for daily tornado frequency exhibit no skill beyond day 10. However, in this extended-range lead window, there are particular model cycles that have exceptionally high forecast skill for tornadoes because of their ability to correctly simulate the future synoptic pattern. Here, model initial conditions that produced a more skillful forecast for tornadoes over the United States were exploited while also highlighting potential causes for low-skill cycles within the Global Ensemble Forecasting System, version 12 (GEFSv12). There were 88 high-skill and 91 low-skill forecasts in which the verifying day-10 synoptic pattern for tornado conditions revealed a western U.S. thermal trough and an eastern U.S. thermal ridge, a favorable configuration for tornadic storm occurrence. Initial conditions for high skill forecasts tended to exhibit warmer sea surface temperatures throughout the tropical Pacific Ocean and Gulf of Mexico, an active Madden–Julian oscillation, and significant modulation of Earth-relative atmospheric angular momentum. Low-skill forecasts were often initialized during La Niña and negative Pacific decadal oscillation conditions. Significant atmospheric blocking over eastern Russia—in which the GEFSv12 overforecast the duration and characteristics of the downstream flow—was a common physical process associated with low-skill forecasts. This work helps to increase our understanding of the common causes of high- or low-skill extended-range tornado forecasts and could serve as a helpful tool for operational forecasters. Significance Statement This research provides a framework for the anticipation of a more (or less) skillful 10-day tornado forecast in an operational numerical weather prediction system. High-skill forecasts were associated with substantial tropical convection and warm sea surface temperature throughout the Pacific Ocean and Gulf of Mexico, whereas the underlying cause of low-skill forecasts were typically associated with a blocking anticyclone over eastern Russia. These findings are important because they permit increased or decreased confidence in a long-range forecast of tornado occurrence based on a dynamical prediction system.

Meteorology & Atmospheric Sciences↗

Complementarities of GEO and LEO IR Sounders for Numerical Weather Prediction in an OSSE Framework

Preparations are underway for the United States’ weather satellite program and the expected advancements in the coming decade. Among these changes is the launch of the proposed NOAA/NASA Geostationary eXtended Observations (GeoXO) Sounder (GXS), coordinated with international counterparts, to form a global ring of hyperspectral infrared (IR) sounders. Coincident with the progress of IR sounders from GEO platforms, the future of low Earth orbit (LEO) sounders is being re-envisioned after decades of beneficial impact in numerical weather prediction (NWP) systems. To address questions about the utility of the two platform types in concert, the Global Modeling and Assimilation Office (GMAO) Observing System Simulation Experiment (OSSE) framework was used to examine the roles of GEO and LEO sounders and their impact on forecast error reduction. Overall, inclusion of both types of platforms in the satellite program produces the largest forecast error reduction with LEO sounders having a strong impact on global skill and GEO sounders providing the most benefit on the scale of targeted regions.

Erica McGrath-Spangler↗

GEOS Atmospheric Model: Challenges at Exascale

The Goddard Earth Observing System (GEOS) model at NASA's Global Modeling and Assimilation Office (GMAO) is used to simulate the multi-scale variability of the Earth's weather and climate, and is used primarily to assimilate conventional and satellite-based observations for weather forecasting and reanalysis. In addition, assimilations coupled to an ocean model are used for longer-term forecasting (e.g., El Nino) on seasonal to interannual times-scales. The GMAO's research activities, including system development, focus on numerous time and space scales, as detailed on the GMAO website, where they are tabbed under five major themes: Weather Analysis and Prediction; Seasonal-Decadal Analysis and Prediction; Reanalysis; Global Mesoscale Modeling, and Observing System Science. A brief description of the GEOS systems can also be found at the GMAO website. GEOS executes as a collection of earth system components connected through the Earth System Modeling Framework (ESMF). The ESMF layer is supplemented with the MAPL (Modeling, Analysis, and Prediction Layer) software toolkit developed at the GMAO, which facilitates the organization of the computational components into a hierarchical architecture. GEOS systems run in parallel using a horizontal decomposition of the Earth's sphere into processing elements (PEs). Communication between PEs is primarily through a message passing framework, using the message passing interface (MPI), and through explicit use of node-level shared memory access via the SHMEM (Symmetric Hierarchical Memory access) protocol. Production GEOS weather prediction systems currently run at 12.5-kilometer horizontal resolution with 72 vertical levels decomposed into PEs associated with 5,400 MPI processes. Research GEOS systems run at resolutions as fine as 1.5 kilometers globally using as many as 30,000 MPI processes. Looking forward, these systems can be expected to see a 2 times increase in horizontal resolution every two to three years, as well as less frequent increases in vertical resolution. Coupling these resolution changes with increases in complexity, the computational demands on the GEOS production and research systems should easily increase 100-fold over the next five years. Currently, our 12.5 kilometer weather prediction system narrowly meets the time-to-solution demands of a near-real-time production system. Work is now in progress to take advantage of a hybrid MPI-OpenMP parallelism strategy, in an attempt to achieve a modest two-fold speed-up to accommodate an immediate demand due to increased scientific complexity and an increase in vertical resolution. Pursuing demands that require a 10- to 100-fold increases or more, however, would require a detailed exploration of the computational profile of GEOS, as well as targeted solutions using more advanced high-performance computing technologies. Increased computing demands of 100-fold will be required within five years based on anticipated changes in the GEOS production systems, increases of 1000-fold can be anticipated over the next ten years.

ESMF↗

Evapotranspiration from nonuniform surfaces - A first approach for short-term numerical weather prediction

Observations of surface heterogeneity of soil moisture from scales of meters to hundreds of kilometers are discussed, and a relationship between grid element size and soil moisture variability is presented. An evapotranspiration model is presented which accounts for the variability of soil moisture, standing surface water, and vegetation internal and stomatal resistance to moisture flow from the soil. The mean values and standard deviations of these parameters are required as input to the model. Tests of this model against field observations are reported, and extensive sensitivity tests are presented which explore the importance of including subgrid-scale variability in an evapotranspiration model.

Wetzel, Peter J.↗

A Case Study of the Impact of AIRS Temperature Retrievals on Numerical Weather Prediction

Large errors in numerical weather prediction are often associated with explosive cyclogenesis. Most studes focus on the under-forecasting error, i.e. cases of rapidly developing cyclones which are poorly predicted in numerical models. However, the over-forecasting error (i.e., to predict an explosively developing cyclone which does not occur in reality) is a very common error that severely impacts the forecasting skill of all models and may also present economic costs if associated with operational forecasting. Unnecessary precautions taken by marine activities can result in severe economic loss. Moreover, frequent occurrence of over-forecasting can undermine the reliance on operational weather forecasting. Therefore, it is important to understand and reduce the prdctions of extreme weather associated with explosive cyclones which do not actually develop. In this study we choose a very prominent case of over-forecasting error in the northwestern Pacific. A 960 hPa cyclone develops in less than 24 hour in the 5-day forecast, with a deepening rate of about 30 hPa in one day. The cyclone is not versed in the analyses and is thus a case of severe over-forecasting. By assimilating AIRS data, the error is largely eliminated. By following the propagation of the anomaly that generates the spurious cyclone, it is found that a small mid-tropospheric geopotential height negative anomaly over the northern part of the Indian subcontinent in the initial conditions, propagates westward, is amplified by orography, and generates a very intense jet streak in the subtropical jet stream, with consequent explosive cyclogenesis over the Pacific. The AIRS assimilation eliminates this anomaly that may have been caused by erroneous upper-air data, and represents the jet stream more correctly. The energy associated with the jet is distributed over a much broader area and as a consequence a multiple, but much more moderate cyclogenesis is observed.

Reale, O.↗

Exploring Flood Predictability in Taiwan through Coupled Atmospheric–Hydrological and High-Performance Hydrodynamic Models

Effective flood simulation capabilities can tremendously support early warning and disaster prevention. To examine the applicability of a fully physics-based and high-performance flood simulation and forecasting modeling framework for a flood-prone region in Taiwan, we conduct a numerical experiment that couples the Weather Research and Forecasting (WRF) Model, WRF-Hydrological modeling system (WRF-Hydro), and the Two-Dimensional Runoff Inundation Toolkit for Operational Needs (TRITON) to perform integrated rainfall, streamflow, and flood simulations. Furthermore, we first use the coupled WRF and WRF-Hydro (WWH) to predict rainfall and streamflow and then drive TRITON with the predicted streamflow hydrographs to simulate flood depth and inundation area. With the refined spatial resolution and parameterization, this framework can better predict rainfall with reasonable spatial patterns. Although WWH could overestimate the amount of rainfall in some areas, the uncertain rainfall–streamflow predictions produce reasonable flood maps able to pinpoint regions at risk of flooding. In terms of model efficiency, the graphics processing unit–based computation can yield a speed-up factor as high as ∼13 compared to the central processing unit–based computation, promoting the efficacy of the coupled modeling framework in practical real-time flood forecasting.

Coupled models↗

Development of the GEOS-5 Atmospheric General Circulation Model: Evolution from MERRA to MERRA2.

The Modern-Era Retrospective Analysis for Research and Applications-2 (MERRA2) version of the GEOS-5 (Goddard Earth Observing System Model - 5) Atmospheric General Circulation Model (AGCM) is currently in use in the NASA Global Modeling and Assimilation Office (GMAO) at a wide range of resolutions for a variety of applications. Details of the changes in parameterizations subsequent to the version in the original MERRA reanalysis are presented here. Results of a series of atmosphere-only sensitivity studies are shown to demonstrate changes in simulated climate associated with specific changes in physical parameterizations, and the impact of the newly implemented resolution-aware behavior on simulations at different resolutions is demonstrated. The GEOS-5 AGCM presented here is the model used as part of the GMAO's MERRA2 reanalysis, the global mesoscale "nature run", the real-time numerical weather prediction system, and for atmosphere-only, coupled ocean-atmosphere and coupled atmosphere-chemistry simulations. The seasonal mean climate of the MERRA2 version of the GEOS-5 AGCM represents a substantial improvement over the simulated climate of the MERRA version at all resolutions and for all applications. Fundamental improvements in simulated climate are associated with the increased re-evaporation of frozen precipitation and cloud condensate, resulting in a wetter atmosphere. Improvements in simulated climate are also shown to be attributable to changes in the background gravity wave drag, and to upgrades in the relationship between the ocean surface stress and the ocean roughness. The series of "resolution aware" parameters related to the moist physics were shown to result in improvements at higher resolutions, and result in AGCM simulations that exhibit seamless behavior across different resolutions and applications.

AGCM↗

Probabilistic Day-Ahead Forecasting Using an Analog Ensemble Approach for Wind Farm Grid Services

Wind resource assessment and wind power forecasting are used in research and industry to anticipate future power output at scales ranging from individual wind turbines to entire wind farms. Probabilistic day-ahead wind forecasting is useful for anticipating how a wind farm could potentially participate in the day-ahead market by providing upper and lower bounds for expected power generation, thus informing grid operators of its uncertainty. Understanding this uncertainty is part of a larger project focused on building a platform that combines efforts in weather forecasting, aerodynamic and economic modeling to create maximum value of a wind plant to better provide services to the grid. This effort is also known as the Atmosphere to Electrons to Grid (A2E2G) project. One method for producing a probabilistic forecast is through the analog ensemble approach (Delle Monache et al., 2011). This method leverages historical forecasts and their corresponding observations as a training data set from which future forecasts can be made. For some future forecast, the most similar historical forecasts (analogs) are identified on a regular time basis such as once per a 3-hour window. The most similar analogs, based on a metric such as root mean square error (RMSE), are recorded and their corresponding verifying observations are used as an ensemble member for this future forecast. Prior work in this area demonstrates improvements over raw Numerical Weather Prediction (NWP) forecasts and shows skill similar to techniques such as logistic regression and machine learning (Delle Monache et al., 2013; Alessandrini et al., 2015). Here, we take the High-Resolution Rapid Refresh model (HRRR) day-ahead forecast (0-36 hours) to create a probabilistic day-ahead forecast using an analog ensemble approach. The HRRR has an hourly temporal resolution, with a spatial resolution of 3 km. The 12 UTC HRRR model run is downloaded every day for one year from August 2019 - July 2020, with the first 11 months serving as a bank of analogs from which the forecasting algorithm can create a probabilistic forecast. Once downloaded, the original HRRR forecast is temporally interpolated to 5-minutes, aligning with both the temporal resolution of the observations as well as the timescale relevant for day-ahead power forecasts. The forecast is validated at the M2 tower at the Flatirons Campus of the National Renewable Energy Laboratory (NREL) at a typical wind turbine height of 80 m. Variables such as wind speed, wind direction, and turbulence intensity are incorporated into the probabilistic forecast model and weighted according to their relative importance to the forecast. Based on metrics such as mean bias error (MBE), mean absolute error (MAE), and root mean square error, the analog ensemble forecast outperforms the raw HRRR forecast during the testing period of July 2020. Figure 1 illustrates an example day-ahead forecast compared against the verifying observations. The general variability and ramps are captured throughout the day, with potential to further improve the analog ensemble model through machine learning techniques.

numerical weather prediction↗