Refine
Year of publication
Document Type
- Article (892) (remove)
Is part of the Bibliography
- yes (892) (remove)
Keywords
- climate change (13)
- Germany (7)
- Climate change (5)
- adaptation (5)
- Europe (4)
- Turkey (4)
- floods (4)
- permafrost (4)
- validation (4)
- Adaptation (3)
Institute
- Institut für Umweltwissenschaften und Geographie (892) (remove)
We present a chronology framework named LegacyAge 1.0 containing harmonized chronologies for 2831 pollen records (downloaded from the Neotoma Paleoecology Database and the supplementary Asian datasets) together with their age control points and metadata in machine-readable data formats.
All chronologies use the Bayesian framework implemented in Bacon version 2.5.3. Optimal parameter settings of priors (accumulation.shape, memory.strength, memory.mean, accumulation.rate, and thickness) were identified based on information in the original publication or iteratively after preliminary model inspection.
The most common control points for the chronologies are radiocarbon dates (86.1 %), calibrated by the latest calibration curves (IntCal20 and SHCal20 for the terrestrial radiocarbon dates in the Northern Hemisphere and Southern Hemisphere and Marine20 for marine materials).
The original publications were consulted when dealing with outliers and inconsistencies. Several major challenges when setting up the chronologies included the waterline issue (18.8% of records), reservoir effect (4.9 %), and sediment deposition discontinuity (4.4 %).
Finally, we numerically compare the LegacyAge 1.0 chronologies to those published in the original publications and show that the reliability of the chronologies of 95.4% of records could be improved according to our assessment.
Our chronology framework and revised chronologies provide the opportunity to make use of the ages and age uncertainties in synthesis studies of, for example, pollen-based vegetation and climate change.
The LegacyAge 1.0 dataset, including metadata, datings, harmonized chronologies, and R code used, is openaccess and available at PANGAEA (https://doi.org/10.1594/PANGAEA.933132; Li et al., 2021) and Zenodo (https://doi.org/10.5281/zenodo.5815192; Li et al., 2022), respectively.
Drought and the availability of mineable phosphorus minerals used for fertilization are two of the important issues agriculture is facing in the future. High phosphorus availability in soils is necessary to maintain high agricultural yields. Drought is one of the major threats for terrestrial ecosystem performance and crop production in future. Among the measures proposed to cope with the upcoming challenges of intensifying drought stress and to decrease the need for phosphorus fertilizer application is the fertilization with silica (Si). Here we tested the importance of soil Si fertilization on wheat phosphorus concentration as well as wheat performance during drought at the field scale. Our data clearly showed a higher soil moisture for the Si fertilized plots. This higher soil moisture contributes to a better plant performance in terms of higher photosynthetic activity and later senescence as well as faster stomata responses ensuring higher productivity during drought periods. The plant phosphorus concentration was also higher in Si fertilized compared to control plots. Overall, Si fertilization or management of the soil Si pools seem to be a promising tool to maintain crop production under predicted longer and more serve droughts in the future and reduces phosphorus fertilizer requirements.
Earthquake site responses or site effects are the modifications of surface geology to seismic waves. How well can we predict the site effects (average over many earthquakes) at individual sites so far? To address this question, we tested and compared the effectiveness of different estimation techniques in predicting the outcrop Fourier site responses separated using the general inversion technique (GIT) from recordings. Techniques being evaluated are (a) the empirical correction to the horizontal-to-vertical spectral ratio of earthquakes (c-HVSR), (b) one-dimensional ground response analysis (GRA), and (c) the square-root-impedance (SRI) method (also called the quarter-wavelength approach). Our results show that c-HVSR can capture significantly more site-specific features in site responses than both GRA and SRI in the aggregate, especially at relatively high frequencies. c-HVSR achieves a "good match" in spectral shape at similar to 80%-90% of 145 testing sites, whereas GRA and SRI fail at most sites. GRA and SRI results have a high level of parametric and/or modeling errors which can be constrained, to some extent, by collecting on-site recordings.
The 2020s are an essential decade for achieving the 2030 Agenda and its Sustainable Development Goals (SDGs). For this, SDG research needs to provide evidence that can be translated into concrete actions. However, studies use different SDG data, resulting in incomparable findings. Researchers primarily use SDG databases provided by the United Nations (UN), the World Bank Group (WBG), and the Bertelsmann Stiftung & Sustainable Development Solutions Network (BE-SDSN). We compile these databases into one unified SDG database and examine the effects of the data selection on our understanding of SDG interactions. Among the databases, we observed more different than similar SDG interactions. Differences in synergies and trade-offs mainly occur for SDGs that are environmentally oriented. Due to the increased data availability, the unified SDG database offers a more nuanced and reliable view of SDG interactions. Thus, the SDG data selection may lead to diverse findings, fostering actions that might neglect or exacerbate trade-offs.
The fluxes of water and solutes in the subsurface compartment of the Critical Zone are temporally dynamic and it is unclear how this impacts microbial mediated nutrient cycling in the spatially heterogeneous subsurface. To investigate this, we undertook numerical modeling, simulating the transport in a wide range of spatially heterogeneous domains, and the biogeochemical transformation of organic carbon and nitrogen compounds using a complex microbial community with four (4) distinct functional groups, in water saturated subsurface compartments. We performed a comprehensive uncertainty analysis accounting for varying residence times and spatial heterogeneity. While the aggregated removal of chemical species in the domains over the entire simulation period was approximately the same as that in steady state conditions, the sub-scale temporal variation of microbial biomass and chemical discharge from a domain depended strongly on the interplay of spatial heterogeneity and temporal dynamics of the forcing. We showed that the travel time and the Damkohler number (Da) can be used to predict the temporally varying chemical discharge from a spatially heterogeneous domain. In homogeneous domains, chemical discharge in temporally dynamic conditions could be double of that in the steady state conditions while microbial biomass varied up to 75% of that in steady state conditions. In heterogeneous domains, the interquartile range of uncertainty in chemical discharge in reaction dominated systems (log(10)Da > 0) was double of that in steady state conditions. However, high heterogeneous domains resulted in outliers where chemical discharge could be as high as 10-20 times of that in steady state conditions in high flow periods. And in transport dominated systems (log(10)Da < 0), the chemical discharge could be half of that in steady state conditions in unusually low flow conditions. In conclusion, ignoring spatio-temporal heterogeneities in a numerical modeling approach may exacerbate inaccurate estimation of nutrient export and microbial biomass. The results are relevant to long-term field monitoring studies, and for homogeneous soil column-scale experiments investigating the role of temporal dynamics on microbial redox dynamics.
Despite the amount of research focussed on the Alpine orogen, different hypotheses still exist regarding varying spatial seismicity distribution patterns throughout the region. Previous measurement-constrained regional 3D models of lithospheric density distribution and thermal field facilitate the generation of a data-based rheological model of the region.
In this study, we compute the long-term lithospheric strength and compare its spatial variation to observed seismicity patterns. We demonstrate how strength maxima within the crust (similar to 1 GPa) and upper mantle (> 2 GPa) occur at temperatures characteristic of the onset of crystal plasticity in those rocks (crust: 200-400 degrees C; mantle: similar to 600 degrees C), with almost all seismicity occurring in these regions. Correlation in the northern and southern forelands between crustal and lithospheric strengths and seismicity show different patterns of event distribution, reflecting their different tectonic settings. Seismicity in the plate boundary setting of the southern foreland corresponds to the integrated lithospheric strength, occurring mainly in the weaker domains surrounding the strong Adriatic plate. In the intraplate setting of the northern foreland, seismicity correlates to modelled crustal strength, and it mainly occurs in the weaker and warmer crust beneath the Upper Rhine Graben.
We, therefore, suggest that seismicity in the upper crust is linked to weak crustal domains, which are more prone to localise deformation promoting failure and, depending on the local properties of the fault, earthquakes at relatively lower levels of accumulated stress than their neighbouring stronger counterparts. Upper mantle seismicity at depths greater than modelled brittle conditions, can be either explained by embrittlement of the mantle due to grain-size sensitive deformation within domains of active or recent slab cooling, or by dissipative weakening mechanisms, such as thermal runaway from shear heating and/or dehydration reactions within an overly ductile mantle.
Results generated in this study are available for open access use to further discussions on the region.
Hundreds of basaltic plateau margins east of the Patagonian Cordillera are undermined by numerous giant slope failures. However, the overall extent of this widespread type of plateau collapse remains unknown and incompletely captured in local maps. To detect giant slope failures consistently throughout the region, we train two convolutional neural networks (CNNs), AlexNet and U-Net, with Sentinel-2 optical data and TanDEM-X topographic data on elevation, surface roughness, and curvature. We validated the performance of these CNNs with independent testing data and found that AlexNet performed better when learned on topographic data, and UNet when learned on optical data. AlexNet predicts a total landslide area of 12,000 km2 in a study area of 450,000 km2, and thus one of Earth's largest clusters of giant landslides. These are mostly lateral spreads and rotational failures in effusive rocks, particularly eroding the margins of basaltic plateaus; some giant landslides occurred along shores of former glacial lakes, but are least prevalent in Quaternary sedimentary rocks. Given the roughly comparable topographic, climatic, and seismic conditions in our study area, we infer that basalts topping weak sedimentary rocks may have elevated potential for large-scale slope failure. Judging from the many newly detected and previously unknown landslides, we conclude that CNNs can be a valuable tool to detect large-scale slope instability at the regional scale. However, visual inspection is still necessary to validate results and correctly outline individual landslide source and deposit areas.
Geostatistics as a subfield of statistics accounts for the spatial correlations encountered in many applications of, for example, earth sciences. Valuable information can be extracted from these correlations, also helping to address the often encountered burden of data scarcity. Despite the value of additional data, the use of geostatistics still falls short of its potential. This problem is often connected to the lack of user-friendly software hampering the use and application of geostatistics. We therefore present GSTools, a Python-based software suite for solving a wide range of geostatistical problems. We chose Python due to its unique balance between usability, flexibility, and efficiency and due to its adoption in the scientific community. GSTools provides methods for generating random fields; it can perform kriging, variogram estimation and much more. We demonstrate its abilities by virtue of a series of example applications detailing their use.
Droughts in São Paulo
(2023)
Literature has suggested that droughts and societies are mutually shaped and, therefore, both require a better understanding of their coevolution on risk reduction and water adaptation. Although the Sao Paulo Metropolitan Region drew attention because of the 2013-2015 drought, this was not the first event. This paper revisits this event and the 1985-1986 drought to compare the evolution of drought risk management aspects. Documents and hydrological records are analyzed to evaluate the hazard intensity, preparedness, exposure, vulnerability, responses, and mitigation aspects of both events. Although the hazard intensity and exposure of the latter event were larger than the former one, the policy implementation delay and the dependency of service areas in a single reservoir exposed the region to higher vulnerability. In addition to the structural and non-structural tools implemented just after the events, this work raises the possibility of rainwater reuse for reducing the stress in reservoirs.
The quantification of spatial propagation of extreme precipitation events is vital in water resources planning and disaster mitigation. However, quantifying these extreme events has always been challenging as many traditional methods are insufficient to capture the nonlinear interrelationships between extreme event time series. Therefore, it is crucial to develop suitable methods for analyzing the dynamics of extreme events over a river basin with a diverse climate and complicated topography. Over the last decade, complex network analysis emerged as a powerful tool to study the intricate spatiotemporal relationship between many variables in a compact way. In this study, we employ two nonlinear concepts of event synchronization and edit distance to investigate the extreme precipitation pattern in the Ganga river basin. We use the network degree to understand the spatial synchronization pattern of extreme rainfall and identify essential sites in the river basin with respect to potential prediction skills. The study also attempts to quantify the influence of precipitation seasonality and topography on extreme events. The findings of the study reveal that (1) the network degree is decreased in the southwest to northwest direction, (2) the timing of 50th percentile precipitation within a year influences the spatial distribution of degree, (3) the timing is inversely related to elevation, and (4) the lower elevation greatly influences connectivity of the sites. The study highlights that edit distance could be a promising alternative to analyze event-like data by incorporating event time and amplitude and constructing complex networks of climate extremes.
Wind erosion of agricultural soils affects their stock of essential elements for plants, like phosphorus (P). It is known that the composition of the eroded sediments varies with height, according to the size and density of the transported substances. Aim of this study was to analyze the concentration and enrichment ratios of P forms in sediments transported by the wind. A wind-tunnel study was performed on a sandy-and a sandy loam soil in order to measure P forms concentrations in the saltating sediments. P concentrations were also measured in the particulate matter (PM) of each soil, gained with the Easy Dust Generator. In both soils, inorganic-(Pi) and organic P (Po) were preferentially transported in PM, with enrichment ratios of 1.8 and 5.5, respectively. Nevertheless, a Pi/Po of 0.9 indicated that the accumulation of the minor Po in PM was more pronounced than Pi. This agrees with P-rich light and easily erodible organic compounds, almost exclusively accumulated in PM, and in relatively heavy and less erodible minerals, like apatites, in lower height sediments. Labile P (Pl) was preferentially transported in saltating sediments of both soils. This was attributed to the selective Bray & Kurtz I's extraction of the abundant inorganic P forms of these sediments. Total P (Pt) copied the transport trends of Pi, the major form. According to the transporting trends, Pi and Po would be re-sedimented at longer distances from the source than Pl. Outcomes become useful for modeling the influence of wind erosion on P cycling.
Land degradation and water availability in semi-arid regions are interdependent challenges for management that are influenced by climatic and anthropogenic changes. Erosion and high sediment loads in rivers cause reservoir siltation and decrease storage capacity, which pose risk on water security for citizens, agriculture, and industry. In regions where resources for management are limited, identifying spatial-temporal variability of sediment sources is crucial to decrease siltation. Despite widespread availability of rigorous methods, approaches simplifying spatial and temporal variability of erosion are often inappropriately applied to very data sparse semi-arid regions. In this work, we review existing approaches for mapping erosional hotspots, and provide an example of spatial-temporal mapping approach in two case study regions. The barriers limiting data availability and their effects on erosion mapping methods, their validation, and resulting prioritization of leverage management areas are discussed.
Understanding the hydrologic connectivity between kettle holes and shallow groundwater, particularly in reaction to the highly variable local meteorological conditions, is of paramount importance for tracing water in a hydro(geo)logically complex landscape and thus for integrated water resource management. This article is aimed at identifying the dominant hydrological processes affecting the kettle holes' water balance and their interactions with the shallow groundwater domain in the Uckermark region, located in the north-east of Germany. For this reason, based on the stable isotopes of oxygen (delta O-18) and hydrogen (delta H-2), an isotopic mass balance model was employed to compute the evaporative loss of water from the kettle holes from February to August 2017. Results demonstrated that shallow groundwater inflow may play the pivotal role in the processes taking part in the hydrology of the kettle holes in the Uckermark region. Based on the calculated evaporation/inflow (E/I) ratios, most of the kettle holes (86.7%) were ascertained to have a partially open, flow-through-dominated system. Moreover, we identified an inverse correlation between E/I ratios and the altitudes of the kettle holes. The same holds for electrical conductivity (EC) and the altitudes of the kettle holes. In accordance with the findings obtained from this study, a conceptual model explaining the interaction between the shallow groundwater and the kettle holes of Uckermark was developed. The model exhibited that across the highest altitudes, the recharge kettle holes are dominant, where a lower ratio of E/I and a lower EC was detected. By contrast, the lowest topographical depressions represent the discharge kettle holes, where a higher ratio of E/I and EC could be identified. The kettle holes existing in between were categorized as flow-through kettle holes through which the recharge takes place from one side and discharge from the other side.
An exceptionally strong stationary planetary wave with Zonal Wavenumber 1 led to a sudden stratospheric warming (SSW) in the Southern Hemisphere in September 2019. Ionospheric data from European Space Agency's Swarm satellite constellation mission show prominent 6-day variations in the dayside low-latitude region at this time, which can be attributed to forcing from the middle atmosphere by the Rossby normal mode "quasi-6-day wave" (Q6DW). Geopotential height measurements by the Microwave Limb Sounder aboard National Aeronautics and Space Administration's Aura satellite reveal a burst of global Q6DW activity in the mesosphere and lower thermosphere during the SSW, which is one of the strongest in the record. The Q6DW is apparently generated in the polar stratosphere at 30-40 km, where the atmosphere is unstable due to strong vertical wind shear connected with planetary wave breaking. These results suggest that an Antarctic SSW can lead to ionospheric variability through wave forcing from the middle atmosphere.
Plain Language Summary: A sudden stratospheric warming (SSW) is an extreme wintertime polar meteorological phenomenon occurring mostly over the Arctic region. Studies have shown that Arctic SSW can influence the entire atmosphere. In September 2019, a rare SSW event occurred in the Antarctic region, providing an opportunity to investigate its broader impact on the whole atmosphere. We present observations from the middle atmosphere and ionosphere during this event, noting unusually strong wave activity throughout this region. Our results suggest that an Antarctic SSW can have a significant impact on the whole atmosphere system similar to those due to Arctic events.
Boreal forests of Siberia play a relevant role in the global carbon cycle. However, global warming threatens the existence of summergreen larch-dominated ecosystems, likely enabling a transition to evergreen tree taxa with deeper active layers. Complex permafrost-vegetation interactions make it uncertain whether these ecosystems could develop into a carbon source rather than continuing atmospheric carbon sequestration under global warming. Consequently, shedding light on the role of current and future active layer dynamics and the feedbacks with the apparent tree species is crucial to predict boreal forest transition dynamics and thus for aboveground forest biomass and carbon stock developments. Hence, we established a coupled model version amalgamating a one-dimensional permafrost multilayer forest land-surface model (CryoGrid) with LAVESI, an individual-based and spatially explicit forest model for larch species (Larix Mill.), extended for this study by including other relevant Siberian forest species and explicit terrain. <br /> Following parameterization, we ran simulations with the coupled version to the near future to 2030 with a mild climate-warming scenario. We focus on three regions covering a gradient of summergreen forests in the east at Spasskaya Pad, mixed summergreen-evergreen forests close to Nyurba, and the warmest area at Lake Khamra in the southeast of Yakutia, Russia. Coupled simulations were run with the newly implemented boreal forest species and compared to runs allowing only one species at a time, as well as to simulations using just LAVESI. Results reveal that the coupled version corrects for overestimation of active layer thickness (ALT) and soil moisture, and large differences in established forests are simulated. We conclude that the coupled version can simulate the complex environment of eastern Siberia by reproducing vegetation patterns, making it an excellent tool to disentangle processes driving boreal forest dynamics.
Previous studies have explored the consequences of flood events for exposed households and companies by focusing on single flood events. Less is known about the consequences of experiencing repeated flood events for the resilience of households and companies. In this paper, we therefore explore how multiple floods experience affects the resilience of exposed households and companies. Resilience was made operational through individual appraisals of households and companies' ability to withstand and recover from material as well as health and psychological impacts of the 2013 flood in Germany. The paper is based on three different datasets including more than 2000 households and 300 companies that were affected by the 2013 flood. The surveys revealed that the resilience of households seems to increase, but only with regard to their subjectively appraised ability to withstand impacts on mobile goods and equipment (e.g., cars, TV, and radios). In regard to the ability of households to withstand overall financial consequences of repetitive floods, evidence for nonlinear (quadratic) trends can be found. With regard to psychological and health-related consequences, the findings are mixed but provide tentative evidence for eroding resilience among households. Companies' resilience increased with respect to material assets but appears to decrease with respect to ability to recover. We conclude by arguing that clear and operational definitions of resilience are required so that evidence-based resilience baselines can be established to assess whether resilience is eroding or improving over time.
Groundwater is the biggest single source of high-quality freshwater worldwide, which is also continuously threatened by the changing climate. In this paper, we investigate the response of the regional groundwater system to climate change under three global warming levels (1.5, 2, and 3 ∘C) in a central German basin (Nägelstedt). This investigation is conducted by deploying an integrated modeling workflow that consists of a mesoscale hydrologic model (mHM) and a fully distributed groundwater model, OpenGeoSys (OGS). mHM is forced with climate simulations of five general circulation models under three representative concentration pathways. The diffuse recharges estimated by mHM are used as boundary forcings to the OGS groundwater model to compute changes in groundwater levels and travel time distributions. Simulation results indicate that groundwater recharges and levels are expected to increase slightly under future climate scenarios. Meanwhile, the mean travel time is expected to decrease compared to the historical average. However, the ensemble simulations do not all agree on the sign of relative change. Changes in mean travel time exhibit a larger variability than those in groundwater levels. The ensemble simulations do not show a systematic relationship between the projected change (in both groundwater levels and travel times) and the warming level, but they indicate an increased variability in projected changes with adjusting the enhanced warming level from 1.5 to 3 ∘C. Correspondingly, it is highly recommended to restrain the trend of global warming.
Submerged sequences of marine terraces potentially provide crucial information of past sea-level positions. However, the distribution and characteristics of drowned marine terrace sequences are poorly known at a global scale. Using bathymetric data and novel mapping and modeling techniques, we studied a submerged sequence of marine terraces in the Bay of Biscay with the objective to identify the distribution and morphologies of submerged marine terraces and the timing and conditions that allowed their formation and preservation. To accomplish the objectives a high-resolution bathymetry (5 m) was analyzed using Geographic Information Systems and TerraceM(R). The successive submerged terraces were identified using a Surface Classification Model, which linearly combines the slope and the roughness of the surface to extract fossil sea-cliffs and fossil rocky shore platforms. For that purpose, contour and hillshaded maps were also analyzed. Then, shoreline angles, a geomorphic marker located at the intersection between the fossil sea-cliff and platform, were mapped analyzing swath profiles perpendicular to the isobaths. Most of the submerged strandlines are irregularly preserved throughout the continental shelf. In summary, 12 submerged terraces with their shoreline angles between approximately: -13 m (T1), -30 and -32 m (T2), -34 and 41 m (T3), -44 and -47 m (T4), -49 and 53 m (T5), -55 and 58 m (T6), -59 and 62 m (T7), -65 and 67 m (T8), -68 and 70 m (T9), -74 and -77 m (T10), -83 and -86 m (T11) and -89 and 92 m (T12). Nevertheless, the ones showing the best lateral continuity and preservation in the central part of the shelf are T3, T4, T5, T7, T8, and T10. The age of the terraces has been estimated using a landscape evolution model. To simulate the formation and preservation of submerged terraces three different scenarios: (i) 20-0 ka; (ii) 128-0 ka; and (iii) 128-20 ka, were compared. The best scenario for terrace generation was between 128 and 20 Ka, where T3, T5, and T7 could have been formed.
Due to the fact that silicon (Si) increases the resistance of plants against diverse abiotic and biotic stresses, Si nowadays is categorized as beneficial substance for plants. However, humans directly influence Si cycling on a global scale. Intensified agriculture and corresponding harvest-related Si exports lead to Si losses in agricultural soils. This anthropogenic desilication might be a big challenge for modern agriculture. However, there is still only little knowledge about Si cycling in agricultural systems of the temperate zone, because most studies focus on rice and sugarcane production in (sub)tropical areas. Furthermore, many studies are performed for a short term only, and thus do not provide the opportunity to analyze slow changes in soil-plant systems (e.g., desilication) over long periods. We analyzed soil and plant samples from an ongoing long-term field experiment (established 1963) in the temperate zone (NE Germany) to evaluate the effects of different nitrogen-phosphoruspotassium (NPK) fertilization rates and crop straw recycling (i.e., straw incorporation) on anthropogenic desilication in the long term. Our results clearly show that crop straw recycling not only prevents anthropogenic desilication (about 43-60% of Si exports can be saved by crop straw recycling in the long term), but also replenishes plant available Si stocks of agricultural soil-plant systems. Furthermore, we found that a reduction of N fertilization rates of about 69% is possible without considerable biomass losses. This economy of the need for N fertilizers potentially can be combined with the benefits of crop straw recycling, i.e., enhancement of carbon sequestration via straw inputs and prevention of anthropogenic desilication of agricultural soil-plant systems. Thus crop straw recycling might have the potential to act as key management practice in sustainable, low fertilization agriculture in the temperate zone in the future.
Pollen records from Siberia are mostly absent in global or Northern Hemisphere synthesis works. Here we present a taxonomically harmonized and temporally standardized pollen dataset that was synthesized using 173 palynological records from Siberia and adjacent areas (northeastern Asia, 42-75 degrees N, 50-180 degrees E). Pollen data were taxonomically harmonized, i.e. the original 437 taxa were assigned to 106 combined pollen taxa. Age-depth models for all records were revised by applying a constant Bayesian age-depth modelling routine. The pollen dataset is available as count data and percentage data in a table format (taxa vs. samples), with age information for each sample. The dataset has relatively few sites covering the last glacial period between 40 and 11.5 ka (calibrated thousands of years before 1950 CE) particularly from the central and western part of the study area. In the Holocene period, the dataset has many sites from most of the area, with the exception of the central part of Siberia. Of the 173 pollen records, 81 % of pollen counts were downloaded from open databases (GPD, EPD, PANGAEA) and 10 % were contributions by the original data gatherers, while a few were digitized from publications. Most of the pollen records originate from peatlands (48 %) and lake sediments (33 %). Most of the records (83 %) have >= 3 dates, allowing the establishment of reliable chronologies. The dataset can be used for various purposes, including pollen data mapping (example maps for Larix at selected time slices are shown) as well as quantitative climate and vegetation reconstructions. The datasets for pollen counts and pollen percentages are available at https://doi.org/10.1594/PANGAEA.898616 (Cao et al., 2019a), also including the site information, data source, original publication, dating data, and the plant functional type for each pollen taxa.