scholarly journals Assessment of a Nonlinear Ensemble Transform Filter for High-Dimensional Data Assimilation

2016 ◽  
Vol 144 (1) ◽  
pp. 409-427 ◽  
Author(s):  
Julian Tödter ◽  
Paul Kirchgessner ◽  
Lars Nerger ◽  
Bodo Ahrens

Abstract This work assesses the large-scale applicability of the recently proposed nonlinear ensemble transform filter (NETF) in data assimilation experiments with the NEMO ocean general circulation model. The new filter constitutes a second-order exact approximation to fully nonlinear particle filtering. Thus, it relaxes the Gaussian assumption contained in ensemble Kalman filters. The NETF applies an update step similar to the local ensemble transform Kalman filter (LETKF), which allows for efficient and simple implementation. Here, simulated observations are assimilated into a simplified ocean configuration that exhibits globally high-dimensional dynamics with a chaotic mesoscale flow. The model climatology is used to initialize an ensemble of 120 members. The number of observations in each local filter update is of the same order resulting from the use of a realistic oceanic observation scenario. Here, an importance sampling particle filter (PF) would require at least 106 members. Despite the relatively small ensemble size, the NETF remains stable and converges to the truth. In this setup, the NETF achieves at least the performance of the LETKF. However, it requires a longer spinup period because the algorithm only relies on the particle weights at the analysis time. These findings show that the NETF can successfully deal with a large-scale assimilation problem in which the local observation dimension is of the same order as the ensemble size. Thus, the second-order exact NETF does not suffer from the PF’s curse of dimensionality, even in a deterministic system.

2016 ◽  
Vol 144 (5) ◽  
pp. 2007-2020 ◽  
Author(s):  
Jonathan Poterjoy ◽  
Jeffrey L. Anderson

This study presents the first application of a localized particle filter (PF) for data assimilation in a high-dimensional geophysical model. Particle filters form Monte Carlo approximations of model probability densities conditioned on observations, while making no assumptions about the underlying error distribution. Unlike standard PFs, the local PF uses a localization function to reduce the influence of distant observations on state variables, which significantly decreases the number of particles required to maintain the filter’s stability. Because the local PF operates effectively using small numbers of particles, it provides a possible alternative to Gaussian filters, such as ensemble Kalman filters, for large geophysical models. In the current study, the local PF is compared with stochastic and deterministic ensemble Kalman filters using a simplified atmospheric general circulation model. The local PF is found to provide stable filtering results over yearlong data assimilation experiments using only 25 particles. The local PF also outperforms the Gaussian filters when observation networks include measurements that have non-Gaussian errors or relate nonlinearly to the model state, like remotely sensed data used frequently in atmospheric analyses. Results from this study encourage further testing of the local PF on more complex geophysical systems, such as weather prediction models.


2013 ◽  
Vol 20 (6) ◽  
pp. 1031-1046 ◽  
Author(s):  
S. G. Penny ◽  
E. Kalnay ◽  
J. A. Carton ◽  
B. R. Hunt ◽  
K. Ide ◽  
...  

Abstract. The most widely used methods of data assimilation in large-scale oceanography, such as the Simple Ocean Data Assimilation (SODA) algorithm, specify the background error covariances and thus are unable to refine the weights in the assimilation as the circulation changes. In contrast, the more computationally expensive Ensemble Kalman Filters (EnKF) such as the Local Ensemble Transform Kalman Filter (LETKF) use an ensemble of model forecasts to predict changes in the background error covariances and thus should produce more accurate analyses. The EnKFs are based on the approximation that ensemble members reflect a Gaussian probability distribution that is transformed linearly during the forecast and analysis cycle. In the presence of nonlinearity, EnKFs can gain from replacing each analysis increment by a sequence of smaller increments obtained by recursively applying the forecast model and data assimilation procedure over a single analysis cycle. This has led to the development of the "running in place" (RIP) algorithm by Kalnay and Yang (2010) and Yang et al. (2012a,b) in which the weights computed at the end of each analysis cycle are used recursively to refine the ensemble at the beginning of the analysis cycle. To date, no studies have been carried out with RIP in a global domain with real observations. This paper provides a comparison of the aforementioned assimilation methods in a set of experiments spanning seven years (1997–2003) using identical forecast models, initial conditions, and observation data. While the emphasis is on understanding the similarities and differences between the assimilation methods, comparisons are also made to independent ocean station temperature, salinity, and velocity time series, as well as ocean transports, providing information about the absolute error of each. Comparisons to independent observations are similar for the assimilation methods but the observation-minus-background temperature differences are distinctly lower for LETKF and RIP. The results support the potential for LETKF to improve the quality of ocean analyses on the space and timescales of interest for seasonal prediction and for RIP to accelerate the spin up of the system.


2017 ◽  
Vol 145 (11) ◽  
pp. 4575-4592 ◽  
Author(s):  
Craig H. Bishop ◽  
Jeffrey S. Whitaker ◽  
Lili Lei

