Super-resolution data assimilation

Ocean Dynamics - Increasing model resolution can improve the performance of a data assimilation system because it reduces model error, the system can more optimally use high-resolution...     

Data assimilation for precipitation nowcasting using Bayesian inference

This work introduces a new variational Bayes data assimilation method for the stochastic estimation of precipitation dynamics using radar observations for short term probabilistic forecasting (nowcasting). A previously developed spatial rainfall model based on the decomposition of the observed precipitation field using a basis function expansion captures the precipitation intensity from radar images as a set of ‘rain cells’. The prior distributions for the basis function parameters are carefully chosen to have a conjugate structure for the precipitation field model to allow a novel variational Bayes method to be applied to estimate the posterior distributions in closed form, based on solving an optimisation problem, in a spirit similar to 3D VAR analysis, but seeking approximations to the posterior distribution rather than simply the most probable state. A hierarchical Kalman filter is used to estimate the advection field based on the assimilated precipitation fields at two times. The model is applied to tracking precipitation dynamics in a realistic setting, using UK Met Office radar data from both a summer convective event and a winter frontal event. The performance of the model is assessed both traditionally and using probabilistic measures of fit based on ROC curves. The model is shown to provide very good assimilation characteristics, and promising forecast skill. Improvements to the forecasting scheme are discussed.     

Use of Scale Recursive Estimation for assimilation of precipitation data from TRMM (PR and TMI) and NEXRAD

The paper shows an application of Scale Recursive Estimation (SRE) used to assimilate rainfall rates estimated during a storm event from three remote sensing devices. These are the TMI radiometer and the PR radar, carried on board of the TRMM satellite and the KNQA Memphis Weather Surveillance radar, belonging to the NEXRAD network, each one providing rain rate estimates at a different spatial scale. The variability of rain rate process in scales is modeled as a multiplicative random cascade, including spatial intermittence. The observational noise in the estimates is modeled according to a multiplicative error. System estimation, including process and observational noise, is carried out using Maximum Likelihood Estimation implemented by a scale recursive Expectation Maximization (EM) algorithm. As a result, new rainfall rate estimates are obtained that feature decreased estimation error as compared to those coming from each device alone. The performance of the SRE-EM approach is compared with that of the latest methods proposed for data fusion of multisensor estimates. The proposed approach improves the current methods adopted for SRE and provides an alternative for data fusion in the field of precipitation.     

Towards a data assimilation system for morphodynamic modeling: bathymetric data assimilation for wave property estimation

Data assimilation is mainly concerned with the proper management of uncertainties. The main objective of the present work is to implement and analyze a data assimilation technique capable of assimilating bathymetric data into a coupled flow, wave, and morphodynamic model. For the case presented here, wave significant height, wave direction of incidence, and wave peak period are being optimized based on bathymetric data taken from a twin experiment. An adjoint-free variational scheme is used. In this approach, a linear reduced order model (ROM) is constructed as an approximation of the full model. The ROM is an autoregressive model of order 1 (AR1) that preserves the parametrization. Since the ROM is linear, the construction of its adjoint is straightforward, making the implementation of 4D variational data assimilation effortless. The scheme is able to update the morphodynamic model satisfactorily despite the fact that the model shows nonlinear behavior even for very small perturbations of all three parameters. The size and direction of the perturbations necessary for constructing the ROM have a significant impact on the performance of the technique.     

Neural networks in auroral data assimilation

Data assimilation is an essential step for improving space weather forecasting by means of a weighted combination between observational data and data from a mathematical model. In the present work data assimilation methods based on Kalman filter (KF) and artificial neural networks are applied to a three-wave model of auroral radio emissions. A novel data assimilation method is presented, whereby a multilayer perceptron neural network is trained to emulate a KF for data assimilation by using cross-validation. The results obtained render support for the use of neural networks as an assimilation technique for space weather prediction.     

Turbulent viscosity optimized by data assimilation

As an alternative approach to classical turbulence modelling using a first or second order closure, the data assimilation method of optimal control is applied to estimate a time and space-dependent turbulent viscosity in a three-dimensional oceanic circulation model. The optimal control method, described for a 3-D primitive equation model, involves the minimization of a cost function that quantifies the discrepancies between the simulations and the observations. An iterative algorithm is obtained via the adjoint model resolution. In a first experiment, a k ± L model is used to simulate the one-dimensional development of inertial oscillations resulting from a wind stress at the sea surface and with the presence of a halocline. These results are used as synthetic observations to be assimilated. The turbulent viscosity is then recovered without the k + L closure, even with sparse and noisy observations. The problems of controllability and of the dimensions of the control are then discussed. A second experiment consists of a two-dimensional schematic simulation. A 2-D turbulent viscosity field is estimated from data on the initial and final states of a coastal upwelling event.     

OpenDA-NEMO framework for ocean data assimilation

