Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “wind forecast error”

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 235 records · Page 13

Bias Characterization, Vertical Interpolation, and Horizontal Interpolation for Distributed Wind Siting Using Mesoscale Wind Resource Estimates

Much like their counterparts in utility-scale wind energy, developers of industrial, small-scale and distributed wind turbine deployments need to understand and accurately characterize the wind resource to properly assess the power generation and financial ramifications during siting and planning. National Renewable Energy Laboratory’s WIND (Wind Integration National Dataset) Toolkit (WTK) provides a best-in-class wind resource dataset generated using the Weather Research and Forecasting (WRF) model. This dataset includes parameters such as the wind speed, wind direction, and temperature at various heights, plus atmospheric stability near the surface. This data is available at 2-km spatial resolution and five-minute temporal resolution across 7 years, from 2007 to 2013 through a publicly accessible API interface. The Tools Assessing Performance (TAP) project seeks to extend this dataset to allow long term resource estimates and leverage it to better equip distributed wind equipment manufacturers, owner-operators, and installation professionals with better tools for practical siting applications. In this report, we present the results from our investigation within the TAP project focused on characterization of bias in WTK-based wind speed estimates and evaluation of vertical and horizontal interpolation techniques. We discuss the tradeoffs between different techniques and their combinations, as well as describe the lower bounds we determine for the studied validation errors. While the specific estimates we present are specific to WTK and the validation dataset we have chosen for this investigation (NREL's Wind Resource Meteorological Database), the overall analysis and the studied techniques are general enough to be applied to a broader set of wind datasets, both simulation-based and observational.

17 WIND ENERGY↗

Uncertainty of Observation Impact Estimation in an Adjoint Model Investigated with an Observing System Simulation Experiment

Adjoint models are often used to estimate the impact of different observations on short-term forecast skill. A common difficulty with the evaluation of short term forecast quality is the choice of verification fields. The use of self-analysis fields for verification is typical but incestuous, and introduces uncertainty due to biases and errors in the analysis field. In this study, an observing system simulation experiment (OSSE) is used to explore the uncertainty in adjoint model estimations of observation impact. The availability of the true state for verification in the OSSE framework in the form of the Nature Run allows calculation of the observation impact without the uncertainties present in self-analysis verification. These impact estimates are compared to estimates calculated using self-analysis verification. The Global Earth Observing System version 5 (GEOS-5) forecast model with Gridpoint Statistical Interpolation (GSI) is used with the National Aeronautics and Space Administration Global Modeling and Assimilation Office (NASA/GMAO) OSSE capability. The adjoint model includes moist processes, with total wet energy selected as the norm for evaluation of observation impacts. The results show that there are measurable but small discrepancies in the adjoint model estimation of observation impact. In general, observations of temperature and winds tend to have overestimated impacts with self-analysis verification, while observations of humidity and moisture-affected observations tend to have underestimated impacts. The small magnitude of the differences in impact estimates supports the robustness of the adjoint method of estimating observation impacts.

N C Prive↗

The relative influences of hydrologic information and dams’ hydropower scheduling decisions on electricity price forecasts

Price dynamics in wholesale electricity markets are driven by supply and demand. In markets with hydroelectric dams, the timing and amount of hydropower offered can influence prices in similar ways to wind and solar power. Unlike variable renewable energy, however, the supply of hydropower in wholesale markets is a function of both water availability and operational decisions at dams. Dam operators maximize revenues in wholesale markets by aligning generation with the periods of highest expected prices, and these scheduling decisions may in turn influence prices. Here, we examine the relative importance of two types of information in predicting forward electricity prices: a) water availability at dams, in the form of short-to-medium-range hydrological forecasts; and b) hourly scheduling decisions at dams. Using softly coupled hydrologic, hydropower scheduling, and power systems models spanning the U.S. Western Interconnection, we quantify the importance of hydrologic forecast accuracy in correctly predicting wholesale electricity prices and compare this with the influence of dam operators’ own hourly scheduling decisions on realized market prices. We find that aligning hydropower generation schedules with the periods of high forecasted prices causes larger, inadvertent price forecast errors than imperfect hydrologic forecasts. This suggests that knowledge of how water is managed by dam operators within the week is more important than weekly inflow forecast errors when predicting forward electricity prices. Our findings have implications for optimal hydropower scheduling by region. Specifically, accounting for price effects is critical in markets dominated by hydropower capacity.

Electricity markets↗

Contrasting effects of lake breeze and urbanization on heat stress in Chicago metropolitan area

