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 721 records · Page 40

A New Tropical Cyclone Dynamic Initialization Technique Using High Temporal and Spatial Density Atmospheric Motion Vectors and Airborne Field Campaign Data

Background: Initialization of tropical cyclones in numerical weather prediction (NWP) systems is a great challenge: Mass-wind field balance; Secondary circulation and heating; Asymmetries. There can be large adjustments in structure and intensity in the first 24 hours if the initial vortex is not in balance: Spurious gravity waves; Spin-up (model and physics). Existing mesoscale NWP model TC (Tropical Cyclone) initialization strategies: Bogus vortex, cold start from global analyses; 3DVAR or 4DVAR, possibly with synthetic observations; EnKF (Ensemble Kalman Filter); Dynamic initialization. Dynamic initialization allows vortex to have improved balance and physics spin-up at the initial time (e.g., Hendricks et al. 2013, 2011; Nguyen and Chen 2011; Fiorino and Warner 1981; Hoke and Anthes 1976). Himawari-8 geostationary satellite has capability of continuous imagery (10-minutes) over the full disk: New GOES-R satellites will have same capability. This will allow for unprecedented observations of tropical cyclones. However, current data assimila1on systems are not capable of ingesting such high temporal observations (Atmospheric Mo1on Vectors - AMVs). Hourly AMVs are produced, and thinned to 100-kilometer spacing in the horizontal. An entirely new data assimilation concept is required to utilize these observations.

remote sensing↗

Uncertainties in Local Intense Precipitation Flood Modeling

In the Task 1 report for the U.S. Nuclear Regulatory Commission (NRC) funded Local Intense Precipitation (LIP) Probabilistic Flood Hazard Assessment (PFHA) Pilot Study, the Pacific Northwest National Laboratory (PNNL) reviewed hydrologic and hydraulic modeling approaches for estimation of LIP flooding and their implementation in readily available simulation software packages (Prasad and Yuan 2020). NRC’s interest is in developing guidance for performing LIP PFHA at nuclear power plant (NPP) sites. This report describes the work done under Task 2 of the study that focused on describing uncertainties associated with LIP flood modeling. Two main types of uncertainties exist in LIP flood modeling: (1) aleatory uncertainties that arise from the inherent natural variability of the hydrometeorological system and (2) epistemic uncertainties that arise from the analysts’ incomplete knowledge of the hydrometeorological and site configuration. Principal sources for aleatory variability in LIP flood simulations at NPP sites include precipitation (i.e., magnitude, duration, and temporal distribution of LIP events), initial conditions (e.g., soil moisture content, stormwater drainage discharge, surface storage), boundary conditions (e.g., upstream discharge, downstream water surface elevations), and long-term temporal trends (e.g., climate change). This report describes approaches used for estimating the aleatory variability in precipitation including precipitation-frequency analysis, numerical weather prediction, and stochastic weather generation. Data sources that can be used to estimate aleatory variability in initial and boundary conditions are also described. Approaches to include effects of long-term trends like those from climate change into estimation of aleatory variability are summarized. Sources of epistemic uncertainty in LIP flood simulations at NPP sites include process representation (e.g., multiple approaches to represent runoff generation, stormwater drainage, and hydraulic routing), site configurations (e.g., site layout, flow features, status of temporary flood protection, blockage of drains), model resolution, and long-term temporal trends (e.g., known/planned site alterations, land-use changes at and in the vicinity of the site). The report describes alternative process representations (methods and models), particularly those implemented in the LIP flood simulation software packages reviewed in Task 1 report, and lists the associated model parameters. Approaches for estimating model parameters when surface and subsurface water exchanges occur are also described.

42 ENGINEERING↗

Resilience Assessment for Distribution Systems during Hurricanes: A Learning-Based Framework

