Assimilation of Meteosat Third Generation (MTG) Lightning Imager (LI) pseudo-observations in AROME-France – proof of concept

This study develops a lightning data assimilation (LDA) scheme for the regional, convection-permitting numerical weather prediction (NWP) model AROME-France. The LDA scheme intends to assimilate total lightning, i.e., cloud-to-ground (CG) and inter- and intra-cloud (IC), of the future Meteosat Third Generation (MTG) Lightning Imager (LI; MTG-LI). MTG-LI proxy data are created, and flash extent density (FED) fields are derived. An FED forward observation operator (FFO) is trained based on modeled, column-integrated graupel mass from 24 storm days in 2018. The FFO is successfully verified for 2 independent storm days. With the FFO, the LDA adapts a 1-dimensional Bayesian (1DBay) retrieval followed by a 3-dimensional variational (3DVar) assimilation approach that is currently run operationally in AROME-France for radar reflectivity data. The 1DBay retrieval derives relative humidity profiles from the background by comparing the FED observations to the FED inferred from the background. Retrieved relative humidity profiles are assimilated as sounding data. The evaluation of the LDA comprises different LDA experiments and four case studies. It is found that all LDA experiments can increase the background integrated water vapor (IWV) in regions where the observed FED exceeds the FED inferred from AROME-France outputs. In addition, IWV can be reduced where spurious FED is modeled. A qualitative analysis of 6 h accumulated rainfall fields reveals that the LDA is capable of locating and initiating some local precipitation fields better than a radar data assimilation (RDA) experiment. However, the LDA also leads to rainfall accumulations that are too high at some locations. Fractions skill scores (FSSs) of 6 h accumulated rainfall are overall similar for the developed LDA and RDA experiments. An approach aiming at mitigating effects due to differences in the optical extents of lightning flashes and the area of the corresponding cloud was developed and included in the LDA; however, it does not always improve the FSS.</p

The AROME-WMED reanalyses of the first special observation period of the Hydrological cycle in the Mediterranean experiment (HyMeX)

To study key processes of the water cycle, two special observation periods (SOPs) of the Hydrological cycle in the Mediterranean experiment (HyMeX) took place during autumn 2012 and winter 2013. The first SOP aimed to study high precipitation systems and flash flooding in the Mediterranean area. The AROME-WMED (western Mediterranean) model (Fourrié et al., 2015) is a dedicated version of the mesoscale Numerical Weather Prediction (NWP) AROME-France model, which covers the western Mediterranean basin providing the HyMeX operational center with daily real-time analyses and forecasts. These products allowed for adequate decision-making for the field campaign observation deployment and the instrument operation. Shortly after the end of the campaign, a first reanalysis with more observations was performed with the first SOP operational software. An ensuing comprehensive second reanalysis of the first SOP, which included field research observations (not assimilated in real time) and some reprocessed observation datasets, was made with AROME-WMED. Moreover, a more recent version of the AROME model was used with updated background error statistics for the assimilation process. This paper depicts the main differences between the real-time version and the benefits brought by HyMeX reanalyses with AROME-WMED. The first reanalysis used 9 % additional data and the second one 24 % more compared to the real-time version. The second reanalysis is found to be closer to observations than the previous AROME-WMED analyses. The second reanalysis forecast errors of surface parameters are reduced up to the 18 and 24 h forecast range. In the middle and upper troposphere, fields are also improved up to the 48 h forecast range when compared to radiosondes. Integrated water vapor comparisons indicate a positive benefit for at least 24 h. Precipitation forecasts are found to be improved with the second reanalysis for a threshold up to 10 mm (24 h)-1. For higher thresholds, the frequency bias is degraded. Finally, improvement brought by the second reanalysis is illustrated with the Intensive Observation Period (IOP8) associated with heavy precipitation over eastern Spain and southern France

Long-term observations minus background monitoring of ground-based brightness temperatures from a microwave radiometer network