This study used the latest Weather Research and Forecasting (WRF) model coupled with multi -layer urban canopy models to investigate contrasting effects from urbanization and lake breeze on summer heat stress over the Chicago metropolitan area (CMA). Comparisons between the model and in situ observations show that this coupled modeling system better captures urban locations' diurnal pattern of surface air temperature, skin temperature and relative humidity, with root mean square error reduced from 1.58 to 1.80 degrees C to 1.14-1.31 degrees C, 3.11-3.55 degrees C to 1.81-2.21 degrees C, and 10.73-11.35% to 7.84-8.60%, respectively, compared to WRF without coupling the urban canopy models. Two sensitivity experiments were conducted to isolate the influence of lake breeze and urbanization: one replaced the urban land use with cropland over the CMA, and the other filled all of Lake Michigan with cropland. Three different heat stress indices were computed to assess the uncertainties of heat stress response to changes in air temperature, relative humidity, and wind conditions. Results show that, when the lake has the largest cooling effect on air temperature, it also increases the relative humidity the most, and vice versa for urban warming and drying effects. Urbanization intensifies heat stress at night, and extends the heat caution period by up to 4 h over inland urban grids; the lake breeze relieves heat stress during afternoon (when the heat stress is the worst), and shortens the heat caution period by 1-3 h over inland urban grids and 3-4 h over coastal urban grids. In conclusion, the intensification of heat stress over the CMA due to urbanization is more than four times greater than the reduction from the lake breeze in the late afternoon and evening.

54 ENVIRONMENTAL SCIENCES↗

Improving Global Analysis and Short-Range Forecast Using Rainfall and Moisture Observations Derived from TRMM and SSM/I Passive Microwave Instruments

The Global Precipitation Mission, a satellite project under consideration as a follow-on to the Tropical Rainfall Measuring Mission (TRMM) by the National Aeronautics and Space Agency (NASA) in the United States, the National Space Development Agency (NASDA) in Japan, and other international partners, comprises an improved TRMM-like satellite and a constellation of 8 satellites carrying passive microwave radiometers to provide global rainfall measurements at 3-hour intervals. The success of this concept relies on the merits of rainfall estimates derived from passive microwave radiometers. This article offers a proof-of-concept demonstration of the benefits of using, rainfall and total precipitable water (TPW) information derived from such instruments in global data assimilation with observations from the TRMM Microwave Imager (TMI) and 2 Special Sensor Microwave/Imager (SSM/I) instruments. Global analyses that optimally combine observations from diverse sources with physical models of atmospheric and land processes can provide a comprehensive description of the climate systems. Currently, such data analyses contain significant errors in primary hydrological fields such as precipitation and evaporation, especially in the tropics. We show that assimilating the 6-h averaged TMI and SSM/I surface rainrate and TPW retrievals improves not only the hydrological cycle but also key climate parameters such as clouds, radiation, and the upper tropospheric moisture in the analysis produced by the Goddard Earth Observing System (GEOS) Data Assimilation System, as verified against radiation measurements by the Clouds and the Earth's Radiant Energy System (CERES) instrument and brightness temperature observations by the TIROS Operational Vertical Sounder (TOVS) instruments. Typically, rainfall assimilation improves clouds and radiation in areas of active convection, as well as the latent heating and large-scale motions in the tropics, while TPW assimilation leads to reduced moisture biases and improved radiative fluxes in clear-sky regions. Ensemble forecasts initialized with analyses that incorporate TMI and SSM/I rainfall and TPW data also yield better short-range predictions of geopotential heights, winds, and precipitation in the tropics. This study offers a compelling illustration of the potential of using rainfall and TPW information derived from passive microwave instruments to significantly improve the quality of 4-dimensional global datasets for climate analysis and weather forecasting applications.

Hou, Arthur Y.↗

Using Wind Observations from Nearby Aircraft to Update the Optimal Descent Trajectory in Real-Time

The ability to meet a controlled time of arrival during a continuous descent operation will enable environmentally friendly and fuel efficient descent operations while simultaneously maintaining airport throughput. However, if the wind forecast used to compute the initial trajectory plan is not accurate enough, the guidance system will need to correct time deviations from the plan during the execution of the descent. Previous work proposed an on-board guidance strategy based on model predictive control, which repeatedly updates the trajectory plan in real-time from the current aircraft state and for the remainder of the descent. However, the wind conditions downstream, at altitudes not explored yet, were difficult to predict due to the lack of data. This paper shows the potential benefits of using wind observations, broadcast by nearby aircraft, to reconstruct the wind profile downstream. The wind profile in the trajectory optimization problem is modeled as a spline, which control points are updated to fit the observations before re-planning the trajectory. Results from simulations using realistic wind data show that the performance of model predictive control significantly improves when including up-to-date wind observations, in terms of time and energy errors at the metering fix and fuel consumption.

