Kalman Filter Data Assimilation Research Articles

Dynamics and Predictability of the Rapid Intensification of Super Typhoon Usagi (2013)

AbstractThis study explores the dynamics and predictability of the rapid intensification (RI) of Super Typhoon Usagi (2013) through a 60‐member convection‐permitting ensemble using the Weather Research and Forecasting (WRF) model and an ensemble Kalman filter (EnKF) data assimilation method. The surface maximum wind speed of Usagi, which was an intense category 4 western North Pacific tropical cyclone (TC), increased by 33 m s−1 over a 24‐hr period. The RI process was captured by the WRF simulation initialized with the global analysis but with a unique forecast challenge of early prediction. We improved the intensity forecasts by assimilating satellite‐derived atmospheric motion vectors into the WRF‐EnKF, which primarily reduced the strength of both the primary and secondary circulations in the TC vortex. Nevertheless, our ensemble forecasts initialized with the EnKF analysis ensemble predicted a significant spread in the intensity with considerable differences in the RI onset timing among individual members. Our analyses show that variation in the RI timing is most sensitive to differences in the initial TC vortex intensity and inner‐core moisture. Ensemble members with similar initial intensities but greater tropospheric moisture content exhibited earlier vortex axisymmetrization and consequently earlier RI. Further sensitivity experiments showed that variations in the inner‐core moisture content have an immediate impact on the structure and strength of inner‐core convection. These variations in inner‐core convection gradually caused differences in intensity between the TC vortices. In this study, we highlight the importance of accurate estimates of the inner‐core moisture content in the modeling and forecasting of TC intensity.

Estimating bamboo forest aboveground biomass using EnKF-assimilated MODIS LAI spatiotemporal data and machine learning algorithms

High-precision LAI (leaf area index) spatiotemporal data obtained from MODIS satellite remote sensing products are important for studying vegetation growth status, biomass carbon reserves, and the spatiotemporal dynamics of carbon cycling. LAI significantly influences biomass accumulation during the growth of bamboo forest in subtropical zones. Therefore, we applied the ensemble Kalman filter (EnKF) data assimilation algorithm to assimilate MODIS LAI products, and used assimilated LAI and the normalized difference vegetation index, enhanced vegetation index, simple ratio index as variables in the random forest model to estimate bamboo forest above ground biomass (AGB) in Zhejiang Province. Assimilated LAI spatiotemporal data using EnKF greatly improve the accuracy of MODIS LAI products, the R2 between assimilated and observed LAI was 0.92, and the RMSE was 0.37. Variations in the assimilated LAI time series were consistent with the seasonal dynamics of bamboo forest growth and had a significant effect on AGB. Moreover, the random forest model had strong predictive capabilities. A comparison of training and testing results produced accuracy (R) values for the random forest model using the assimilated LAI time series of 0.71 and 0.73, respectively. Using the assimilated LAI achieved a more accurate AGB estimate than using MODIS LAI time series products, as the R values were 54.3% and 58.7% higher, and the RMSE values were 19.2% and 19.1% lower for training and testing results, respectively. The calculated spatial distribution of bamboo forest AGB in Zhejiang province was consistent with the observed values. By combining assimilation technology of the MODIS LAI time series with the random forest model to more accurately estimate bamboo forest AGB in Zhejiang province, this study provided a new method for estimating large scale forest AGB based on low-resolution time series data.

Evaluation on surface current observing network of high frequency ground wave radars in the Gulf of Thailand