Data assimilation methods provide a means to handle the modeling errors and uncertainties in sophisticated ocean models. In this study, we have created an OpenDA-NEMO framework unlocking the data assimilation tools available in OpenDA for use with NEMO models. This includes data assimilation methods, automatic parallelization, and a recently implemented automatic localization algorithm that removes spurious correlations in the model based on uncertainties in the computed Kalman gain matrix. We have set up a twin experiment where we assimilate sea surface height (SSH) satellite measurements. From the experiments, we can conclude that the OpenDA-NEMO framework performs as expected and that the automatic localization significantly improves the performance of the data assimilation algorithm by successfully removing spurious correlations. Based on these results, it looks promising to extend the framework with new kinds of observations and work on improving the computational speed of the automatic localization technique such that it becomes feasible to include large number of observations.     

Review of the Kalman-type hydrological data assimilation

ABSTRACT

There is great potential in Data Assimilation (DA) for the purposes of uncertainty identification, reduction and real-time correction of hydrological models. This paper reviews the latest developments in Kalman filters (KFs), particularly the Extended KF (EKF) and the Ensemble KF (EnKF) in hydrological DA. The hydrological DA targets, methodologies and their applicability are examined. The recent applications of the EKF and EnKF in hydrological DA are summarized and assessed critically. Furthermore, this review highlights the existing challenges in the implementation of the EKF and EnKF, especially error determination and joint parameter estimation. A detailed review of these issues would benefit not only the Kalman-type DA but also provide an important reference to other hydrological DA types.

Editor D. Koutsoyiannis; Associate editor F. Pappenberger     

On the assimilation of total-ozone satellite data

A two-dimensional model for advection and data assimilation of total-ozone data has been developed. The Assimilation Model KNMI (AMK) is a global model describing the transport of the column amounts of ozone, by a wind field at a single pressure level, assuming that total ozone behaves as a passive tracer. In this study, ozone column amounts measured by the TIROS Operational Vertical Sounder (TOVS) instrument on the National Oceanic and Atmospheric Administration (NOAA) polar satellites and wind fields from the Meteorological Archive and Retrieval System (MARS) archives at ECMWF have been used. By means of the AMK, the incomplete space-time distribution of the TOVS measurements is filled in and global total-ozone maps at any given time can be obtained. The choice of wind field to be used for transporting column amounts of ozone is extensively discussed. It is shown that the 200-hPa wind field is the optimal single-pressure-level wind field for advecting total ozone. Assimilated ozone fields are the basic information for research on atmospheric chemistry and dynamics, but are also important for the validation of ozone measurements.     

Toward a global ocean data assimilation system based on ensemble optimum interpolation: altimetry data assimilation experiment

A global ocean data assimilation system based on the ensemble optimum interpolation (EnOI) has been under development as the Chinese contribution to the Global Ocean Data Assimilation Experiment. The system uses a global ocean general circulation model, which is eddy permitting, developed by the Institute of Atmospheric Physics of the Chinese Academy of Sciences. In this paper, the implementation of the system is described in detail. We describe the sampling strategy to generate the stationary ensembles for EnOI. In addition, technical methods are introduced to deal with the requirement of massive memory space to hold the stationary ensembles of the global ocean. The system can assimilate observations such as satellite altimetry, sea surface temperature (SST), in situ temperature and salinity from Argo, XBT, Tropical Atmosphere Ocean (TAO), and other sources in a straightforward way. As a first step, an assimilation experiment from 1997 to 2001 is carried out by assimilating the sea level anomaly (SLA) data from TOPEX/Poseidon. We evaluate the performance of the system by comparing the results with various types of observations. We find that SLA assimilation shows very positive impact on the modeled fields. The SST and sea surface height fields are clearly improved in terms of both the standard deviation and the root mean square difference. In addition, the assimilation produces some improvements in regions where mesoscale processes cannot be resolved with the horizontal resolution of this model. Comparisons with TAO profiles in the Pacific show that the temperature and salinity fields have been improved to varying degrees in the upper ocean. The biases with respect to the independent TAO profiles are reduced with a maximum magnitude of about 0.25°C and 0.1 psu for the time-averaged temperature and salinity. The improvements on temperature and salinity also lead to positive impact on the subsurface currents. The equatorial under current is enhanced in the Pacific although it is still underestimated after the assimilation.     

Error covariance modeling in sequential data assimilation

 The efficiency of a sequential data assimilation scheme relies on the capability to describe the error covariance. This aspect is all the more relevant if one needs accurate statistics on the estimation error. Frequently an ad hoc function depending on a few parameters is proposed, and these parameters are tuned, estimated or updated. This usually requires that the covariance is second-order stationary (i.e. depends only on the distance between two points). In this paper, we discuss this feature and show that even in simple applications (such as one-dimensional hydrodynamics), this assumption does not hold and may lead to poorly described estimation errors. We propose a method relying on the analysis of the error term and the use of the hydrodynamical model to generate one part of the covariance function, the other part being modeled using a second-order stationary approach. This method is discussed using a twin experiment in the case where a physical parameter is erroneous, and improves significantly the results: the model bias is strongly reduced and the estimation error is well described. Moreover, it enables a better adaptation of the Kalman gain to the actual estimation error.     