This paper presents a proactive strategy for hurricane-resilient distribution systems. It proposes a Bayesian Neural Network-based outage prediction model considering various parameters, including electrical components, and weather and environmental factors. Addressing challenges in imbalanced outage datasets, a Bias-Variance Tradeoff method is proposed. A resilience assessment model quantifies resilience indices, providing insights into system weaknesses. The approach identifies weak points and serves as a planning benchmark. Numerical results on the modified IEEE 123-node test system demonstrate effectiveness in realistic hurricane scenarios.

Vahedi, Soroush↗

Simulating Virtual Terminal Area Weather Data Bases for Use in the Wake Vortex Avoidance System (Wake VAS) Prediction Algorithm

During the research project, sounding datasets were generated for the region surrounding 9 major airports, including Dallas, TX, Boston, MA, New York, NY, Chicago, IL, St. Louis, MO, Atlanta, GA, Miami, FL, San Francico, CA, and Los Angeles, CA. The numerical simulation of winter and summer environments during which no instrument flight rule impact was occurring at these 9 terminals was performed using the most contemporary version of the Terminal Area PBL Prediction System (TAPPS) model nested from 36 km to 6 km to 1 km horizontal resolution and very detailed vertical resolution in the planetary boundary layer. The soundings from the 1 km model were archived at 30 minute time intervals for a 24 hour period and the vertical dependent variables as well as derived quantities, i.e., 3-dimensional wind components, temperatures, pressures, mixing ratios, turbulence kinetic energy and eddy dissipation rates were then interpolated to 5 m vertical resolution up to 1000 m elevation above ground level. After partial validation against field experiment datasets for Dallas as well as larger scale and much coarser resolution observations at the other 8 airports, these sounding datasets were sent to NASA for use in the Virtual Air Space and Modeling program. The application of these datasets being to determine representative airport weather environments to diagnose the response of simulated wake vortices to realistic atmospheric environments. These virtual datasets are based on large scale observed atmospheric initial conditions that are dynamically interpolated in space and time. The 1 km nested-grid simulated datasets providing a very coarse and highly smoothed representation of airport environment meteorological conditions. Details concerning the airport surface forcing are virtually absent from these simulated datasets although the observed background atmospheric processes have been compared to the simulated fields and the fields were found to accurately replicate the flows surrounding the airport where coarse verification data were available as well as where airport scale datasets were available.

Kaplan, Michael L.↗

Atmospheric Delay Reduction Using KARAT for GPS Analysis and Implications for VLBI

We have been developing a state-of-the-art tool to estimate the atmospheric path delays by raytracing through mesoscale analysis (MANAL) data, which is operationally used for numerical weather prediction by the Japan Meteorological Agency (JMA). The tools, which we have named KAshima RAytracing Tools (KARAT)', are capable of calculating total slant delays and ray-bending angles considering real atmospheric phenomena. The KARAT can estimate atmospheric slant delays by an analytical 2-D ray-propagation model by Thayer and a 3-D Eikonal solver. We compared PPP solutions using KARAT with that using the Global Mapping Function (GMF) and Vienna Mapping Function 1 (VMF1) for GPS sites of the GEONET (GPS Earth Observation Network System) operated by Geographical Survey Institute (GSI). In our comparison 57 stations of GEONET during the year of 2008 were processed. The KARAT solutions are slightly better than the solutions using VMF1 and GMF with linear gradient model for horizontal and height positions. Our results imply that KARAT is a useful tool for an efficient reduction of atmospheric path delays in radio-based space geodetic techniques such as GNSS and VLBI.

Ichikawa, Ryuichi↗

Simulation of 1986 South China Sea Monsoon with a Regional Climate Model