Due to the high cost of ocean observation system, the scientific design of observation network becomes much important. The current network of the high frequency radar system in the Gulf of Thailand has been studied using a three-dimensional coastal ocean model. At first, the observations from current radars have been assimilated into this coastal model and the forecast results have improved due to the data assimilation. But the results also show that further optimization of the observing network is necessary. And then, a series of experiments were carried out to assess the performance of the existing high frequency ground wave radar surface current observation system. The simulated surface current data in three regions were assimilated sequentially using an efficient ensemble Kalman filter data assimilation scheme. The experimental results showed that the coastal surface current observation system plays a positive role in improving the numerical simulation of the currents. Compared with the control experiment without assimilation, the simulation precision of surface and subsurface current had been improved after assimilated the surface currents observed at current networks. However, the improvement for three observing regions was quite different and current observing network in the Gulf of Thailand is not effective and a further optimization is required. Based on these evaluations, a manual scheme has been designed by discarding the redundant and inefficient locations and adding new stations where the performance after data assimilation is still low. For comparison, an objective scheme based on the idea of data assimilation has been obtained. Results show that all the two schemes of observing network perform better than the original network and optimal scheme-based data assimilation is much superior to the manual scheme that based on the evaluation of original observing network in the Gulf of Thailand. The distributions of the optimal network of radars could be a useful guidance for future design of observing system in this region.

Assimilation of GOES-13 Imager Clear-Sky Water Vapor (6.5 μm) Radiances into a Warn-on-Forecast System

A prototype convection-allowing system using the Advanced Research version of the Weather Research and Forecasting (WRF-ARW) Model and employing an ensemble Kalman filter (EnKF) data assimilation technique has been developed and used during the spring 2016 and 2017 Hazardous Weather Testbeds. This system assimilates WSR-88D reflectivity and radial velocity, geostationary satellite cloud water path (CWP) retrievals, and available surface observations over a regional domain with a 3-km horizontal resolution at 15-min intervals, with 3-km initial conditions provided by an experimental High-Resolution Rapid Refresh ensemble (HRRR-e). However, no information on upper-level thermodynamic conditions in cloud-free regions is currently assimilated, as few timely observations exist. One potential solution is to also assimilate clear-sky satellite radiances, which provide information on mid- and upper-tropospheric temperature and moisture conditions. This research assimilates GOES-13 imager water vapor band (6.5 μm) radiances using the GSI-EnKF system to take advantage of the Community Radiative Transfer Model (CRTM) integration. Results using four cases from May 2016 showed that assimilating radiances generally had a neutral-to-positive impact on the model analysis, reducing humidity bias and/or errors at the appropriate model levels where verification observations were present. The effects on high-impact weather forecasts, as verified against forecast reflectivity and updraft helicity, were mixed. Three cases (9, 22, and 24 May) showed some improvement in skill, while the other (25 May) performed worse, despite the improved environment. This research represents the first step in designing a high-resolution ensemble data assimilation system to use GOES-16 Advanced Baseline Imager data, which provides additional water vapor bands and increased spatial and temporal resolution.

Leveraging Google Earth Engine for Drought Assessment using Global Soil Moisture Data.

Soil moisture is considered a key variable to assess crop and drought conditions. However, readily available soil moisture datasets developed for monitoring agricultural drought conditions are uncommon. The aim of this work is to examine two global soil moisture data sets and a set of soil moisture web-based processing tools developed to demonstrate the value of the soil moisture data for drought monitoring and crop forecasting using Google Earth Engine (GEE). The two global soil moisture data sets discussed in the paper are generated by integrating Soil Moisture Ocean Salinity (SMOS) and Soil Moisture Active Passive (SMAP) satellite-derived observations into the modified two-layer Palmer model using a 1-D Ensemble Kalman Filter (EnKF) data assimilation approach. The web-based tools are designed to explore soil moisture variability as a function of land cover change and to easily estimate drought characteristics such as drought duration and intensity using soil moisture anomalies, and to inter-compare them against alternative drought indicators. To demonstrate the utility of these tools for agricultural drought monitoring, the soil moisture products, vegetation- and precipitation-based products are assessed over drought prone regions in South Africa and Ethiopia. Overall, the 3-month scale Standardized Precipitation Index (SPI) and Normalized Vegetation Index (NDVI) showed higher agreement with the root zone soil moisture anomalies. Soil moisture anomalies exhibited lower drought duration but higher intensity compare to SPIs. Inclusion of the global soil moisture data into GEE data catalog and the development of the web-based tools described in the paper enable a vast diversity of users to quickly and easily assess the impact of drought and improve planning related to drought risk assessment and early warning. GEE also improves the accessibility and usability of the earth observation data and related tools by making them available to a wide range of researchers and the public. In particular, the cloud-based nature of GEE is useful for providing access to the soil moisture data and scripts to users in developing countries that lack adequate observational soil moisture data or the necessary computational resources required to develop them.

