Search NASA⌕ Search

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.

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

At least 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

The Potential of GEO Sounders for Numerical Weather Prediction

The future of Earth observations from space is currently being decided upon in the United States. To better understand the impacts of various configurations, the Global Modeling and Assimilation Office (GMAO) observing system simulation experiment (OSSE) framework has been utilized to examine the impact of various hyperspectral infrared sounder constellations from the perspective of global numerical weather prediction (NWP). With an emphasis on the proposed NOAA/NASA Geostationary eXtended Observations (GeoXO) Sounder (GXS), planned to launch in the mid-2030s, the NWP benefit of a global ring of geostationary sounders, consistent with the WMO’s 2040 vision is evaluated. In addition, considerations extend to the future of the Low Earth Orbit (LEO) sounder program with the loss of several existing instruments and plans for future missions ongoing. The GMAO OSSE framework has examined the potential complementarities of GEO and LEO sounders in weather prediction improvement. Overall, the inclusion of both GEO and LEO sounders produces the most beneficial impact on NWP for regional and global scales and the most accurate global weather forecasts as evaluated using several metrics, including hurricanes and the forecast sensitivity observation impact (FSOI).

Erica Mcgrath-Spangler↗
DOE OSTI2024

Learning nonlinear operators in latent spaces for real-time predictions of complex dynamics in physical systems

Abstract Predicting complex dynamics in physical applications governed by partial differential equations in real-time is nearly impossible with traditional numerical simulations due to high computational cost. Neural operators offer a solution by approximating mappings between infinite-dimensional Banach spaces, yet their performance degrades with system size and complexity. We propose an approach for learning neural operators in latent spaces, facilitating real-time predictions for highly nonlinear and multiscale systems on high-dimensional domains. Our method utilizes the deep operator network architecture on a low-dimensional latent space to efficiently approximate underlying operators. Demonstrations on material fracture, fluid flow prediction, and climate modeling highlight superior prediction accuracy and computational efficiency compared to existing methods. Notably, our approach enables approximating large-scale atmospheric flows with millions of degrees, enhancing weather and climate forecasts. Here we show that the proposed approach enables real-time predictions that can facilitate decision-making for a wide range of applications in science and engineering.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1985

A mesoscale numerical forecast of an intense convective snowburst along the East Coast

Mesoscale numerical forecasts utilizing the Mesoscale Atmospheric Simulation System (MASS) are documented for a convective snowburst in the Washington, D.C.-Baltimore, Maryland, region on March 8, 1984 that was associated with secondary cyclogenesis along the East Coast. The forecasts are presented to demonstrate the ability of a mesoscale model to simulate dynamical interactions and diabatic process for a wintertime convective event that was inadequately predicted by local forecasters and to note some of the possible benefits of using mesoscale models for day-to-day forecasting. The results from this and other recent mesoscale modeling studies indicate that three-hourly output of key model fields, when combined with other data sources, can be a valuable aid to forecasters concerned with predicting weather events that are mesoscale in character.

Kocin, P. J.↗
DOE OSTI2022

Airborne Lidar Measurements of XCO 2 in Synoptically Active Environment and Associated Comparisons With Numerical Simulations

Frontal boundaries have been shown to cause large changes in CO 2 mole-fractions, but clouds and the complex vertical structure of fronts make these gradients difficult to observe. It remains unclear how the column average CO 2 dry air mole-fraction (XCO 2 ) changes spatially across fronts, and how well airborne lidar observations, data assimilation systems, and numerical models without assimilation capture XCO 2 frontal contrasts (ΔXCO 2 , i.e., warm minus cold sector average of XCO 2 ). We demonstrated the potential of airborne Multifunctional Fiber Laser Lidar (MFLL) measurements in heterogeneous weather conditions (i.e., frontal environment) to investigate the ΔXCO 2 during four seasonal field campaigns of the Atmospheric Carbon and Transport-America (ACT-America) mission. Most frontal cases in summer (winter) reveal higher (lower) XCO 2 in the warm (cold) sector than in the cold (warm) sector. During the transitional seasons (spring and fall), no clear signal in ΔXCO 2 was observed. Intercomparison among the MFLL, assimilated fields from NASA's Global Modeling and Assimilation Office (GMAO), and simulations from the Weather Research and Forecasting-—Chemistry (WRF-Chem) showed that (a) all products had a similar sign of ΔXCO 2 though with different levels of agreement in ΔXCO 2 magnitudes among seasons; (b) ΔXCO 2 in summer decreases with altitude; and (c) significant challenges remain in observing and simulating XCO 2 frontal contrasts. A linear regression analyses between ΔXCO 2 for MFLL versus GMAO, and MFLL versus WRF-Chem for summer-2016 cases yielded a correlation coefficient of 0.95 and 0.88, respectively. The reported ΔXCO 2 variability among four seasons provide guidance to the spatial structures of XCO 2 transport errors in models and satellite measurements of XCO 2 in synoptically-active weather systems.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2003

