ALBERT

Description: Using data assimilation (DA) to improve model forecast accuracy is a powerful approach that requires available observations. Infrared satellite measurements of volcanic ash mass loadings are often used as input observations for the assimilation scheme. However, because these primary satellite-retrieved data are often two-dimensional (2-D) and the ash plume is usually vertically located in a narrow band, directly assimilating the 2-D ash mass loadings in a three-dimensional (3-D) volcanic ash model (with an integral observational operator) can usually introduce large artificial/spurious vertical correlations.In this study, we look at an approach to avoid the artificial vertical correlations by not involving the integral operator. By integrating available data of ash mass loadings and cloud top heights, as well as data-based assumptions on thickness, we propose a satellite observational operator (SOO) that translates satellite-retrieved 2-D volcanic ash mass loadings to 3-D concentrations. The 3-D SOO makes the analysis step of assimilation comparable in the 3-D model space.Ensemble-based DA is used to assimilate the extracted measurements of ash concentrations. The results show that satellite DA with SOO can improve the estimate of volcanic ash state and the forecast. Comparison with both satellite-retrieved data and aircraft in situ measurements shows that the effective duration of the improved volcanic ash forecasts for the distal part of the Eyjafjallajökull volcano is about 6 h.

Description: In this study, we investigate a strategy to accelerate the data assimilation (DA) algorithm. Based on evaluations of the computational time, the analysis step of the assimilation turns out to be the most expensive part. After a study of the characteristics of the ensemble ash state, we propose a mask-state algorithm which records the sparsity information of the full ensemble state matrix and transforms the full matrix into a relatively small one. This will reduce the computational cost in the analysis step. Experimental results show the mask-state algorithm significantly speeds up the analysis step. Subsequently, the total amount of computing time for volcanic ash DA is reduced to an acceptable level. The mask-state algorithm is generic and thus can be embedded in any ensemble-based DA framework. Moreover, ensemble-based DA with the mask-state algorithm is promising and flexible, because it implements exactly the standard DA without any approximation and it realizes the satisfying performance without any change in the full model.

Description: The development and application of chemistry transport models has a long tradition. Within the Netherlands the LOTOS–EUROS model has been developed by a consortium of institutes, after combining its independently developed predecessors in 2005. Recently, version 2.0 of the model was released as an open-source version. This paper presents the curriculum vitae of the model system, describing the model's history, model philosophy, basic features and a validation with EMEP stations for the new benchmark year 2012, and presents cases with the model's most recent and key developments. By setting the model developments in context and providing an outlook for directions for further development, the paper goes beyond the common model description.With an origin in ozone and sulfur modelling for the models LOTOS and EUROS, the application areas were gradually extended with persistent organic pollutants, reactive nitrogen, and primary and secondary particulate matter. After the combination of the models to LOTOS–EUROS in 2005, the model was further developed to include new source parametrizations (e.g. road resuspension, desert dust, wildfires), applied for operational smog forecasts in the Netherlands and Europe, and has been used for emission scenarios, source apportionment, and long-term hindcast and climate change scenarios. LOTOS–EUROS has been a front-runner in data assimilation of ground-based and satellite observations and has participated in many model intercomparison studies. The model is no longer confined to applications over Europe but is also applied to other regions of the world, e.g. China. The increasing interaction with emission experts has also contributed to the improvement of the model's performance. The philosophy for model development has always been to use knowledge that is state of the art and proven, to keep a good balance in the level of detail of process description and accuracy of input and output, and to keep a good record on the effect of model changes using benchmarking and validation. The performance of v2.0 with respect to EMEP observations is good, with spatial correlations around 0.8 or higher for concentrations and wet deposition. Temporal correlations are around 0.5 or higher. Recent innovative applications include source apportionment and data assimilation, particle number modelling, and energy transition scenarios including corresponding land use changes as well as Saharan dust forecasting. Future developments would enable more flexibility with respect to model horizontal and vertical resolution and further detailing of model input data. This includes the use of different sources of land use characterization (roughness length and vegetation), detailing of emissions in space and time, and efficient coupling to meteorology from different meteorological models.