Dalmau, Ramon↗

Evaluating Model Performance of an Ensemble-based Chemical Data Assimilation System During INTEX-B Field Mission

We present a global chemical data assimilation system using a global atmosphere model, the Community Atmosphere Model (CAM3) with simplified chemistry and the Data Assimilation Research Testbed (DART) assimilation package. DART is a community software facility for assimilation studies using the ensemble Kalman filter approach. Here, we apply the assimilation system to constrain global tropospheric carbon monoxide (CO) by assimilating meteorological observations of temperature and horizontal wind velocity and satellite CO retrievals from the Measurement of Pollution in the Troposphere (MOPITT) satellite instrument. We verify the system performance using independent CO observations taken on board the NSFINCAR C-130 and NASA DC-8 aircrafts during the April 2006 part of the Intercontinental Chemical Transport Experiment (INTEX-B). Our evaluations show that MOPITT data assimilation provides significant improvements in terms of capturing the observed CO variability relative to no MOPITT assimilation (i.e. the correlation improves from 0.62 to 0.71, significant at 99% confidence). The assimilation provides evidence of median CO loading of about 150 ppbv at 700 hPa over the NE Pacific during April 2006. This is marginally higher than the modeled CO with no MOPITT assimilation (-140 ppbv). Our ensemble-based estimates of model uncertainty also show model overprediction over the source region (i.e. China) and underprediction over the NE Pacific, suggesting model errors that cannot be readily explained by emissions alone. These results have important implications for improving regional chemical forecasts and for inverse modeling of CO sources and further demonstrate the utility of the assimilation system in comparing non-coincident measurements, e.g. comparing satellite retrievals of CO with in-situ aircraft measurements. The work described above also brought to light several short-comings of the data assimilation approach for CO profiles. Because of the limited vertical resolution of the measurement, the retrievals at different altitudes are correlated which can lead to problems with numerical error and overall efficiency. This has resulted in a manuscript that is about to be submitted to JGR:

Arellano, A. F., Jr.↗

The Vertical Error Characteristics of GOES-derived Winds: Description and Impact on Numerical Weather Prediction

Errors in the height assignment of some satellite-derived winds exist because the satellites sense radiation emitted from a finite layer of the atmosphere rather than a specific level. Potential problems in data assimilation may arise because the motion of a measured layer is often represented by a single-level value. In this research, cloud and water vapor motion winds that are derived from the Geostationary Operational Environmental Satellites (GOES winds) are compared to collocated rawinsonde observations (RAOBs). An important aspect of this work is that in addition to comparisons at each assigned height, the GOES winds are compared to the entire profile of the collocated RAOB data to determine the vertical error characteristics of the GOES winds. The impact of these results on numerical weather prediction is then investigated. The comparisons at individual vector height assignments indicate that the error of the GOES winds range from approx. 3 to 10 m/s and generally increase with height. However, if taken as a percentage of the total wind speed, accuracy is better at upper levels. As expected, comparisons with the entire profile of the collocated RAOBs indicate that clear-air water vapor winds represent deeper layers than do either infrared or water vapor cloud-tracked winds. This is because in cloud-free regions the signal from water vapor features may result from emittance over a thicker layer. To further investigate characteristics of the clear-air water vapor winds, they are stratified into two categories that are dependent on the depth of the layer represented by the vector. It is found that if the vertical gradient of moisture is smooth and uniform from near the height assignment upwards, the clear-air water vapor wind tends to represent a relatively deep layer. The information from the comparisons is then used in numerical model simulations of two separate events to determine the forecast impacts. Four simulations are performed for each case: 1) A control simulation that assimilates no satellite wind data, 2) assimilation of all GOES winds according to their assigned single level height, 3) assimilation of all GOES winds spread over multiple levels, and 4) assimilation of all GOES winds spread over multiple levels, but with variations in the vertical influence of clear-air water vapor winds based on the moisture profile in the model. In the first case, a strong mid-latitude cyclone is present and the use of the satellite data results in improved storm tracks during the initial approx. 36 h forecast period. This is because the satellite data improves the analysis of the environment into which the storm progresses. Statistics for mean wind vector and height differences show that, with the exception of the height field at later times in the first case, the use of GOES winds improves the simulation with time. The simulation results suggest that it is beneficial to spread the GOES wind information over multiple levels, particularly when the moisture profile is used to define the vertical influence.

Rao, P. Anil↗

Development of an OSSE Framework for a Global Atmospheric Data Assimilation System

