Search

Air-pollution modelling at the local scale requires accurate meteorological inputs such as from the velocity field. These meteorological fields are generally simulated with microscale models (here Code_Saturne), which are forced with boundary conditions provided by larger scale models or observations. Local atmospheric simulations are very sensitive to the boundary conditions, whose accurate estimation is difficult but crucial. When observations of the wind speed and turbulence or pollutant concentration are available inside the domain, they provide supplementary information via data assimilation, to enhance the simulation accuracy by modifying the boundary conditions. Among the existing data assimilation methods, the iterative ensemble Kalman smoother (IEnKS) is adapted to urban-scale simulations. This method has already been found to increase the accuracy of wind-resource assessment. Here we assess the ability of the IEnKS method to improve scalar-dispersion modelling—an important component of air-quality modelling—by assimilating perturbed measurements inside the urban canopy. To test the data assimilation method in urban conditions, we use the observations provided by the Mock Urban Setting Test field campaign and consider cases with neutral and stable conditions, and the boundary conditions consisting of the horizontal velocity components and turbulence. We prove the capacity of the IEnKS method to assimilate observations of velocity as well as pollutant concentration. In both cases, the accuracy of pollutant concentration estimates is enhanced by 40–60%. We also show that assimilating both types of observations allows further improvements of turbulence predictions by the model. [ABSTRACT FROM AUTHOR]

Data assimilation for systems possessing many scales of motions is a substantial methodological and technological challenge. Systems with these features are found in many areas of computational physics and are becoming common thanks to increased computational power allowing to resolve finer scales and to couple together several sub-components. Coupled data assimilation (CDA) distinctively appears as a main concern in numerical weather and climate prediction with major efforts put forward by meteo services worldwide. The core issue is the scale separation acting as a barrier that hampers the propagation of the information across model components (e.g. ocean and atmosphere). We provide a brief survey of CDA, and then focus on CDA using the ensemble Kalman filter (EnKF), a widely used Monte Carlo Gaussian method. Our goal is to elucidate the mechanisms behind information propagation across model components. We consider first a coupled system of equations with temporal scale difference, and deduce that: (i) cross components effects are strong from the slow to the fast scale, but, (ii) intra-component effects are much stronger in the fast scale. While observing the slow scale is desirable and benefits the fast, the latter must be observed with high frequency otherwise the error will grow up to affect the slow scale. Numerical experiments are performed using the atmosphere-ocean model, MAOOAM. Six configurations are considered, differing for the strength of the atmosphere-ocean coupling and/or the number of model modes. The performance of the EnKF depends on the model configuration, i.e. on its dynamical features. A comprehensive dynamical characterisation of the model configurations is provided by examining the Lyapunov spectrum, Kolmogorov entropy and Kaplan–Yorke attractor dimension. We also compute the covariant Lyapunov vectors and use them to explain how model instabilities act on different model's modes according to the coupling strength. The experiments confirm the importance of observing the fast scale, but show also that, despite its slow temporal scale, frequent observations in the ocean are beneficial. The relation between the ensemble size, N, and the unstable subspace dimension, n 0 , has been studied. Results largely ratify what known for uncoupled system: the condition N ≥ n 0 is necessary for the EnKF to work satisfactorily. Nevertheless the quasi-degeneracy of the Lyapunov spectrum of MAOOAM, with many near-zero exponents, is potentially the cause of the smooth gradual reduction of the analysis error observed for some model configurations, even when N > n 0 . Future prospects for the EnKF in the context of coupled ocean-atmosphere systems are finally discussed. [ABSTRACT FROM AUTHOR]

We consider localisation problems belonging to the chiral symmetry classes, in which sublattice symmetry is responsible for singular behaviour at a band centre. As an account for the talk given in TH2002, we focus here on the chiral time-reversal invariant symmetry class. We formulate a model which has the relevant symmetries and which is a generalisation of the network model introduced previously in the context of the integer quantum Hall plateau transition. We show that this generalisation can be re-expressed as corresponding to the introduction of absorption and amplification into the original network model. This model represents a convenient starting point for analytic discussions and computational studies. It exhibits both localised and critical phases, and band-centre singularities in the critical phase approach more closely in small systems the expected asymptotic form than in other known realisations of the symmetry class. In addition, we demonstrate that by imposing appropriate constraints on disorder, a lattice version of the Dirac equation with a random vector potential can be obtained. [ABSTRACT FROM AUTHOR]