An Assessment of the Impact of Land Thermal Infrared Observation on Regional Weather Forecasts Using Two Different Data Assimilation Approaches.

Recent studies have shown the unique value of satellite-observed land surface thermal infrared (TIR) information (e.g., skin temperature) and the feasibility of assimilating land surface temperature (LST) into land surface models (LSMs) to improve the simulation of land-atmosphere water and energy exchanges. In this study, two different types of LST assimilation techniques are implemented and the benefits from the techniques are compared. One of the techniques is to directly assimilate LST using ensemble Kalman filter (EnKF) data assimilation (DA) utilities. The other is to use the Atmosphere-Land Exchange Inversion model (ALEXI) as an “observation operator” that converts LST retrievals into the soil moisture (SM) proxy based on the ratio of actual to potential evapotranspiration (fPET), which is then assimilated into an LSM. While most current studies have shown some success in both directly the assimilating LST and assimilating ALEXI SM proxy into offline LSMs, the potential impact of the assimilation of TIR information through coupled numerical weather prediction (NWP) models is unclear. In this study, a semi-coupled Land Information System (LIS) and Weather Research and Forecast (WRF) system is employed to assess the impact of the two different techniques for assimilating the TIR observations from NOAA GOES satellites on WRF model forecasts. The NASA LIS, equipped with a variety of LSMs and advanced data assimilation tools (e.g., the ensemble Kalman Filter (EnKF)), takes atmospheric forcing data from the WRF model run, generates updated initial land surface conditions with the assimilation of either LST- or TIR-based SM and returns them to WRF for initializing the forecasts. The WRF forecasts using the daily updated initializations with the TIR data assimilation are evaluated against ground weather observations and re-analysis products. It is found that WRF forecasts with the LST-based SM assimilation have better agreement with the ground weather observations than those with the direct LST assimilation or without the land TIR data assimilation.

Ensemble Kalman Filter Data Assimilation for the Model for Prediction Across Scales (MPAS)

A global atmospheric analysis and forecast system is constructed based on the atmospheric component of the Model for Prediction Across Scales (MPAS-A) and the Data Assimilation Research Testbed (DART) ensemble Kalman filter. The system is constructed using the unstructured MPAS-A Voronoi (nominally hexagonal) mesh and thus facilitates multiscale analysis and forecasting without the need for developing new covariance models at different scales. Cycling experiments with the assimilation of real observations show that the global ensemble system is robust and reliable throughout a one-month period for both quasi-uniform and variable-resolution meshes. The variable-mesh assimilation system consistently provides higher-quality analyses than those from the coarse uniform mesh, in addition to the benefits of the higher-resolution forecasts, which leads to substantial improvements in 5-day forecasts. Using the fractions skill score, the spatial scale for skillful precipitation forecasts is evaluated over the high-resolution area of the variable-resolution mesh. Skill decreases more rapidly at smaller scales, but the variable mesh consistently outperforms the coarse uniform mesh in precipitation forecasts at all times and thresholds. Use of incremental analysis updates (IAU) greatly decreases high-frequency noise overall and improves the quality of EnKF analyses, particularly in the tropics. Important aspects of the system design related to the unstructured Voronoi mesh are also investigated, including algorithms for handling the C-grid staggered horizontal velocities.

Development of an ensemble Kalman filter data assimilation system for the Venusian atmosphere

