Browsing by Subject "data assimilation"
Now showing 1 - 3 of 3
Results Per Page
Sort Options
Item Ensemble Statistics and Error Covariance of a Rapidly Intensifying Hurricane(2010-01-16) Rigney, Matthew C.This thesis presents an investigation of ensemble Gaussianity, the effect of non- Gaussianity on covariance structures, storm-centered data assimilation techniques, and the relationship between commonly used data assimilation variables and the underlying dynamics for the case of Hurricane Humberto. Using an Ensemble Kalman Filter (EnKF), a comparison of data assimilation results in Storm-centered and Eulerian coordinate systems is made. In addition, the extent of the non-Gaussianity of the model ensemble is investigated and quantified. The effect of this non-Gaussianity on covariance structures, which play an integral role in the EnKF data assimilation scheme, is then explored. Finally, the correlation structures calculated from a Weather Research Forecast (WRF) ensemble forecast of several state variables are investigated in order to better understand the dynamics of this rapidly intensifying cyclone. Hurricane Humberto rapidly intensified in the northwestern Gulf of Mexico from a tropical disturbance to a strong category one hurricane with 90 mph winds in 24 hours. Numerical models did not capture the intensification of Humberto well. This could be due in large part to initial condition error, which can be addressed by data assimilation schemes. Because the EnKF scheme is a linear theory developed on the assumption of the normality of the ensemble distribution, non-Gaussianity in the ensemble distribution used could affect the EnKF update. It is shown that multiple state variables do indeed show significant non-Gaussianity through an inspection of statistical moments. In addition, storm-centered data assimilation schemes present an alternative to traditional Eulerian schemes by emphasizing the centrality of the cyclone to the assimilation window. This allows for an update that is most effective in the vicinity of the storm center, which is of most concern in mesoscale events such as Humberto. Finally, the effect of non-Gaussian distributions on covariance structures is examined through data transformations of normal distributions. Various standard transformations of two Gaussian distributions are made. Skewness, kurtosis, and correlation between the two distributions are taken before and after the transformations. It can be seen that there is a relationship between a change in skewness and kurtosis and the correlation between the distributions. These effects are then taken into consideration as the dynamics contributing to the rapid intensification of Humberto are explored through correlation structures.Item Mesoscale ensemble-based data assimilation and parameter estimation(Texas A&M University, 2005-11-01) Aksoy, AltugThe performance of the ensemble Kalman filter (EnKF) in forced, dissipative flow under imperfect model conditions is investigated through simultaneous state and parameter estimation where the source of model error is the uncertainty in the model parameters. Two numerical models with increasing complexity are used with simulated observations. For lower complexity, a two-dimensional, nonlinear, hydrostatic, non-rotating, and incompressible sea breeze model is developed with buoyancy and vorticity as the prognostic variables. Model resolution is 4 km horizontally and 50 m vertically. The ensemble size is set at 40. Forcing is maintained through an explicit heating function with additive stochastic noise. Simulated buoyancy observations on land surface with 40-km spacing are assimilated every 3 hours. Up to six model parameters are successfully subjected to estimation attempts in various experiments. The overall EnKF performance in terms of the error statistics is found to be superior to the worst-case scenario (when there is parameter error but no parameter estimation is performed) with an average error reduction in buoyancy and vorticity of 40% and 46%, respectively, for the simultaneous estimation of six parameters. The model chosen to represent the complexity of operational weather forecasting is the Pennsylvania State University-National Center for Atmospheric Research MM5 model with a 36-km horizontal resolution and 43 vertical layers. The ensemble size for all experiments is chosen as 40 and a 41st member is generated as the truth with the same ensemble statistics. Assimilations are performed with a 12-hour interval with simulated sounding and surface observations of horizontal winds and temperature. Only single-parameter experiments are performed focusing on a constant inserted into the code as the multiplier of the vertical eddy mixing coefficient. Estimation experiments produce very encouraging results and the mean estimated parameter value nicely converges to the true value exhibiting a satisfactory level of variability.Item Tropical Cyclone Data Assimilation: Experiments with a Coupled Global-Limited-Area Analysis System(2014-04-22) Holt, ChristinaThis study investigates the benefits of employing a limited-area data assimilation (DA) system to enhance lower-resolution global analyses in the Northwest Pacific tropical cyclone (TC) basin. Numerical experiments are carried out with a global analysis system at horizontal resolution T62 and a limited-area analysis system at resolutions from 200 km to 36 km. The global and limited-area DA systems, which are both based on the Local Ensemble Transform Kalman Filter algorithm, are implemented using a unique configuration, in which the global DA system provides information about the large-scale analysis and background uncertainty to the limited-area DA system. In experiments that address the global-to-limited-area resolution ratio, the limited-area analyses of the storm locations for experiments in which the ratio is 1:2 are, on average, more accurate than those from the global analyses. Increasing the resolution of the limited-area system beyond 100 km adds little direct benefit to the analysis of position or intensity, although 48 km analyses reduce boundary effects of coupling the models and may benefit analyses in which observations with larger representativeness error are assimilated. Two factors contribute to the higher accuracy of the limited-area analyses. First, the limited-area system improves the accuracy of the location estimates for strong storms, which is introduced when the background is updated by the global assimilation. Second, it improves the accuracy of the background estimate of the storm locations for moderate and weak storms. Improvements in the steering flow analysis due to increased resolution are modest and short-lived in the forecasts. Limited-area track forecasts are more accurate, on average, than global forecasts, independently of the strength of the storms up to five days. This forecast improvement is due to the more accurate analysis of the initial position of storms and the better representation of the interactions between the storms and their immediate environment. Experiments that test the treatment and quality control (QC) methods of TC observations show that significant gainful improvements can be achieved in the analyses and forecasts of TCs when observations with large representativeness error are not discarded in the online QC procedure. These experiments examine the impact of assimilating TCVitals SLP, QuikSCAT 10 m wind components, and reconnaissance dropsondes alongside the conventional observations assimilated by NCEP in real time. Implementing a Combined method that clips the special TC observations via Huberization when multiple observation types are unavailable, and keeping the TCVital observation when other special observations are present, showed significant systematic improvements for strong and moderate storm analyses and forecasts.