Dynamical Structures of Cross-Domain Forecast Error Covariance of a Simulated Tropical Cyclone in a Convection-Permitting Coupled Atmosphere–Ocean Model

2021 ◽  
Vol 149 (1) ◽  
pp. 41-63
Author(s):  
Xingchao Chen ◽  
Robert G. Nystrom ◽  
Christopher A. Davis ◽  
Colin M. Zarzycki
Keyword(s):  
Tropical Cyclone ◽  
Data Assimilation ◽  
Forecast Error ◽  
State Variables ◽  
Ensemble Spread ◽  
Strongly Coupled ◽  
Error Covariance ◽  
Cross Domain ◽  
Forecast Lead Time ◽  
Coupled Data Assimilation

AbstractUnderstanding the dynamics of the flow-dependent forecast error covariance across the air–sea interface is beneficial toward revealing the potential influences of strongly coupled data assimilation on tropical cyclone (TC) initialization in coupled models, and the fundamental dynamics associated with TC air–sea interactions. A 200-member ensemble of convection-permitting forecasts from a coupled atmosphere–ocean regional model is used to investigate the forecast error covariance across the oceanic and atmospheric domains during the rapid intensification of Hurricane Florence (2018). Forecast uncertainties in both atmospheric and oceanic domains, from an Eulerian perspective, increase with forecast lead time, mainly from TC displacement errors. In a storm-relative framework, the ensemble forecast uncertainties in both domains are predominantly caused by differences in the simulated storm intensity and structure. The largest ensemble spread in the atmospheric pressure, temperature, and wind fields can be found within the TC inner-core region. Alternatively, the largest ensemble spread in the upper-ocean currents and temperature fields are located along the cold wake behind the storm. Cross-domain ensemble correlations between simulated atmospheric (oceanic) observations and oceanic (atmospheric) state variables in the storm-relative coordinates are highly anisotropic, variable dependent, and ultimately driven by the dynamics of TC air–sea interactions. Meaningful and dynamically consistent cross-domain ensemble correlations suggest that it is possible to use atmospheric and oceanic observations to simultaneously update state variables associated with the coupled ocean–atmosphere prediction of TCs using strongly coupled data assimilation. Sensitivity experiments demonstrate that at least 60–80 ensemble members are required to represent physically consistent cross-domain correlations and minimize sampling errors.

Air-sea strongly coupled data assimilation for tropical cyclone prediction

2021 ◽  
Author(s):  
Xingchao Chen
Keyword(s):  
Tropical Cyclone ◽  
Data Assimilation ◽  
Forecast Error ◽  
Regional Scale ◽  
Forecast Errors ◽  
State Variables ◽  
Strongly Coupled ◽  
Weakly Coupled ◽  
Coupled Data Assimilation ◽  
Tc Prediction

<p>Air-sea interactions are critical to tropical cyclone (TC) energetics. However, oceanic state variables are still poorly initialized, and are inconsistent with atmospheric initial fields in most operational coupled TC forecast models. In this study, we first investigate the forecast error covariance across the oceanic and atmospheric domains during the rapid intensification of Hurricane Florence (2018) using a 200-member ensemble of convection-permitting forecasts from a coupled atmosphere-ocean regional model. Meaningful and dynamically consistent cross domain ensemble error correlations suggest that it is possible to use atmospheric and oceanic observations to simultaneously update model state variables associated with the coupled ocean-atmosphere prediction of TCs using strongly coupled data assimilation (DA). A regional-scale strongly coupled DA system based on the ensemble Kalman filter (EnKF) is then developed for TC prediction. The potential impacts of different atmospheric and oceanic observations on TC analysis and prediction are examined through observing system simulation experiments (OSSEs) of Hurricane Florence (2018). Results show that strongly coupled DA resulted in better analysis and forecast of both the oceanic and atmospheric variables than weakly coupled DA. Compared to weakly coupled DA in which the analysis update is performed separately for the atmospheric and oceanic domains, strongly coupled DA reduces the forecast errors of TC track and intensity. Results show promise in potential further improvement in TC prediction through assimilation of both atmospheric and oceanic observations using the ensemble-based strongly coupled DA system.</p>

