SEARCH · Search NASA

Results for “numerical weather 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.

SourceResearch connectionFrom yearThrough year

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

At least 361 records · Page 20

All sources NASA NTRS DOE OSTI

NASA NTRS1984

Simulation studies of the impact of advanced observing systems on numerical weather prediction

To study the potential impact of advanced passive sounders and lidar temperature, pressure, humidity, and wind observing systems on large-scale numerical weather prediction, a series of realistic simulation studies between the European Center for medium-range weather forecasts, the National Meteorological Center, and the Goddard Laboratory for Atmospheric Sciences is conducted. The project attempts to avoid the unrealistic character of earlier simulation studies. The previous simulation studies and real-data impact tests are reviewed and the design of the current simulation system is described. Consideration is given to the simulation of observations of space-based sounding systems.

Atlas, R.↗

NASA NTRS2010

Improving the Representation of Snow Crystal Properties with a Single-Moment Mircophysics Scheme

Single-moment microphysics schemes are utilized in an increasing number of applications and are widely available within numerical modeling packages, often executed in near real-time to aid in the issuance of weather forecasts and advisories. In order to simulate cloud microphysical and precipitation processes, a number of assumptions are made within these schemes. Snow crystals are often assumed to be spherical and of uniform density, and their size distribution intercept may be fixed to simplify calculation of the remaining parameters. Recently, the Canadian CloudSat/CALIPSO Validation Project (C3VP) provided aircraft observations of snow crystal size distributions and environmental state variables, sampling widespread snowfall associated with a passing extratropical cyclone on 22 January 2007. Aircraft instrumentation was supplemented by comparable surface estimations and sampling by two radars: the C-band, dual-polarimetric radar in King City, Ontario and the NASA CloudSat 94 GHz Cloud Profiling Radar. As radar systems respond to both hydrometeor mass and size distribution, they provide value when assessing the accuracy of cloud characteristics as simulated by a forecast model. However, simulation of the 94 GHz radar signal requires special attention, as radar backscatter is sensitive to the assumed crystal shape. Observations obtained during the 22 January 2007 event are used to validate assumptions of density and size distribution within the NASA Goddard six-class single-moment microphysics scheme. Two high resolution forecasts are performed on a 9-3-1 km grid, with C3VP-based alternative parameterizations incorporated and examined for improvement. In order to apply the CloudSat 94 GHz radar to model validation, the single scattering characteristics of various crystal types are used and demonstrate that the assumption of Mie spheres is insufficient for representing CloudSat reflectivity derived from winter precipitation. Furthermore, snow density and size distribution characteristics are allowed to vary with height, based upon direct aircraft estimates obtained from C3VP data. These combinations improve the representation of modeled clouds versus their radar-observed counterparts, based on profiles and vertical distributions of reflectivity. These meteorological events are commonplace within the mid-latitude cold season and present a challenge to operational forecasters. This study focuses on one event, likely representative of others during the winter season, and aims to improve the representation of snow for use in future operational forecasts.

Molthan, Andrew L.↗

DOE OSTI2022

“Gray Zone” Simulations Using a Three-Dimensional Planetary Boundary Layer Parameterization in the Weather Research and Forecasting Model