Description: Random walk models are a powerful tool for the investigation of transport processes in turbulent flows. However, standard random walk methods are applicable only when the flow velocities and diffusivity are sufficiently smooth functions. In practice there are some regions where the rapid but continuous change in diffusivity may be represented by a discontinuity. The random walk model based on backward Îto calculus can be used for these problems. This model was proposed by LaBolle et al. (2000). The latter is best suited to the problems under consideration. It is then applied to two test cases with discontinuous diffusivity, highlighting the advantages of this method.

Description: Random walk models are a powerful tool for the investigation of transport processes in turbulent flows. However, standard random walk methods are applicable only when the flow velocities and diffusivity are sufficiently smooth functions. In practice there are some regions where the rapid but continuous change in diffusivity may be represented by a discontinuity. The random walk model based on backward Îto calculus can be used for these problems. This model was proposed by LaBolle et al. (2000). The latter is best suited to the problems under consideration. It is then applied for two test cases with discontinuous diffusivity, highlighting the advantages of this method.

Description: The development and application of chemistry transport models has a long tradition. Within the Netherlands the LOTOS-EUROS model has been developed by a consortium of institutes, after combination of its independently developed predecessors in 2005. Recently, version 2.0 of the model was released as an open source version. This paper presents the curriculum vitae of the model system, describing the model’s history, model philosophy, basic features, a validation with EMEP stations for the new benchmark year 2012, and presents cases with the model's most recent and key developments. By setting the model developments in context and providing an outlook for directions for further development, the paper goes beyond the common model description. With an origin in ozone and sulphur modelling for the models LOTOS and EUROS, the application areas were gradually extended with POPs, reactive nitrogen and primary and secondary particulate matter. After the combination of the models to LOTOS-EUROS in 2005, the model was further developed to include new source parametrizations (e.g. road resuspension, desert dust, wildfires), applied for operational smog forecasts in the Netherlands and Europe, and has been used for emission scenarios, source apportionment and long-term hindcast and climate change scenarios. LOTOS-EUROS has been a front-runner in data assimilation of ground-based and satellite observations and has participated in many model intercomparison studies. The model is no longer confined to applications over Europe but is also applied to other regions of the world, e.g. China. Also the increasing interaction with emission experts has contributed to the improvement of the model’s performance. The philosophy for model development has always been to use knowledge that is state of the art and proven, to keep good balance in the level of detail of process description and accuracy of input and output, and to keep a good track on the effect of model changes using benchmarking and validation. The performance of v2.0 with respect to EMEP observations is good, with spatial correlations around 0.8 or higher for concentrations and wet deposition. Temporal correlations are around 0.5 or higher. Recent innovative applications include source apportionment and data assimilation, particle number modelling, energy transition scenarios including corresponding land use changes as well as Saharan dust forecasting. Future developments would enable more flexibility with respect to model horizontal and vertical resolution and further detailing of model input data. This includes use of different sources of land use characterization (roughness length and vegetation), detailing of emissions in space and time, and efficient coupling to meteorology from different meteorological models.

Description: Data assimilation algorithms rely on a basic assumption of an unbiased observation error. However, the presence of inconsistent measurements with nontrivial biases or inseparable baselines is unavoidable in practice. Assimilation analysis might diverge from reality, since the data assimilation itself cannot distinguish whether the differences between model simulations and observations are due to the biased observations or model deficiencies. Unfortunately, modeling of observation biases or baselines which show strong spatiotemporal variability is a challenging task. In this study, we report how data-driven machine learning can be used to perform observation bias correction for data assimilation through a real application, which is the dust emission inversion using PM10 observations. PM10 observations are considered as unbiased, however, a bias correction is necessary if they are used as a proxy for dust during dust storms since they actually represent a sum of dust particles and non-dust aerosols. Two observation bias correction methods have been designed in order to use PM10 measurements as proxy for the dust storm loads under severe dust conditions. The first one is the conventional chemical transport (CTM) model that simulates life cycles of non-dust aerosols. The other one is the machine learning model that describes the relations between the regular PM10 and other air quality measurement. The latter is trained by learning two-year's historical samples. The machine learning based non-dust model is shown to be in better agreements with observations compared to the CTM. The dust emission inversion tests have been performed, either through assimilating the raw measurements, or the bias-corrected dust observations using either the CTM or machine learning model. The emission field, surface dust concentration and forecast skill are evaluated. The worst case is when we directly assimilate the original observations. The forecasts driven by the posterior emission in this case even results in larger errors than the reference prediction. This shows the necessities of bias correction in data assimilation. The best results are obtained when using a machine learning model for bias correction, with the existing measurements used more precisely and the resulting forecasts close to reality.