The size and mass of Venus is similar to those of the Earth; however, its atmospheric dynamics are considerably different and they are poorly understood due to limited observations and computational difficulties. Here, we developed a data assimilation system based on the local ensemble transform Kalman filter (LETKF) for a Venusian Atmospheric GCM for the Earth Simulator (VAFES), to make full use of the observational data. To examine the validity of the system, two datasets were assimilated separately into the VAFES forecasts forced with solar heating that excludes the diurnal component Qz; one was created from a VAFES run forced with solar heating that includes the diurnal component Qt, whereas the other was based on observations made by the Venus Monitoring Camera (VMC) onboard the Venus Express. The VAFES-LETKF system rapidly reduced the errors between the analysis and forecasts. In addition, the VAFES-LETKF system successfully reproduced the thermal tide excited by the diurnal component of solar heating, even though the second datasets only included horizontal winds at a single altitude on the dayside with a long interval of approximately one Earth day. This advanced system could be useful in the analysis of future datasets from the Venus Climate Orbiter ‘Akatsuki’.

A modified iterative ensemble Kalman filter data assimilation method

High nonlinearity is a typical characteristic associated with data assimilation systems. Additionally, iterative ensemble based methods have attracted a large amount of research attention, which has been focused on dealing with nonlinearity problems. To solve the local convergence problem of the iterative ensemble Kalman filter, a modified iterative ensemble Kalman filter algorithm was put forward, which was based on a global convergence strategy from the perspective of a Gauss-Newton iteration. Through self-adaption, the step factor was adjusted to enable every iteration to approach expected values during the process of the data assimilation. A sensitivity experiment was carried out in a low dimensional Lorenz-63 chaotic system, as well as a Lorenz-96 model. The new method was tested via ensemble size, observation variance, and inflation factor changes, along with other aspects. Meanwhile, comparative research was conducted with both a traditional ensemble Kalman filter and an iterative ensemble Kalman filter. The results showed that the modified iterative ensemble Kalman filter algorithm was a data assimilation method that was able to effectively estimate a strongly nonlinear system state.

Unscented Kalman filter assimilation of time-lapse self-potential data for monitoring solute transport

Monitoring the extent and evolution of contaminant plumes in local and regional groundwater systems from existing landfills is critical in contamination control and remediation. The self-potential survey is an efficient and economical nondestructive geophysical technique that can be used to investigate underground contaminant plumes. Based on the unscented transform, we have built a Kalman filtering cycle to conduct time-lapse data assimilation for monitoring the transport of solute based on the solute transport experiment using a bench-scale physical model. The data assimilation was formed by modeling the evolution based on the random walk model and observation correcting based on the self-potential forward. Thus, monitoring self-potential data can be inverted by the data assimilation technique. As a result, we can reconstruct the dynamic process of the contaminant plume instead of using traditional frame-to-frame static inversion, which may cause inversion artifacts. The data assimilation inversion algorithm was evaluated through noise-added synthetic time-lapse self-potential data. The result of the numerical experiment shows validity, accuracy and tolerance to the noise of the dynamic inversion. To validate the proposed algorithm, we conducted a scaled-down sandbox self-potential observation experiment to generate time-lapse data that closely mimics the real-world contaminant monitoring setup. The results of physical experiments support the idea that the data assimilation method is a potentially useful approach for characterizing the transport of contamination plumes using the unscented Kalman filter (UKF) data assimilation technique applied to field time-lapse self-potential data.

Data-assimilated computational fluid dynamics modeling of convection-diffusion-reaction problems

This study focuses on the development and application of a data-assimilated multi-algorithm which combines a fourth-order finite-volume computational fluid dynamics (CFD) algorithm and the ensemble Kalman filter (EnKF) data assimilation algorithm for solving the problems involving important physics relevant to wave convection, molecular diffusion, and reaction of interest to engineering science. This data-assimilated CFD algorithm is applied to estimate an uncertain model parameter. By restricting the problems to one- and two-dimensional spatial space, this study allows us to develop a greater understanding of the integrated algorithm without the complexity or computational cost associated with three dimensions. Although scalar partial differential equations are used, fundamental issues in the algorithm development and application of EnKF to the physical processes occurring in a domain on engineering scales are sufficiently illuminated. Results of the one-dimensional convection-diffusion-reaction problem and the two-dimensional flame propagation demonstrate the validity of the data-assimilated CFD modeling system.

A model integrating longshore and cross‐shore processes for predicting long‐term shoreline response to climate change