Abstract Generating accurate weather forecasts of planetary boundary layer (PBL) properties is challenging in many geographical regions, oftentimes due to complex topography or horizontal variability in, for example, land characteristics. While recent advances in high-performance computing platforms have led to an increase in the spatial resolution of numerical weather prediction (NWP) models, the horizontal gridcell spacing (Δ x ) of many regional-scale NWP models currently fall within or are beginning to approach the gray zone (i.e., Δ x ≈ 100–1000 m). At these gridcell spacings, three-dimensional (3D) effects are important, as the most energetic turbulent eddies are neither fully parameterized (as in traditional mesoscale simulations) nor fully resolved [as in traditional large-eddy simulations (LES)]. In light of this modeling challenge, we have implemented a 3D PBL parameterization for high-resolution mesoscale simulations using the Weather Research and Forecasting Model. The PBL scheme, which is based on the algebraic model developed by Mellor and Yamada, accounts for the 3D effects of turbulence by calculating explicitly the momentum, heat, and moisture flux divergences in addition to the turbulent kinetic energy. In this study, we present results from idealized simulations in the gray zone that illustrate the benefit of using a fully consistent turbulence closure framework under convective conditions. While the 3D PBL scheme reproduces the evolution of convective features more appropriately than the traditional 1D PBL scheme, we highlight the need to improve the turbulent length scale formulation. Significance Statement The spatial resolution of weather models continues to increase at a rapid rate in accordance with the enhancement of computing power. As a result, smaller-scale atmospheric features become more explicitly resolved. However, most numerical models still ignore the impact of horizontal weather variations on boundary layer flows, which becomes more important at these smaller spatial scales. To address this issue, we have implemented a new modeling approach, using fundamental principles, which accounts for horizontal variability. Our results show that including three-dimensional effects of turbulence is necessary to achieve realistic boundary layer characteristics. This novel technique may be useful for many applications including complex terrain flows, pollutant dispersion, and surface–atmosphere interaction studies.

Juliano, Timothy W.↗

NASA NTRS2008

CALIPSO Satellite Lidar Identification Of Elevated Dust Over Australia Compared With Air Quality Model PM60 Forecasts

Global measurements of the vertical distribution of clouds and aerosols have been recorded by the lidar on board the CALIPSO (Cloud Aerosol Lidar Infrared Pathfinder Satellite Observations) satellite since June 2006. Such extensive, height-resolved measurements provide a rare and valuable opportunity for developing, testing and validating various atmospheric models, including global climate, numerical weather prediction, chemical transport and air quality models. Here we report on the initial results of an investigation into the performance of the Australian Air Quality Forecast System (AAQFS) model in forecasting the distribution of elevated dust over the Australian region. The model forecasts of PM60 dust distribution are compared with the CALIPSO lidar Vertical Feature Mask (VFM) data product. The VFM classifies contiguous atmospheric regions of enhanced backscatter as either cloud or aerosols. Aerosols are further classified into six subtypes. By comparing forecast PM60 concentration profiles to the spatial distribution of dust reported in the CALIPSO VFM, we can assess the model s ability to predict the occurrence and the vertical and horizontal extents of dust events within the study area.

Young, Stuart A.↗

NASA NTRSDate not supplied

Building a Real-Time Flood Prediction Model for Improving Early Warning Systems in Ellicott City, Maryland

As flood events in the United States grow in frequency and intensity, the use of applied remote sensing analyses is increasingly necessary for effective flood monitoring and warning systems. The NASA DEVELOP National Program partnered with the local government of Howard County, Maryland, to investigate the use of machine learning for advanced flood risk detection, and to test the feasibility of integrating this approach into the county’s flood early warning system. To strengthen the efforts of the Howard County Office of Emergency Management (OEM), the project developed a statistical model capable of hindcasting the two severe flash flood events that devastated Ellicott City and transitioned to a ‘Long Short-Term Memory’ based sequence-to-sequence deep learning model with 8-hour forecast capability. The team combined data inputs from public sources including river and precipitation gauges, NASA and NOAA Earth observations, and numerical weather model products using scripts written in the Google Colaboratory Python scripting environment. In addition to designing the deep learning architecture, the team trained and tested the model, and evaluated its performance using Nash-Sutcliffe Efficiency. The final product, the Sequentially Trained Real-time EstimAted Model (STREAM) predicts stage height for the Hudson Branch gauge in Ellicott City using data products available in near real-time, including the High-Resolution Rapid Refresh model’s accumulated precipitation forecasts supplemented by stream gauge data from the OEM and the U.S. Geological Survey. STREAM was incorporated into an online dashboard in a user-friendly interface capable of triggering the alarms that initiate emergency response protocols up to 8 hours in advance of a predicted severe flood event. The project demonstrated the potential for the integration of open data and Earth observations into a flood risk forecasting tool capable of informing near real-time decision making.

NASA DEVELOP↗

DOE OSTI2024

naturf: a package for generating urban parameters for numerical weather modeling

