Search

Methane is the second-most important greenhouse gas after carbon dioxide and accounts for around 10 % of total European Union greenhouse gas emissions. Given that the atmospheric methane budget over a region depends on its terrestrial and aquatic methane sources, inverse modelling techniques appear as powerful tools for identifying critical areas that can later be submitted to emission mitigation strategies. In this regard, an inverse modelling system of methane emissions for Europe is being implemented based on the Weather Research and Forecasting (WRF) model: the Aarhus University Methane Inversion Algorithm (AUMIA) v1.0. The forward modelling component of AUMIA consists of the WRF model coupled to a multipurpose global database of methane anthropogenic emissions. To assure transport consistency during the inversion process, the backward modelling component will be based on the WRF model coupled to a Lagrangian particle dispersion module. A description of the modelling tools, input data sets, and 1-year forward modelling evaluation from 1 April 2018 to 31 March 2019 is provided in this paper. The a posteriori methane emission estimates, including a more focused inverse modelling for Denmark, will be provided in a second paper. A good general agreement is found between the modelling results and observations based on the TROPOspheric Monitoring Instrument (TROPOMI) onboard the Sentinel-5 Precursor satellite. Model–observation discrepancies for the summer peak season are in line with previous studies conducted over urban areas in central Europe, with relative differences between simulated concentrations and observational data in this study ranging from 1 % to 2 %. Domain-wide correlation coefficients and root-mean-square errors for summer months ranged from 0.4 to 0.5 and from 27 to 30 ppb, respectively. On the other hand, model–observation discrepancies for winter months show a significant overestimation of anthropogenic emissions over the study region, with relative differences ranging from 2 % to 3 %. Domain-wide correlation coefficients and root-mean-square errors in this case ranged from 0.1 to 0.4 and from 33 to 50 ppb, respectively, indicating that a more refined inverse analysis assessment will be required for this season. According to modelling results, the methane enhancement above the background concentrations came almost entirely from anthropogenic sources; however, these sources contributed with only up to 2 % to the methane total-column concentration. Contributions from natural sources (wetlands and termites) and biomass burning were not relevant during the study period. The results found in this study contribute with a new model evaluation of methane concentrations over Europe and demonstrate a huge potential for methane inverse modelling using improved TROPOMI products in large-scale applications. [ABSTRACT FROM AUTHOR]

In this study, we assess the skill of a stochastic weather generator (SWG) to forecast precipitation in several cities in western Europe. The SWG is based on a random sampling of analogs of the geopotential height at 500 hPa (Z500). The SWG is evaluated for two reanalyses (NCEP and ERA5). We simulate 100-member ensemble forecasts on a daily time increment. We evaluate the performance of SWG with forecast skill scores and we compare it to ECMWF forecasts. Results show significant positive skill score (continuous rank probability skill score and correlation) compared with persistence and climatology forecasts for lead times of 5 and 10 d for different areas in Europe. We find that the low predictability episodes of our model are related to specific weather regimes, depending on the European region. Comparing the SWG forecasts to ECMWF forecasts, we find that the SWG shows a good performance for 5 d. This performance varies from one region to another. This paper is a proof of concept for a stochastic regional ensemble precipitation forecast. Its parameters (e.g., region for analogs) must be tuned for each region in order to optimize its performance. [ABSTRACT FROM AUTHOR]

This is the second of two papers that document the creation of the New European Wind Atlas (NEWA). In Part 1, we described the sensitivity experiments and accompanying evaluation done to arrive at the final mesoscale model setup used to produce the mesoscale wind atlas. In this paper, Part 2, we document how we made the final wind atlas product, covering both the production of the mesoscale climatology generated with the Weather Research and Forecasting (WRF) model and the microscale climatology generated with the Wind Atlas Analysis and Applications Program (WAsP). The paper includes a detailed description of the technical and practical aspects that went into running the mesoscale simulations and the downscaling using WAsP. We show the main results from the final wind atlas and present a comprehensive evaluation of each component of the NEWA model chain using observations from a large set of tall masts located all over Europe. The added value of the WRF and WAsP downscaling of wind climatologies is evaluated relative to the performance of the driving ERA5 reanalysis and shows that the WRF downscaling reduces the mean wind speed bias and spread relative to that of ERA5 from -1.50±1.30 to 0.02±0.78 m s -1. The WAsP downscaling has an added positive impact relative to that of the WRF model in simple terrain. In complex terrain, where the assumptions of the linearized flow model break down, both the mean bias and spread in wind speed are worse than those from the raw mesoscale results. [ABSTRACT FROM AUTHOR]

