Agencia Estatal de Meteorología

Abstract Global climate change impacts can already be tracked in many physical and biological systems; in particular, terrestrial ecosystems provide a consistent picture of observed changes. One of the preferred indicators is phenology, the science of natural recurring events, as their recorded dates provide a high‐temporal resolution of ongoing changes. Thus, numerous analyses have demonstrated an earlier onset of spring events for mid and higher latitudes and a lengthening of the growing season. However, published single‐site or single‐species studies are particularly open to suspicion of being biased towards predominantly reporting climate change‐induced impacts. No comprehensive study or meta‐analysis has so far examined the possible lack of evidence for changes or shifts at sites where no temperature change is observed. We used an enormous systematic phenological network data set of more than 125 000 observational series of 542 plant and 19 animal species in 21 European countries (1971–2000). Our results showed that 78% of all leafing, flowering and fruiting records advanced (30% significantly) and only 3% were significantly delayed, whereas the signal of leaf colouring/fall is ambiguous. We conclude that previously published results of phenological changes were not biased by reporting or publication predisposition: the average advance of spring/summer was 2.5 days decade −1 in Europe. Our analysis of 254 mean national time series undoubtedly demonstrates that species' phenology is responsive to temperature of the preceding months (mean advance of spring/summer by 2.5 days°C −1 , delay of leaf colouring and fall by 1.0 day°C −1 ). The pattern of observed change in spring efficiently matches measured national warming across 19 European countries (correlation coefficient r =−0.69, P <0.001).