The Neighborhood Adaptive Tissues for Urban Resilience Futures tool (NATURF) is a Python workflow that generates files readable by the Weather Research and Forecasting (WRF) model. NATURF uses geopandas and hamilton to calculate 132 building parameters from shapefiles with building footprint and height information. These parameters can be collected and used in many formats, and the primary output is a binary file configured for input to WRF. This workflow is a flexible adaptation of the National/World Urban Database and Access Portal Tool (NUDAPT/WUDAPT) that can be used with any study area at any spatial resolution. The climate modeling community and urban planners can identify the effects of building/neighborhood morphology on the microclimate using the urban parameters and WRF-readable files produced by NATURF. More information on the urban parameters calculated can be found in the documentation.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRSDate not supplied

Building a Real-Time Predictive Flood Model for Improving Early Warning Systems in Ellicott City, Maryland

As flood events in the United States grow in frequency and intensity, the use of applied remote sensing analyses is increasingly necessary for effective flood monitoring and warning systems. The NASA DEVELOP National Program partnered with the Howard County government in Maryland to investigate the use of machine learning for advanced flood risk detection, and to test the feasibility of integrating this approach into the county’s flood early warning system. To strengthen the efforts of the Howard County Office of Emergency Management (OEM), the project developed a prediction model capable of hindcasting the two severe flash flood events that devastated Ellicott City, and transitioned to an Long Short-Term Memory (LSTM) based sequence-to-sequence deep learning model with 8-hour forecast capability. The team combined data inputs from public sources including river and precipitation gauges, NASA and NOAA Earth observations, and numerical weather model products using scripts written in the Google Colaboratory Python scripting environment. In addition to designing the deep learning architecture, the team trained and tested the model, and evaluated its performance using the Nash-Sutcliffe Efficiency (NSE). The final product, called the Sequentially Trained Real-time EstimAted Model (STREAM), predicts stage height for the Hudson Branch gauge in Ellicott City using data products available in near real-time, including the High-Resolution Rapid Refresh (HRRR) model’s accumulated precipitation forecasts supplemented by stream gauge data from the OEM and the U.S. Geological Survey. STREAM was incorporated into an online dashboard in a user-friendly interface capable of triggering the alarms that initiate the OEM’s emergency response protocols up to 8 hours in advance of a predicted severe flood event. The project demonstrated the potential for the integration of open data and Earth observations into a flood risk forecasting tool capable of informing near real-time decision making.

Ryan Hammock↗

NASA NTRS1983

Initialization of mesoscale models - The possible impact of remotely sensed data

Little or no improvement appears to have been achieved in mesoscale numerical prediction. This has been blamed on inherent unpredictability, insufficient spatial resolution, inadequate parameterization of important physical processes, and errors in the numerical discretization of the governing hydrodynamic equations. Attention is presently given to the more fundamental difficulty of forecasting when present conditions are not known with sufficient accuracy, and to the matching of remotely sensed sounding-capability satellite and all-weather Doppler radar data with mesoscale numerical models. This matching is not straightforward; simple and variational forms of four-dimensional assimilation, and Kalman filtering, are suggested approaches but their practical testing remains to be accomplished. The most severe problem faced is that of the temperature and humidity profiles' vertical resolution.

Gal-Chen, T.↗

DOE OSTI2020

PyDA: A Hands-On Introduction to Dynamical Data Assimilation with Python

Dynamic data assimilation offers a suite of algorithms that merge measurement data with numerical simulations to predict accurate state trajectories. Meteorological centers rely heavily on data assimilation to achieve trustworthy weather forecast. With the advance in measurement systems, as well as the reduction in sensor prices, data assimilation (DA) techniques are applicable to various fields, other than meteorology. However, beginners usually face hardships digesting the core ideas from the available sophisticated resources requiring a steep learning curve. In this tutorial, we lay out the mathematical principles behind DA with easy-to-follow Python module implementations so that this group of newcomers can quickly feel the essence of DA algorithms. We explore a series of common variational, and sequential techniques, and highlight major differences and potential extensions. We demonstrate the presented approaches using an array of fluid flow applications with varying levels of complexity.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS2000