AbstractWe present a shoreline change model for coastal hazard assessment and management planning. The model, CoSMoS‐COAST (Coastal One‐line Assimilated Simulation Tool), is a transect‐based, one‐line model that predicts short‐term and long‐term shoreline response to climate change in the 21st century. The proposed model represents a novel, modular synthesis of process‐based models of coastline evolution due to longshore and cross‐shore transport by waves and sea level rise. Additionally, the model uses an extended Kalman filter for data assimilation of historical shoreline positions to improve estimates of model parameters and thereby improve confidence in long‐term predictions. We apply CoSMoS‐COAST to simulate sandy shoreline evolution along 500 km of coastline in Southern California, which hosts complex mixtures of beach settings variably backed by dunes, bluffs, cliffs, estuaries, river mouths, and urban infrastructure, providing applicability of the model to virtually any coastal setting. Aided by data assimilation, the model is able to reproduce the observed signal of seasonal shoreline change for the hindcast period of 1995–2010, showing excellent agreement between modeled and observed beach states. The skill of the model during the hindcast period improves confidence in the model's predictive capability when applied to the forecast period (2010–2100) driven by GCM‐projected wave and sea level conditions. Predictions of shoreline change with limited human intervention indicate that 31% to 67% of Southern California beaches may become completely eroded by 2100 under sea level rise scenarios of 0.93 to 2.0 m.

Quality assessment of the TOPAZ4 reanalysis in the Arctic over the period 1991–2013

Abstract. Long dynamical atmospheric reanalyses are widely used for climate studies, but data-assimilative reanalyses of ocean and sea ice in the Arctic are less common. TOPAZ4 is a coupled ocean and sea ice data assimilation system for the North Atlantic and the Arctic that is based on the HYCOM ocean model and the ensemble Kalman filter data assimilation method using 100 dynamical members. A 23-year reanalysis has been completed for the period 1991–2013 and is the multi-year physical product in the Copernicus Marine Environment Monitoring Service (CMEMS) Arctic Marine Forecasting Center (ARC MFC). This study presents its quantitative quality assessment, compared to both assimilated and unassimilated observations available in the whole Arctic region, in order to document the strengths and weaknesses of the system for potential users. It is found that TOPAZ4 performs well with respect to near-surface ocean variables, but some limitations appear in the interior of the ocean and for ice thickness, where observations are sparse. In the course of the reanalysis, the skills of the system are improving as the observation network becomes denser, in particular during the International Polar Year. The online bias estimation successfully maintains a low bias in our system. In addition, statistics of the reduced centered random variables (RCRVs) confirm the reliability of the ensemble for most of the assimilated variables. Occasional discontinuities of these statistics are caused by the changes of the input data sets or the data assimilation settings, but the statistics remain otherwise stable throughout the reanalysis, regardless of the density of observations. Furthermore, no data type is severely less dispersed than the others, even though the lack of consistently reprocessed observation time series at the beginning of the reanalysis has proven challenging.

Impact of coastal radar observability on the forecast of the track and rainfall of Typhoon Morakot (2009) using WRF-based ensemble Kalman filter data assimilation

This study explored the impact of coastal radar observability on the forecast of the track and rainfall of Typhoon Morakot (2009) using a WRF-based ensemble Kalman filter (EnKF) data assimilation (DA) system. The results showed that the performance of radar EnKF DA was quite sensitive to the number of radars being assimilated and the DA timing relative to the landfall of the tropical cyclone (TC). It was found that assimilating radial velocity (Vr) data from all the four operational radars during the 6 h immediately before TC landfall was quite important for the track and rainfall forecasts after the TC made landfall. The TC track forecast error could be decreased by about 43% and the 24-h rainfall forecast skill could be almost tripled. Assimilating Vr data from a single radar outperformed the experiment without DA, though with less improvement compared to the multiple-radar DA experiment. Different forecast performances were obtained by assimilating different radars, which was closely related to the first-time wind analysis increment, the location of moisture transport, the quasi-stationary rainband, and the local convergence line. However, only assimilating Vr data when the TC was farther away from making landfall might worsen TC track and rainfall forecasts. Besides, this work also demonstrated that Vr data from multiple radars, instead of a single radar, should be used for verification to obtain a more reliable assessment of the EnKF performance.