Abstract. Ground-based microwave radiometers (MWRs) offer the capability to provide continuous, high-temporal-resolution observations of the atmospheric thermodynamic state in the planetary boundary layer (PBL) with low maintenance. This makes MWR an ideal instrument to supplement radiosonde and satellite observations when initializing numerical weather prediction (NWP) models through data assimilation. State-of-the-art data assimilation systems (e.g. variational schemes) require an accurate representation of the differences between model (background) and observations, which are then weighted by their respective errors to provide the best analysis of the true atmospheric state. In this perspective, one source of information is contained in the statistics of the differences between observations and their background counterparts (O–B). Monitoring of O–B statistics is crucial to detect and remove systematic errors coming from the measurements, the observation operator, and/or the NWP model. This work illustrates a 1-year O–B analysis for MWR observations in clear-sky conditions for an European-wide network of six MWRs. Observations include MWR brightness temperatures (TB) measured by the two most common types of MWR instruments. Background profiles are extracted from the French convective-scale model AROME-France before being converted into TB. The observation operator used to map atmospheric profiles into TB is the fast radiative transfer model RTTOV-gb. It is shown that O–B monitoring can effectively detect instrument malfunctions. O–B statistics (bias, standard deviation, and root mean square) for water vapour channels (22.24–30.0 GHz) are quite consistent for all the instrumental sites, decreasing from the 22.24 GHz line centre ( ∼  2–2.5 K) towards the high-frequency wing ( ∼  0.8–1.3 K). Statistics for zenith and lower-elevation observations show a similar trend, though values increase with increasing air mass. O–B statistics for temperature channels show different behaviour for relatively transparent (51–53 GHz) and opaque channels (54–58 GHz). Opaque channels show lower uncertainties (< 0.8–0.9 K) and little variation with elevation angle. Transparent channels show larger biases ( ∼  2–3 K) with relatively low standard deviations ( ∼  1–1.5 K). The observations minus analysis TB statistics are similar to the O–B statistics, suggesting a possible improvement to be expected by assimilating MWR TB into NWP models. Lastly, the O–B TB differences have been evaluated to verify the normal-distribution hypothesis underlying variational and ensemble Kalman filter-based DA systems. Absolute values of excess kurtosis and skewness are generally within 1 and 0.5, respectively, for all instrumental sites, demonstrating O–B normal distribution for most of the channels and elevations angles

Hydrometeorological multi-model ensemble simulations of the 4 November 2011 flash flood event in Genoa, Italy, in the framework of the DRIHM Project

The e-Science environment developed in the framework of the EU-funded DRIHM project was used to demonstrate its ability to provide relevant, meaningful hydrometeorological forecasts. This was illustrated for the tragic case of 4 November 2011, when Genoa, Italy, was flooded as the result of heavy, convective precipitation that inundated the Bisagno catchment. The Meteorological Model Bridge (MMB), an innovative software component developed within the DRIHM project for the interoperability of meteorological and hydrological models, is a key component of the DRIHM e-Science environment. The MMB allowed three different rainfall-discharge models (DRiFt, RIBS and HBV) to be driven by four mesoscale limited-area atmospheric models (WRF-NMM, WRF-ARW, Meso-NH and AROME) and a downscaling algorithm (RainFARM) in a seamless fashion. In addition to this multi-model configuration, some of the models were run in probabilistic mode, thus giving a comprehensive account of modelling errors and a very large amount of likely hydrometeorological scenarios (> 1500). The multi-model approach proved to be necessary because, whilst various aspects of the event were successfully simulated by different models, none of the models reproduced all of these aspects correctly. It was shown that the resulting set of simulations helped identify key atmospheric processes responsible for the large rainfall accumulations over the Bisagno basin. The DRIHM e-Science environment facilitated an evaluation of the sensitivity to atmospheric and hydrological modelling errors. This showed that both had a significant impact on predicted discharges, the former being larger than the latter. Finally, the usefulness of the set of hydrometeorological simulations was assessed from a flash flood early-warning perspective

Molecular Characterization of HIV-1 CRF01_AE in Mekong Delta, Vietnam, and Impact of T-Cell Epitope Mutations on HLA Recognition (ANRS 12159)

To date, 11 HIV-1 subtypes and 48 circulating recombinant forms have been described worldwide. The underlying reason why their distribution is so heterogeneous is not clear. Host genetic factors could partly explain this distribution. The aim of this study was to describe HIV-1 strains circulating in an unexplored area of Mekong Delta, Vietnam, and to assess the impact of optimal epitope mutations on HLA binding.We recruited 125 chronically antiretroviral-naive HIV-1-infected subjects from five cities in the Mekong Delta. We performed high-resolution DNA typing of HLA class I alleles, sequencing of Gag and RT-Prot genes and phylogenetic analysis of the strains. Epitope mutations were analyzed in patients bearing the HLA allele restricting the studied epitope. Optimal wild-type epitopes from the Los Alamos database were used as reference. T-cell epitope recognition was predicted using the immune epitope database tool according to three different scores involved in antigen processing (TAP and proteasome scores) and HLA binding (MHC score). with a Vietnamese specificity held by two different haplotypes. The percentage of homology between Mekong and B consensus HIV-1 sequences was above 85%. Divergent epitopes had TAP and proteasome scores comparable with wild-type epitopes. MHC scores were significantly lower in divergent epitopes with a mean of 2.4 (±0.9) versus 2 (±0.7) in non-divergent ones (p<0.0001).Our study confirms the wide predominance of CRF01_AE in the Mekong Delta where patients harbor a specific HLA pattern. Moreover, it demonstrates the lower MHC binding affinity among divergent epitopes. This weak immune pressure combined with a narrow genetic diversity favors immune escape and could explain why CRF01_AE is still predominant in Vietnam, particularly in the Mekong area

The HIV-1 Integrase α4-Helix Involved in LTR-DNA Recognition Is also a Highly Antigenic Peptide Element