This is the first of two papers that document the creation of the New European Wind Atlas (NEWA). It describes the sensitivity analysis and evaluation procedures that formed the basis for choosing the final setup of the mesoscale model simulations of the wind atlas. The suitable combination of model setup and parameterizations, bound by practical constraints, was found for simulating the climatology of the wind field at turbine-relevant heights with the Weather Research and Forecasting (WRF) model. Initial WRF model sensitivity experiments compared the wind climate generated by using two commonly used planetary boundary layer schemes and were carried out over several regions in Europe. They confirmed that the most significant differences in annual mean wind speed at 100 m a.g.l. (above ground level) mostly coincide with areas of high surface roughness length and not with the location of the domains or maximum wind speed. Then an ensemble of more than 50 simulations with different setups for a single year was carried out for one domain covering northern Europe for which tall mast observations were available. We varied many different parameters across the simulations, e.g. model version, forcing data, various physical parameterizations, and the size of the model domain. These simulations showed that although virtually every parameter change affects the results in some way, significant changes in the wind climate in the boundary layer are mostly due to using different physical parameterizations, especially the planetary boundary layer scheme, the representation of the land surface, and the prescribed surface roughness length. Also, the setup of the simulations, such as the integration length and the domain size, can considerably influence the results. We assessed the degree of similarity between winds simulated by the WRF ensemble members and the observations using a suite of metrics, including the Earth Mover's Distance (EMD), a statistic that measures the distance between two probability distributions. The EMD was used to diagnose the performance of each ensemble member using the full wind speed and direction distribution, which is essential for wind resource assessment. We identified the most realistic ensemble members to determine the most suitable configuration to be used in the final production run, which is fully described and evaluated in the second part of this study. [ABSTRACT FROM AUTHOR]

The prediction of climate on time scales of years to decades is attracting the interest of both climate researchers and stakeholders. The German Ministry for Education and Research (BMBF) has launched a major research programme on decadal climate prediction called MiKlip (Mittelfristige Klimaprognosen, Decadal Climate Prediction) in order to investigate the prediction potential of global and regional climate models (RCMs). In this paper we describe a regional predictive hindcast ensemble, its validation, and the added value of regional downscaling. Global predictions are obtained from an ensemble of simulations by the MPI-ESM-LR model (baseline 0 runs), which were downscaled for Europe using the COSMO-CLM regional model. Decadal hindcasts were produced for the 5 decades starting in 1961 until 2001. Observations were taken from the E-OBS data set. To identify decadal variability and predictability, we removed the long-term mean, as well as the long-term linear trend from the data. We split the resulting anomaly time series into two parts, the first including lead times of 1-5 years, reflecting the skill which originates mainly from the initialisation, and the second including lead times from 6-10 years, which are more related to the representation of low frequency climate variability and the effects of external forcing. We investigated temperature averages and precipitation sums for the summer and winter half-year. Skill assessment was based on correlation coefficient and reliability. We found that regional downscaling preserves, but mostly does not improve the skill and the reliability of the global predictions for summer half-year temperature anomalies. In contrast, regionalisation improves global decadal predictions of half-year precipitation sums in most parts of Europe. The added value results from an increased predictive skill on grid-point basis together with an improvement of the ensemble spread, i.e. the reliability. [ABSTRACT FROM AUTHOR]

Biogenic volatile organic compounds (BVOCs) emitted from the natural ecosystem are highly reactive and can thus impact air quality and aerosol radiative forcing. BVOC emission models (e.g., Model of Emissions of Gases and Aerosols from Nature – MEGAN) in global and regional chemical transport models still have large uncertainties in estimating biogenic trace gases because of uncertainties in emission activity factors, specification of vegetation type, and plant emission factors. This study evaluates a set of updates made to MEGAN v2.04 in the Weather Research and Forecasting model coupled with chemistry (WRF-Chem version 3.9). Our study considers four simulations for each update made to MEGAN v2.04: (i) a control run with no changes to MEGAN, (ii) a simulation with the emission activity factors modified following MEGAN v2.10, (iii) a simulation considering the changes to the plant functional type (PFT) emission factor, and (iv) a simulation with the isoprene emission factor calculated within the MEGAN module instead of being prescribed by the input database. We evaluate two regions, Europe and the southeastern United States, by comparing WRF-Chem results to ground-based monitoring observations in Europe (i.e., AirBase database) and aircraft observations obtained during the NOMADSS field campaign. We find that the updates to MEGAN v2.04 in WRF-Chem caused overpredictions in ground-based ozone concentrations in Europe and in isoprene mixing ratios compared to aircraft observations in the southeastern US. The update in emission activity factors caused the largest biases. These results suggest that further experimental and modeling studies should be conducted to address potential shortcomings in BVOC emission models. [ABSTRACT FROM AUTHOR]

