Search

The rank histogram filter (RHF) and the ensemble Kalman filter (EnKF) are assessed for soil moisture estimation using perfect model (identical twin) synthetic data assimilation experiments. The primary motivation is to gauge the impact on analysis quality attributable to the consideration of non-Gaussian forecast error distributions. Using the NASA Catchment land surface model, the two filters are compared at 18 globally-distributed single-catchment locations for a 10-year experiment period. It is shown that both filters yield adequate estimates of soil moisture, with the RHF having a small but significant performance advantage. Most notably, the RHF systematically increases the normalized information contribution (NIC) score of the mean absolute bias by 0.05 over that of the EnKF for surface, root-zone and profile soil moisture. The RHF also increases the NIC score for the anomaly correlation of surface soil moisture by 0.02 over that of the EnKF (at a 5% significance level). Results also demonstrate that the performance of both filters is somewhat improved when the ensemble priors are adaptively inflated to offset the negative effects of systematic errors.

Much of solar activity within a sunspot cycle occurs as bursts, or 'seasons' of strong activity over several months, separated by periods of much less activity. The most important space weather effects occur during these bursts. Previous modeling and forecasting efforts have focused on time-scales of hours-to-days and decades-to-centuries. The recent discovery of Rossby waves in the Sun, together with recently developed global models of solar MHD Rossby waves and their interactions with differential rotation and spot-producing magnetic fields, reveal the opportunity to simulate and predict the occurrence, strength and location of enhanced activity bursts a few weeks up to several months in advance. We now have a golden opportunity to fill in this gap in time-scales of forecasting space weather. This requires a) continuous observations of solar Rossby waves by various techniques; b) development of coupled nonlinear MHD models that simulate both global Rossby waves and the much smaller spatial scale emergence of new active regions; c) application of advanced data assimilation techniques to couple surface observations to update the model-system to integrate forward in time for creating forecasts months ahead. Then it will be possible to build operational prediction models to meet the needs of customers and stakeholders, including support of future NASA missions, regarding what kind and level of space weather to expect a few weeks to several months ahead.

Lidar backscatter and wind retrievals of the planetary boundary layer height (PBLH) are assimilated into 22 hourly forecasts from the NASA Unified - Weather and Research Forecast (NU-WRF) model during the Plains Elevated Convection Convection at Night (PECAN) campaign on July 11, 2015 in Greensburg, Kansas, using error statistics collected from the model profiles to compute the necessary covariance matrices. Two separate forecast runs using different PBL physics schemes were employed, and comparisons with 6 independent radiosonde profiles were made for each run. Both of the forecast runs accurately predicted the PBLH and the state variable profiles within the planetary boundary layer during the early morning, and the assimilation had a small impact during this time. In the late afternoon, the forecast runs showed decreased accuracy as the convective boundary layer developed. However, assimilation of the Doppler lidar PBLH observations were found to improve the temperature and V velocity profiles relative to independent radiosonde profiles. Water vapor was overcorrected, leading to increased differences with independent data. Errors in the U velocity were made slightly larger. The computed forecast error covariances between the PBLH and state variables were found to rise in the late afternoon, leading to the larger improvements in the afternoon. This work represents the first effort to assimilate PBLH into forecast states using ensemble methods.