The Madden-Julian Oscillation and its Impact on Northern Hemisphere Weather Predictability during Wintertime

The Madden-Julian Oscillation (MJO) is known as the dominant mode of tropical intraseasonal variability and has an important role in the coupled-atmosphere system. This study used twin numerical model experiments to investigate the influence of the MJO activity on weather predictability in the midlatitudes of the Northern Hemisphere during boreal winter. The National Aeronautics and Space Administration (NASA) Goddard laboratory for the Atmospheres (GLA) general circulation model was first used in a 10-yr simulation with fixed climatological SSTs to generate a validation data set as well as to select initial conditions for active MJO periods and Null cases. Two perturbation numerical experiments were performed for the 75 cases selected [(4 MJO phases + Null phase) _ 15 initial conditions in each]. For each alternative initial condition, the model was integrated for 90 days. Mean anomaly correlations in the midlatitudes of the Northern Hemisphere (2O deg N_60 deg.N) and standardized root-mean-square errors were computed to validate forecasts and control run. The analyses of 500-hPa geopotential height, 200-hPa Streamfunction and 850-hPa zonal wind component systematically show larger predictability during periods of active MJO as opposed to quiescent episodes of the oscillation.

Jones, Charles↗
NASA NTRSDate not supplied

Assimilation of Soil Moisture Observations Over Land Improves Analysis and Prediction of Tropical Cyclone Idai

Soil moisture conditions can impact the circulation and structure of a tropical cyclone (TC) when part or all of the circulation is over land. Dry land surface conditions may lead to faster dissipation of a TC over land, whereas very wet conditions may lead to a prolonged maintenance of its intensity. While this relationship is relatively well understood in theory, applications of these findings in the context of numerical weather prediction (NWP) have been limited. Here we present a case study that explores the potential of improving TC predictions through an improved soil moisture initialization in an NWP framework. Specifically, we examine the impact of assimilating observations from the NASA Soil Moisture Active Passive (SMAP) mission into the NASA Goddard Earth Observing System (GEOS) global weather model on the prediction of South-West Indian Ocean TC Idai (2019). SMAP provides accurate L-band (1.4 GHz) brightness temperatures (Tb) observations that are sensitive to soil moisture globally and at high revisit times of 2-3 days. It has previously been shown that the assimilation of SMAP Tbs significantly improves modeled land surface states. Here we evaluate: (i) forecasts initialized from an analysis that is comparable to the GEOS operational analysis (without SMAP Tb assimilation) and (ii) forecasts initialized from an analysis that additionally assimilates SMAP Tb observations. We find that in the analysis with SMAP assimilation, the TC has a better-defined, more aligned vertical structure over land relative to the control run; moreover, the analyzed TC size, as measured by the wind speed radius, better matches the observed TC size. We further find significant reductions in the forecast intensity error and the forecast along-track error, measured against observations. The largest error reductions occur at lead times of 36 to 72 hours, suggesting that the land with its longer memory gains in importance as a source of predictability at this timescale. An investigation of the underlying mechanisms leading to the skill improvements from SMAP data assimilation revealed that the assimilation of SMAP leads to wetter soil moisture conditions and an increased latent heat flux in the SMAP analysis, which results in a TC with higher column-integrated total moisture content and total energy compared to the control analysis.

Jana Kolassa↗
NASA NTRS1990

An evaluation of the structure of tropical intraseasonal oscillations in three general circulation models

The attributes for the simulation of tropical intraseasonal oscillations in the 30 to 60 day range by three GCMs are compared with analyses of the European Center for Medium-Range Weather Forecasts. Models under consideration include the Goddard Laboratory of Atmospheric sciences (GLAS), Goddard Laboratory for Atmospheres (GLA), and UCLA GCMs. They are characterized by the same resolutions of four degrees latitude by five degrees longitude with 9 levels. Observations reveal that the eastward traveling planetary scale structure becomes more stationary over the Indonesian region and accelerates over the Central Pacific. In the tropics a Kelvin wave-type structure is found to be dominant near the center of the oscillation. The simulated winds are realistic although the meridional component is too strong, especially in the GLA GCM. The differences in the structure of the oscillation in the GLAS GCM and GLA GCM are considered to be a consequence of the different numerical schemes used. The GCMs are characterized by preferred zones for diabatic heating, with a turn-on heating occurring when the rising branch of the intraseasonal oscillation passes over these convective regions.

