Balance conditions in variational data assimilation for a high-resolution forecast model

This paper explores the role of balance relationships for background error covariance modelling as the model's grid box decreases to convective-scales. Data assimilation (DA) analyses are examined from a simplified convective-scale model and DA system (called ABC-DA) with a grid box size of 1.5km in a 2D 540km (longitude), 15km (height) domain. The DA experiments are performed with background error covariance matrices B modelled and calibrated by switching on/off linear balance (LB) and hydrostatic balance (HB), and by observing a subset of the ABC variables, namely v, meridional wind, r', scaled density (a pressure-like variable), and b', buoyancy (a temperature-like variable). Calibration data are sourced from two methods of generating proxies of forecast errors. One uses forecasts from different latitude slices of a 3D parent model (here called the `latitude slice method'), and the other uses sets of differences between forecasts of different lengths but valid at the same time (the National Meteorological Center method). Root-mean-squared errors computed over the domain from identical twin DA experiments suggest that there is no combination of LB/HB switches that give the best analysis for all model quantities. It is frequently found though that the B-matrices modelled with both LB and HB do perform the best. A clearer picture emerges when the errors are examined at different spatial scales. In particular it is shown that switching on HB in B mostly has a neutral/positive effect on the DA accuracy at `large' scales, and switching off the HB has a neutral/positive effect at `small' scales. The division between `large' and `small' scales is between 10 and 100km. Furthermore, one hour forecast error correlations computed between control parameters find that correlations are small at large scales when balances are enforced, and at small scales when balances are not enforced (ideal control parameters have zero cross correlations). This points the way to modelling B with scale-dependent balances

River water level height measurements obtained from river cameras near Tewkesbury

This dataset contains extracted water level observations (WLOs) from four river cameras located on rivers Severn and Avon near Tewkesbury, UK, and owned by Farson Digital Ltd The data is extracted from hourly daylight river camera images between 21st Nov – 5th Dec 2012. The dataset consists of spreadsheets containing extracted water level data for each camera along with necessary metadata, all available Farson Digital Ltd river camera images between 21st Nov – 5th Dec 2012, and 3D point measurements of the location of each camera and locations of all measured points in the field-of-view for each camera

Ensemble data assimilation in the presence of cloud

In numerical weather prediction (NWP) data assimilation (DA) methods are used to combine available observations with numerical model estimates. This is done by minimising measures of error on both observations and model estimates with more weight given to data that can be more trusted. For any DA method an estimate of the initial forecast error covariance matrix is required. For convective scale data assimilation, however, the properties of the error covariances are not well understood. An effective way to investigate covariance properties in the presence of convection is to use an ensemble-based method for which an estimate of the error covariance is readily available at each time step. In this work, we investigate the performance of the ensemble square root filter (EnSRF) in the presence of cloud growth applied to an idealised 1D convective column model of the atmosphere. We show that the EnSRF performs well in capturing cloud growth, but the ensemble does not cope well with discontinuities introduced into the system by parameterised rain. The state estimates lose accuracy, and more importantly the ensemble is unable to capture the spread (variance) of the estimates correctly. We also find, counter-intuitively, that by reducing the spatial frequency of observations and/or the accuracy of the observations, the ensemble is able to capture the states and their variability successfully across all regimes

Advances in the study of boundary value problems for nonlinear integrable PDEs

In this review I summarise some of the most significant advances of the last decade in the analysis and solution of boundary value problems for integrable partial differential equations in two independent variables. These equations arise widely in mathematical physics, and in order to model realistic applications, it is essential to consider bounded domain and inhomogeneous boundary conditions. I focus specifically on a general and widely applicable approach, usually referred to as the Unified Transform or Fokas Transform, that provides a substantial generalisation of the classical Inverse Scattering Transform. This approach preserves the conceptual efficiency and aesthetic appeal of the more classical transform approaches, but presents a distinctive and important difference. While the Inverse Scattering Transform follows the "separation of variables" philosophy, albeit in a nonlinear setting, the Unified Transform is a based on the idea of synthesis, rather than separation, of variables. I will outline the main ideas in the case of linear evolution equations, and then illustrate their generalisation to certain nonlinear cases of particular significance