1DVAR Analysis of Temperature and Humidity Using GPS Radio Occultation Data

The Global Positioning System enables positioning in 3 dimensions about our planet. It has been operational since 1994. Twenty-four satellites are used to aclile\,e this performance. The signals sent by these satellites are electromagnetic waves travelling through our atmosphere down to the small receivers used by the civilian community and the military. Because of varying meteorological conditions (namely, temperature and humidity changes along the ray path), the rays do not travel in a straight line. They bend towards the surface. As a consequence, the ray path between two points is longer than a straight line, and the time it takes for a signal to travel this distance is longer. In 1995, a small GPS receiver was launched on a satellite (GPS/MET). It become possible to perform radio occultations around the Earth: the source - one of the 24 GPS satellites - is seen by the receiver as it rises or sets around the other side of the Earth. When the source disappears, the receiver progressively loses the signals. By measuring accurately the time delay between the emission and the reception of the signal, it is possible to infer which part of the delay is due to the atmosphere. We use GPS/MET data to retrieve temperature and humidity profiles simultaneously. A specific method is implemented: it combines information from numerical forecasts and GPS observations in an optimal way. Comparing the result with an independent source of observations (weather balloons), we demonstrate that GPS data have the potential to improve weather analyses. We also show that improved temperature and humidity profiles can be obtained using information from a forecast model. This confirms results obtained in this study using simulated data.

Poli, Paul↗

NASA NTRS2012

Dual-Doppler Feasibility Study

When two or more Doppler weather radar systems are monitoring the same region, the Doppler velocities can be combined to form a three-dimensional (3-D) wind vector field thus providing for a more intuitive analysis of the wind field. A real-time display of the 3-D winds can assist forecasters in predicting the onset of convection and severe weather. The data can also be used to initialize local numerical weather prediction models. Two operational Doppler Radar systems are in the vicinity of Kennedy Space Center (KSC) and Cape Canaveral Air Force Station (CCAFS); these systems are operated by the 45th Space Wing (45 SW) and the National Weather Service Melbourne, Fla. (NWS MLB). Dual-Doppler applications were considered by the 45 SW in choosing the site for the new radar. Accordingly, the 45th Weather Squadron (45 WS), NWS MLB and the National Aeronautics and Space Administration tasked the Applied Meteorology Unit (AMU) to investigate the feasibility of establishing dual-Doppler capability using the two existing systems. This study investigated technical, hardware, and software requirements necessary to enable the establishment of a dual-Doppler capability. Review of the available literature pertaining to the dual-Doppler technique and consultation with experts revealed that the physical locations and resulting beam crossing angles of the 45 SW and NWS MLB radars make them ideally suited for a dual-Doppler capability. The dual-Doppler equations were derived to facilitate complete understanding of dual-Doppler synthesis; to determine the technical information requirements; and to determine the components of wind velocity from the equation of continuity and radial velocity data collected by the two Doppler radars. Analysis confirmed the suitability of the existing systems to provide the desired capability. In addition, it is possible that both 45 SW radar data and Terminal Doppler Weather Radar data from Orlando International Airport could be used to alleviate any radar geometry issues at the NWS MLB radar, such as the "cone of silence" or beam blockage. In the event of a radar outage at one of the sites, the multi-radar algorithms would provide continuing coverage of the area through use of the data from the remaining operational radar sites. There are several options to collect, edit, synthesize and display dual-Doppler data sets. These options include commercial packages available for purchase and a variety of freeware packages available from the National Center for Atmospheric Research (NCAR) for processing raw radar data. However, evaluation of the freeware packages revealed that they do not have sufficient documentation and configuration control to be certified for 45 SW use. Additionally, a TI data line must be installed/leased from the NWS MLB office and CCAFS to enable the receipt of NWS MLB raw radar data to use in the dual-Doppler synthesis. Integration of the TI data line into the Eastern Range infrastructure that will meet the security requirements necessary for 45 SW use is time-consuming and costly. Overall evaluation indicates that establishment of the dual-Doppler capability using the existing operational radar systems is desirable and feasible with no technical concerns. Installation of such a system represents a significant enhancement to forecasting capabilities at the 45 WS and at NWS MLB. However, data security and cost considerations must be evaluated in light of current budgetary constraints. In any case, gaining the dual-Doppler capability will provide opportunities for better visualization of the wind field and better forecasting of the onset of convection and severe weather events to support space launch operations at KSC and CCAFS.