scholarly journals Application of a local attractor dimension to reduced space strongly coupled data assimilation for chaotic multiscale systems

2020 ◽  
Vol 27 (1) ◽  
pp. 51-74 ◽  
Author(s):  
Courtney Quinn ◽  
Terence J. O'Kane ◽  
Vassili Kitsios
Keyword(s):  
Data Assimilation ◽  
Accurate Estimate ◽  
Variable Number ◽  
Frobenius Norm ◽  
Accurate Analysis ◽  
Attractor Dimension ◽  
Strongly Coupled ◽  
Cross Domain ◽  
Local Attractor ◽  
Coupled Data Assimilation

Abstract. The basis and challenge of strongly coupled data assimilation (CDA) is the accurate representation of cross-domain covariances between various coupled subsystems with disparate spatio-temporal scales, where often one or more subsystems are unobserved. In this study, we explore strong CDA using ensemble Kalman filtering methods applied to a conceptual multiscale chaotic model consisting of three coupled Lorenz attractors. We introduce the use of the local attractor dimension (i.e. the Kaplan–Yorke dimension, dimKY) to prescribe the rank of the background covariance matrix which we construct using a variable number of weighted covariant Lyapunov vectors (CLVs). Specifically, we consider the ability to track the nonlinear trajectory of each of the subsystems with different variants of sparse observations, relying only on the cross-domain covariance to determine an accurate analysis for tracking the trajectory of the unobserved subdomain. We find that spanning the global unstable and neutral subspaces is not sufficient at times where the nonlinear dynamics and intermittent linear error growth along a stable direction combine. At such times a subset of the local stable subspace is also needed to be represented in the ensemble. In this regard the local dimKY provides an accurate estimate of the required rank. Additionally, we show that spanning the full space does not improve performance significantly relative to spanning only the subspace determined by the local dimension. Where weak coupling between subsystems leads to covariance collapse in one or more of the unobserved subsystems, we apply a novel modified Kalman gain where the background covariances are scaled by their Frobenius norm. This modified gain increases the magnitude of the innovations and the effective dimension of the unobserved domains relative to the strength of the coupling and timescale separation. We conclude with a discussion on the implications for higher-dimensional systems.

scholarly journals Multivariate localization functions for strongly coupled data assimilation in the bivariate Lorenz ’96 system

2021 ◽  
Author(s):  
Zofia Stanley ◽  
Ian Grooms ◽  
William Kleiber
Keyword(s):  
Data Assimilation ◽  
Spatial Scales ◽  
Background Error ◽  
Strongly Coupled ◽  
Localization Function ◽  
Background Error Covariance ◽  
Error Covariance ◽  
Coupled Data Assimilation ◽  
Lorenz 96 ◽  
The Impact

Abstract. Localization is widely used in data assimilation schemes to mitigate the impact of sampling errors on ensemble-derived background error covariance matrices. Strongly coupled data assimilation allows observations in one component of a coupled model to directly impact another component through inclusion of cross-domain terms in the background error covariance matrix. When different components have disparate dominant spatial scales, localization between model domains must properly account for the multiple length scales at play. In this work we develop two new multivariate localization functions, one of which is a multivariate extension of the fifth-order piecewise rational Gaspari-Cohn localization function; the within-component localization functions are standard Gaspari-Cohn with different localization radii while the cross-localization function is newly constructed. The functions produce non-negative definite localization matrices, which are suitable for use in variational data assimilation schemes. We compare the performance of our two new multivariate localization functions to two other multivariate localization functions and to the univariate analogs of all four functions in a simple experiment with the bivariate Lorenz '96 system. In our experiment the multivariate Gaspari-Cohn function leads to better performance than any of the other localization functions.

scholarly journals Correlation-Cutoff Method for Covariance Localization in Strongly Coupled Data Assimilation