Park, Chung-Kyu↗
NASA NTRS1984

VAS data assimilation into a mesoscale model

In numerical weather prediction, one of the arising problems has always been related to insufficient data, especially in cases involving scales smaller than the synoptic rawinsonde network. The Visible Infrared Spin-Scan Radiometer (VISSR) Atmospheric Sounder (VAS) instrument on the Geostationary Operational Environmental Satellite (GOES) can provide mesoscale vertical profiles of temperature and moisture several times over a twelve hour period. However, as this data is only available over a relatively small area at any one time, the VAS data is more suited to data assimilation. A variational method has been developed to assimilate VAS data into a mesoscale model. In the present paper, the assimilation method and VAS data from one case are discussed. It is found that the greatest improvements to the forecast resulted when the VAS data resolved the mesoscale structure of the temperature and relative humidity fields.

Cram, J. M.↗
NASA NTRS2010

Tropical Cyclones, Hurricanes, and Climate: NASA's Global Cloud-Scale Simulations and New Observations that Characterize the Lifecycle of Hurricanes

One of the primary interests of Global Change research is the impact of climate changes and climate variability on extreme weather events, such as intense tropical storms and hurricanes. Atmospheric climate models run at resolutions of global weather models have been used to study the impact of climate variability, as seen in sea surface temperatures, on the frequency and intensity of tropical cyclones. NASA's Goddard Earth Observing System Model, version 5 (GEOS-5) in ensembles run at 50 km resolution has been able to reproduce the interannual variations of tropical cyclone frequency seen in nature. This, and other global models, have found it much more difficult to reproduce the interannual changes in intensity, a result that reflects the inability of the models to simulate the intensities of the most extreme storms. Better representation of the structures of cyclones requires much higher resolution models. Such improved representation is also fundamental to making best use of satellite observations. In collaboration with NOAA's Geophysical Fluid Dynamics Laboratory, GEOS-5 now has the capability of running at much higher resolution to better represent cloud-scale resolutions. Global simulations at cloud-permitting resolutions (10- to 3.5-km) allows for the development of realistic tropical cyclones from tropical storm 119 km/hr winds) to category 5 (>249km1hr winds) intensities. GEOS-5 has produced realistic rain-band and eye-wall structures in tropical cyclones that can be directly analyzed against satellite observations. For the first time a global climate model is capable of representing realistic intensity and track variability on a seasonal scale across basins. GEOS-5 is also used in assimilation mode to test the impact of NASA's observations on tropical cyclone forecasts. One such test, for tropical cyclone Nargis in the Indian Ocean in May 2008, showed that observations from Atmospheric Infrared Sounder (AIRS) and the Advanced Microwave Sounding Unit (AMSU-A) on Aqua substantially reduced forecast track errors. Tropical cyclones in the northern Indian Ocean pose serious challenges to operational weather forecasting systems, partly due to their shorter lifespan and more erratic track, compared to those in the Atlantic and the Pacific. SA is also bringing several state of the art instruments in recent field campaigns to peer under the clouds and study the inner workings of the tropical storms. With the Genesis and Rapid Intensification Processes (GRIP) experiment, a NASA Earth science field experiment in 2010 that includes the Global Hawk Unmanned Airborne System (UAS) configured with a suite of in situ and remote sensing instruments that are observing and characterizing the lifecycle of hurricanes, we expect significant improvement in our understanding of the track and intensification processes with the assimilation of the satellite and field campaign observations of meteorological parameters in the numerical prediction models.

Putman, William M.↗
NASA NTRS2009

Forecasting Lightning Threat using Cloud-resolving Model Simulations