Flow-dependent assimilation of sea surface temperature in isopycnal coordinates with the Norwegian Climate Prediction Model

We document a pilot stochastic re-analysis computed by assimilating sea surface temperature (SST) anomalies into the ocean component of the coupled Norwegian Climate Prediction Model (NorCPM) for the period 1950–2010 (doi: 10.11582/2016.00002). NorCPM is based on the Norwegian Earth System Model and uses the ensemble Kalman filter for data assimilation (DA). Here, we assimilate SST from the stochastic HadISST2 historical reconstruction. The accuracy, reliability and drift are investigated using both assimilated and independent observations. NorCPM is slightly overdispersive against assimilated observations but shows stable performance through the analysis period. It demonstrates skills against independent measurements: sea surface height, heat and salt content, in particular in the Equatorial and North Pacific, the North Atlantic Subpolar Gyre (SPG) region and the Nordic Seas. Furthermore, NorCPM provides a reliable monitoring of the SPG index and represents the vertical temperature variability there, in good agreement with observations. The monitoring of the Atlantic meridional overturning circulation is also encouraging. The benefit of using a flow-dependent assimilation method and constructing the covariance in isopycnal coordinates are investigated in the SPG region. Isopycnal coordinates discretisation is found to better capture the vertical structure than standard depth-coordinate discretisation, because it leads to a deeper influence of the assimilated surface observations. The vertical covariance shows a pronounced seasonal and decadal variability that highlights the benefit of flow-dependent DA method. This study demonstrates the potential of NorCPM to compute an ocean re-analysis for the 19th and 20th centuries when SST observations are available.

Impact of assimilating Taiwan’s coastal radar radial velocity on forecasting Typhoon Morakot (2009) in southeastern China using a WRF-based EnKF

This study explores for the first time the impact of assimilating radial velocity (Vr) observations from a single or multiple Taiwan’s coastal radars on tropical cyclone (TC) forecasting after landfall in the Chinese mainland by using a Weather Research and Forecasting model (WRF)-based ensemble Kalman filter (EnKF) data assimilation system. Typhoon Morakot (2009), which caused widespread damage in the southeastern coastal regions of the mainland after devastating Taiwan, was chosen as a case study. The results showed that assimilating Taiwan’s radar Vr data improved environmental field and steering flow and produced a more realistic TC position and structure in the final EnKF cycling analysis. Thus, the subsequent TC track and rainfall forecasts in southeastern China were improved. In addition, better observations of the TC inner core by Taiwan’s radar was a primary factor in improving TC rainfall forecast in the Chinese mainland.

Multivariate hydrological data assimilation of soil moisture and groundwater head

Abstract. Observed groundwater head and soil moisture profiles are assimilated into an integrated hydrological model. The study uses the ensemble transform Kalman filter (ETKF) data assimilation method with the MIKE SHE hydrological model code. The method was firstly tested on synthetic data in a catchment of less complexity (the Karup catchment in Denmark), and later implemented using data from real observations in a larger and more complex catchment (the Ahlergaarde catchment in Denmark). In the Karup model, several experiments were designed with respect to different observation types, ensemble sizes and localization schemes, to investigate the assimilation performance. The results showed the necessity of using localization, especially when assimilating both groundwater head and soil moisture. The proposed scheme with both distance localization and variable localization was shown to be more robust and provide better results. Using the same assimilation scheme in the Ahlergaarde model, groundwater head and soil moisture were successfully assimilated into the model. The hydrological model with assimilation showed an overall improved performance compared to the model without assimilation.

Mesoscale Hybrid Data Assimilation System based on JMA Nonhydrostatic Model