A Regional Land-Atmosphere Climate Simulation System (RELACS) project is being developed at NASA Goddard Space Flight Center. One of the major goals of RELACS is to use a regional scale model with improved physical processes and in particular land-related processes, to understand the role of the land surface and its interaction with convection and radiation as well as the water/energy cycles in the IndoChina/South China Sea (SCS) region. The Penn State/NCAR MM5 atmospheric modeling system, a state of the art atmospheric numerical model designed to simulate regional weather and climate, has been successfully coupled to the Parameterization for Land-Atmosphere-Cloud Exchange (PLACE) land surface model. The original MM5 model (without PLACE) includes the option for either a simple slab soil model or a five-layer soil model (MRF) in which the soil moisture availability evolves over time. However, the MM5 soil models do not include the effects of vegetation, and thus important physical processes such as evapotranspiration and interception are precluded. The PLACE model incorporates vegetation type and has been shown in international comparisons to accurately predict evapotranspiration and runoff over a wide variety of land surfaces. The coupling of MM5 and PLACE creates a numerical modeling system with the potential to more realistically simulate atmosphere and land surface processes including land-sea interaction, regional circulations such as monsoons, and flash flood events. In addition, the Penn State/NCAR MM5 atmospheric modeling system has been: (1) coupled to the Goddard Ice Microphysical scheme; (2) coupled to a turbulent kinetic energy (TKE) scheme; (3) modified to ensure cloud budget balance; and (4) incorporated initialization with the Goddard EOS data sets at NASA/Goddard Laboratory for Atmospheres. The improved MM5 with two nested domains (60 and 20 km horizontal resolution) was used to simulate convective activity over IndoChina and the South China Sea, during the monsoon season, from May 6 to May 20, 1986. The model results captured several dominant observed features, such as twin cyclones, a depression system over the Bay of Bengal, strong south-westerly winds over IndoChina before and during the on-set of convection over the SCS, and a vortex over the SCS. Two additional MM5 runs with different land process models, Blackadar and MRF, were performed, and their results are compared to the run with PLACE. The preliminary results indicate that the MM5 results using PLACE and Blackadar are in very good agreement, but the results using MRF do not contain the south-westerly wind over IndoChina prior to the on-set of convection over the SCS.

Tao, W. -K.↗

Comparison of Optimum Interpolation and Cressman Analyses

The objective of this investigation is to develop a state-of-the-art optimum interpolation (O/I) objective analysis procedure for use in numerical weather prediction studies. A three-dimensional multivariate O/I analysis scheme has been developed. Some characteristics of the GLAS O/I compared with those of the NMC and ECMWF systems are summarized. Some recent enhancements of the GLAS scheme include a univariate analysis of water vapor mixing ratio, a geographically dependent model prediction error correlation function and a multivariate oceanic surface analysis.

Baker, W. E.↗

Comparison of Optimum Interpolation and Cressman Analyses

The development of a state of the art optimum interpolation (O/I) objective analysis procedure for use in numerical weather prediction studies was investigated. A three dimensional multivariate O/I analysis scheme was developed. Some characteristics of the GLAS O/I compared with those of the NMC and ECMWF systems are summarized. Some recent enhancements of the GLAS scheme include a univariate analysis of water vapor mixing ratio, a geographically dependent model prediction error correlation function and a multivariate oceanic surface analysis.

Baker, W. E.↗

Using MODIS estimates of fractional snow cover area to improve streamflow forecasts in interior Alaska