The size distribution of atmospheric aerosols plays a key role for understanding and quantifying the uncertainties related to aerosol–radiation and aerosol–cloud interactions. These interactions ultimately depend on the size distribution through optical properties (such as aerosol optical depth, AOD) or cloud microphysical properties. Hence, the main objective of this contribution is to disentangle the impact of the representation of aerosol size distribution on aerosol optical properties over central Europe, particularly over the Mediterranean Basin, during a summertime aerosol episode. To fulfill this objective, a sensitivity test has been conducted using the coupled chemistry–meteorology model WRF-Chem (Weather Research Forecast model coupled with Chemistry). The test modified the parameters defining a lognormal size distribution (geometric diameter and standard deviation) by 10 %, 20 %, and 50 %. Results reveal that the reduction in the standard deviation of the accumulation mode leads to the largest impacts on AOD due to a transfer of particles from the accumulation mode to the coarse mode. A reduction in the geometric diameter of the accumulation mode also has an influence on AOD representation since particles in this mode are assumed to be smaller. In addition, an increase in the geometric diameter of the coarse mode produces a redistribution through the total size distribution by relocating particles from the finer modes to the coarse. [ABSTRACT FROM AUTHOR]

In this work we present downscaling experiments with the Weather Research and Forecasting model (WRF) to test the sensitivity to resolving aerosol–radiation and aerosol–cloud interactions on simulated regional climate for the EURO-CORDEX domain. The sensitivities mainly focus on the aerosol–radiation interactions (direct and semi-direct effects) with four different aerosol optical depth datasets (Tegen, MAC-v1, MACC, GOCART) being used and changes to the aerosol absorptivity (single scattering albedo) being examined. Moreover, part of the sensitivities also investigates aerosol–cloud interactions (indirect effect). Simulations have a resolution of 0.44 ∘ and are forced by the ERA-Interim reanalysis. A basic evaluation is performed in the context of seasonal-mean comparisons to ground-based (E-OBS) and satellite-based (CM SAF SARAH, CLARA) benchmark observational datasets. The impact of aerosols is calculated by comparing it against a simulation that has no aerosol effects. The implementation of aerosol–radiation interactions reduces the direct component of the incoming surface solar radiation by 20 %–30 % in all seasons, due to enhanced aerosol scattering and absorption. Moreover the aerosol–radiation interactions increase the diffuse component of surface solar radiation in both summer (30 %–40 %) and winter (5 %–8 %), whereas the overall downward solar radiation at the surface is attenuated by 3 %–8 %. The resulting aerosol radiative effect is negative and is comprised of the net effect from the combination of the highly negative direct aerosol effect (-17 to -5 W m -2) and the small positive changes in the cloud radiative effect (+5 W m -2), attributed to the semi-direct effect. The aerosol radiative effect is also stronger in summer (-12 W m -2) than in winter (-2 W m -2). We also show that modelling aerosol–radiation and aerosol–cloud interactions can lead to small changes in cloudiness, mainly regarding low-level clouds, and circulation anomalies in the lower and mid-troposphere, which in some cases, mainly close to the Black Sea in autumn, can be of statistical significance. Precipitation is not affected in a consistent pattern throughout the year by the aerosol implementation, and changes do not exceed ± 5 % except for the case of unrealistically absorbing aerosol. Temperature, on the other hand, systematically decreases by -0.1 to -0.5 ∘ C due to aerosol–radiation interactions with regional changes that can be up to -1.5 ∘ C. [ABSTRACT FROM AUTHOR]

In this study we compare monthly gross primary productivity (GPP) time series (2000-2007), computed for Europe with the Biosphere Energy Transfer Hydrology (BETHY/DLR) model with monthly data from the eddy covariance measurements network FLUXNET. BETHY/DLR with a spatial resolution of 1 km2 is designed for regional and continental applications (here Europe) and operated at the German Aerospace Center (DLR). It was adapted from the BETHY scheme to be driven by remote sensing data (leaf area index (LAI) and land cover information) and meteorology. Time series of LAI obtained from the CYCLOPES database are used to control the phenology of vegetation. Meteorological time series from the European Centre for Medium-Range Weather Forecasts (ECMWF) are used as driver. These comprise daily information on temperature, precipitation, wind speed and radiation. Additionally, static maps such as land cover, elevation, and soil type are used. To validate our model results we used eddy covariance measurements from the FLUXNET network of 74 towers across Europe. For forest sites we found that our model predicts between 20 and 40% higher annual GPP sums. In contrast, for cropland sites BETHY/DLR results show about 18% less GPP than eddy covariance measurements. For grassland sites, between 10% more and 16% less GPP was calculated with BETHY/DLR. A mean total carbon uptake of 2.5 PgC a-1 (±0.17 PgC a-1) was found for Europe. In addition, this study reports on risks that arise from the comparison of modelled data to FLUXNET measurements and their interpretation width. Furthermore we investigate reasons for uncertainties in model results and focus here on Vmax values, and finally embed our results into a broader context of model validation studies published during the last years in order to evaluate differences or similarities in analysed error sources. [ABSTRACT FROM AUTHOR]