An explicit four-dimensional variational data assimilation method

A new data assimilation method called the explicit four-dimensional variational (4DVAR) method is proposed. In this method, the singular value decomposition (SVD) is used to construct the orthogonal basis vectors from a forecast ensemble in a 4D space. The basis vectors represent not only the spatial structure of the analysis variables but also the temporal evolution. After the analysis variables are ex-pressed by a truncated expansion of the basis vectors in the 4D space, the control variables in the cost function appear explicitly, so that the adjoint model, which is used to derive the gradient of cost func-tion with respect to the control variables, is no longer needed. The new technique significantly simpli-fies the data assimilation process. The advantage of the proposed method is demonstrated by several experiments using a shallow water numerical model and the results are compared with those of the conventional 4DVAR. It is shown that when the observation points are very dense, the conventional 4DVAR is better than the proposed method. However, when the observation points are sparse, the proposed method performs better. The sensitivity of the proposed method with respect to errors in the observations and the numerical model is lower than that of the conventional method.     

Nonlinear balance constraints in 3DVAR data assimilation

Consideration of the so-called balance properties/ constraints is of great importance in the data assimila- tion. First, it is desirable to separate the slow modes (e.g., Rossby waves) from the fast modes (e.g., gravity waves and deep convection). Second, consideration of the balance constraints enables us to infer information about all variables that are balanced with the observed variable and hereby improve the quality of the analy- sis. The balance properties are usually represented by some…     

Application of altimetry data assimilation on mesoscale eddies simulation

Mesoscale eddy plays an important role in the ocean circulation. In order to improve the simulation accuracy of the mesoscale eddies, a three-dimensional variation (3DVAR) data assimilation system called Ocean Variational Analysis System (OVALS) is coupled with a POM model to simulate the mesoscale eddies in the Northwest Pacific Ocean. In this system, the sea surface height anomaly (SSHA) data by satellite altimeters are assimilated and translated into pseudo temperature and salinity (T-S) profile data. Then, these profile data are taken as observation data to be assimilated again and produce the three-dimensional analysis T-S field. According to the characteristics of mesoscale eddy, the most appropriate assimilation parameters are set up and testified in this system. A ten years mesoscale eddies simulation and comparison experiment is made, which includes two schemes: assimilation and non-assimilation. The results of comparison between two schemes and the observation show that the simulation accuracy of the assimilation scheme is much better than that of non-assimilation, which verified that the altimetry data assimilation method can improve the simulation accuracy of the mesoscale dramatically and indicates that it is possible to use this system on the forecast of mesoscale eddies in the future.     

Effects of sample density on the assimilation performance of an explicit four-dimensional variational data assimilation method

The concepts of sample sphere radius and sample density are proposed in this paper to help illustrate that different vector transformations result in diverse sample density with the same sample ensemble, which finally affects their assimilation performance. Several numerical experiments using a one-dimensional (1-D) soil water equation and synthetic observations are conducted to evaluate this new theory in land data assimilation.     

Temperature assimilation into a coastal ocean-biogeochemical model: assessment of weakly and strongly coupled data assimilation

Ocean Dynamics - Satellite data of both physical properties as well as ocean colour can be assimilated into coupled ocean-biogeochemical models with the aim to improve the model state. The physical...     

Trend analysis in Turkish precipitation data

This study aims to determine trends in the long‐term annual mean and monthly total precipitation series using non‐parametric methods (i.e. the Mann–Kendall and Sen's T tests). The change per unit time in a time series having a linear trend was estimated by applying a simple non‐parametric procedure, namely Sen's estimator of slope. Serial correlation structure in the data was accounted for determining the significance level of the results of the Mann–Kendall test. The data network used in this study, which is assumed to reflect regional hydroclimatic conditions, consists of 96 precipitation stations across Turkey. Monthly totals and annual means of the monthly totals are formed for each individual station, spanning from 1929 to 1993. In this case, a total of 13 precipitation variables at each station are subjected to trend detection analysis. In addition, regional average precipitation series are established for the same analysis purpose. The application of a trend detection framework resulted in the identification of some significant trends, especially in January, February, and September precipitations and in the annual means. A noticeable decrease in the annual mean precipitation was observed mostly in western and southern Turkey, as well as along the coasts of the Black Sea. Regional average series also displayed trends similar to those for individual stations. Copyright © 2005 John Wiley & Sons, Ltd.