Description: Data assimilation is a powerful tool that requires available observations to improve model forecast accuracy. Infrared satellite measurements of volcanic ash mass loadings are often used as input observations into the assimilation scheme. However, these satellite-retrieved data are often two-dimensional (2D), and cannot be easily combined with a three-dimensional (3D) volcanic ash model to continuously improve the volcanic ash state in a data assimilation system. By integrating available data including ash mass loadings, cloud top heights and thickness information, we propose a satellite observational operator (SOO) that translates satellite-retrieved 2D volcanic ash mass loadings to 3D concentrations at the top layer of the ash cloud. Ensemble-based data assimilation is used to continuously assimilate the extracted measurements of ash concentrations. The results show that satellite data assimilation can force the volcanic ash state to match the satellite observations, and that it improves the forecast of the ash state. Comparison with highly accurate aircraft in-situ measurements shows that the effective duration of the improved volcanic ash forecasts is about a half day. It is shown that an effective half-day ash forecast significantly improves the quality of the advice given to aviation over continental Europe.

Description: The forecast accuracy of distal volcanic ash clouds is important for providing valid aviation advice during volcanic ash eruption. However, because the distal part of volcanic ash plume is far from the volcano, the influence of eruption information on this part becomes rather indirect and uncertain, resulting in inaccurate volcanic ash forecasts in these distal areas. In our approach, we use real-life aircraft in-situ observations, measured in the North-West part of Germany during the 2010 Eyjafjallajokull eruption, in an ensemble-based data assimilation system combined with a volcanic ash transport model to investigate the potential improvement on the forecast accuracy with regard to the distal volcanic ash plume. We show that the error of the analyzed volcanic ash state can be significantly reduced through assimilating real-life in-situ measurements. After a continuous assimilation, it is shown that the aviation advice for Germany, the Netherlands and Belgium can be significantly improved. We suggest that with suitable aircrafts measuring once per day across the distal volcanic ash plume, the description and prediction of volcanic ash clouds in these areas can be greatly improved.

Description: Data assimilation algorithms rely on a basic assumption of an unbiased observation error. However, the presence of inconsistent measurements with nontrivial biases or inseparable baselines is unavoidable in practice. Assimilation analysis might diverge from reality since the data assimilation itself cannot distinguish whether the differences between model simulations and observations are due to the biased observations or model deficiencies. Unfortunately, modeling of observation biases or baselines which show strong spatiotemporal variability is a challenging task. In this study, we report how data-driven machine learning can be used to perform observation bias correction for data assimilation through a real application, which is the dust emission inversion using PM10 observations. PM10 observations are considered unbiased; however, a bias correction is necessary if they are used as a proxy for dust during dust storms since they actually represent a sum of dust particles and non-dust aerosols. Two observation bias correction methods have been designed in order to use PM10 measurements as proxy for the dust storm loads under severe dust conditions. The first one is the conventional chemistry transport model (CTM) that simulates life cycles of non-dust aerosols. The other one is the machine-learning model that describes the relations between the regular PM10 and other air quality measurements. The latter is trained by learning using 2 years of historical samples. The machine-learning-based non-dust model is shown to be in better agreement with observations compared to the CTM. The dust emission inversion tests have been performed, through assimilating either the raw measurements or the bias-corrected dust observations using either the CTM or machine-learning model. The emission field, surface dust concentration, and forecast skill are evaluated. The worst case is when we directly assimilate the original observations. The forecasts driven by the a posteriori emission in this case even result in larger errors than the reference prediction. This shows the necessities of bias correction in data assimilation. The best results are obtained when using the machine-learning model for bias correction, with the existing measurements used more precisely and the resulting forecasts close to reality.