To ameliorate suboptimality in ensemble data assimilation, methods have been introduced that involve expanding the ensemble size. Such expansions can incorporate model space covariance localization and/or estimates of climatological or model error covariances. Model space covariance localization in the vertical overcomes problematic aspects of ensemble-based satellite data assimilation. In the case of the ensemble transform Kalman filter (ETKF), the expanded ensemble size associated with vertical covariance localization would also enable the simultaneous update of entire vertical columns of model variables from hyperspectral and multispectral satellite sounders. However, if the original formulation of the ETKF were applied to an expanded ensemble, it would produce an analysis ensemble that was the same size as the expanded forecast ensemble. This article describes a variation on the ETKF called the gain ETKF (GETKF) that takes advantage of covariances from the expanded ensemble, while producing an analysis ensemble that has the required size of the unexpanded forecast ensemble. The approach also yields an inflation factor that depends on the localization length scale that causes the GETKF to perform differently to an ensemble square root filter (EnSRF) using the same expanded ensemble. Experimentation described herein shows that the GETKF outperforms a range of alternative ETKF-based solutions to the aforementioned problems. In cycling data assimilation experiments with a newly developed storm-track version of the Lorenz-96 model, the GETKF analysis root-mean-square error (RMSE) matches the EnSRF RMSE at shorter than optimal localization length scales but is superior in that it yields smaller RMSEs for longer localization length scales.


Ocean Science ◽  
2012 ◽  
Vol 8 (2) ◽  
pp. 143-159 ◽  
Author(s):  
S. Cailleau ◽  
J. Chanut ◽  
J.-M. Lellouche ◽  
B. Levier ◽  
C. Maraldi ◽  
...  

Abstract. The regional ocean operational system remains a key element in downscaling from large scale (global or basin scale) systems to coastal ones. It enables the transition between systems in which the resolution and the resolved physics are quite different. Indeed, coastal applications need a system to predict local high frequency events (inferior to the day) such as storm surges, while deep sea applications need a system to predict large scale lower frequency ocean features. In the framework of the ECOOP project, a regional system for the Iberia-Biscay-Ireland area has been upgraded from an existing V0 version to a V2. This paper focuses on the improvements from the V1 system, for which the physics are close to a large scale basin system, to the V2 for which the physics are more adapted to shelf and coastal issues. Strong developments such as higher regional physics resolution in the NEMO Ocean General Circulation Model for tides, non linear free surface and adapted vertical mixing schemes among others have been implemented in the V2 version. Thus, regional thermal fronts due to tidal mixing now appear in the latest version solution and are quite well positioned. Moreover, simulation of the stratification in shelf areas is also improved in the V2.


2007 ◽  
Vol 64 (11) ◽  
pp. 3766-3784 ◽  
Author(s):  
Philippe Lopez

Abstract This paper first reviews the current status, issues, and limitations of the parameterizations of atmospheric large-scale and convective moist processes that are used in numerical weather prediction and climate general circulation models. Both large-scale (resolved) and convective (subgrid scale) moist processes are dealt with. Then, the general question of the inclusion of diabatic processes in variational data assimilation systems is addressed. The focus is put on linearity and resolution issues, the specification of model and observation error statistics, the formulation of the control vector, and the problems specific to the assimilation of observations directly affected by clouds and precipitation.


2017 ◽  
Vol 24 (4) ◽  
pp. 681-694 ◽  
Author(s):  
Yuxin Zhao ◽  
Xiong Deng ◽  
Shaoqing Zhang ◽  
Zhengyu Liu ◽  
Chang Liu ◽  
...  

Abstract. Climate signals are the results of interactions of multiple timescale media such as the atmosphere and ocean in the coupled earth system. Coupled data assimilation (CDA) pursues balanced and coherent climate analysis and prediction initialization by incorporating observations from multiple media into a coupled model. In practice, an observational time window (OTW) is usually used to collect measured data for an assimilation cycle to increase observational samples that are sequentially assimilated with their original error scales. Given different timescales of characteristic variability in different media, what are the optimal OTWs for the coupled media so that climate signals can be most accurately recovered by CDA? With a simple coupled model that simulates typical scale interactions in the climate system and twin CDA experiments, we address this issue here. Results show that in each coupled medium, an optimal OTW can provide maximal observational information that best fits the characteristic variability of the medium during the data blending process. Maintaining correct scale interactions, the resulting CDA improves the analysis of climate signals greatly. These simple model results provide a guideline for when the real observations are assimilated into a coupled general circulation model for improving climate analysis and prediction initialization by accurately recovering important characteristic variability such as sub-diurnal in the atmosphere and diurnal in the ocean.


2007 ◽  
Vol 4 (5) ◽  
pp. 3413-3440 ◽  
Author(s):  
E. P. Maurer ◽  
H. G. Hidalgo