Observing system simulation experiments (OSSEs) are powerful tools for estimating the usefulness of various configurations of envisioned observing systems and data assimilation techniques. Their utility stems from their being conducted in an entirely simulated context, utilizing simulated observations having simulated errors and drawn from a simulation of the earth's environment. Observations are generated by applying physically based algorithms to the simulated state, such as performed during data assimilation or using other appropriate algorithms. Adding realistic instrument plus representativeness errors, including their biases and correlations, can be critical for obtaining realistic assessments of the impact of a proposed observing system or analysis technique. If estimates of the expected accuracy of proposed observations are realistic, then the OSSE can be also used to learn how best to utilize the new information, accelerating its transition to operations once the real data are available. As with any inferences from simulations, however, it is first imperative that some baseline OSSEs are performed and well validated against corresponding results obtained with a real observing system. This talk provides an overview of, and highlights critical issues related to, the development of an OSSE framework for the tropospheric weather prediction component of the NASA GEOS-5 global atmospheric data assimilation system. The framework includes all existing observations having significant impact on short-term forecast skill. Its validity has been carefully assessed using a range of metrics that can be evaluated in both the OSSE and real contexts, including adjoint-based estimates of observation impact. A preliminary application to the Aeolus Doppler wind lidar mission, scheduled for launch by the European Space Agency in 2014, has also been investigated.

Gelaro, Ronald↗

Impact of Bias in the Marine Air Temperature Observation Set on Atmospheric Reanalyses

Marine air temperature observations are collected by volunteer observing ships carrying automatic meteorological sensor packages. Close proximity of the package mount to a dark colored ship's deck is known to introduce a warm bias on sunny afternoons in summer under low wind conditions. This observation set is used as an influential constraint in the MERRA-2 reanalysis, as well as several other reanalyses, in the poorly sampled surface marine boundary layer. Here we explore the consequences of bias in the marine air temperature observation set during the 36 year period 1980-2015. Marine air temperature observations have mainly been collected along major shipping routes, with the densest coverage in the North Atlantic and North Pacific. In the first part of this presentation we review the historical observation coverage and compare air temperature and contemporaneous SST measurements to identify systematic differences. The comparison suggests the presence of diurnal bias in the air temperature observations peaking in subtropics and midlatitudes latitudes in mid-afternoon local time. In the second part we explore the impact of the marine air temperature observations on MERRA-2 through analysis of the terms in the temperature forecast tendency equation, the equation which describes the contributions to changes in temperature during each analysis cycle. Examination of the temperature forecast tendency equation allows us to separate systematic errors due to such processes and radiation, mixing, and advection and allows us to identify the major source of systematic forecast tendencies. The results show that the marine air temperature measurements introduce an error that can reach 0.4K/dy which also peaks in local afternoon at subtropical latitudes. This bias varies with season and has gradually reduced from the beginning of the MERRA-2 reanalysis in 1980. In the final part we conduct numerical experiments with a version of the MERRA-2 analysis system, comparing experiments when the marine air temperature observations are included or withheld from the analysis system.

Akella, Santharam↗

Mesoscale to Microscale Simulations over Complex Terrain with the Immersed Boundary Method in the Weather Research and Forecasting Model

Improvements to the Weather Research and Forecasting (WRF) Model are made to enable multiscale simulations over highly complex terrain with dynamically downscaled boundary conditions from the mesoscale to the microscale. Over steep terrain, the WRF Model develops numerical errors that are due to grid deformation of the terrain-following coordinates. An alternative coordinate system, the immersed boundary method (IBM), has been implemented into WRF, allowing for simulations over highly complex terrain; however, the new coordinate system precluded nesting within mesoscale simulations using WRF’s native terrain-following coordinates. Here, the immersed boundary method and WRF’s grid-nesting framework are modified to seamlessly work together. This improved framework for the first time allows for large-eddy simulation over complex (urban) terrain with IBM to be nested within a typical mesoscale WRF simulation. Simulations of the Joint Urban 2003 field campaign in Oklahoma City, Oklahoma, are performed using a multiscale five-domain nested configuration, spanning horizontal grid resolutions from 6 km to 2 m. These are compared with microscale-only simulations with idealized lateral boundary conditions and with observations of wind speed/direction and SF 6 concentrations from a controlled release from intensive observation period 3. The multiscale simulation, which is configured independent of local observations, shows similar model skill predicting wind speed/direction and improved skill predicting SF 6 concentrations when compared with the idealized simulations, which require use of observations to set mean flow conditions. Use of this improved multiscale framework shows promise for enabling large-eddy simulation over highly complex terrain with dynamically downscaled boundary conditions from mesoscale models.