2018 ◽  
Vol 146 (9) ◽  
pp. 2881-2889 ◽  
Author(s):  
Takuma Yoshida ◽  
Eugenia Kalnay
Keyword(s):  
Kalman Filter ◽  
Data Assimilation ◽  
Coupled Model ◽  
Error Variance ◽  
State Variables ◽  
Background Error ◽  
Strongly Coupled ◽  
Error Correlation ◽  
Analysis Error ◽  
Coupled Data Assimilation

Abstract Strongly coupled data assimilation (SCDA), where observations of one component of a coupled model are allowed to directly impact the analysis of other components, sometimes fails to improve the analysis accuracy with an ensemble Kalman filter (EnKF) as compared with weakly coupled data assimilation (WCDA). It is well known that an observation’s area of influence should be localized in EnKFs since the assimilation of distant observations often degrades the analysis because of spurious correlations. This study derives a method to estimate the reduction of the analysis error variance by using estimates of the cross covariances between the background errors of the state variables in an idealized situation. It is shown that the reduction of analysis error variance is proportional to the squared background error correlation between the analyzed and observed variables. From this, the authors propose an offline method to systematically select which observations should be assimilated into which model state variable by cutting off the assimilation of observations when the squared background error correlation between the observed and analyzed variables is small. The proposed method is tested with the local ensemble transform Kalman filter (LETKF) and a nine-variable coupled model, in which three Lorenz models with different time scales are coupled with each other. The covariance localization with the correlation-cutoff method achieves an analysis more accurate than either the full SCDA or the WCDA methods, especially with smaller ensemble sizes.

scholarly journals Multivariate localization functions for strongly coupled data assimilation in the bivariate Lorenz 96 system

2021 ◽  
Vol 28 (4) ◽  
pp. 565-583
Author(s):  
Zofia Stanley ◽  
Ian Grooms ◽  
William Kleiber
Keyword(s):  
Data Assimilation ◽  
Spatial Scales ◽  
Background Error ◽  
Strongly Coupled ◽  
Localization Function ◽  
Background Error Covariance ◽  
Error Covariance ◽  
Coupled Data Assimilation ◽  
Lorenz 96 ◽  
The Impact

Abstract. Localization is widely used in data assimilation schemes to mitigate the impact of sampling errors on ensemble-derived background error covariance matrices. Strongly coupled data assimilation allows observations in one component of a coupled model to directly impact another component through the inclusion of cross-domain terms in the background error covariance matrix. When different components have disparate dominant spatial scales, localization between model domains must properly account for the multiple length scales at play. In this work, we develop two new multivariate localization functions, one of which is a multivariate extension of the fifth-order piecewise rational Gaspari–Cohn localization function; the within-component localization functions are standard Gaspari–Cohn with different localization radii, while the cross-localization function is newly constructed. The functions produce positive semidefinite localization matrices which are suitable for use in both Kalman filters and variational data assimilation schemes. We compare the performance of our two new multivariate localization functions to two other multivariate localization functions and to the univariate and weakly coupled analogs of all four functions in a simple experiment with the bivariate Lorenz 96 system. In our experiments, the multivariate Gaspari–Cohn function leads to better performance than any of the other multivariate localization functions.

scholarly journals Application of local attractor dimension to reduced space strongly coupled data assimilation for chaotic multiscale systems

2019 ◽  
Author(s):  
Courtney Quinn ◽  
Terence J. O'Kane ◽  
Vassili Kitsios
Keyword(s):  
Data Assimilation ◽  
Accurate Estimate ◽  
Variable Number ◽  
Frobenius Norm ◽  
Accurate Analysis ◽  
Attractor Dimension ◽  
Strongly Coupled ◽  
Cross Domain ◽  
Local Attractor ◽  
Coupled Data Assimilation