Abstract. Downscaling of climate model data is essential to most impact analysis. We compare two methods of statistical downscaling to produce continuous, gridded time series of precipitation and surface air temperature at a 1/8-degree (approximately 140 km² per grid cell) resolution over the western U.S. We use NCEP/NCAR Reanalysis data from 1950–1999 as a surrogate General Circulation Model (GCM). The two methods included are constructed analogues (CA) and a bias correction and spatial downscaling (BCSD), both of which have been shown to be skillful in different settings, and BCSD has been used extensively in hydrologic impact analysis. Both methods use the coarse scale Reanalysis fields of precipitation and temperature as predictors of the corresponding fine scale fields. CA downscales daily large-scale data directly and BCSD downscales monthly data, with a random resampling technique to generate daily values. The methods produce comparable skill in producing downscaled, gridded fields of precipitation and temperatures at a monthly and seasonal level. For daily precipitation, both methods exhibit some skill in reproducing both observed wet and dry extremes and the difference between the methods is not significant, reflecting the general low skill in daily precipitation variability in the reanalysis data. For low temperature extremes, the CA method produces greater downscaling skill than BCSD for fall and winter seasons. For high temperature extremes, CA demonstrates higher skill than BCSD in summer. We find that the choice of most appropriate downscaling technique depends on the variables, seasons, and regions of interest, on the availability of daily data, and whether the day to day correspondence of weather from the GCM needs to be reproduced for some applications. The ability to produce skillful downscaled daily data depends primarily on the ability of the climate model to show daily skill.


2006 ◽  
Vol 24 (8) ◽  
pp. 2075-2089 ◽  
Author(s):  
A. Chakraborty ◽  
R. S. Nanjundiah ◽  
J. Srinivasan

Abstract. A theory is proposed to determine the onset of the Indian Summer Monsoon (ISM) in an Atmospheric General Circulation Model (AGCM). The onset of ISM is delayed substantially in the absence of global orography. The impact of orography over different parts of the Earth on the onset of ISM has also been investigated using five additional perturbed simulations. The large difference in the date of onset of ISM in these simulations has been explained by a new theory based on the Surface Moist Static Energy (SMSE) and vertical velocity at the mid-troposphere. It is found that onset occurs only after SMSE crosses a threshold value and the large-scale vertical motion in the middle troposphere becomes upward. This study shows that both dynamics and thermodynamics play profound roles in the onset of the monsoon.


2021 ◽  
Author(s):  
Leonardo Lima ◽  
Stefania Angela Ciliberti ◽  
Ali Aydogdu ◽  
Romain Escudier ◽  
Simona Masina ◽  
...  

<p>Ocean reanalyses are becoming increasingly important to reconstruct and provide an overview of the ocean state from the past to the present-day. These products require advanced scientific methods and techniques to produce a more accurate ocean representation. In the scope of the Copernicus Marine Environment Monitoring Service (CMEMS), a new Black Sea (BS) reanalysis, BS-REA (BSE3R1 system), has been produced by using an advanced variational data assimilation method to combine the best available observations with a state-of-the-art ocean general circulation model. The hydrodynamical model is based on Nucleus for European Modeling of the Ocean (NEMO, v3.6), implemented for the BS domain with horizontal resolution of 1/27° x 1/36°, and 31 unevenly distributed vertical levels. NEMO is forced by atmospheric surface fluxes computed via bulk formulation and forced by ECMWF ERA5 atmospheric reanalysis product. At the surface, the model temperature is relaxed to daily objective analysis fields of sea surface temperature from CMEMS SST TAC. The exchange with Mediterranean Sea is simulated through relaxation of the temperature and salinity near Bosporus toward a monthly climatology computed from a high-resolution multi-year simulation, and the barotropic Bosporus Strait transport is corrected to balance the variations of the freshwater flux and the sea surface height measured by multi-satellite altimetry observations. A 3D-Var ocean data assimilation scheme (OceanVar) is used to assimilate sea level anomaly along-track observations from CMEMS SL TAC and available in situ vertical profiles of temperature and salinity from both SeaDataNet and CMEMS INS TAC products. Comparisons against the previous Black Sea reanalysis (BSE2R2 system) show important improvements for temperature and salinity, such that errors have significantly decreased (about 50%). Temperature fields present a continuous warming in the layer between 25-150 m, within which there is the presence of the Black Sea Cold Intermediate Layer (CIL). SST exhibits a positive bias and relatively higher root mean square error (RMSE) values are present in the summer season. Spatial maps of sea level anomaly reveal the largest RMSE close to the shelf areas, which are related to the mesoscale activity along the Rim current. The BS-REA catalogue includes daily and monthly means for 3D temperature, salinity, and currents and 2D sea surface height, bottom temperature, mixed layer fields, from Jan 1993 to Dec 2019.  The BSE3R1 system has produced very accurate estimates which makes it very suitable for assessing more realistic climate trends and indicators for important ocean properties.</p>


Sign in / Sign up

Export Citation Format

Share Document