42 ENGINEERING↗

Understanding Shock Dynamics in the Inner Heliosphere with Modeling and Type II Radio Data: the 2010-04-03 Event

The 2010 April 03 solar event was studied using observations from STEREO SECCHI, SOHO LASCO, and Wind kilometric Type II data (kmTII) combined with WSA-Cone-ENLIL model simulations performed at the Community Coordinated Modeling Center (CCMC). In particular, we identified the origin of the coronal mass ejection (CME) using STEREO EUVI and SOHO EIT images. A flux-rope model was fit to the SECCHI A and B, and LASCO images to determine the CMEs direction, size, and actual speed. J-maps from STEREO COR2HI-1HI-2 and simulations fromCCMC were used to study the formation and evolution of the shock in the inner heliosphere. In addition, we also studied the time-distance profile of the shock propagation from kmTII radio burst observations. The J-maps together with in-situ datafrom the Wind spacecraft provided an opportunity to validate the simulation results andthe kmTII prediction. Here we report on a comparison of two methods of predictinginterplanetary shock arrival time: the ENLIL model and the kmTII method; andinvestigate whether or not using the ENLIL model density improves the kmTIIprediction. We found that the ENLIL model predicted the kinematics of shock evolutionwell. The shock arrival times (SAT) and linear-fit shock velocities in the ENLILmodel agreed well with those measurements in the J-maps along both the CME leading edge and the Sun-Earth line. The ENLIL model also reproduced most of the largescale structures of the shock propagation and gave the SAT prediction at Earth with an error of 17 hours. The kmTII method predicted the SAT at Earth with an error of 15 hours when using n0 4.16 cm3, the ENLIL model plasma density near Earth; but itimproved to 2 hours when using n0 6.64 cm3, the model density near the CMEleading edge at 1 AU.

Space weather forecast↗

Collaborative Research: Advancing Arctic Climate Projection Capability at Seasonal to Decadal Scales (Final Technical Report)