Abstract. The basis and challenge of strongly coupled data assimilation (CDA) is the accurate representation of cross-domain covariances between various coupled subsystems with disparate spatio-temporal scales, where often one or more subsystems are unobserved. In this study, we explore strong CDA using ensemble Kalman filtering methods applied to a conceptual multiscale chaotic model consisting of three coupled Lorenz attractors. We introduce the use of the local attractor dimension (i.e. the Kaplan-Yorke dimension, dimKY) to determine the rank of the background covariance matrix which we construct using a variable number of weighted covariant Lyapunov vectors (CLVs). Specifically, we consider the ability to track the nonlinear trajectory of each of the subsystems with different variants of sparse observations, relying only on the cross-domain covariance to determine an accurate analysis for tracking the trajectory of the unobserved subdomain. We find that spanning the global unstable and neutral subspaces is not sufficient at times where the nonlinear dynamics and intermittent linear error growth along a stable direction combine. At such times a subset of the local stable subspace is also needed to be represented in the ensemble. In this regard the local dimKY provides an accurate estimate of the required rank. Additionally, we show that spanning the full space does not improve performance significantly relative to spanning only the subspace determined by the local dimension. Where weak coupling between subsystems leads to covariance collapse in one or more of the unobserved subsystems, we apply a novel modified Kalman gain where the background covariances are scaled by their Frobenius norm. This modified gain increases the magnitude of the innovations and the effective dimension of the unobserved domains relative to the strength of the coupling and time-scale separation. We conclude with a discussion on the implications for higher dimensional systems.

Theoretical Application of Ensemble Kalman Filter to Adsorptive Solute Cr(VI) Transfer from Soil into Surface Runoff

2014 ◽  
Vol 919-921 ◽  
pp. 1257-1261
Author(s):  
Chao Qun Tan ◽  
Ju Xiu Tong ◽  
Bill X. Hu ◽  
Jin Zhong Yang
Keyword(s):  
Kalman Filter ◽  
Data Assimilation ◽  
Ensemble Kalman Filter ◽  
Surface Runoff ◽  
Amplification Factor ◽  
Forecast Error ◽  
Assimilation Effect ◽  
Solute Transfer ◽  
Error Covariance ◽  
Theoretical Application

This paper mainly discusses some details when applying data assimilation method via an ensemble Kalman filter (EnKF) to improve prediction of adsorptive solute Cr(VI) transfer from soil into runoff. Based on this work, we could make better use of our theoretical model to predict adsorptive solute transfer from soil into surface runoff in practice. The results show that the ensemble number of 100 is reasonable, considering assimilation effect and efficiency after selecting its number from 25 to 225 at an interval of 25. While the initial ensemble value makes little difference to data assimilation (DA) results. Besides, DA results could be improved by multiplying an amplification factor to forecast error covariance matrix due to underestimation of forecast error.

scholarly journals Facilitating Strongly Coupled Ocean–Atmosphere Data Assimilation with an Interface Solver

2015 ◽  
Vol 144 (1) ◽  
pp. 3-20 ◽  
Author(s):  
Sergey Frolov ◽  
Craig H. Bishop ◽  
Teddy Holt ◽  
James Cummings ◽  
David Kuhl
Keyword(s):  
Data Assimilation ◽  
Impact Analysis ◽  
Coupled Model ◽  
Deep Ocean ◽  
Accurate Solution ◽  
Fluid Interface ◽  
Strongly Coupled ◽  
Grid Points ◽  
Coupled Data Assimilation ◽  
Observation Vector

Abstract In a strongly coupled data assimilation (DA), a cross-fluid covariance is specified that allows measurements from a coupled fluid (e.g., atmosphere) to directly impact analysis increments in a target fluid (e.g., ocean). The exhaustive solution to this coupled DA problem calls for a covariance where all available measurements can influence all grid points in all fluids. Solution of such a large algebraic problem is computationally expensive, often calls for a substantial rewrite of existing fluid-specific DA systems, and, as shown in this paper, can be avoided. The proposed interface solver assumes that covariances between coupled measurements and target fluid are often close to null (e.g., between stratospheric observations and the deep ocean within a 6-h forecast cycle). In the interface solver, two separate DA solvers are run in parallel: one that produces an analysis solution in the atmosphere, and one in the ocean. Each system uses a coupled observation vector where in addition to resident measurements in the target fluid it also includes nonresident measurements in the coupled fluid that are likely to have significant influence on the analysis in the target fluid (interface measurements). An ensemble-based method is employed and a localization function for coupled ensembles is proposed. Using a coupled model for the Mediterranean Sea (in a twin setting), it is demonstrated that (i) the solution of the interface solver converges to the exhaustive solution and (ii) that in presence of poorly known error covariances, the interface solver can be configured to produce a more accurate solution than an exhaustive solver.