As numerical forecasts capable of resolving individual convective clouds become more common, it is of interest to see if quantitative forecasts of lightning flash rate density are possible, based on fields computed by the numerical model. Previous observational research has shown robust relationships between observed lightning flash rates and inferred updraft and large precipitation ice fields in the mixed phase regions of storms, and that these relationships might allow simulated fields to serve as proxies for lightning flash rate density. It is shown in this paper that two simple proxy fields do indeed provide reasonable and cost-effective bases for creating time-evolving maps of predicted lightning flash rate density, judging from a series of diverse simulation case study events in North Alabama for which Lightning Mapping Array data provide ground truth. One method is based on the product of upward velocity and the mixing ratio of precipitating ice hydrometeors, modeled as graupel only, in the mixed phase region of storms at the -15\dgc\ level, while the second method is based on the vertically integrated amounts of ice hydrometeors in each model grid column. Each method can be calibrated by comparing domainwide statistics of the peak values of simulated flash rate proxy fields against domainwide peak total lightning flash rate density data from observations. Tests show that the first method is able to capture much of the temporal variability of the lightning threat, while the second method does a better job of depicting the areal coverage of the threat. A blended solution is designed to retain most of the temporal sensitivity of the first method, while adding the improved spatial coverage of the second. Weather Research and Forecast Model simulations of selected North Alabama cases show that this model can distinguish the general character and intensity of most convective events, and that the proposed methods show promise as a means of generating quantitatively realistic fields of lightning threat. However, because models tend to have more difficulty in correctly predicting the instantaneous placement of storms, forecasts of the detailed location of the lightning threat based on single simulations can be in error. Although these model shortcomings presently limit the precision of lightning threat forecasts from individual runs of current generation models, the techniques proposed herein should continue to be applicable as newer and more accurate physically-based model versions, physical parameterizations, initialization techniques and ensembles of cloud-allowing forecasts become available.

McCaul, E. W., Jr.↗
NASA NTRS1986

Comparison of vertical velocities analyzed by a numerical model and measured by a VHF wind profiler

The use of wind profilers for measuring vertical velocities in the troposphere and lower stratosphere is potentially of great interest for verification of forecasts, diagnosis of mesoscale circulations, and studies of wave motions. The studies of profiler vertical velocities to date have shown that the observed patterns of ascent and subsidence are reasonable when compared to the synoptic conditions. However, difficulties arise when a direct verification of the profiler vertical winds is sought. Since no other technique can measure the vertical velocities over the same height range and with the same claimed accuracy as the profilers, direct comparisons are impossible. The only alternative is to compare the measurements to analyzed vertical velocity fields. Here, researchers compare vertical measurements made with the SOUSY VHF radar over a period of 11 days at the beginning of November 1981 to the analyzed vertical velocities produced by the European Center for Medium-range Weather Forecasting (ECMWF) model for grid points near the radar site.

Larsen, M. F.↗
DOE OSTI2022

Physics-Informed AI for Climate and Weather Risk Prediction (Final Report)

As part of the work Terrafuse developed 1) a model for wildfire risk in California, and 2) a model for downscaled wind fields from Numerical Weather Prediction (NWP) numerical models. The wildfire model is based on 20 years of historical data and captures the dependence of wildfire incidence and spread on climatic, weather and land-use variables by training a machine learning model. Nonlinear relationships between input features are learned and expressed by the model and model transparency allows features to be ranked and interpreted. The fire model is of use for accurately predicting real-time daily and long-term wildfire risk for use cases in energy and insurance. The downscaled wind model is a spatio-temporal deep learning model that emulates the influence of high-resolution variables on wind speed, allowing coarse-resolution operational NWP models to be accurately expressed on fine grids at high resolution, with application to wind energy and weather forecasting.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Evaluation of the Four-Dimensional Ensemble-Variational Hybrid Data Assimilation with Self-Consistent Regional Background Error Covariance for Improved Hurricane Intensity Forecasts