This paper discusses the benefits of using a hybrid ensemble Kalman filter and four-dimensional variational (4D-Var) data assimilation (DA) system rather than a 4D-Var system employing the National Meteorological Center (NMC, now known as NCEP) method (4D-Var-Bnmc) to predict severe weather events. An adjoint-based 4D-Var system was employed with a background error covariance matrix [Formula: see text] constructed from the NMC method and perturbations in a local ensemble transform Kalman filter system. The DA systems are based on the Japan Meteorological Agency’s nonhydrostatic model. To reduce the sampling noise, three types of implementation (the spatial localization, spectral localization, and neighboring ensemble approaches) were tested. The assimilation of a pseudosingle observation of sea level pressure located at a tropical cyclone (TC) center yielded analysis increments physically consistent with what is expected of a mature TC in the hybrid systems at the beginning of the assimilation window, whereas analogous experiments performed using the 4D-Var-Bnmc system did not. At the end, the structures of the 4D-Var-based increments became similar to one another, while the analysis increment by the 4D-Var-Bnmc system was broad in the horizontal direction. Realistic DA experiments showed that all of the hybrid systems provided initial conditions that yielded more accurate TC track and intensity forecasts than those achievable by the 4D-Var-Bnmc system. The hybrid systems also yielded some statistically significant improvements in forecasting local heavy rainfall events in terms of fraction skill scores when a 160 km × 160 km window size was used. The overall skills of the hybrid systems were relatively independent of the choice of implementation.

An integrated wind-forecast system based on the weather research and forecasting model, Kalman filter, and data assimilation with nacelle-wind observation

To forecast wind speed at hub-height is a challenging task for wind energy application, especially in Japan where the terrain feature is very complex and large fluctuations are observed in surface wind field. In this study, an integrated system to predict the hub-height wind speed has been developed by combining data assimilation and Kalman filter with the high resolution Weather Research and Forecasting (WRF) model. Assimilating the nacelle wind data (quality-controlled) and the Kalman filter algorithm effectively improves accuracy of the WRF model forecast by optimizing initial condition and post-processing the model output, respectively. It is found that the WRF model forecasts can be markedly improved after assimilating the nacelle wind data through the Gridpoint Statistical Interpolation analysis system, with the relative improvements of 34.3%, 23.9%, and 8.8% in ME (mean error), RMSE (root mean square error), and IA (index of agreement), respectively. The implementation of the Kalman filter can significantly reduce ME and RMSE while increases the value of IA as well. Further improvement can be achieved if the Kalman filter and nacelle wind data assimilation are implemented simultaneously. It is observed that the role of the Kalman filter is more dominant for the wind band of rated out speeds, while data assimilation is effective in reducing the random errors and becomes more important in rare or extreme weather conditions. Both data assimilation and Kalman filter modules apply the nacelle wind data which is routinely available, so the system can be easily adopted in different wind farm sites for operational use.

Scale-Dependent Representation of the Information Content of Observations in the Global Ensemble Kalman Filter Data Assimilation

Abstract Global data assimilation systems for numerical weather prediction (NWP) are characterized by significant uncertainties in tropical analysis fields. Furthermore, the largest spread of global ensemble forecasts in the short range on all scales is in the tropics. The presented results suggest that these properties hold even in the perfect-model framework and the ensemble Kalman filter data assimilation with a globally homogeneous network of wind and temperature profiles. The reasons for this are discussed by using the modal analysis, which provides information about the scale dependency of analysis and forecast uncertainties and information about the efficiency of data assimilation to reduce the prior uncertainties in the balanced and inertio-gravity dynamics. The scale-dependent representation of variance reduction of the prior ensemble by the data assimilation shows that the peak efficiency of data assimilation is on the synoptic scales in the midlatitudes that are associated with quasigeostrophic dynamics. In contrast, the variance associated with the inertia–gravity modes is less successfully reduced on all scales. A smaller information content of observations on planetary scales with respect to the synoptic scales is discussed in relation to the large-scale tropical uncertainties that current data assimilation methodologies do not address successfully. In addition, it is shown that a smaller reduction of the large-scale uncertainties in the prior state for NWP in the tropics than in the midlatitudes is influenced by the applied radius for the covariance localization.