Remotely sensed snow cover observations give an opportunity to improve operational snow melt and stream flow forecasting in remote regions. This is particularly true in Alaska, where remote basins and a spatially and temporally sparse gaging network plague efforts to understand and forecast the hydrology of subarctic boreal basins and where climate change is leading to rapid shifts in basin function. In this study,the operational framework employed by the United States (US) National Weather Service, including the Alaska Pacific River Forecast Center, is adapted to integrate Moderate Resolution Imaging Spectroradiometer (MODIS) remotely sensed observations of fractional snow cover area (fSCA) to determine if these data improve streamflow forecasts in interior Alaska river basins. Two versions of MODIS fSCA are tested against a base case extent of snow cover derived by aerial depletion curves: the MODIS 10A1(MOD10A1) and the MODIS Snow Cover Area and Grain size (MODSCAG) product over the period 2000–2010. Observed runoff is compared to simulated runoff to calibrate both iterations of the model. MODIS-forced simulations haveimproved snow depletion timing compared with snow telemetry sites in the basins, with discernable increases in skill for the streamflow simulations.The MODSCAG fSCA version provides moderate increases in skill but is similar to the MOD10A1 results. The basins with the largest improvement in streamflow simulations have the sparsest streamflow observations.Considering the numerous low-quality gages (discontinuous, short, or unreliable) and ungauged systems throughout the high-latitude regions of the globe, this result is valuable and indicates the utility of the MODIS fSCAdata in these regions. Additionally, while improvements in predicted discharge values are subtle, the snow model better represents the physical conditions of the snow pack and therefore provides more robust simulations,which are consistent with the US National Weather Service's move toward a physically based National Water Model. Physically based models may also be more capable of adapting to changing climates than statistical models corrected to past regimes. This work provides direction for both the Alaska Pacific River Forecast Center and other forecast centers across the US to implement remote-sensing observations within their operational framework, to refine the representation of snow, and to improve streamflow forecasting skill in basins with few or poor-quality observations.

54 ENVIRONMENTAL SCIENCES↗

A GLOBAL TWO-TEMPERATURE CORONA AND INNER HELIOSPHERE MODEL: A COMPREHENSIVE VALIDATION STUDY

The recent solar minimum with very low activity provides us a unique opportunity for validating solar wind models. During CR2077 (2008 November 20 through December 17), the number of sunspots was near the absolute minimum of solar cycle 23. For this solar rotation, we perform a multi-spacecraft validation study for the recently developed three-dimensional, two-temperature, Alfven-wave-driven global solar wind model (a component within the Space Weather Modeling Framework). By using in situ observations from the Solar Terrestrial Relations Observatory (STEREO) A and B, Advanced Composition Explorer (ACE), and Venus Express, we compare the observed proton state (density, temperature, and velocity) and magnetic field of the heliosphere with that predicted by the model. Near the Sun, we validate the numerical model with the electron density obtained from the solar rotational tomography of Solar and Heliospheric Observatory/Large Angle and Spectrometric Coronagraph C2 data in the range of 2.4 to 6 solar radii. Electron temperature and density are determined from differential emission measure tomography (DEMT) of STEREO A and B Extreme Ultraviolet Imager data in the range of 1.035 to 1.225 solar radii. The electron density and temperature derived from the Hinode/Extreme Ultraviolet Imaging Spectrometer data are also used to compare with the DEMT as well as the model output. Moreover, for the first time, we compare ionic charge states of carbon, oxygen, silicon, and iron observed in situ with the ACE/Solar Wind Ion Composition Spectrometer with those predicted by our model. The validation results suggest that most of the model outputs for CR2077 can fit the observations very well. Based on this encouraging result, we therefore expect great improvement for the future modeling of coronal mass ejections (CMEs) and CME-driven shocks.

79 ASTRONOMY AND ASTROPHYSICS↗

Evaluating GXS Impact in the Context of International Coordination

The proposed NOAA/NASA Geostationary eXtended Observations (GeoXO) program plans to include a hyperspectral infrared (IR) sounder on its central satellite. Expected to launch in the mid-2030s, the GeoXO Sounder (GXS) will join international counterparts in a geostationary orbit. Ahead of launch, the NASA Global Modeling and Assimilation Office (GMAO) assessed the potential effectiveness of GXS both as a single GEO IR sounder and as part of a global ring of such instruments, including those already being built by international agencies. Using an observing system simulation experiment (OSSE) framework, GXS was assessed from a global numerical weather prediction (NWP) perspective. The ability of GXS, both alone and as part of a global ring of GEO sounders, to improve weather prediction of thermodynamic variables was evaluated globally and regionally. Compared to a control, GXS dominated regional analysis and forecast improvements, and contributed significantly to global increases in forecast skill. However, more sustained global improvements on the order of 4 days rely on international partnerships. Using the FSOI metric over CONUS, the GXS observations provide the strongest radiance impact on the moist energy error norm reduction. Additionally, GXS shows the capability to improve hurricane forecast track errors, resulting in improved forecast warnings. Overall, the persistent atmospheric profile information from GXS over much of the western hemisphere provide an opportunity to improve the representation of weather systems and their forecasts.