The feasibility of a hurricane initialization framework based on the Gridpoint Statistical Interpolation (GSI)-based four-dimensional ensemble-variational (GSI-4DEnVar) hybrid data assimilation system for the Hurricane Weather Research and Forecasting model (HWRF) model is evaluated in this study. The system considers the temporal evolution of error covariances via the use of four-dimensional ensemble perturbations that are provided by high-resolution, self-consistent HWRF ensemble forecasts. It is different from the configuration of the GSI-based three-dimensional ensemble-variational (GSI-3DEnVar) hybrid data assimilation system, similar to that used in the operational HWRF, which employs background error covariances provided by coarser-resolution global ensembles from the National Centers for Environmental Prediction (NCEP) Global Forecast System (GFS) ensemble Kalman filtering data assimilation system. In addition, our proposed initialization framework discards the empirical intensity correction in the vortex initialization package that is employed by the GSI-3DEnVar initialization framework in operational HWRF. Data assimilation and numerical simulation experiments for Hurricanes Joaquin (2015), Patricia (2015), and Matthew (2016) are conducted during their intensity changes. The impacts of two initialization frameworks on the HWRF analyses and forecasts are compared. It is found that GSI-4DEnVar leads to a reduction in track, minimum sea level pressure (MSLP), and maximum surface wind (MSW) forecast errors in all of the HWRF simulations, compared with the GSI-3DEnVar initialization framework. With assimilating high-resolution observations within the hurricane inner-core region, GSI-4DEnVar can produce the initial hurricane intensity reasonably well without the empirical vortex intensity correction. Further diagnoses with Hurricane Joaquin indicate that GSI-4DEnVar can significantly alleviate the imbalances in the initial conditions and enhance the performance of the data assimilation and subsequent hurricane intensity and precipitation forecasts

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2010

Extracting Independent Local Oscillatory Geophysical Signals by Geodetic Tropospheric Delay

Zenith Tropospheric Delay (ZTD) due to water vapor derived from space geodetic techniques and numerical weather prediction simulated-reanalysis data exhibits non-linear and non-stationary properties akin to those in the crucial geophysical signals of interest to the research community. These time series, once decomposed into additive (and stochastic) components, have information about the long term global change (the trend) and other interpretable (quasi-) periodic components such as seasonal cycles and noise. Such stochastic component(s) could be a function that exhibits at most one extremum within a data span or a monotonic function within a certain temporal span. In this contribution, we examine the use of the combined Ensemble Empirical Mode Decomposition (EEMD) and Independent Component Analysis (ICA): the EEMD-ICA algorithm to extract the independent local oscillatory stochastic components in the tropospheric delay derived from the European Centre for Medium-Range Weather Forecasts (ECMWF) over six geodetic sites (HartRAO, Hobart26, Wettzell, Gilcreek, Westford, and Tsukub32). The proposed methodology allows independent geophysical processes to be extracted and assessed. Analysis of the quality index of the Independent Components (ICs) derived for each cluster of local oscillatory components (also called the Intrinsic Mode Functions (IMFs)) for all the geodetic stations considered in the study demonstrate that they are strongly site dependent. Such strong dependency seems to suggest that the localized geophysical signals embedded in the ZTD over the geodetic sites are not correlated. Further, from the viewpoint of non-linear dynamical systems, four geophysical signals the Quasi-Biennial Oscillation (QBO) index derived from the NCEP/NCAR reanalysis, the Southern Oscillation Index (SOI) anomaly from NCEP, the SIDC monthly Sun Spot Number (SSN), and the Length of Day (LoD) are linked to the extracted signal components from ZTD. Results from the synchronization analysis show that ZTD and the geophysical signals exhibit (albeit subtle) site dependent phase synchronization index.

Botai, O. J.↗
NASA NTRSDate not supplied

Estimates of Future Numerical Weather Prediction Impacts from Hyperspectral Sounders

The United States’ satellite programs are currently in the process of evaluating the future of satellite observations of the Earth’s atmosphere, both as an individual nation and within the context of international partnerships. In an effort to better understand the impacts of various alternatives, the Global Modeling and Assimilation Office (GMAO) observing system simulation experiment (OSSE) framework has been utilized to evaluate several configurations of the international hyperspectral infrared sounder constellation from the perspective of global numerical weather prediction (NWP). This was done with an emphasis on the proposed NOAA/NASA Geostationary eXtended Observations (GeoXO) Sounder (GXS), planned to launch in the mid-2030s, and similar planned missions from international partners. GXS, in addition to contributions from EUMETSAT and JMA, will form a global ring of geostationary sounders consistent with the WMO’s 2040 vision. In addition to a novel geostationary IR sounder program, considerations extend to the future of the Low Earth Orbit (LEO) sounder program with the upcoming demise of several existing instruments and plans for future missions ongoing. The GMAO OSSE framework has examined the potential roles of GEO and LEO sounders in weather prediction improvement. Overall, the inclusion of both GEO sounders in a global ring and LEO sounders produces the most beneficial observation impact and the most accurate global weather forecasts as evaluated using several metrics, including hurricanes and the FSOI.

Erica Mcgrath-Spangler↗
DOE OSTI2025

Quantifying UAS Observation Error Variance Used in Data Assimilation Systems and Its Impact on Predictive Skill