Abstract We present a dataset of daily resolution climatic time series that has been compiled for the European Climate Assessment (ECA). As of December 2001, this ECA dataset comprises 199 series of minimum, maximum and/or daily mean temperature and 195 series of daily precipitation amount observed at meteorological stations in Europe and the Middle East. Almost all series cover the standard normal period 1961–90, and about 50% extends back to at least 1925. Part of the dataset (90%) is made available for climate research on CDROM and through the Internet (at http://www.knmi.nl/samenw/eca ). A comparison of the ECA dataset with existing gridded datasets, having monthly resolution, shows that correlation coefficients between ECA stations and nearest land grid boxes between 1946 and 1999 are higher than 0.8 for 93% of the temperature series and for 51% of the precipitation series. The overall trends in the ECA dataset are of comparable magnitude to those in the gridded datasets. The potential of the ECA dataset for climate studies is demonstrated in two examples. In the first example, it is shown that the winter (October–March) warming in Europe in the 1976–99 period is accompanied by a positive trend in the number of warm‐spell days at most stations, but not by a negative trend in the number of cold‐spell days. Instead, the number of cold‐spell days increases over Europe. In the second example, it is shown for winter precipitation between 1946 and 1999 that positive trends in the mean amount per wet day prevail in areas that are getting drier and wetter. Because of its daily resolution, the ECA dataset enables a variety of empirical climate studies, including detailed analyses of changes in the occurrence of extremes in relation to changes in mean temperature and total precipitation. Copyright © 2002 Royal Meteorological Society.

The summer of 2010 was exceptionally warm in eastern Europe and large parts of Russia. We provide evidence that the anomalous 2010 warmth that caused adverse impacts exceeded the amplitude and spatial extent of the previous hottest summer of 2003. "Mega-heatwaves" such as the 2003 and 2010 events likely broke the 500-year-long seasonal temperature records over approximately 50% of Europe. According to regional multi-model experiments, the probability of a summer experiencing mega-heatwaves will increase by a factor of 5 to 10 within the next 40 years. However, the magnitude of the 2010 event was so extreme that despite this increase, the likelihood of an analog over the same region remains fairly low until the second half of the 21st century.

A multi-model ensemble-based system for seasonal-to-interannual prediction has been developed in a joint European project known as DEMETER (Development of a European Multimodel Ensemble Prediction System for Seasonal to Interannual Prediction). The DEMETER system comprises seven global atmosphere–ocean coupled models, each running from an ensemble of initial conditions. Comprehensive hindcast evaluation demonstrates the enhanced reliability and skill of the multimodel ensemble over a more conventional single-model ensemble approach. In addition, innovative examples of the application of seasonal ensemble forecasts in malaria and crop yield prediction are discussed. The strategy followed in DEMETER deals with important problems such as communication across disciplines, downscaling of climate simulations, and use of probabilistic forecast information in the applications sector, illustrating the economic value of seasonal-to-interannual prediction for society as a whole.

The recently initiated Arecibo Legacy Fast ALFA (ALFALFA) survey aims to map 7000 square degrees of the high galactic latitude sky visible from Arecibo, providing a HI line spectral database covering the redshift range between -1600 km/s and 18,000 km/s with 5 km/s resolution. Exploiting Arecibo's large collecting area and small beam size, ALFALFA is specifically designed to probe the faint end of the HI mass function in the local universe and will provide a census of HI in the surveyed sky area to faint flux limits, making it especially useful in synergy with wide area surveys conducted at other wavelengths. ALFALFA will also provide the basis for studies of the dynamics of galaxies within the Local and nearby superclusters, will allow measurement of the HI diameter function, and enable a first wide-area blind search for local HI tidal features, HI absorbers at z < 0.06 and OH megamasers in the redshift range 0.16 < z < 0.25. Although completion of the survey will require some five years, public access to the ALFALFA data and data products will be provided in a timely manner, thus allowing its application for studies beyond those targeted by the ALFALFA collaboration. ALFALFA adopts a two-pass, minimum intrusion, drift scan observing technique which samples the same region of sky at two separate epochs to aid in the discrimination of cosmic signals from noise and terrestrial interference. Survey simulations, which take into account large scale structure in the mass distribution and incorporate experience with the ALFA system gained from tests conducted during its commissioning phase, suggest that ALFALFA will detect on the order of 20,000 extragalactic HI line sources out to z=0.06, including several hundred with HI masses of less than 10^{7.5} msun.

Abstract The paper describes the turbulence scheme implemented in the Meso‐NH community research model, and reports on some validation studies. Since the model is intended to perform both large‐eddy and mesoscale simulations, we have developed a full three‐dimensional scheme, based on the original method of Redelsperger and Sommeria. A prognostic equation for the turbulent kinetic energy is used, together with conservative variables for moist non‐precipitating processes. A particularity of the scheme is the use of variable turbulent Prandtl and Schmidt numbers, consistently derived from the complete set of second‐order turbulent‐moment equations. The results of three idealized boundary‐layer simulations allowing detailed comparisons with other large‐eddy simulation (LES) models are discussed, and lead to the conclusion that the model is performing satisfactorily. The vertical flux and gradient computation can be run in isolation from the rest of the scheme, providing an efficient single‐column parametrization for the mesoscale configuration of the model, if an appropriate parametrization of the eddy length‐scale is used. The mixing‐length specification is then the only aspect of the scheme which differs from the LES to the mesoscale configuration, and the numerical constants used for the closure terms are the same in both configurations. The scheme is run in single‐column mode for the same three cases as above, and a comparison of single‐column and LES results again leads to satisfactory results. It is believed that this result is original, and is due to the proper formulation of the parametrized mixing length and of the turbulent Prandtl and Schmidt numbers. In fact, a comparison of the parametrized mixing length with the length‐scale of the energy‐containing eddies deduced by spectral analysis of the LES shows interesting similarity.

9 Pags., 6 Figs. The Content from this work may be used under the terms of the Creative Commons Attribution 3.0 licence. Any further
\ndistribution of this work must maintain attribution to the author(s) and the
\nof the work, journal citation and DOI.

The Tropospheric Ozone Assessment Report (TOAR) is an activity of the International Global Atmospheric Chemistry Project. This paper is a component of the report, focusing on the present-day distribution and trends of tropospheric ozone relevant to climate and global atmospheric chemistry model evaluation. Utilizing the TOAR surface ozone database, several figures present the global distribution and trends of daytime average ozone at 2702 non-urban monitoring sites, highlighting the regions and seasons of the world with the greatest ozone levels. Similarly, ozonesonde and commercial aircraft observations reveal ozone’s distribution throughout the depth of the free troposphere. Long-term surface observations are limited in their global spatial coverage, but data from remote locations indicate that ozone in the 21st century is greater than during the 1970s and 1980s. While some remote sites and many sites in the heavily polluted regions of East Asia show ozone increases since 2000, many others show decreases and there is no clear global pattern for surface ozone changes since 2000. Two new satellite products provide detailed views of ozone in the lower troposphere across East Asia and Europe, revealing the full spatial extent of the spring and summer ozone enhancements across eastern China that cannot be assessed from limited surface observations. Sufficient data are now available (ozonesondes, satellite, aircraft) across the tropics from South America eastwards to the western Pacific Ocean, to indicate a likely tropospheric column ozone increase since the 1990s. The 2014–2016 mean tropospheric ozone burden (TOB) between 60°N–60°S from five satellite products is 300 Tg ± 4%. While this agreement is excellent, the products differ in their quantification of TOB trends and further work is required to reconcile the differences. Satellites can now estimate ozone’s global long-wave radiative effect, but evaluation is difficult due to limited in situ observations where the radiative effect is greatest.

This paper documents a global Bayesian variational inversion of CO 2 surface fluxes during the period 1988–2008. Weekly fluxes are estimated on a 3.75° × 2.5° (longitude‐latitude) grid throughout the 21 years. The assimilated observations include 128 station records from three large data sets of surface CO 2 mixing ratio measurements. A Monte Carlo approach rigorously quantifies the theoretical uncertainty of the inverted fluxes at various space and time scales, which is particularly important for proper interpretation of the inverted fluxes. Fluxes are evaluated indirectly against two independent CO 2 vertical profile data sets constructed from aircraft measurements in the boundary layer and in the free troposphere. The skill of the inversion is evaluated by the improvement brought over a simple benchmark flux estimation based on the observed atmospheric growth rate. Our error analysis indicates that the carbon budget from the inversion should be more accurate than the a priori carbon budget by 20% to 60% for terrestrial fluxes aggregated at the scale of subcontinental regions in the Northern Hemisphere and over a year, but the inversion cannot clearly distinguish between the regional carbon budgets within a continent. On the basis of the independent observations, the inversion is seen to improve the fluxes compared to the benchmark: the atmospheric simulation of CO 2 with the Bayesian inversion method is better by about 1 ppm than the benchmark in the free troposphere, despite possible systematic transport errors. The inversion achieves this improvement by changing the regional fluxes over land at the seasonal and at the interannual time scales.

Abstract The aim of this article is to describe the reference configuration of the convection-permitting numerical weather prediction (NWP) model HARMONIE-AROME, which is used for operational short-range weather forecasts in Denmark, Estonia, Finland, Iceland, Ireland, Lithuania, the Netherlands, Norway, Spain, and Sweden. It is developed, maintained, and validated as part of the shared ALADIN–HIRLAM system by a collaboration of 26 countries in Europe and northern Africa on short-range mesoscale NWP. HARMONIE–AROME is based on the model AROME developed within the ALADIN consortium. Along with the joint modeling framework, AROME was implemented and utilized in both northern and southern European conditions by the above listed countries, and this activity has led to extensive updates to the model’s physical parameterizations. In this paper the authors present the differences in model dynamics and physical parameterizations compared with AROME, as well as important configuration choices of the reference, such as lateral boundary conditions, model levels, horizontal resolution, model time step, as well as topography, physiography, and aerosol databases used. Separate documentation will be provided for the atmospheric and surface data-assimilation algorithms and observation types used, as well as a separate description of the ensemble prediction system based on HARMONIE–AROME, which is called HarmonEPS.

Abstract The European CORDEX (EURO-CORDEX) initiative is a large voluntary effort that seeks to advance regional climate and Earth system science in Europe. As part of the World Climate Research Programme (WCRP) - Coordinated Regional Downscaling Experiment (CORDEX), it shares the broader goals of providing a model evaluation and climate projection framework and improving communication with both the General Circulation Model (GCM) and climate data user communities. EURO-CORDEX oversees the design and coordination of ongoing ensembles of regional climate projections of unprecedented size and resolution (0.11° EUR-11 and 0.44° EUR-44 domains). Additionally, the inclusion of empirical-statistical downscaling allows investigation of much larger multi-model ensembles. These complementary approaches provide a foundation for scientific studies within the climate research community and others. The value of the EURO-CORDEX ensemble is shown via numerous peer-reviewed studies and its use in the development of climate services. Evaluations of the EUR-44 and EUR-11 ensembles also show the benefits of higher resolution. However, significant challenges remain. To further advance scientific understanding, two flagship pilot studies (FPS) were initiated. The first investigates local-regional phenomena at convection-permitting scales over central Europe and the Mediterranean in collaboration with the Med-CORDEX community. The second investigates the impacts of land cover changes on European climate across spatial and temporal scales. Over the coming years, the EURO-CORDEX community looks forward to closer collaboration with other communities, new advances, supporting international initiatives such as the IPCC reports, and continuing to provide the basis for research on regional climate impacts and adaptation in Europe.

Existing descriptions of bi-directional ammonia (NH3) land-atmosphere exchange incorporate temperature and moisture controls, and are beginning to be used in regional chemical transport models. However, such models have typically applied simpler emission factors to upscale the main NH3 emission terms. While this approach has successfully simulated the main spatial patterns on local to global scales, it fails to address the environment- and climate-dependence of emissions. To handle these issues, we outline the basis for a new modelling paradigm where both NH3 emissions and deposition are calculated online according to diurnal, seasonal and spatial differences in meteorology. We show how measurements reveal a strong, but complex pattern of climatic dependence, which is increasingly being characterized using ground-based NH3 monitoring and satellite observations, while advances in process-based modelling are illustrated for agricultural and natural sources, including a global application for seabird colonies. A future architecture for NH3 emission-deposition modelling is proposed that integrates the spatio-temporal interactions, and provides the necessary foundation to assess the consequences of climate change. Based on available measurements, a first empirical estimate suggests that 5°C warming would increase emissions by 42 per cent (28-67%). Together with increased anthropogenic activity, global NH3 emissions may increase from 65 (45-85) Tg N in 2008 to reach 132 (89-179) Tg by 2100.

Abstract. NASA's Orbiting Carbon Observatory-2 (OCO-2) has been measuring carbon dioxide column-averaged dry-air mole fraction, XCO2, in the Earth's atmosphere for over 2 years. In this paper, we describe the comparisons between the first major release of the OCO-2 retrieval algorithm (B7r) and XCO2 from OCO-2's primary ground-based validation network: the Total Carbon Column Observing Network (TCCON). The OCO-2 XCO2 retrievals, after filtering and bias correction, agree well when aggregated around and coincident with TCCON data in nadir, glint, and target observation modes, with absolute median differences less than 0.4 ppm and RMS differences less than 1.5 ppm. After bias correction, residual biases remain. These biases appear to depend on latitude, surface properties, and scattering by aerosols. It is thus crucial to continue measurement comparisons with TCCON to monitor and evaluate the OCO-2 XCO2 data quality throughout its mission.

Abstract. Small changes in the radiation budget at the earth's surface can lead to large climatological responses when persistent over time. With the increasing debate on anthropogenic influences on climatic processes during the 1980s the need for accurate radiometric measurements with higher temporal resolution was identified, and it was determined that the existing measurement networks did not have the resolution or accuracy required to meet this need. In 1988 the WMO therefore proposed the establishment of a new international Baseline Surface Radiation Network (BSRN), which should collect and centrally archive high-quality ground-based radiation measurements in 1 min resolution. BSRN began its work in 1992 with 9 stations; currently (status 2018-01-01), the network comprises 59 stations (delivering data to the archive) and 9 candidates (stations recently accepted into the network with data forthcoming to the archive) distributed over all continents and oceanic environments. The BSRN database is the World Radiation Monitoring Center (WRMC). It is hosted at the Alfred Wegener Institute (AWI) in Bremerhaven, Germany, and now offers more than 10 300 months of data from the years 1992 to 2017. All data are available at https://doi.org/10.1594/PANGAEA.880000 free of charge.

Abstract. Since September 2014, NASA's Orbiting Carbon Observatory-2 (OCO-2) satellite has been taking measurements of reflected solar spectra and using them to infer atmospheric carbon dioxide levels. This work provides details of the OCO-2 retrieval algorithm, versions 7 and 8, used to derive the column-averaged dry air mole fraction of atmospheric CO2 (XCO2) for the roughly 100 000 cloud-free measurements recorded by OCO-2 each day. The algorithm is based on the Atmospheric Carbon Observations from Space (ACOS) algorithm which has been applied to observations from the Greenhouse Gases Observing SATellite (GOSAT) since 2009, with modifications necessary for OCO-2. Because high accuracy, better than 0.25 %, is required in order to accurately infer carbon sources and sinks from XCO2, significant errors and regional-scale biases in the measurements must be minimized. We discuss efforts to filter out poor-quality measurements, and correct the remaining good-quality measurements to minimize regional-scale biases. Updates to the radiance calibration and retrieval forward model in version 8 have improved many aspects of the retrieved data products. The version 8 data appear to have reduced regional-scale biases overall, and demonstrate a clear improvement over the version 7 data. In particular, error variance with respect to TCCON was reduced by 20 % over land and 40 % over ocean between versions 7 and 8, and nadir and glint observations over land are now more consistent. While this paper documents the significant improvements in the ACOS algorithm, it will continue to evolve and improve as the CO2 data record continues to expand.

To date, projections of European crop yields under climate change have been based almost entirely on the outputs of crop-growth models. While this strategy can provide good estimates of the effects of climatic factors, soil conditions and management on crop yield, these models usually do not capture all of the important aspects related to crop management, or the relevant environmental factors. Moreover, crop-simulation studies often have severe limitations with respect to the number of crops covered or the spatial extent. The present study, based on agroclimatic indices, provides a general picture of agroclimatic conditions in western and central Europe (study area lays between 8.5°W–27°E and 37–63.5°N), which allows for a more general assessment of climate-change impacts. The results obtained from the analysis of data from 86 different sites were clustered according to an environmental stratification of Europe. The analysis was carried for the baseline (1971–2000) and future climate conditions (time horizons of 2030, 2050 and with a global temperature increase of 5 °C) based on outputs of three global circulation models. For many environmental zones, there were clear signs of deteriorating agroclimatic condition in terms of increased drought stress and shortening of the active growing season, which in some regions become increasingly squeezed between a cold winter and a hot summer. For most zones the projections show a marked need for adaptive measures to either increase soil water availability or drought resistance of crops. This study concludes that rainfed agriculture is likely to face more climate-related risks, although the analyzed agroclimatic indicators will probably remain at a level that should permit rainfed production. However, results suggests that there is a risk of increasing number of extremely unfavorable years in many climate zones, which might result in higher interannual yield variability and constitute a challenge for proper crop management.

We analyze century‐long daily temperature and precipitation records for stations in Europe west of 60°E. A set of climatic indices derived from the daily series, mainly focusing on extremes, is defined. Linear trends in these indices are assessed over the period 1901–2000. Average trends, for 75 stations mostly representing Europe west of 20°E, show a warming for all temperature indices. Winter has, on average, warmed more (∼1.0°C/100 yr) than summer (∼0.8°C), both for daily maximum (TX) and minimum (TN) temperatures. Overall, the warming of TX in winter was stronger in the warm tail than in the cold tail (1.6 and 1.5°C for 98th and 95th, but ∼1.0°C for 2nd, 5th and 10th percentiles). There are, however, large regional differences in temperature trend patterns. For summer, there is a tendency for stronger warming, both for TX and TN, in the warm than in the cold tail only in parts of central Europe. Winter precipitation totals, averaged over 121 European stations north of 40°N, have increased significantly by ∼12% per 100 years. Trends in 90th, 95th and 98th percentiles of daily winter precipitation have been similar. No overall long‐term trend occurred in summer precipitation totals, but there is an overall weak (statistically insignificant and regionally dependent) tendency for summer precipitation to have become slightly more intense but less common. Data inhomogeneities and relative sparseness of station density in many parts of Europe preclude more robust conclusions. It is of importance that new methods are developed for homogenizing daily data.

The Galician coast (NW Spain) is a region that is strongly influenced by the presence of low pressure systems in the mid-Atlantic Ocean and the periodic passage of storms that give rise to severe sea states. Since its wave climate is one of the most energetic in Europe, the objectives of this paper were twofold. The first objective was to characterize the most extreme wave height events in Galicia over the wintertime of a two-year period (2015–2016) by using reliable high-frequency radar wave parameters in concert with predictions from a regional wave (WAV) forecasting system running operationally in the Iberia-Biscay-Ireland (IBI) area, denominatedIBI-WAV. The second objective was to showcase the application of satellite wave altimetry (in particular, remote-sensed three-hourly wave height estimations) for the daily skill assessment of the IBI-WAV model product. Special attention was focused on monitoring Ophelia—one of the major hurricanes on record in the easternmost Atlantic—during its 3-day track over Ireland and the UK (15–17 October 2017). Overall, the results reveal the significant accuracy of IBI-WAV forecasts and prove that a combined observational and modeling approach can provide a comprehensive characterization of severe wave conditions in coastal areas and shows the benefits from the complementary nature of both systems.

Abstract. Absorption photometers for real time application have been available since the 1980s, but the use of filter-based instruments to derive information on aerosol properties (absorption coefficient and black carbon, BC) is still a matter of debate. Several workshops have been conducted to investigate the performance of individual instruments over the intervening years. Two workshops with large sets of aerosol absorption photometers were conducted in 2005 and 2007. The data from these instruments were corrected using existing methods before further analysis. The inter-comparison shows a large variation between the responses to absorbing aerosol particles for different types of instruments. The unit to unit variability between instruments can be up to 30% for Particle Soot Absorption Photometers (PSAPs) and Aethalometers. Multi Angle Absorption Photometers (MAAPs) showed a variability of less than 5%. Reasons for the high variability were identified to be variations in sample flow and spot size. It was observed that different flow rates influence system performance with respect to response to absorption and instrumental noise. Measurements with non absorbing particles showed that the current corrections of a cross sensitivity to particle scattering are not sufficient. Remaining cross sensitivities were found to be a function of the total particle load on the filter. The large variation between the response to absorbing aerosol particles for different types of instruments indicates that current correction functions for absorption photometers are not adequate.

In support of the first Tropospheric Ozone Assessment Report (TOAR) a relational database of global surface ozone observations has been developed and populated with hourly measurement data and enhanced metadata. A comprehensive suite of ozone data products including standard statistics, health and vegetation impact metrics, and trend information, are made available through a common data portal and a web interface. These data form the basis of the TOAR analyses focusing on human health, vegetation, and climate relevant ozone issues, which are part of this special feature. Cooperation among many data centers and individual researchers worldwide made it possible to build the world’s largest collection of in-situ hourly surface ozone data covering the period from 1970 to 2015. By combining the data from almost 10,000 measurement sites around the world with global metadata information, new analyses of surface ozone have become possible, such as the first globally consistent characterisations of measurement sites as either urban or rural/remote. Exploitation of these global metadata allows for new insights into the global distribution, and seasonal and long-term changes of tropospheric ozone and they enable TOAR to perform the first, globally consistent analysis of present-day ozone concentrations and recent ozone changes with relevance to health, agriculture, and climate. Considerable effort was made to harmonize and synthesize data formats and metadata information from various networks and individual data submissions. Extensive quality control was applied to identify questionable and erroneous data, including changes in apparent instrument offsets or calibrations. Such data were excluded from TOAR data products. Limitations of a posteriori data quality assurance are discussed. As a result of the work presented here, global coverage of surface ozone data for scientific analysis has been significantly extended. Yet, large gaps remain in the surface observation network both in terms of regions without monitoring, and in terms of regions that have monitoring programs but no public access to the data archive. Therefore future improvements to the database will require not only improved data harmonization, but also expanded data sharing and increased monitoring in data-sparse regions.