Erica McGrath-Spangler↗

Evaluating GXS Impact in the Context of International Coordination

The proposed NOAA/NASA Geostationary eXtended Observations (GeoXO) program plans to include a hyperspectral infrared (IR) sounder on its central satellite. Expected to launch in the mid-2030s, the GeoXO Sounder (GXS) will join international counterparts in a geostationary orbit. Ahead of launch, the NASA Global Modeling and Assimilation Office (GMAO) assessed the potential effectiveness of GXS both as a single GEO IR sounder and as part of a global ring of such instruments, including those already being built by international agencies. Using an observing system simulation experiment (OSSE) framework, GXS was assessed from a global numerical weather prediction (NWP) perspective. The ability of GXS, both alone and as part of a global ring of GEO sounders, to improve weather prediction of thermodynamic variables was evaluated globally and regionally. Compared to a control, GXS dominated regional analysis and forecast improvements, and contributed significantly to global increases in forecast skill. However, more sustained global improvements on the order of 4 days rely on international partnerships. Using the FSOI metric over CONUS, the GXS observations provide the strongest radiance impact on the moist energy error norm reduction. Additionally, GXS shows the capability to improve hurricane forecast track errors, resulting in improved forecast warnings. Overall, the persistent atmospheric profile information from GXS over much of the western hemisphere provide an opportunity to improve the representation of weather systems and their forecasts.

Erica McGrath-Spangler↗

Laboratory for Atmospheres 2008 Technical Highlights

The 2008 Technical Highlights describes the efforts of all members of the Laboratory for Atmospheres. Their dedication to advancing Earth Science through conducting research, developing and running models, designing instruments, managing projects, running field campaigns, and numerous other activities, is highlighted in this report. The Laboratory for Atmospheres (Code 613) is part of the Earth Sciences Division (Code 610), formerly the Earth Sun Exploration Division, under the Sciences and Exploration Directorate (Code 600) based at NASA s Goddard Space Flight Center in Greenbelt, Maryland. In line with NASA s Exploration Initiative, the Laboratory executes a comprehensive research and technology development program dedicated to advancing knowledge and understanding of the atmospheres of Earth and other planets. The research program is aimed at understanding the influence of solar variability on the Earth s climate; predicting the weather and climate of Earth; understanding the structure, dynamics, and radiative properties of precipitation, clouds, and aerosols; understanding atmospheric chemistry, especially the role of natural and anthropogenic trace species on the ozone balance in the stratosphere and the troposphere; and advancing our understanding of physical properties of Earth s atmosphere. The research program identifies problems and requirements for atmospheric observations via satellite missions. Laboratory scientists conceive, design, develop, and implement ultraviolet, infrared, optical, radar, laser, and lidar technology for remote sensing of the atmosphere. Laboratory members conduct field measurements for satellite data calibration and validation, and carry out numerous modeling activities. These modeling activities include climate model simulations, modeling the chemistry and transport of trace species on regional-to-global scales, cloud-resolving models, and development of next-generation Earth system models. Interdisciplinary research is carried out in collaboration with other laboratories and research groups within the Earth Sciences Division, across the Sciences and Exploration Directorate, and with partners in universities and other Government agencies. The Laboratory for Atmospheres is a vital participant in NASA s research agenda. Our Laboratory often has relatively large programs, sizable satellite missions, and observational campaigns that require the cooperative and collaborative efforts of many scientists. We ensure an appropriate balance between our scientists responsibility for these large collaborative projects and their need for an active individual research agenda. This balance allows members of the Laboratory to continuously improve their scientific credentials. Members of the Laboratory interact with the general public to support a wide range of interests in the atmospheric sciences. Among other activities, the Laboratory raises the public s awareness of atmospheric science by presenting public lectures and demonstrations, by making scientific data available to wide audiences, by teaching, and by mentoring students and teachers. The Laboratory makes substantial efforts to attract new scientists to the various areas of atmospheric research. We strongly encourage the establishment of partnerships with Federal and state agencies that have operational responsibilities to promote the societal application of our science products. This report describes our role in NASA s mission, gives a broad description of our research, and summarizes our scientists major accomplishments during calendar year 2008. The report also contains useful information on human resources, scientific interactions, and outreach activities.