Observation error determines the weights of the observations and background state used in data assimilation to generate analyses. Quantifying observation error is critical for the optimal assimilation of observational data sets. Uncrewed Aircraft System (UAS) observations have shown potential benefits in filling observational gaps in the lower atmosphere; however, characterization of their error characteristics has been limited. To optimize the use of UAS observations in numerical weather prediction, UAS observation error is estimated based on the 3‐cornered hat diagnostic approach which uses three independent estimates of the atmospheric state. This approach is applied to data from the 2018 Lower Atmospheric Profiling Studies at Elevation‐a Remotely‐piloted Aircraft Team Experiment field campaign using collocated UAS and rawinsonde observations along with output from a set of convection‐permitting model simulations. The estimated observation error values for UAS temperature, wind, and relative humidity measurements were found to be only weakly dependent on height AGL with mean values equal to 0.5°C, 0.8 m s −1 , and 3%, respectively. Only the newly estimated observation error for temperature differed from that previously used to assimilate commercial aircraft observations into global models (1.0°C). However, using this reduced temperature observation error produced more accurate mesoscale analyses and forecasts of both terrain‐driven flows and convection initiation generated by colliding outflow boundaries within the San Luis Valley of Colorado.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2001

Assimilation of GOES Land Surface Data: Benefits to Numerical Weather Prediction

A technique has been developed for assimilating GOES-derived skin temperature tendencies and insolation into the surface energy budget equation of a mesoscale model so that the simulated rate of temperature change closely agrees with the satellite observations. A critical assumption of the technique is that the availability of moisture (either from the soil or vegetation) is the least known term in the model's surface energy budget. Therefore, the simulated latent heat flux, which is a function of surface moisture availability, is adjusted based upon differences between the modeled and satellite observed skin temperature tendencies. An advantage of this technique is that satellite temperature tendencies are assimilated in an energetically consistent manner that avoids energy imbalances and surface stability problems that arise from direct assimilation of surface shelter temperatures. The fact that the rate of change of the satellite skin temperature is used rather than the absolute temperature means that sensor calibration is not as critical. The technique has been employed on a semi-operational basis at the GHCC within the PSU/CAR MM5 since 1 November 1998. A one-way nested grid configuration was employed with a 75 kin CONUS domain and a 25 km grid over the southeastern United States. Initial conditions were obtained from the 12 UTC Early Eta Data Assimilation System analyses and lateral boundary conditions from the Early Eta forecast available at 3 hour intervals. The satellite-derived land surface temperature tendencies and insolation were assimilated between two and five hours (1400 and 1650 UTC) of the forecast. We performed the assimilation on the Southeastern domain only. In addition, a control run without assimilation was performed to provide insight into the performance of the assimilation technique.

Lapenta, William M.↗
DOE OSTI2025

Quantifying Uncertainties in Modeling Wind Resource Data from Different PBL Schemes in the WRF Model: A Case Study Over the Puerto Rico Region

This study examines the modeling uncertainty of wind resource data stemming from the use of various planetary boundary layer (PBL) parameterizations available in the Weather Research and Forecasting (WRF) model. WRF-based wind simulations spanning 20 years at 3-km resolution using 11 different PBL schemes are used to objectively investigate the uncertainty in modeling wind speed for land-based wind (LBW) and offshore wind (OSW) locations in Puerto Rico. The uncertainty in the wind modeling for the 20-year dataset is quantified using the spread index (SI) and standard deviation (SD). For virtual LBW and OSW sites, the SI and SD values are analyzed as calculated across various spatial and temporal scales. Because the PBL's atmospheric stability conditions can be characterized into two dominant categories, the study focuses on analyzing the SI and SD for daytime (mainly unstable PBL conditions) and nighttime (mainly stable PBL conditions). For wind shear (10 m-200 m) at the OSW and LBW sites, WRF-based numerical experiments indicate the following SI (or SD) ranges: 39%-94% (0.74 m/s-1.44 m/s) during the daytime for OSW, 50%-75% (0.68 m/s-1.19 m/s) during the daytime for LBW, 37%-60% (0.73 m/s-1.12 m/s) during the nighttime for OSW, and 57%-143 % (0.65 m/s-1.43 m/s) during the nighttime for LBW. While a high SI is observed when modeling LBW during the nighttime, there are notable modeling uncertainties during the daytime on the leeward side of the orographic barriers for Puerto Rico.

17 WIND ENERGY↗