scholarly journals The Development of a Hybrid EnKF-3DVAR Algorithm for Storm-Scale Data Assimilation

2013 ◽  
Vol 2013 ◽  
pp. 1-12 ◽  
Author(s):  
Jidong Gao ◽  
Ming Xue ◽  
David J. Stensrud
Keyword(s):  
Data Assimilation ◽  
Forecast Error ◽  
Control Variable ◽  
Radar Data ◽  
Ensemble Forecasts ◽  
Advanced Regional Prediction System ◽  
Error Covariance ◽  
Variable Approach ◽  
Regional Prediction ◽  
Rms Errors

A hybrid 3DVAR-EnKF data assimilation algorithm is developed based on 3DVAR and ensemble Kalman filter (EnKF) programs within the Advanced Regional Prediction System (ARPS). The hybrid algorithm uses the extended alpha control variable approach to combine the static and ensemble-derived flow-dependent forecast error covariances. The hybrid variational analysis is performed using an equal weighting of static and flow-dependent error covariance as derived from ensemble forecasts. The method is first applied to the assimilation of simulated radar data for a supercell storm. Results obtained using 3DVAR (with static covariance entirely), hybrid 3DVAR-EnKF, and the EnKF are compared. When data from a single radar are used, the EnKF method provides the best results for the model dynamic variables, while the hybrid method provides the best results for hydrometeor related variables in term of rms errors. Although storm structures can be established reasonably well using 3DVAR, the rms errors are generally worse than seen from the other two methods. With two radars, the results from 3DVAR are closer to those from EnKF. Our tests indicate that the hybrid scheme can reduce the storm spin-up time because it fits the observations, especially the reflectivity observations, better than the EnKF and the 3DVAR at the beginning of the assimilation cycles.

scholarly journals Localized Ensemble-Based Tangent Linear Models and Their Use in Propagating Hybrid Error Covariance Models

2016 ◽  
Vol 144 (4) ◽  
pp. 1383-1405 ◽  
Author(s):  
Sergey Frolov ◽  
Craig H. Bishop
Keyword(s):  
Data Assimilation ◽  
Short Duration ◽  
Linear Models ◽  
Forecast Error ◽  
Weather Forecast ◽  
Step Function ◽  
Rank Deficiency ◽  
One Dimensional ◽  
Error Covariance ◽  
Covariance Models

Abstract Hybrid error covariance models that blend climatological estimates of forecast error covariances with ensemble-based, flow-dependent forecast error covariances have led to significant reductions in forecast error when employed in 4DVAR data assimilation schemes. Tangent linear models (TLMs) designed to predict the differences between perturbed and unperturbed simulations of the weather forecast are a key component of such 4DVAR schemes. However, many forecasting centers have found that TLMs and their adjoints do not scale well computationally and are difficult to create and maintain—particularly for coupled ocean–wave–ice–atmosphere models. In this paper, the authors create ensemble-based TLMs (ETLMs) and test their ability to propagate both climatological and flow-dependent parts of hybrid error covariance models. These tests demonstrate that rank deficiency limits the utility of unlocalized ETLMs. High-rank, time-evolving, flow-adaptive localization functions are constructed and tested using recursive application of short-duration ETLMs, each of which is localized using a static localization. Since TLM operators do not need to be semipositive definite, the authors experiment with a variety of localization approaches including step function localization. The step function localization leads to a local formulation that was found to be highly effective. In tests using simple one-dimensional models with both dispersive and nondispersive dynamics, it is shown that practical ETLM configurations were effective at propagating covariances as far as four error correlation scales.

Sign in / Sign up

Export Citation Format

Share Document