Cote, Charles E.↗

The GOES-R Lightning Mapper Sensor

The Lightning Mapper Sensor on GOES-R builds on previous measurements of lightning from low earth orbit by the OTD (Optical Transient Detector) and LIS (Lightning Imaging Sensor) sensors. Unlike observations from low earth orbit, the GOES-R platform will allow continuous monitoring of lightning activity over the Continental United States and southern Canada, Central and South America, and portions of the Atlantic and Pacific Oceans. The LMS will detect total (cloud-to-ground and intracloud) lightning at storm scale resolution (approx. 8 km) using a highly sensitive Charge Coupled Device (CCD) detector array. Discrimination between lightning optical transients and a bright sunlit background scene is accomplished by employing spectral, spatial, and temporal filtering along with a background subtraction technique. The result is 24 hour detection capability of total lightning. These total lightning observations can be made available to users within about 20 seconds. Research indicates a number of ways that total lightning observations from LMS could benefit operational activities, including 1) potential increases in lead times and reduced false alarms for severe thunderstorm and tornado Warnings, 2) improved routing of &rail around thunderstorms, 3) support for spacecraft launches and landings, 4) improved ability to monitor tropical cyclone intensity, 5) ability to monitor thunderstorm intensification/weakening during radar outages or where radar coverage is poor, 6) better identification of deep convection for the initialization of numerical prediction models, 7) improved forest fire forecasts, 8) identification of convective initiation, 9) identification of heavy convective snowfall, and 10) enhanced temporal resolution of storm evolution (1 minute) than is available from radar observations. Total lightning data has been used in an operational environment since July 2003 at the Huntsville, Alabama National Weather Service office. Total lightning measurements are obtained by the North Alabama Lightning Mapping Array (LMA) and have successfully been used in warning decisions. Every 2 minutes, total lightning counts in 2 km by 2 km horizontal, 1 km vertical grids are available to forecasters on an AWIPS (Advanced Weather Interactive Processing System) workstation. Experience with the LMA total lightning data is used to illustrate the potential use of LMS data that would be available to forecasters across the US. This abstract is for submission as a presentation to the National Weather Association Annual Meeting to be held 16-21 October 2004 in Portland, OR. This abstract will be published in the conference proceedings.

Buechler, Dennis↗

The Super Tuesday Outbreak: Forecast Sensitivities to Single-Moment Microphysics Schemes

Forecast precipitation and radar characteristics are used by operational centers to guide the issuance of advisory products. As operational numerical weather prediction is performed at increasingly finer spatial resolution, convective precipitation traditionally represented by sub-grid scale parameterization schemes is now being determined explicitly through single- or multi-moment bulk water microphysics routines. Gains in forecasting skill are expected through improved simulation of clouds and their microphysical processes. High resolution model grids and advanced parameterizations are now available through steady increases in computer resources. As with any parameterization, their reliability must be measured through performance metrics, with errors noted and targeted for improvement. Furthermore, the use of these schemes within an operational framework requires an understanding of limitations and an estimate of biases so that forecasters and model development teams can be aware of potential errors. The National Severe Storms Laboratory (NSSL) Spring Experiments have produced daily, high resolution forecasts used to evaluate forecast skill among an ensemble with varied physical parameterizations and data assimilation techniques. In this research, high resolution forecasts of the 5-6 February 2008 Super Tuesday Outbreak are replicated using the NSSL configuration in order to evaluate two components of simulated convection on a large domain: sensitivities of quantitative precipitation forecasts to assumptions within a single-moment bulk water microphysics scheme, and to determine if these schemes accurately depict the reflectivity characteristics of well-simulated, organized, cold frontal convection. As radar returns are sensitive to the amount of hydrometeor mass and the distribution of mass among variably sized targets, radar comparisons may guide potential improvements to a single-moment scheme. In addition, object-based verification metrics are evaluated for their utility in gauging model performance and QPF variability.

