12 results
Search Results
Now showing 1 - 10 of 12
- ItemTarget categorization of aerosol and clouds by continuous multiwavelength-polarization lidar measurements(Katlenburg-Lindau : Copernicus, 2017) Baars, Holger; Seifert, Patric; Engelmann, Ronny; Wandinger, UllaAbsolute calibrated signals at 532 and 1064 nm and the depolarization ratio from a multiwavelength lidar are used to categorize primary aerosol but also clouds in high temporal and spatial resolution. Automatically derived particle backscatter coefficient profiles in low temporal resolution (30 min) are applied to calibrate the lidar signals. From these calibrated lidar signals, new atmospheric parameters in temporally high resolution (quasi-particle-backscatter coefficients) are derived. By using thresholds obtained from multiyear, multisite EARLINET (European Aerosol Research Lidar Network) measurements, four aerosol classes (small; large, spherical; large, non-spherical; mixed, partly nonspherical) and several cloud classes (liquid, ice) are defined. Thus, particles are classified by their physical features (shape and size) instead of by source. The methodology is applied to 2 months of continuous observations (24 h a day, 7 days a week) with the multiwavelength-Raman-polarization lidar PollyXT during the High-Definition Clouds and Precipitation for advancing Climate Prediction (HD(CP)2) Observational Prototype Experiment (HOPE) in spring 2013. Cloudnet equipment was operated continuously directly next to the lidar and is used for comparison. By discussing three 24 h case studies, it is shown that the aerosol discrimination is very feasible and informative and gives a good complement to the Cloudnet target categorization. Performing the categorization for the 2-month data set of the entire HOPE campaign, almost 1 million pixel (5 min×30 m) could be analysed with the newly developed tool. We find that the majority of the aerosol trapped in the planetary boundary layer (PBL) was composed of small particles as expected for a heavily populated and industrialized area. Large, spherical aerosol was observed mostly at the top of the PBL and close to the identified cloud bases, indicating the importance of hygroscopic growth of the particles at high relative humidity. Interestingly, it is found that on several days non-spherical particles were dispersed from the ground into the atmosphere.
- ItemHelicopter-borne observations of the continental background aerosol in combination with remote sensing and ground-based measurements(Katlenburg-Lindau : EGU, 2018) Düsing, Sebastian; Wehner, Birgit; Seifert, Patric; Ansmann, Albert; Baars, Holger; Ditas, Florian; Henning, Silvia; Ma, Nan; Poulain, Laurent; Siebert, Holger; Wiedensohler, Alfred; Macke, AndreasThis paper examines the representativeness of ground-based in situ measurements for the planetary boundary layer (PBL) and conducts a closure study between airborne in situ and ground-based lidar measurements up to an altitude of 2300 m. The related measurements were carried out in a field campaign within the framework of the High-Definition Clouds and Precipitation for Advancing Climate Prediction (HD(CP)2) Observational Prototype Experiment (HOPE) in September 2013 in a rural background area of central Europe. The helicopter-borne probe ACTOS (Airborne Cloud and Turbulence Observation System) provided measurements of the aerosol particle number size distribution (PNSD), the aerosol particle number concentration (PNC), the number concentration of cloud condensation nuclei (CCN-NC), and meteorological atmospheric parameters (e.g., temperature and relative humidity). These measurements were supported by the ground-based 3+2 wavelength polarization lidar system PollyXT, which provided profiles of the particle backscatter coefficient (σbsc) for three wavelengths (355, 532, and 1064 nm). Particle extinction coefficient (σext) profiles were obtained by using a fixed backscatter-to-extinction ratio (also lidar ratio, LR). A new approach was used to determine profiles of CCN-NC for continental aerosol. The results of this new approach were consistent with the airborne in situ measurements within the uncertainties. In terms of representativeness, the PNSD measurements on the ground showed a good agreement with the measurements provided with ACTOS for lower altitudes. The ground-based measurements of PNC and CCN-NC are representative of the PBL when the PBL is well mixed. Locally isolated new particle formation events on the ground or at the top of the PBL led to vertical variability in the cases presented here and ground-based measurements are not entirely representative of the PBL. Based on Mie theory (Mie, 1908), optical aerosol properties under ambient conditions for different altitudes were determined using the airborne in situ measurements and were compared with the lidar measurements. The investigation of the optical properties shows that on average the airborne-based particle light backscatter coefficient is 50.1 % smaller for 1064 nm, 27.4 % smaller for 532 nm, and 29.5 % smaller for 355 nm than the measurements of the lidar system. These results are quite promising, since in situ measurement-based Mie calculations of the particle light backscattering are scarce and the modeling is quite challenging. In contrast, for the particle light extinction coefficient we found a good agreement. The airborne-based particle light extinction coefficient was just 8.2 % larger for 532 nm and 3 % smaller for 355 nm, for an assumed LR of 55 sr. The particle light extinction coefficient for 1064 nm was derived with a LR of 30 sr. For this wavelength, the airborne-based particle light extinction coefficient is 5.2 % smaller than the lidar measurements. For the first time, the lidar ratio of 30 sr for 1064 nm was determined on the basis of in situ measurements and the LR of 55 sr for 355 and 532 nm wavelength was reproduced for European continental aerosol on the basis of this comparison. Lidar observations and the in situ based aerosol optical properties agree within the uncertainties. However, our observations indicate that a determination of the PNSD for a large size range is important for a reliable modeling of aerosol particle backscattering.
- ItemThe HD(CP)2 Observational Prototype Experiment (HOPE) - An overview(Katlenburg-Lindau : EGU, 2017) Macke, Andreas; Seifert, Patric; Baars, Holger; Barthlott, Christian; Beekmans, Christoph; Behrendt, Andreas; Bohn, Birger; Brueck, Matthias; Bühl, Johannes; Crewell, Susanne; Damian, Thomas; Deneke, Hartwig; Düsing, Sebastian; Foth, Andreas; Di Girolamo, Paolo; Hammann, Eva; Heinze, Rieke; Hirsikko, Anne; Kalisch, John; Kalthoff, Norbert; Kinne, Stefan; Kohler, Martin; Löhnert, Ulrich; Madhavan, Bomidi Lakshmi; Maurer, Vera; Muppa, Shravan Kumar; Schween, Jan; Serikov, Ilya; Siebert, Holger; Simmer, Clemens; Späth, Florian; Steinke, Sandra; Träumner, Katja; Trömel, Silke; Wehner, Birgit; Wieser, Andreas; Wulfmeyer, Volker; Xie, XinxinThe HD(CP)2 Observational Prototype Experiment (HOPE) was performed as a major 2-month field experiment in Jülich, Germany, in April and May 2013, followed by a smaller campaign in Melpitz, Germany, in September 2013. HOPE has been designed to provide an observational dataset for a critical evaluation of the new German community atmospheric icosahedral non-hydrostatic (ICON) model at the scale of the model simulations and further to provide information on land-surface-atmospheric boundary layer exchange, cloud and precipitation processes, as well as sub-grid variability and microphysical properties that are subject to parameterizations. HOPE focuses on the onset of clouds and precipitation in the convective atmospheric boundary layer. This paper summarizes the instrument set-ups, the intensive observation periods, and example results from both campaigns.
HOPE-Jülich instrumentation included a radio sounding station, 4 Doppler lidars, 4 Raman lidars (3 of them provide temperature, 3 of them water vapour, and all of them particle backscatter data), 1 water vapour differential absorption lidar, 3 cloud radars, 5 microwave radiometers, 3 rain radars, 6 sky imagers, 99 pyranometers, and 5 sun photometers operated at different sites, some of them in synergy. The HOPE-Melpitz campaign combined ground-based remote sensing of aerosols and clouds with helicopter- and balloon-based in situ observations in the atmospheric column and at the surface.
HOPE provided an unprecedented collection of atmospheric dynamical, thermodynamical, and micro- and macrophysical properties of aerosols, clouds, and precipitation with high spatial and temporal resolution within a cube of approximately 10 × 10 × 10km3. HOPE data will significantly contribute to our understanding of boundary layer dynamics and the formation of clouds and precipitation. The datasets have been made available through a dedicated data portal.
First applications of HOPE data for model evaluation have shown a general agreement between observed and modelled boundary layer height, turbulence characteristics, and cloud coverage, but they also point to significant differences that deserve further investigations from both the observational and the modelling perspective. - ItemProfiling of Saharan dust from the Caribbean to western Africa - Part 2: Shipborne lidar measurements versus forecasts(Katlenburg-Lindau : EGU, 2017) Ansmann, Albert; Rittmeister, Franziska; Engelmann, Ronny; Basart, Sara; Jorba, Oriol; Spyrou, Christos; Remy, Samuel; Skupin, Annett; Baars, Holger; Seifert, Patric; Senf, Fabian; Kanitz, ThomasA unique 4-week ship cruise from Guadeloupe to Cabo Verde in April-May 2013 see part 1, Rittmeister et al. (2017) is used for an in-depth comparison of dust profiles observed with a polarization/Raman lidar aboard the German research vessel Meteor over the remote tropical Atlantic and respective dust forecasts of a regional (SKIRON) and two global atmospheric (dust) transport models (NMMB/BSC-Dust, MACC/CAMS). New options of model-observation comparisons are presented. We analyze how well the modeled fine dust (submicrometer particles) and coarse dust contributions to light extinction and mass concentration match respective lidar observations, and to what extent models, adjusted to aerosol optical thickness observations, are able to reproduce the observed layering and mixing of dust and non-dust (mostly marine) aerosol components over the remote tropical Atlantic. Based on the coherent set of dust profiles at well-defined distances from Africa (without any disturbance by anthropogenic aerosol sources over the ocean), we investigate how accurately the models handle dust removal at distances of 1500g km to more than 5000g km west of the Saharan dust source regions. It was found that (a) dust predictions are of acceptable quality for the first several days after dust emission up to 2000g km west of the African continent, (b) the removal of dust from the atmosphere is too strong for large transport paths in the global models, and (c) the simulated fine-to-coarse dust ratio (in terms of mass concentration and light extinction) is too high in the models compared to the observations. This deviation occurs initially close to the dust sources and then increases with distance from Africa and thus points to an overestimation of fine dust emission in the models.
- ItemOptical and geometrical properties of cirrus clouds in Amazonia derived from 1 year of ground-based lidar measurements(Katlenburg-Lindau : EGU, 2017) Gouveia, Diego A.; Barja, Boris; Barbosa, Henrique M. J.; Seifert, Patric; Baars, Holger; Pauliquevis, Theotonio; Artaxo, PauloCirrus clouds cover a large fraction of tropical latitudes and play an important role in Earth's radiation budget. Their optical properties, altitude, vertical and horizontal coverage control their radiative forcing, and hence detailed cirrus measurements at different geographical locations are of utmost importance. Studies reporting cirrus properties over tropical rain forests like the Amazon, however, are scarce. Studies with satellite profilers do not give information on the diurnal cycle, and the satellite imagers do not report on the cloud vertical structure. At the same time, ground-based lidar studies are restricted to a few case studies. In this paper, we derive the first comprehensive statistics of optical and geometrical properties of upper-tropospheric cirrus clouds in Amazonia. We used 1 year (July 2011 to June 2012) of ground-based lidar atmospheric observations north of Manaus, Brazil. This dataset was processed by an automatic cloud detection and optical properties retrieval algorithm. Upper-tropospheric cirrus clouds were observed more frequently than reported previously for tropical regions. The frequency of occurrence was found to be as high as 88 % during the wet season and not lower than 50 % during the dry season. The diurnal cycle shows a minimum around local noon and maximum during late afternoon, associated with the diurnal cycle of precipitation. The mean values of cirrus cloud top and base heights, cloud thickness, and cloud optical depth were 14.3 ± 1.9 (SD) km, 12.9 ± 2.2 km, 1.4 ± 1.1 km, and 0.25 ± 0.46, respectively. Cirrus clouds were found at temperatures down to ĝ'90 °C. Frequently cirrus were observed within the tropical tropopause layer (TTL), which are likely associated to slow mesoscale uplifting or to the remnants of overshooting convection. The vertical distribution was not uniform, and thin and subvisible cirrus occurred more frequently closer to the tropopause. The mean lidar ratio was 23.3 ± 8.0 sr. However, for subvisible cirrus clouds a bimodal distribution with a secondary peak at about 44 sr was found suggesting a mixed composition. A dependence of the lidar ratio with cloud temperature (altitude) was not found, indicating that the clouds are vertically well mixed. The frequency of occurrence of cirrus clouds classified as subvisible (τ < 0. 03) were 41.6 %, whilst 37.8 % were thin cirrus (0. 03;lt; τ < 0. 3) and 20.5 % opaque cirrus (τ > 0. 3). Hence, in central Amazonia not only a high frequency of cirrus clouds occurs, but also a large fraction of subvisible cirrus clouds. This high frequency of subvisible cirrus clouds may contaminate aerosol optical depth measured by sun photometers and satellite sensors to an unknown extent.
- ItemExtreme levels of Canadian wildfire smoke in the stratosphere over central Europe on 21-22 August 2017(Katlenburg-Lindau : EGU, 2018) Ansmann, Albert; Baars, Holger; Chudnovsky, Alexandra; Mattis, Ina; Veselovskii, Igor; Haarig, Moritz; Seifert, Patric; Engelmann, Ronny; Wandinger, UllaLight extinction coefficients of 500 Mm1, about 20 times higher than after the Pinatubo volcanic eruptions in 1991, were observed by European Aerosol Research Lidar Network (EARLINET) lidars in the stratosphere over central Europe on 21-22 August 2017. Pronounced smoke layers with a 1-2 km vertical extent were found 2-5 km above the local tropopause. Optically dense layers of Canadian wildfire smoke reached central Europe 10 days after their injection into the upper troposphere and lower stratosphere which was caused by rather strong pyrocumulonimbus activity over western Canada. The smoke-related aerosol optical thickness (AOT) identified by lidar was close to 1.0 at 532 nm over Leipzig during the noon hours on 22 August 2017. Smoke particles were found throughout the free troposphere (AOT of 0.3) and in the pronounced 2 km thick stratospheric smoke layer at an altitude of 14-16 km (AOT of 0.6). The lidar observations indicated peak mass concentrations of 70-100 μgm-3 in the stratosphere. In addition to the lidar profiles, we analyzed Moderate Resolution Imaging Spectroradiometer (MODIS) fire radiative power (FRP) over Canada, and the distribution of MODIS AOT and Ozone Monitoring Instrument (OMI) aerosol index across the North Atlantic. These instruments showed a similar pattern and a clear link between the western Canadian fires and the aerosol load over Europe. In this paper, we also present Aerosol Robotic Network (AERONET) sun photometer observations, compare photometer and lidar-derived AOT, and discuss an obvious bias (the smoke AOT is too low) in the photometer observations. Finally, we compare the strength of this recordbreaking smoke event (in terms of the particle extinction coefficient and AOT) with major and moderate volcanic events observed over the northern midlatitudes.
- ItemIce-nucleating particle versus ice crystal number concentrationin altocumulus and cirrus layers embedded in Saharan dust: a closure study(Katlenburg-Lindau : EGU, 2019) Ansmann, Albert; Mamouri, Rodanthi-Elisavet; Bühl, Johannes; Seifert, Patric; Engelmann, Ronny; Hofer, Julian; Nisantzi, Argyro; Atkinson, James D.; Kanji, Zamin A.; Sierau, Berko; Vrekoussis, Mihalis; Sciare, JeanFor the first time, a closure study of the relationship between the ice-nucleating particle concentration (INP; INPC) and ice crystal number concentration (ICNC) in altocumulus and cirrus layers, solely based on groundbased active remote sensing, is presented. Such aerosol- cloud closure experiments are required (a) to better understand aerosol-cloud interaction in the case of mixed-phase clouds, (b) to explore to what extent heterogeneous ice nucleation can contribute to cirrus formation, which is usually controlled by homogeneous freezing, and (c) to check the usefulness of available INPC parameterization schemes, applied to lidar profiles of aerosol optical and microphysical properties up to the tropopause level. The INPC-ICNC closure studies were conducted in Cyprus (Limassol and Nicosia) during a 6-week field campaign in March-April 2015 and during the 17-month CyCARE (Cyprus Clouds Aerosol and Rain Experiment) campaign. The focus was on altocumulus and cirrus layers which developed in pronounced Saharan dust layers at heights from 5 to 11 km. As a highlight, a long-lasting cirrus event was studied which was linked to the development of a very strong dust-infused baroclinic storm (DIBS) over Algeria. The DIBS was associated with strong convective cloud development and lifted large amounts of Saharan dust into the upper troposphere, where the dust influenced the evolution of an unusually large anvil cirrus shield and the subsequent transformation into an cirrus uncinus cloud system extending from the eastern Mediterranean to central Asia, and thus over more than 3500 km. Cloud top temperatures of the three discussed closure study cases ranged from - 20 to -57 °C. The INPC was estimated from polarization/Raman lidar observations in combination with published INPC parameterization schemes, whereas the ICNC was retrieved from combined Doppler lidar, aerosol lidar, and cloud radar observations of the terminal velocity of falling ice crystals, radar reflectivity, and lidar backscatter in combination with the modeling of backscattering at the 532 and 8.5 mm wavelengths. A good-to-acceptable agreement between INPC (observed before and after the occurrence of the cloud layer under investigation) and ICNC values was found in the discussed three proof-of-concept closure experiments. In these case studies, INPC and ICNC values matched within an order of magnitude (i.e., within the uncertainty ranges of the INPC and ICNC estimates), and they ranged from 0.1 to 10 L-1 in the altocumulus layers and 1 to 50 L-1 in the cirrus layers observed between 8 and 11 km height. The successful closure experiments corroborate the important role of heterogeneous ice nucleation in atmospheric ice formation processes when mineral dust is present. The observed longlasting cirrus event could be fully explained by the presence of dust, i.e., without the need for homogeneous ice nucleation processes. © 2019 Author(s).
- ItemThe Arctic Cloud Puzzle: Using ACLOUD/PASCAL Multiplatform Observations to Unravel the Role of Clouds and Aerosol Particles in Arctic Amplification(Boston, Mass. : ASM, 2019) Wendisch, Manfred; Macke, Andreas; Ehrlich, André; Lüpkes, Christof; Mech, Mario; Chechin, Dmitry; Dethloff, Klaus; Velasco, Carola Barrientos; Bozem, Heiko; Brückner, Marlen; Clemen, Hans-Christian; Crewell, Susanne; Donth, Tobias; Dupuy, Regis; Ebell, Kerstin; Egerer, Ulrike; Engelmann, Ronny; Engler, Christa; Eppers, Oliver; Gehrmann, Martin; Gong, Xianda; Gottschalk, Matthias; Gourbeyre, Christophe; Griesche, Hannes; Hartmann, Jörg; Hartmann, Markus; Heinold, Bernd; Herber, Andreas; Herrmann, Hartmut; Heygster, Georg; Hoor, Peter; Jafariserajehlou, Soheila; Jäkel, Evelyn; Järvinen, Emma; Jourdan, Olivier; Kästner, Udo; Kecorius, Simonas; Knudsen, Erlend M.; Köllner, Franziska; Kretzschmar, Jan; Lelli, Luca; Leroy, Delphine; Maturilli, Marion; Mei, Linlu; Mertes, Stephan; Mioche, Guillaume; Neuber, Roland; Nicolaus, Marcel; Nomokonova, Tatiana; Notholt, Justus; Palm, Mathias; van Pinxteren, Manuela; Quaas, Johannes; Richter, Philipp; Ruiz-Donoso, Elena; Schäfer, Michael; Schmieder, Katja; Schnaiter, Martin; Schneider, Johannes; Schwarzenböck, Alfons; Seifert, Patric; Shupe, Matthew D.; Siebert, Holger; Spreen, Gunnar; Stapf, Johannes; Stratmann, Frank; Vogl, Teresa; Welti, André; Wex, Heike; Wiedensohler, Alfred; Zanatta, Marco; Zeppenfeld, SebastianClouds play an important role in Arctic amplification. This term represents the recently observed enhanced warming of the Arctic relative to the global increase of near-surface air temperature. However, there are still important knowledge gaps regarding the interplay between Arctic clouds and aerosol particles, and surface properties, as well as turbulent and radiative fluxes that inhibit accurate model simulations of clouds in the Arctic climate system. In an attempt to resolve this so-called Arctic cloud puzzle, two comprehensive and closely coordinated field studies were conducted: the Arctic Cloud Observations Using Airborne Measurements during Polar Day (ACLOUD) aircraft campaign and the Physical Feedbacks of Arctic Boundary Layer, Sea Ice, Cloud and Aerosol (PASCAL) ice breaker expedition. Both observational studies were performed in the framework of the German Arctic Amplification: Climate Relevant Atmospheric and Surface Processes, and Feedback Mechanisms (AC) project. They took place in the vicinity of Svalbard, Norway, in May and June 2017. ACLOUD and PASCAL explored four pieces of the Arctic cloud puzzle: cloud properties, aerosol impact on clouds, atmospheric radiation, and turbulent dynamical processes. The two instrumented Polar 5 and Polar 6 aircraft; the icebreaker Research Vessel (R/V) Polarstern; an ice floe camp including an instrumented tethered balloon; and the permanent ground-based measurement station at Ny-Ålesund, Svalbard, were employed to observe Arctic low- and mid-level mixed-phase clouds and to investigate related atmospheric and surface processes. The Polar 5 aircraft served as a remote sensing observatory examining the clouds from above by downward-looking sensors; the Polar 6 aircraft operated as a flying in situ measurement laboratory sampling inside and below the clouds. Most of the collocated Polar 5/6 flights were conducted either above the R/V Polarstern or over the Ny-Ålesund station, both of which monitored the clouds from below using similar but upward-looking remote sensing techniques as the Polar 5 aircraft. Several of the flights were carried out underneath collocated satellite tracks. The paper motivates the scientific objectives of the ACLOUD/PASCAL observations and describes the measured quantities, retrieved parameters, and the applied complementary instrumentation. Furthermore, it discusses selected measurement results and poses critical research questions to be answered in future papers analyzing the data from the two field campaigns.
- ItemIce crystal number concentration from lidar, cloud radar and radar wind profiler measurements(Katlenburg-Lindau : Copernicus, 2019) Bühl, Johannes; Seifert, Patric; Radenz, Martin; Ansmann, AlbertA new method for the retrieval of ice crystal number concentration (ICNC) from combined active remote-sensing measurements of Raman lidar, cloud radar and radar wind profiler is presented. We exploit – for the first time – measurements of terminal fall velocity together with the radar reflectivity factor and/or the lidar-derived particle extinction coefficient in clouds for retrieving the number concentration of pristine ice particles with presumed particle shapes. A lookup table approach for the retrieval of the properties of the particle size distribution from observed parameters is presented. Analysis of methodological uncertainties and error propagation is performed, which shows that a retrieval of ice particle number concentration based on terminal fall velocity is possible within 1 order of magnitude. Comparison between a retrieval of the number concentration based on terminal fall velocity on the one hand and lidar and cloud radar on the other shows agreement within the uncertainties of the retrieval.
- ItemVertical aerosol distribution in the southern hemispheric midlatitudes as observed with lidar in Punta Arenas, Chile (53.2° and 70.9° W), during ALPACA(Katlenburg-Lindau : EGU, 2019) Foth, Andreas; Kanitz, Thomas; Engelmann, Ronny; Baars, Holger; Radenz, Martin; Seifert, Patric; Barja, Boris; Fromm, Michael; Kalesse, Heike; Ansmann, AlbertWithin this publication, lidar observations of the vertical aerosol distribution above Punta Arenas, Chile (53.2 S and 70.9 W), which have been performed with the Raman lidar PollyXT from December 2009 to April 2010, are presented. Pristine marine aerosol conditions related to the prevailing westerly circulation dominated the measurements. Lofted aerosol layers could only be observed eight times during the whole measurement period. Two case studies are presented showing long-range transport of smoke from biomass burning in Australia and regionally transported dust from the Patagonian Desert, respectively. The aerosol sources are identified by trajectory analyses with the Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) and FLEXible PARTicle dispersion model (FLEXPART). However, seven of the eight analysed cases with lofted layers show an aerosol optical thickness of less than 0.05. From the lidar observations, a mean planetary boundary layer (PBL) top height of 1150 350m was determined. An analysis of particle backscatter coefficients confirms that the majority of the aerosol is attributed to the PBL, while the free troposphere is characterized by a very low background aerosol concentration. The ground-based lidar observations at 532 and 1064 nm are supplemented by the Aerosol Robotic Network (AERONET) Sun photometers and the space-borne Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) aboard the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO). The averaged aerosol optical thickness (AOT) determined by CALIOP was 0:02 0:01 in Punta Arenas from 2009 to 2010. © Author(s) 2019. This work is distributed under the Creative Commons Attribution 4.0 License.