Monoclonal antibodies (MAbas) constitute remarkable tools to analyze the relationship between the structure and the function of a protein. By immunizing a mouse with a 29mer peptide (K159) formed by residues 147 to 175 of the HIV-1 integrase (IN), we obtained a monoclonal antibody (MAba4) recognizing an epitope lying in the N-terminal portion of K159 (residues 147–166 of IN). The boundaries of the epitope were determined in ELISA assays using peptide truncation and amino acid substitutions. The epitope in K159 or as a free peptide (pep-a4) was mostly a random coil in solution, while in the CCD (catalytic core domain) crystal, the homologous segment displayed an amphipathic helix structure (α4-helix) at the protein surface. Despite this conformational difference, a strong antigenic crossreactivity was observed between pep-a4 and the protein segment, as well as K156, a stabilized analogue of pep-a4 constrained into helix by seven helicogenic mutations, most of them involving hydrophobic residues. We concluded that the epitope is freely accessible to the antibody inside the protein and that its recognition by the antibody is not influenced by the conformation of its backbone and the chemistry of amino acids submitted to helicogenic mutations. In contrast, the AA →Glu mutations of the hydrophilic residues Gln148, Lys156 and Lys159, known for their interactions with LTRs (long terminal repeats) and inhibitors (

Modelling of the effect of ELMs on fuel retention at the bulk W divertor of JET

Effect of ELMs on fuel retention at the bulk W target of JET ITER-Like Wall was studied with multi-scale calculations. Plasma input parameters were taken from ELMy H-mode plasma experiment. The energetic intra-ELM fuel particles get implanted and create near-surface defects up to depths of few tens of nm, which act as the main fuel trapping sites during ELMs. Clustering of implantation-induced vacancies were found to take place. The incoming flux of inter-ELM plasma particles increases the different filling levels of trapped fuel in defects. The temperature increase of the W target during the pulse increases the fuel detrapping rate. The inter-ELM fuel particle flux refills the partially emptied trapping sites and fills new sites. This leads to a competing effect on the retention and release rates of the implanted particles. At high temperatures the main retention appeared in larger vacancy clusters due to increased clustering rate

The effect of beryllium oxide on retention in JET ITER-like wall tiles

Preliminary results investigating the microstructure, bonding and effect of beryllium oxide formation on retention in the JET ITER-like wall beryllium tiles, are presented. The tiles have been investigated by several techniques: Scanning Electron Microscopy (SEM) equipped with Energy Dispersive X-ray (EDX), Transmission Electron microscopy (TEM) equipped with EDX and Electron Energy Loss Spectroscopy (EELS), Raman Spectroscopy and Thermal Desorption Spectroscopy (TDS). This paper focuses on results from melted materials of the dump plate tiles in JET. From our results and the literature, it is concluded, beryllium can form micron deep oxide islands contrary to the nanometric oxides predicted under vacuum conditions. The deepest oxides analyzed were up to 2-micron thicknesses. The beryllium Deuteroxide (BeOxDy) bond was found with Raman Spectroscopy. Application of EELS confirmed the oxide presence and stoichiometry. Literature suggests these oxides form at temperatures greater than 700 °C where self-diffusion of beryllium ions through the surface oxide layer can occur. Further oxidation is made possible between oxygen plasma impurities and the beryllium ions now present at the wall surface. Under Ultra High Vacuum (UHV) nanometric Beryllium oxide layers are formed and passivate at room temperature. After continual cyclic heating (to the point of melt formation) in the presence of oxygen impurities from the plasma, oxide growth to the levels seen experimentally (approximately two microns) is proposed. This retention mechanism is not considered to contribute dramatically to overall retention in JET, due to low levels of melt formation. However, this mechanism, thought the result of operation environment and melt formation, could be of wider concern to ITER, dependent on wall temperatures

Modelling of tungsten erosion and deposition in the divertor of JET-ILW in comparison to experimental findings

The erosion, transport and deposition of tungsten in the outer divertor of JET-ILW has been studied for an HMode discharge with low frequency ELMs. For this specific case with an inter-ELM electron temperature at the strike point of about 20 eV, tungsten sputtering between ELMs is almost exclusively due to beryllium impurity and self-sputtering. However, during ELMs tungsten sputtering due to deuterium becomes important and even dominates. The amount of simulated local deposition of tungsten relative to the amount of sputtered tungsten in between ELMs is very high and reaches values of 99% for an electron density of 5E13 cm$^{-3}$ at the strike point and electron temperatures between 10 and 30 eV. Smaller deposition values are simulated with reduced electron density. The direction of the B-field significantly influences the local deposition and leads to a reduction if the E×B drift directs towards the scrape-off-layer. Also, the thermal force can reduce the tungsten deposition, however, an ion temperature gradient of about 0.1 eV/mm or larger is needed for a significant effect. The tungsten deposition simulated during ELMs reaches values of about 98% assuming ELM parameters according to free-streaming model. The measured WI emission profiles in between and within ELMs have been reproduced by the simulation. The contribution to the overall net tungsten erosion during ELMs is about 5 times larger than the one in between ELMs for the studied case. However, this is due to the rather low electron temperature in between ELMs, which leads to deuterium impact energies below the sputtering threshold for tungsten