Huddleston, Lisa L.↗

NASA NTRS1984

Observing system simulation experiments related to space-borne LIDAR wind profiling. I - Forecast impacts of highly idealized observing systems. II - Sensitivity to atmospheric and instrumental influences

It is pointed out that extended-range numerical weather prediction requires complete global and accurate three-dimensional fields of temperature, moisture, and winds as a function of pressure throughout the troposphere, lower stratosphere, and at the earth's surface. Recent developments in the use of incoherent optical radar (LIDAR) systems indicate the feasibility of obtaining pressure as a function of height with presumably very high vertical resolution and accuracy. In evaluating the feasibility of an employment of spaceborne LIDAR systems as an aid for obtaining data required for weather prediction, simulation studies could provide important information. The present investigation is concerned with the conduction of such simulation experiments. It is found that the results of a composite system of temperature and pressure show significant improvements over those of temperature only. The effects of simulating several sources of LIDAR wind error for two different representations of nature are also examined.

Dlouhy, R.↗

DOE OSTI2024

Sensitivity of Fine‐Resolution Urban Heat Island Simulations to Soil Moisture Parameterization

ABSTRACT Urban areas experience the impact of natural disasters, such as heatwaves and flash floods, disparately in different neighbourhoods across a city. The demand for precise urban hydrometeorological and hydroclimatological modelling to examine this disparity, and the interacting challenges posed by climate change and urbanisation, has thus surged. The Weather Research and Forecasting (WRF) model has served such operational and research purposes for decades. Recent advancements in WRF, including enhanced numerical schemes and sophisticated urban atmospheric‐hydrological parameterizations, have empowered the simulation of urban geophysical processes at high resolution (~1 km), but even this resolution misses significant urban microclimate variability. This study applies the large‐eddy simulations (LES) mode within WRF, coupled with single‐layer urban canopy models (SLUCM), to enable even finer‐scale modelling (150 m) of the Urban Heat Island (UHI) effect in the Baltimore metropolitan area. We run nine scenarios to evaluate various methods of initializing soil moisture and various spinup lead times, and to assess the impact of WRF's Mosaic approach in depicting subgrid‐scale processes. We evaluate the scenarios by comparing the WRF simulated land surface temperature (LST) against Landsat LST and the WRF simulated hourly 2‐m air temperatures (AT) with observations from eight weather stations across the domain. Results underscore the paramount influence of the lead spinup time on the spatiotemporal distribution of simulated soil moisture, consequently shaping WRF's efficacy in predicting the UHI. Furthermore, interpolating soil moisture‐related parameters from the parent for child domain initialization yields a notable reduction in mean and root‐mean‐squared errors. This improvement was particularly evident in simulations with the longest spinup time, affirming the importance of carefully designing the initialization of soil moisture for improved urban temperature predictions.

Talebpour, Mahdad↗

DOE OSTI2022

Smoke from 2020 United States wildfires responsible for substantial solar energy forecast errors

Abstract The 2020 wildfire season (May through December) in the United States was exceptionally active, with the National Interagency Fire Center reporting over 10 million acres ( > 40 000 km 2 ) burned. During the September 2020 wildfire events, large concentrations of smoke particulates were emitted into the atmosphere. As a result, smoke was responsible for ∼10%–30% reduction in solar power production during peak hours as recorded by the California Independent System Operator (CAISO) sites. In this study, we focus on a 9 d period in September when wildfire smoke had a profound impact on solar energy production. During the smoke episodes, hour-ahead forecasts utilized by CAISO did not include the effects of smoke and therefore overestimated the expected power production by ∼10%–50%. Here we use multiple observational networks and a numerical weather prediction (NWP) model to show that the wildfire events of 2020 had a significantly detrimental influence on solar energy production due to high aerosol loading. We find that including the contribution of biomass burning particles greatly improves the day-ahead solar energy bias forecast of both global horizontal irradiance and direct normal irradiance by nearly ∼50%. Our results suggest that a more comprehensive treatment of aerosols, including biomass burning aerosols, in NWP models may be an important consideration for energy grid balancing, in addition to solar resource assessment, as solar power reliance increases.