The Regional Arctic System Model (RASM) at process resolving configurations has been used to (i) advance understanding of physical processes and feedbacks involved in Arctic amplification and (ii) understand and potentially reduce uncertainty in prediction of arctic climate change at seasonal to decadal scales. RASM consists the atmosphere (Weather and Research Forecasting model, WRF), ocean (Parallel Ocean Program, POP), sea ice (CICE), land hydrology (Variable Infiltration Capacity model, VIC), river routing scheme (RVIC), marine biogeochemistry components and the coupling framework (CPL7). Its domain is pan-Arctic, with the atmosphere and land components configured on a 50-km or 25-km grid and four configurations of the ocean and sea ice components: 1/12°(~9.3km) or 1/48°(~2.4km) and 45 or 60 vertical layers. These RASM configurations have been motivated by the emerging exascale capability for high performance computing to improve model fidelity. The dynamical downscaling of reanalysis allows comparison of RASM results with observations in place and time to: (i) advance system level understanding of physical processes and coupling involved in an event, (ii) optimize model parameter space, (iii) diagnose and reduce model biases and (iv) produce realistic and consistent across all the components initial conditions for predictions and predictability studies, which are all unique capabilities not available in global Earth System Models (ESMs). An evaluation of RASM 1.0 (Cassano et al. 2017) revealed that it had a cold bias over the oceans and a warm bias over land areas due largely to cloud and radiation biases in the model, with too little cloud cover simulated over land and too much cloud cover simulated over sub-polar oceans. This study has motivated an upgrade to WRF version 3.7.1 in RASM and allowed for the inclusion of the radiative impact of convective clouds. A variety of atmospheric physics parameterizations were evaluated against observations (e.g. data from the Arctic Clouds in Summer Experiment (ACSE); Sedlar et al. 2020) to identify an optimal suite of WRF physics options in RASM. The RASM with the optimized WRF physics were used to study the impact of strong mesoscale winds over the ocean around the southern tip of Greenland (DuVivier and Cassano 2016) and their impact on oceanic convection (DuVivier et al. 2017a). Data from the PolarWinds field campaign were used to evaluate WRF boundary layer physics and resolution impacts on the simulation of a Greenland barrier wind event (DuVivier et al. 2017b). The RVIC streamflow routing model has been implemented in RASM to realistically represent high-resolution streamflow processes (Hamman et al. 2017) and to couple the land buoyancy fluxes to the ocean. The RASM-RVIC high-resolution data set of all coastal freshwater fluxes for the Arctic drainage basin and surrounding areas for 1979-2014 was published as a separate product (https://doi.org/10.5281/zenodo.293037). The fidelity of atmospheric momentum transfer to and the response of polar marine Ekman layer in RASM and Community Earth System Model (CESM) was investigated by Roberts et al. (2015). The increased frequency of oceanic flux exchange in CESM, following the RASM guidance, caused a considerable increase in the median inertial ice speed across the Southern Ocean and parts of the Arctic. A comprehensive evaluation of the RASM1.0 atmosphere-ocean-sea ice-land interface was completed by Brunke et al. (2018). RASM was also demonstrated for its capability to simulate extreme events in agreement with observations in space and time (Lee et al. submitted). In particular, the development of three open water events, known as polynyas, have been simulated north of Greenland in February of 2011, 2017 and 2018, in agreement with satellite observations for the past four decades. The optimized RASM sea ice results have been favorably evaluated against satellite observations and a subset of eleven CMIP6 models (Watts et al. submitted). In a complementary project, Jin et al. (2018) have shown that RASM with higher-resolution and new sea-ice processes contributed to lower model errors in sea-ice conditions, concentrations of nutrients and ice algae, in comparison to results from the coarse-resolution (1°) CESM. In two other complementary studies, RASM results were used (i) to explain areas of concentrated use by bowhead whales, the seasonal progression in the use, and the physical environment within those areas (Citta et al. 2015) and (ii) for a synthesis of fall bowhead whales distribution and migration in the Bering-Chukchi-Beaufort (BCB) Sea to investigate whale movements and feeding to the local ocean hydrography and currents (Citta et al. 2018). However, the multi-decadal output from the CESM Large Ensemble yielded unrealistic forcing. Instead, the shorter NCEP CFSv2 9-month forecasts were successfully tested and afforded an increased ensemble size (~30) to demonstrate gains of dynamical downscaling at sub-seasonal to intra-annual time scales. The improved model physics and coupling among RASM model components have yielded more realistic representation of the sea ice cover and consistent across all model components initial conditions. Consequently, RASM demonstrates significant gains compared to simulation of sea ice in the NCEP reanalysis. In addition, RASM 6-month ensemble forecasts yield very realistic sea ice distribution, which demonstrates both significant gains of dynamical downscaling and the continued impact of the initial conditions on forecasts out to 6 months (https://nps.edu/web/rasm/predictions). A follow up study using RASM for dynamical downscaling of the more realistic CESM initialized Decadal Prediction Large Ensemble output is currently ongoing as part of the DOE RGMA HiLAT-RASM project.

54 ENVIRONMENTAL SCIENCES↗

Daily grass reference evapotranspiration with Meteosat Second Generation shortwave radiation and reference ET products

Here we assess the accuracy of estimating daily grass reference evapotranspiration (PM-ET o ) using daily shortwave radiation (R s ) and reference evapotranspiration (ET REF ) products provided by the Meteosat Second Generation (MSG) geostationary satellite delivered by the Satellite Applications Facility on Land Surface Analysis (LSA-SAF) framework. The accuracy of using reanalysis ERA5 shortwave radiation data (R s ERA5 ) provided by the European Center for Medium-Range Weather Forecasts (ECMWF) is also evaluated. The assessments were performed using observed weather variables at 37 weather stations distributed across continental Portugal, where climate conditions range from semi-arid to humid, and 12 weather stations located in Azores islands, characterized by humid, windy and often cloudy conditions. This study&#x27;s use of data from a variety of climate conditions contributed to a unique and innovative assessment of the usability of LSA-SAF and ERA5 products for ET o estimation. The first assessment focused on comparing LSA-SAF estimates of R s (R s LSA-SAF ) against ground stations (R s ground ). The results showed a good matching between the two R s data sets for continental Portugal but a tendency for R s LSA-SAF to under-estimate R s ground in the cloudy islands of Azores. ET o values computed using R s LSA-SAF data and observed temperature, humidity and wind speed (ET o LSA-SAF ) were then compared with PM-ET o estimates with ground-based data, which were used as benchmark; input data of temperature and humidity needed for PM-ET o were quality checked for surface aridity effects. It was observed that ET o LSA-SAF is strongly correlated with PM-ET o (R 2 > 0.97) for most locations in continental Portugal, with regression coefficient of a linear regression forced to the origin ranging between 0.95 and 1.05, mean root mean square error (RMSE) of 0.13 mm d –1 , and Nash and Sutcliff efficiency of modeling (EF) above 0.95. For most Azores locations, ET o LSA-SAF over-estimated PM-ET o . This is likely a consequence of the high spatio-temporal heterogeneity of weather conditions that occur in these oceanic islands together with the different footprints of satellite (averaged over the pixel) and station observations. Reanalysis ERA5 shortwave radiation data presented similar behavior to the LSA-SAF products, however with slightly lower accuracy. The daily LSA-SAF ET REF product (ET REF LSA-SAF ) was assessed and results have shown a good accuracy of this product, with acceptable RMSE and high EF values, for continental Portugal but a low accuracy for the Azores islands. A simplified bias correction approach was shown to improve both ET o derived from the LSA-SAF products, namely for Azores stations, which seem to be representative of smaller areas. The use of the FAO-PM temperature approach (PMT) was also assessed using the R s LSA-SAF and R s ERA5 data, which showed a superiority of the LSA-SAF product for ET o estimations (ET o PMT LSA-SAF ). No significant differences (p < 0.05) were observed in terms of the median value of the RMSE when adopting ET o PMT and ET REF LSA-SAF . Differently, results showed that using the R s LSA-SAF in the PMT approach (ET o PMT LSA-SAF ) produces significantly better RMSE results than ET o PMT and ET REF LSA-SAF . Overall, the performed assessment allows concluding that the use of R s LSA-SAF, and to a lesser extent the use of the R s ERA5 , highly improves the accuracy of computation of ET o when R s observations are not available, including when only temperature data are accessible. The use of the ET REF LSA-SAF product is a good alternative when observed weather data are not available.

54 ENVIRONMENTAL SCIENCES↗

Radar sensitivity and antenna scan pattern study for a satellite-based Radar Wind Sounder (RAWS)

Modeling global atmospheric circulations and forecasting the weather would improve greatly if worldwide information on winds aloft were available. Recognition of this led to the inclusion of the LAser Wind Sounder (LAWS) system to measure Doppler shifts from aerosols in the planned for Earth Observation System (EOS). However, gaps will exist in LAWS coverage where heavy clouds are present. The RAdar Wind Sensor (RAWS) is an instrument that could fill these gaps by measuring Doppler shifts from clouds and rain. Previous studies conducted at the University of Kansas show RAWS as a feasible instrument. This thesis pertains to the signal-to-noise ratio (SNR) sensitivity, transmit waveform, and limitations to the antenna scan pattern of the RAWS system. A dop-size distribution model is selected and applied to the radar range equation for the sensitivity analysis. Six frequencies are used in computing the SNR for several cloud types to determine the optimal transmit frequency. the results show the use of two frequencies, one higher (94 GHz) to obtain sensitivity for thinner cloud, and a lower frequency (24 GHz) to obtain sensitivity for thinner cloud, and a lower frequency (24 GHz) for better penetration in rain, provide ample SNR. The waveform design supports covariance estimation processing. This estimator eliminates the Doppler ambiguities compounded by the selection of such high transmit frequencies, while providing an estimate of the mean frequency. the unambiguous range and velocity computation shows them to be within acceptable limits. The design goal for the RAWS system is to limit the wind-speed error to less than 1 ms(exp -1). Due to linear dependence between vectors for a three-vector scan pattern, a reasonable wind-speed error is unattainable. Only the two-vector scan pattern falls within the wind-error limits for azimuth angles between 16 deg to 70 deg. However, this scan only allows two components of the wind to be determined. As a result, a technique is then shown, based on the Z-R-V relationships, that permit the vertical component (i.e., rain) to be computed. Thus the horizontal wind components may be obtained form the covariance estimator and the vertical component from the reflectivity factor. Finally, a new candidate system is introduced which summarizes the parameters taken from previous RAWS studies, or those modified in this thesis.

Stuart, Michael A.↗

Verification of Real-Time WSA-ENLIL+Cone Simulations of CME Arrival-Time at the CCMC from 2010 to 2016

The Wang-Sheeley-Arge (WSA)-ENLIL+Cone model is used extensively in space weather operations world-wide to model coronal mass ejection (CME) propagation. As such, it is important to assess its performance. We present validation results of the WSA-ENLIL+Cone model installed at the Community Coordinated Modeling Center (CCMC) and executed in real-time by the CCMC space weather team. CCMC uses the WSA-ENLIL+Cone model to predict CME arrivals at NASA missions throughout the inner heliosphere. In this work we compare model predicted CME arrival-times to in situ interplanetary coronal mass ejection leading edge measurements at Solar TErrestrial RElations Observatory-Ahead (STEREO-A), Solar TErrestrial RElations Observatory-Behind (STEREO-B), and Earth (Wind and ACE) for simulations completed between March 2010 and December 2016 (over 1,800 CMEs). We report hit, miss, false alarm, and correct rejection statistics for all three locations. For all predicted CME arrivals, the hit rate is 0.5, and the false alarm rate is 0.1. For the 273 events where the CME was predicted to arrive at Earth, STEREO-A, or STEREO-B, and was actually observed (hit event), the mean absolute arrival-time prediction error was 10.4 +/- 0.9 h, with a tendency to early prediction error of -4.0 h. We show the dependence of the arrival-time error on CME input parameters. We also explore the impact of the multi-spacecraft observations used to initialize the model CME inputs by comparing model verification results before and after the STEREO-B communication loss (since September 2014) and STEREO-A sidelobe operations (August 2014- December 2015). There is an increase of 1.7 h in the CME arrival time error during single, or limited two-viewpoint periods, compared to the three-spacecraft viewpoint period. This trend would apply to a future space weather mission at L5 or L4 as another coronagraph viewpoint to reduce CME arrival time errors compared to a single L1 viewpoint.

MHD / modeling / validation / forecasting / corona↗

Assimilation of Satellite-Derived Soil Moisture for Improved Forecasts of the Great Plains Low-Level Jet

In the context of forecasting societally impactful Great Plains low-level jets (GPLLJs), the potential added value of satellite soil moisture (SM) data assimilation (DA) is high. GPLLJs are both sensitive to regional soil moisture gradients and frequent drivers of severe weather, including mesoscale convective systems. An untested hypothesis is that SM DA is more effective in forecasts of weakly synoptically forced, or uncoupled GPLLJs, than in forecasts of cyclone-induced coupled GPLLJs. Using the NASA Unified Weather Research and Forecasting (NU-WRF) Model, 75 GPLLJs are simulated at 9-km resolution both with and without NASA Soil Moisture Active Passive SM DA. Differences in modeled SM, surface sensible (SH) and latent heat (LH) fluxes, 2-m temperature (T2), 2-m humidity (Q2), PBL height (PBLH), and 850-hPa wind speed (W850) are quantified for individual jets and jet-type event subsets over the south-central Great Plains, as well as separately for each GPLLJ sector (entrance, core, and exit). At the GPLLJ core, DA-related changes of up to 5.4 kg m–2 in SM can result in T2, Q2, LH, SH, PBLH, and W850 differences of 0.68°C, 0.71 g kg –2 , 59.9 W m –2 , 52.4 W m –2 , 240 m, and 4 m s –1 , respectively. W850 differences focus along the jet axis and tend to increase from south to north. Further, jet-type differences are most evident at the GPLLJ exit where DA increases and decreases W850 in uncoupled and coupled GPLLJs, respectively. Data assimilation marginally reduces negative wind speed bias for all jets, but the correction is greater for uncoupled GPLLJs, as hypothesized.

54 ENVIRONMENTAL SCIENCES↗

Gap Flows Along the Columbia River Observed During the WFIP2 Field Campaign

Data collected in the Columbia River Gorge and Basin Area during the second Wind Forecasting Improvement Project (WFIP2) was used to study gap flow events in the region. In particular, the relationship between 200 m zonal wind speed and the zonal mesoscale pressure gradient at Troutdale (West of the gorge), Wasco (east of the gorge), and Boardman (farther east of the gorge) sites was examined. An objective criterion was used to identify 169 (Troutdale), 161 (Wasco), and 113 (Boardman) gap flow days that had weak synoptic forcing. The 200 m zonal wind speed exhibited a linear relationship with the zonal pressure gradient at Wasco during both easterly and westerly gap flow days at daily timescales. Such a relationship was only observed during easterly gap flow days at Troutdale and was not observed at Boardman. The zonal wind speed exhibited a diurnal cycle at Wasco and Boardman, with greater changes in the nighttime winds than the daytime winds due to changes in pressure gradient. The root mean square differences (RMSDs) between the observed pressure gradient at Wasco and that simulated by the High Resolution Rapid Refresh (HRRR) model in 3 km (HRRR-3km) and 750 m (HRRR-750m) horizontal resolution were 3.14 hPa and 0.57 hPa, respectively. However, the RMSDs for the 200 m zonal wind speed were 1.40 m s -1 for HRRR-3km and 1.20 m s -1 for HRRR- 750m. Improvement in simulating the zonal pressure gradient accompanied by negligible changes in simulating zonal wind speed due to higher horizontal resolution points to a complex interplay between different components of the model physics when simulating the meteorology. The HRRR-750m resolution model accurately mimicked the observed covariability (phase diagram) between the zonal pressure gradient and the 200 m zonal wind speed during the nighttime, but not during the daytime, suggesting model errors might be related to the representation of boundary layer processes.

54 ENVIRONMENTAL SCIENCES↗