Molthan, Andrew L.↗

A new scheme for the parameterization of the turbulent planetary boundary layer in the GLAS fourth order GCM

Methods being used to increase the horizontal and vertical resolution and to implement more sophisticated parameterization schemes for general circulation models (GCM) run on newer, more powerful computers are described. Attention is focused on the NASA-Goddard Laboratory for Atmospherics fourth order GCM. A new planetary boundary layer (PBL) model has been developed which features explicit resolution of two or more layers. Numerical models are presented for parameterizing the turbulent vertical heat, momentum and moisture fluxes at the earth's surface and between the layers in the PBL model. An extended Monin-Obhukov similarity scheme is applied to express the relationships between the lowest levels of the GCM and the surface fluxes. On-line weather prediction experiments are to be run to test the effects of the higher resolution thereby obtained for dynamic atmospheric processes.

Helfand, H. M.↗

Acceleration of the IMplicit–EXplicit nonhydrostatic unified model of the atmosphere on manycore processors

We present the acceleration of an IMplicit-EXplicit (IMEX) nonhydrostatic atmospheric model on manycore processors such as graphic processing units (GPUs) and Intel's Many Integrated Core (MIC) architecture. IMEX time integration methods sidestep the constraint imposed by the Courant-Friedrichs-Lewy condition on explicit methods through corrective implicit solves within each time step. In this work, we implement and evaluate the performance of IMEX on manycore processors relative to explicit methods. Using 3D-IMEX at Courant number C = 15, we obtained a speedup of about 4x relative to an explicit time stepping method run with the maximum allowable C = 1. Moreover, the unconditional stability of IMEX with respect to the fast waves means the speedup can increase significantly with the Courant number as long as the accuracy of the resulting solution is acceptable. We show a speedup of 100x at C = 150 using 1D-IMEX to demonstrate this point. Several improvements on the IMEX procedure were necessary in order to outperform our results with explicit methods: (a) reducing the number of degrees of freedom of the IMEX formulation by forming the Schur complement, (b) formulating a horizontally explicit vertically implicit 1D-IMEX scheme that has a lower workload and better scalability than 3D-IMEX, (c) using high-order polynomial preconditioners to reduce the condition number of the resulting system, and (d) using a direct solver for the 1D-IMEX method by performing and storing LU factorizations once to obtain a constant cost for any Courant number. Without all of these improvements, explicit time integration methods turned out to be difficult to beat. We discuss in detail the IMEX infrastructure required for formulating and implementing efficient methods on manycore processors. Several parametric studies are conducted to demonstrate the gain from each of the abovementioned improvements. Lastly, we validate our results with standard benchmark problems in numerical weather prediction and evaluate the performance and scalability of the IMEX method using up to 4192 GPUs and 16 Knights Landing processors.

97 MATHEMATICS AND COMPUTING↗

An Introduction to Observing System Simulation Experiments

Observing System Simulation Experiments (OSSEs) are used to estimate the potential impact of proposed new instruments and data on numerical weather prediction. OSSEs can also be used to help design new observing platforms and to investigate the behavior of data assimilation systems. A basic overview of how to design and perform an OSSE will be given, as well as best practices and pitfalls. Some examples using the OSSE framework developed at the NASA Global Modeling and Assimilation Office will be shown.

OSSE↗