14 SOLAR ENERGY↗

NASA NTRS1988

Doppler lidar wind measurements on the EOS - LAWS

This paper discusses the efforts currently underway to prepare for the deployment of the Laser Atmospheric Wind Sounder (LAWS) instrument on an Earth Observing System (EOS) polar orbiter and on the Space Station. Attention is given to the measurement techiques, which include ground-based and airborne measurements in the atmosphere, and to the scientific and the technical issues being addressed. Global wind profiles obtained with the LAWS will effect dramatic improvement in numerical weather prediction and will be useful in determining the midtroposphere steering currents that are important in hurricane track forecasts. The global wind data will also provide information on large-scale atmospheric circulation and climate dynamics and on biogeochemical and hydrologic cycles.

Fitzjarrald, Daniel E.↗

NASA NTRS2009

Simulating the Effects of Irrigation over the U.S. in a Land Surface Model Based on Satellite Derived Agricultural Data

A novel method is introduced for integrating satellite derived irrigation data and high-resolution crop type information into a land surface model (LSM). The objective is to improve the simulation of land surface states and fluxes through better representation of agricultural land use. Ultimately, this scheme could enable numerical weather prediction (NWP) models to capture land-atmosphere feedbacks in managed lands more accurately and thus improve forecast skill. Here we show that application of the new irrigation scheme over the continental US significantly influences the surface water and energy balances by modulating the partitioning of water between the surface and the atmosphere. In our experiment, irrigation caused a 12% increase in evapotranspiration (QLE) and an equivalent reduction in the sensible heat flux (QH) averaged over all irrigated areas in the continental US during the 2003 growing season. Local effects were more extreme: irrigation shifted more than 100 W/m from QH to QLE in many locations in California, eastern Idaho, southern Washington, and southern Colorado during peak crop growth. In these cases, the changes in ground heat flux (QG), net radiation (RNET), evapotranspiration (ET), runoff (R), and soil moisture (SM) were more than 3 W/m(sup 2), 20 W/m(sup 2), 5 mm/day, 0.3 mm/day, and 100 mm, respectively. These results are highly relevant to continental- to global-scale water and energy cycle studies that, to date, have struggled to quantify the effects of agricultural management practices such as irrigation. Based on the results presented here, we expect that better representation of managed lands will lead to improved weather and climate forecasting skill when the new irrigation scheme is incorporated into NWP models such as NOAA's Global Forecast System (GFS).

Ozdogan, Mutlu↗

NASA NTRS2009

Sensitivity of Short-Term Weather Forecasts to Assimilated AIRS Data: Implications for NPOESS Applications

The Atmospheric Infrared Sounder (AIRS) is acting as a heritage and risk reduction instrument for the Cross-track lnfrared Sounder (CrIS) to be flown aboard the NPP and NPOESS satellites. The hyperspectral nature of AIRS and CrIS provides high-quality soundings that, along with their asynoptic observation time over North America, make them attractive sources to fill the spatial and temporal data voids in upper air temperature and moisture measurements for use in data assimilation and numerical weather prediction. Observations from AlRS can be assimilated either as direct radiances or retrieved thermodynamic profiles, and the Short-Term Prediction Research and Transition (SPORT) Center at NASA's Marshall Space Flight Center has used both data types to improve short-term (0-48h), regional forecasts. The purpose of this paper is to share SPORT'S experiences using AlRS radiances and retrieved profiles in regional data assimilation activities by showing that proper handling of issues-including cloud contamination and land emissivity characterization-are necessary to produce optimal analyses and forecasts.

Zavodsky, Bradley↗

NASA NTRS2017

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↗