ALBERT

Description: Abstract. Clouds are assumed to play an important role in the Arctic amplification process. This motivated a detailed investigation of cloud processes, including radiative and turbulent fluxes. Data from the aircraft campaign ACLOUD were analyzed with a focus on the mean and turbulent structure of the cloudy boundary layer over the Fram Strait marginal sea ice zone in late spring and early summer 2017. Vertical profiles of turbulence moments are presented from contrasting atmospheric boundary layers (ABLs) from 4 d. They differ by the magnitude of wind speed, boundary-layer height, stability, the strength of the cloud-top radiative cooling and the number of cloud layers. Turbulence statistics up to third-order moments are presented, which were obtained from horizontal-level flights and from slanted profiles. It is shown that both of these flight patterns complement each other and form a data set that resolves the vertical structure of the ABL turbulence well. The comparison of the 4 d shows that especially during weak wind, even in shallow Arctic ABLs with mixing ratios below 3 g kg-1, cloud-top cooling can serve as a main source of turbulent kinetic energy (TKE).Well-mixed ABLs are generated where TKE is increased and vertical velocity variance shows pronounced maxima in the cloud layer. Negative vertical velocity skewness points then to upside-down convection. Turbulent heat fluxes are directed upward in the cloud layer as a result of cold downdrafts. In two cases with single-layer stratocumulus, turbulent transport of heat flux and of temperature variance are both negative in the cloud layer, suggesting an important role of large eddies. In contrast, in a case with weak cloud-top cooling, these quantities are positive in the ABL due to the heating from the surface. Based on observations and results of a mixed-layer model it is shown that the maxima of turbulent fluxes are, however, smaller than the jump of the net terrestrial radiation flux across the upper part of a cloud due to the (i) shallowness of the mixed layer and (ii) the presence of a downward entrainment heat flux. The mixed-layer model also shows that the buoyancy production of TKE is substantially smaller in stratocumulus over the Arctic sea ice compared to subtropics due to a smaller surface moisture flux and smaller decrease in specific humidity (or even humidity inversions) right above the cloud top. In a case of strong wind, wind shear shapes the ABL turbulent structure, especially over rough sea ice, despite the presence of a strong cloud-top cooling. In the presence of mid-level clouds, cloud-top radiative cooling and thus also TKE in the lowermost cloud layer are strongly reduced, and the ABL turbulent structure becomes governed by stability, i.e., by the surface–air temperature difference and wind speed. A comparison of slightly unstable and weakly stable cases shows a strong reduction of TKE due to increased stability even though the absolute value of wind speed was similar. In summary, the presented study documents vertical profiles of the ABL turbulence with a high resolution in a wide range of conditions. It can serve as a basis for turbulence closure evaluation and process studies in Arctic clouds.

Description: The presence of reactive bromine in polar regions is a widespread phenomenon that plays an important role in the photochemistry of the Arctic and Antarctic lower troposphere, including the destruction of ozone, the disturbance of radical cycles, and the oxidation of gaseous elemental mercury. The chemical mechanisms leading to the heterogeneous release of gaseous bromine compounds from saline surfaces are in principle well understood. There are, however, substantial uncertainties about the contribution of different potential sources to the release of reactive bromine, such as sea ice, brine, aerosols, and the snow surface, as well as about the seasonal and diurnal variation and the vertical distribution of reactive bromine. Here we use continuous long-term measurements of the vertical distribution of bromine monoxide (BrO) and aerosols at the two Antarctic sites Neumayer (NM) and Arrival Heights (AH), covering the periods of 2003–2021 and 2012–2021, respectively, to investigate how chemical and physical parameters affect the abundance of BrO. We find the strongest correlation between BrO and aerosol extinction (R=0.56 for NM and R=0.28 for AH during spring), suggesting that the heterogeneous release of Br2 from saline airborne particles (blowing snow and aerosols) is a dominant source for reactive bromine. Positive correlations between BrO and contact time of air masses, both with sea ice and the Antarctic ice sheet, suggest that reactive bromine is not only emitted by the sea ice surface but by the snowpack on the ice shelf and in the coastal regions of Antarctica. In addition, the open ocean appears to represent a source for reactive bromine during late summer and autumn when the sea ice extent is at its minimum. A source–receptor analysis based on back trajectories and sea ice maps shows that main source regions for BrO at NM are the Weddell Sea and the Filchner–Ronne Ice Shelf, as well as coastal polynyas where sea ice is newly formed. A strong morning peak in BrO frequently occurring during summer and that is particularly strong during autumn suggests a night-time build-up of Br2 by heterogeneous reaction of ozone on the saline snowpack in the vicinity of the measurement sites. We furthermore show that BrO can be sustained for at least 3 d while travelling across the Antarctic continent in the absence of any saline surfaces that could serve as a source for reactive bromine.

Description: Springtime Arctic mixed-phase convection over open water in the Fram Strait as observed during the recent ACLOUD (Arctic CLoud Observations Using airborne measurements during polar Day) field campaign is simulated at turbulence-resolving resolutions. The first objective is to assess the skill of large-eddy simulation (LES) in reproducing the observed mixed-phase convection. The second goal is to then use the model to investigate how aerosol modulates the way in which turbulent mixing and clouds transform the low-level air mass. The focus lies on the low-level thermal structure and lapse rate, the heating efficiency of turbulent entrainment, and the low-level energy budget. A composite case is constructed based on data collected by two research aircraft on 18 June 2017. Simulations are evaluated against independent datasets, showing that the observed thermodynamic, cloudy, and turbulent states are well reproduced. Sensitivity tests on cloud condensation nuclei (CCN) concentration are then performed, covering a broad range between pristine polar and polluted continental values. We find a significant response in the resolved mixed-phase convection, which is in line with previous LES studies. An increased CCN substantially enhances the depth of convection and liquid cloud amount, accompanied by reduced surface precipitation. Initializing with the in situ CCN data yields the best agreement with the cloud and turbulence observations, a result that prioritizes its measurement during field campaigns for supporting high-resolution modeling efforts. A deeper analysis reveals that CCN significantly increases the efficiency of radiatively driven entrainment in warming the boundary layer. The marked strengthening of the thermal inversion plays a key role in this effect. The low-level heat budget shifts from surface driven to radiatively driven. This response is accompanied by a substantial reduction in the surface energy budget, featuring a weakened flow of solar radiation into the ocean. Results are interpreted in the context of air–sea interactions, air mass transformations, and climate feedbacks at high latitudes.

Description: The thermokarst lakes of permafrost regions play a major role in the global carbon cycle. These lakes are sources of methane to the atmosphere although the methane flux is restricted by an ice cover for most of the year. How methane concentrations and fluxes in these waters are affected by the presence of an ice cover is poorly understood. To relate water body morphology, ice formation and methane to each other, we studied the ice of three different water bodies in locations typical of the transition of permafrost from land to ocean in a continuous permafrost coastal region in Siberia. In total, 11 ice cores were analyzed as records of the freezing process and methane composition during the winter season. The three water bodies differed in terms of connectivity to the sea, which affected fall freezing. The first was a bay underlain by submarine permafrost (Tiksi Bay, BY), the second a shallow thermokarst lagoon cut off from the sea in winter (Polar Fox Lagoon, LG) and the third a land-locked freshwater thermokarst lake (Goltsovoye Lake, LK). Ice on all water bodies was mostly methane-supersaturated with respect to atmospheric equilibrium concentration, except for three cores from the isolated lake. In the isolated thermokarst lake, ebullition from actively thawing basin slopes resulted in the localized integration of methane into winter ice. Stable δ13C-CH4 isotope signatures indicated that methane in the lagoon ice was oxidized to concentrations close to or below the calculated atmospheric equilibrium concentration. Increasing salinity during winter freezing led to a micro-environment on the lower ice surface where methane oxidation occurred and the lagoon ice functioned as a methane sink. In contrast, the ice of the coastal marine environment was slightly supersaturated with methane, consistent with the brackish water below. Our interdisciplinary process study shows how water body morphology affects ice formation which mitigates methane fluxes to the atmosphere.

Description: Methane emissions from boreal and arctic wetlands, lakes, and rivers are expected to increase in response to warming and associated permafrost thaw. However, the lack of appropriate land cover datasets for scaling field-measured methane emissions to circumpolar scales has contributed to a large uncertainty for our understanding of present-day and future methane emissions. Here we present the Boreal–Arctic Wetland and Lake Dataset (BAWLD), a land cover dataset based on an expert assessment, extrapolated using random forest modelling from available spatial datasets of climate, topography, soils, permafrost conditions, vegetation, wetlands, and surface water extents and dynamics. In BAWLD, we estimate the fractional coverage of five wetland, seven lake, and three river classes within 0.5 × 0.5∘ grid cells that cover the northern boreal and tundra biomes (17 % of the global land surface). Land cover classes were defined using criteria that ensured distinct methane emissions among classes, as indicated by a co-developed comprehensive dataset of methane flux observations. In BAWLD, wetlands occupied 3.2 × 106 km2 (14 % of domain) with a 95 % confidence interval between 2.8 and 3.8 × 106 km2. Bog, fen, and permafrost bog were the most abundant wetland classes, covering ∼ 28 % each of the total wetland area, while the highest-methane-emitting marsh and tundra wetland classes occupied 5 % and 12 %, respectively. Lakes, defined to include all lentic open-water ecosystems regardless of size, covered 1.4 × 106 km2 (6 % of domain). Low-methane-emitting large lakes (〉10 km2) and glacial lakes jointly represented 78 % of the total lake area, while high-emitting peatland and yedoma lakes covered 18 % and 4 %, respectively. Small (〈0.1 km2) glacial, peatland, and yedoma lakes combined covered 17 % of the total lake area but contributed disproportionally to the overall spatial uncertainty in lake area with a 95 % confidence interval between 0.15 and 0.38 × 106 km2. Rivers and streams were estimated to cover 0.12  × 106 km2 (0.5 % of domain), of which 8 % was associated with high-methane-emitting headwaters that drain organic-rich landscapes. Distinct combinations of spatially co-occurring wetland and lake classes were identified across the BAWLD domain, allowing for the mapping of “wetscapes” that have characteristic methane emission magnitudes and sensitivities to climate change at regional scales. With BAWLD, we provide a dataset which avoids double-accounting of wetland, lake, and river extents and which includes confidence intervals for each land cover class. As such, BAWLD will be suitable for many hydrological and biogeochemical modelling and upscaling efforts for the northern boreal and arctic region, in particular those aimed at improving assessments of current and future methane emissions. Data are freely available at https://doi.org/10.18739/A2C824F9X (Olefeldt et al., 2021).

Description: Lakes in permafrost regions are dynamiclandscape components and play an important role for climatechange feedbacks. Lake processes such as mineralizationand flocculation of dissolved organic carbon (DOC), oneof the main carbon fractions in lakes, contribute to thegreenhouse effect and are part of the global carbon cycle.These processes are in the focus of climate research, butstudies so far are limited to specific study regions. Inour synthesis, we analyzed 2167 water samples from 1833lakes across the Arctic in permafrost regions of Alaska,Canada, Greenland, and Siberia to provide first pan-Arcticinsights for linkages between DOC concentrations andthe environment. Using published data and unpublisheddatasets from the author team, we report regional DOCdifferences linked to latitude, permafrost zones, ecoregions,geology, near-surface soil organic carbon contents, andground ice classification of each lake region. The lakeDOC concentrations in our dataset range from 0 to1130 mg L−1(10.8 mg L−1median DOC concentration).Regarding the permafrost regions of our synthesis, wefound median lake DOC concentrations of 12.4 mg L−1(Siberia), 12.3 mg L−1(Alaska), 10.3 mg L−1(Greenland),and 4.5 mg L−1(Canada). Our synthesis shows a significantrelationship between lake DOC concentration and lakeecoregion. We found higher lake DOC concentrationsat boreal permafrost sites compared to tundra sites. Wefound significantly higher DOC concentrations in lakesin regions with ice-rich syngenetic permafrost deposits(yedoma) compared to non-yedoma lakes and a weak butsignificant relationship between soil organic carbon contentand lake DOC concentration as well as between ground icecontent and lake DOC. Our pan-Arctic dataset shows that theDOC concentration of a lake depends on its environmentalproperties, especially on permafrost extent and ecoregion, aswell as vegetation, which is the most important driver of lakeDOC in this study. This new dataset will be fundamental toquantify a pan-Arctic lake DOC pool for estimations of theimpact of lake DOC on the global carbon cycle and climatechange.

Description: Arctic river deltas and deltaic near-shore zones represent important land–ocean transition zones influencing sediment dynamics and nutrient fluxes from permafrost-affected terrestrial ecosystems into the coastal Arctic Ocean. To accurately model fluvial carbon and freshwater export from rapidly changing river catchments as well as assess impacts of future change on the Arctic shelf and coastal ecosystems, we need to understand the sea floor characteristics and topographic variety of the coastal zones. To date, digital bathymetrical data from the poorly accessible, shallow, and large areas of the eastern Siberian Arctic shelves are sparse. We have digitized bathymetrical information for nearly 75 000 locations from large-scale (1:25 000–1:500 000) current and historical nautical maps of the Lena Delta and the Kolyma Gulf region in northeastern Siberia. We present the first detailed and seamless digital models of coastal zone bathymetry for both delta and gulf regions in 50 and 200 m spatial resolution. We validated the resulting bathymetry layers using a combination of our own water depth measurements and a collection of available depth measurements, which showed a strong correlation (r〉0.9). Our bathymetrical models will serve as an input for a high-resolution coupled hydrodynamic–ecosystem model to better quantify fluvial and coastal carbon fluxes to the Arctic Ocean, but they may be useful for a range of other studies related to Arctic delta and near-shore dynamics such as modeling of submarine permafrost, near-shore sea ice, or shelf sediment transport. The new digital high-resolution bathymetry products are available on the PANGAEA data set repository for the Lena Delta (https://doi.org/10.1594/PANGAEA.934045; Fuchs et al., 2021a) and Kolyma Gulf region (https://doi.org/10.1594/PANGAEA.934049; Fuchs et al., 2021b), respectively. Likewise, the depth validation data are available on PANGAEA as well (https://doi.org/10.1594/PANGAEA.933187; Fuchs et al., 2021c).

Description: TR17-08, a marine sedimentary core (14.6 m), was collected during 2017 from the Edisto Inlet (Ross Sea, Antarctica), a small fjord near Cape Hallett. The core is characterized by expanded laminated sedimentary sequences making it suitable for studying submillennial processes during the Early Holocene. By studying different well-known foraminifera species (Globocassidulina biora, G. subglobosa, Trifarina angulosa, Nonionella iridea, Epistominella exigua, Stainforthia feylingi, Miliammina arenacea, Paratrochammina bartrami and Portatrochammina antarctica), we were able to identify five different foraminiferal assemblages over the last ∼ 2000 years BP. Comparison with diatom assemblages and other geochemical proxies retrieved from nearby sediment cores in the Edisto Inlet (BAY05-20 and HLF17-1) made it possible to distinguish three different phases characterized by different environmental settings: (1) a seasonal phase (from 2012 to 1486 years BP) characterized by the dominance of calcareous species, indicating a seasonal opening of the inlet by more frequent events of melting of the sea-ice cover during the austral summer and, in general, a higher-productivity, more open and energetic environment; (2) a transitional phase (from 1486 to 696 years BP) during which the fjord experienced less extensive sea-ice melting, enhanced oxygen-poor conditions and carbonate dissolution conditions, indicated by the shifts from calcareous-dominated association to agglutinated-dominated association probably due to a freshwater input from the retreat of three local glaciers at the start of this period; and (3) a cooler phase (from 696 years BP to present) during which the sedimentation rate decreased and few to no foraminiferal specimens were present, indicating ephemeral openings or a more prolonged cover of the sea ice during the austral summer, affecting the nutrient supply and the sedimentation regime.

Description: Published

Description: 95–115

Description: OSA2: Evoluzione climatica: effetti e loro mitigazione

Description: JCR Journal

Description: 〈jats:p〉Abstract. The risk of carbon emissions from permafrost ground is linked to ground temperature and thus in particular to thermal insulation by vegetation and organic soil layers in summer and snow cover in winter. This ground insulation is strongly influenced by the presence of large herbivorous animals browsing for food. In this study, we examine the potential impact of large herbivore presence on the ground carbon storage in thermokarst landscapes of northeastern Siberia. Our aim is to understand how intensive animal grazing may affect permafrost thaw and hence organic matter decomposition, leading to different ground carbon storage, which is significant in the active layer. Therefore, we analysed sites with differing large herbivore grazing intensity in the Pleistocene Park near Chersky and measured maximum thaw depth, total organic carbon content and decomposition state by δ13C isotope analysis. In addition, we determined sediment grain size composition as well as ice and water content. We found the thaw depth to be shallower and carbon storage to be higher in intensively grazed areas compared to extensively and non-grazed sites in the same thermokarst basin. The intensive grazing presumably leads to a more stable thermal ground regime and thus to increased carbon storage in the thermokarst deposits and active layer. However, the high carbon content found within the upper 20 cm on intensively grazed sites could also indicate higher carbon input rather than reduced decomposition, which requires further studies. We connect our findings to more animal trampling in winter, which causes snow disturbance and cooler winter ground temperatures during the average annual 225 days below freezing. This winter cooling overcompensates ground warming due to the lower insulation associated with shorter heavily grazed vegetation during the average annual 140 thaw days. We conclude that intensive grazing influences the carbon storage capacities of permafrost areas and hence might be an actively manageable instrument to reduce net carbon emission from these sites. 〈/jats:p〉

Description: We combine satellite data products to provide a first and general overview of the physical sea ice conditions along the drift of the international Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition and a comparison with previous years (2005–2006 to 2018–2019). We find that the MOSAiC drift was around 20 % faster than the climatological mean drift, as a consequence of large-scale low-pressure anomalies prevailing around the Barents–Kara–Laptev sea region between January and March. In winter (October–April), satellite observations show that the sea ice in the vicinity of the Central Observatory (CO; 50 km radius) was rather thin compared to the previous years along the same trajectory. Unlike ice thickness, satellite-derived sea ice concentration, lead frequency and snow thickness during winter months were close to the long-term mean with little variability. With the onset of spring and decreasing distance to the Fram Strait, variability in ice concentration and lead activity increased. In addition, the frequency and strength of deformation events (divergence, convergence and shear) were higher during summer than during winter. Overall, we find that sea ice conditions observed within 5 km distance of the CO are representative for the wider (50 and 100 km) surroundings. An exception is the ice thickness; here we find that sea ice within 50 km radius of the CO was thinner than sea ice within a 100 km radius by a small but consistent factor (4 %) for successive monthly averages. Moreover, satellite acquisitions indicate that the formation of large melt ponds began earlier on the MOSAiC floe than on neighbouring floes.

Description: 〈jats:p〉Abstract. The combined effect of hot and dry extremes can have disastrous consequences for the society, economy, and the environment. While a significant number of studies have been conducted regarding the variability of the individual hot or dry extremes in Romania, the evaluation of the combined effect of these extremes (e.g. compound effect) is still lacking for this region. Thus, in this study we have assessed the spatio-temporal variability and trends of hot and dry summers in the eastern part of Europe, focusing on Romania, between 1950 and 2020 and we have analyzed the relationship between the frequency of hot summers and the prevailing large-scale atmospheric circulation. The length, spatial extent and frequency of HWs in Romania has increased significantly over the last 70 years, while for the drought conditions no significant changes have been observed. The rate of increase in the frequency and spatial extent of HWs has accelerated significantly after the 1990’s, while the smallest number of HWs was observed between 1970 and 1985. The hottest years, in terms of heatwave duration and frequency, were 2007, 2012, 2015, and 2019. One of the key drivers of hot summers, over our analyzed region, is the prevailing large-scale circulation, featuring an anticyclonic circulation over the central and eastern parts of Europe and enhanced atmospheric blocking activity associated with positive temperature anomalies underneath. We conclude that our study can help improve our understanding of the spatio-temporal variability of hot and dry summers, especially at the regional scale, as well as their driving mechanisms which might lead to a better predictability of these extreme events. 〈/jats:p〉

Description: In this study, the first fully continuous monitoring of water vapour isotopic composition at Neumayer Station III, Antarctica, during the 2-year period from February 2017 to January 2019 is presented. Seasonal and synoptic-scale variations in both stable water isotopes H182O and HDO are reported, and their links to variations in key meteorological variables are analysed. In addition, the diurnal cycle of isotope variations during the summer months (December and January 2017/18 and 2018/19) has been examined. Changes in local temperature and specific humidity are the main drivers for the variability in δ18O and δD in vapour at Neumayer Station III, on both seasonal and shorter timescales. In contrast to the measured δ18O and δD variations, no seasonal cycle in the Deuterium excess signal (d) in vapour is detected. However, a rather high uncertainty in measured d values especially in austral winter limits the confidence of this finding. Overall, the d signal shows a stronger inverse correlation with specific humidity than with temperature, and this inverse correlation between d and specific humidity is stronger for the cloudy-sky conditions than for clear-sky conditions during summertime. Back-trajectory simulations performed with the FLEXPART model show that seasonal and synoptic variations in δ18O and δD in vapour coincide with changes in the main sources of water vapour transported to Neumayer Station III. In general, moisture transport pathways from the east lead to higher temperatures and more enriched δ18O values in vapour, while weather situations with southerly winds lead to lower temperatures and more depleted δ18O values. However, on several occasions, δ18O variations linked to wind direction changes were observed, which were not accompanied by a corresponding temperature change. Comparing isotopic compositions of water vapour at Neumayer Station III and snow samples taken in the vicinity of the station reveals almost identical slopes, both for the δ18O–δD relation and for the temperature–δ18O relation.

Description: Stable water isotopologues of snow, firn and ice cores provide valuable information on past climate variations. Yet single profiles are generally not suitable for robust climate reconstructions. Stratigraphic noise, introduced by the irregular deposition, wind-driven erosion and redistribution of snow, impacts the utility of high-resolution isotope records, especially in low-Accumulation areas. However, it is currently unknown how stratigraphic noise differs across the East Antarctic Plateau and how it is affected by local environmental conditions. Here, we assess the amount and structure of stratigraphic noise at seven sites along a 120 km transect on the plateau of Dronning Maud Land, East Antarctica. Replicated oxygen isotope records of 1 m length were used to estimate signal-To-noise ratios as a measure of stratigraphic noise at sites characterised by different accumulation rates (43-64 mm w.e. a-1), snow surface roughnesses and slope inclinations. While we found a high level of stratigraphic noise at all sites, there was also considerable variation between sites. At sastrugi-dominated sites, greater stratigraphic noise coincided with stronger surface roughnesses, steeper slopes and lower accumulation rates, probably related to increased wind speeds. These results provide a first step to modelling stratigraphic noise and might guide site selection and sampling strategies for future expeditions to improve high-resolution climate reconstructions from low-Accumulation regions.

Description: Despite the importance of surface energy budgets (SEBs) for land-climate interactions in the Arctic, uncertainties in their prediction persist. In situ observational data of SEB components - useful for research and model validation - are collected at relatively few sites across the terrestrial Arctic, and not all available datasets are readily interoperable. Furthermore, the terrestrial Arctic consists of a diversity of vegetation types, which are generally not well represented in land surface schemes of current Earth system models. This dataset describes the data generated in a literature synthesis, covering 358 study sites on vegetation or glacier (〉=60°N latitude), which contained surface energy budget observations. The literature synthesis comprised 148 publications searched on the ISI Web of Science Core Collection.

Description: Organic carbon (OC) stored in Arctic permafrost represents one of Earth's largest and most vulnerable terrestrial carbon pools. Amplified climate warming across the Arctic results in widespread permafrost thaw. Permafrost deposits exposed at river cliffs and coasts are particularly susceptible to thawing processes. Accelerating erosion of terrestrial permafrost along shorelines leads to increased transfer of organic matter (OM) to nearshore waters. However, the amount of terrestrial permafrost carbon and nitrogen as well as the OM quality in these deposits are still poorly quantified. Here, we characterise the sources and the quality of OM supplied to the Lena River at a rapidly eroding permafrost river shoreline cliff in the eastern part of the delta (Sobo-Sise Island). Our multi-proxy approach captures bulk elemental, molecular geochemical and carbon isotopic analyses of late Pleistocene Yedoma permafrost and Holocene cover deposits, discontinuously spanning the last ~52 ka. We show that the ancient permafrost exposed in the Sobo-Sise cliff has a high organic carbon content (mean of about 5 wt%).We found that the OM quality, which we define as the intrinsic potential to further transformation, decomposition, and mineralization, is also high as inferred by the lipid biomarker inventory. The oldest sediments stem from Marine Isotope Stage (MIS) 3 interstadial deposits (dated to 52 to 28 cal kyr BP) and is overlaid by Last Glacial MIS 2 (dated to 28 to 15 cal ka BP) and Holocene MIS 1 (dated to 7–0 cal ka BP) deposits. The relatively high average chain length (ACL) index of n-alkanes along the cliff profile indicates a predominant contribution of vascular plants to the OM composition. The elevated ratio of iso and anteiso-branched FAs relative to long chain (C ≥ 20) n-FAs in the interstadial MIS 3 and the interglacial MIS 1 deposits, suggests stronger microbial activity and consequently higher input of bacterial biomass during these climatically warmer periods. The overall high carbon preference index (CPI) and higher plant fatty acid (HPFA) values as well as high C / N ratios point to a good quality of the preserved OM and thus to a high potential of the OM for decomposition upon thaw. A decrease of HPFA values downwards along the profile probably indicates a relatively stronger OM decomposition in the oldest (MIS 3) deposits of the cliff.

Description: Surface air temperature measurements obtained from different sensors are used to construct a unique time series with one minute time-interval. Apart from differences in design and environmental exposition, periods of missing data also exist in the data series of each sensor. A primary data set was selected in terms of quality and temporal extension. A combination of two different techniques is applied to complete this data set: one is based on the autocorrelation of the series and the other on measurements taken from other sensors. The resulting values constitute a complete series of surface air temperature at AGGO.

Description: The rewetting of peatlands is a promising measure to mitigate greenhouse gas (GHG) emissions by preventing the further mineralization of the peat soil through aeration. In coastal peatland, the rewetting with brackish water can increase the GHG mitigation potential by the introduction of sulfate, a terminal electron acceptor (TEA). Sulfate is known to lower the CH4 production and thus, its emission by favoring the growth of sulfate-reducers, which outcompete methanogens for substrate. The data contain porewater variables such as pH, electrical conductivity (EC) and sulfate, chloride, dissolved CO2 and CH4 concentrations, as well as absolute abundances of methane- and sulfate-cycling microbial communities. The data were collected in spring and autumn 2019 after a storm surge with brackish water inflow in January 2019. Field sampling was conducted in the nature reserve Heiligensee and Hütelmoor in North-East Germany, close to the Southern Baltic Sea coast. We took peat cores using a Russian peat corer in addition to pore water diffusion samplers and plastic liners (length: 60cm; inner diameter 10 cm) at four locations along a transect from further inland towards the Baltic Sea. We wanted to compare the soil and pore water geochemistry as well as the microbial communities after the brackish water inflow to the common freshwater rewetting state. Pore water was extracted using pore water suction samplers in the lab and environmental variables were quantified with an ICP. Microbial samples were sampled from the peat core using sterile equipment. We used quantitative polymerase chain reaction (qPCR) to characterize pools of DNA and cDNA targeting total and putatively active bacteria and archaea. qPCR was performed on key functional genes of methane production (mcrA), aerobic methane oxidation (pmoA) and sulfate reduction (dsrB) in addition to the 16S rRNA gene for the absolute abundance of total prokaryotes. Furthermore, we retrieved soil plugs to determine the concentrations and isotopic signatures of dissolved trace gases (CO2/DIC and CH4) in the pore water.

Description: This dataset provides annually resolved microfacies data from ICDP core 5017-1-A, retrieved from the deep northern Dead Sea basin in 2010/11, for the last glacial-interglacial transition (ca. 14-13 ka BP). Sediments of the Lisan Formation were investigated between ~94.7 and 91.8 m sediment depth below lake floor (lithozone C2) by continuous thin section microscopy. Thin sections were prepared following the standard procedure by Brauer and Casanova (2001) that was adjusted for salty sediments. Thin section analyses were performed on overlapping large-scale thin sections using a Zeiss Axiolab pol microscope at magnifications of 50-400x. Microfacies analyses included varve counting and measurements of varve and sublayer thickness. The amount of varves in erosional gaps was interpolated and the position of mass flow deposits (MFD) is marked.

Description: We present sea ice temperature and salinity data from first-year ice (FYI) and second-year ice (SYI) relevant to the temporal development of sea ice permeability and brine drainage efficiency from the early growth phase in October 2019 to the onset of spring warming in May 2020. Our dataset was collected in the central Arctic Ocean during the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) Expedition in 2019 to 2020. MOSAiC was an international transpolar drift expedition in which the German icebreaker RV Polarstern anchored into an ice floe to gain new insights into Arctic climate over a full annual cycle. In October 2019, RV Polarstern moored to an ice floe in the Siberian sector of the Arctic at 85 degrees north and 137 degrees east to begin the drift towards the North Pole and the Fram Strait via the Transpolar Drift Stream. The data presented here were collected during the first three legs of the expedition, so all the coring activities took place on the same floe. The end dates of legs 1, 2, and 3 were 13 December, 24 February, and 4 June, respectively. The dataset contributed to a baseline study entitled, Deciphering the properties of different Arctic ice types during the growth phase of the MOSAiC floes: Implications for future studies. The study highlights downward directed gas pathways in FYI and SYI by inferring sea ice permeability and potential brine release from several time series of temperature and salinity measurements. The physical properties presented in this paper lay the foundation for subsequent analyses on actual gas contents measured in the ice cores, as well as air-ice and ice-ocean gas fluxes. Sea ice cores were collected with a Kovacs Mark II 9 cm diameter corer. To measure ice temperatures, about 4.5 cm deep holes were drilled into the core (intervals varied by site and leg) . The temperatures were measured by a digital thermometer within minutes after the cores were retrieved. The ice cores were placed into pre-labelled plastic sleeves sealed at the bottom end. The ice cores were transported to RV Polarstern and stored in a -20 degrees Celsius freezer. Each of the cores was sub-sampled, melted at room temperature, and processed for salinity within one or two days. The practical salinity was estimated by measuring the electrical conductivity and temperature of the melted samples using a WTW Cond 3151 salinometer equipped with a Tetra-Con 325 four-electrode conductivity cell. The practical salinity represents the the salinity estimated from the electrical conductivity of the solution. The dataset also contains derived variables, including sea ice density, brine volume fraction, and the Rayleigh number.

Description: Organic carbon (OC) stored in Arctic permafrost represents one of Earth’s largest and most vulnerable terrestrial carbon pools. Amplified climate warming across the Arctic results in widespread permafrost thaw. Permafrost deposits exposed at river cliffs and coasts are particularly susceptible to thawing processes. Accelerating erosion of terrestrial permafrost along shorelines leads to increased transfer of organic matter (OM) to nearshore waters. However, the amount of terrestrial permafrost carbon and nitrogen as well as the OM quality in these deposits are still poorly quantified. Here, we characterise the sources and the quality of OM supplied to the Lena River at a rapidly eroding permafrost river shoreline cliff in the eastern part of the delta (Sobo-Sise Island). Our multi-proxy approach captures bulk elemental, molecular geochemical and carbon isotopic analyses of late Pleistocene Yedoma permafrost and Holocene cover deposits, discontinuously spanning the last ~52 ka. We show that the ancient permafrost exposed in the Sobo-Sise cliff has a high organic carbon content (mean of about 5 wt%).We found that the OM quality, which we define as the intrinsic potential to further transformation, decomposition, and mineralization, is also high as inferred by the lipid biomarker inventory. The oldest sediments stem from Marine Isotope Stage (MIS) 3 interstadial deposits (dated to 52 to 28 cal kyr BP) and is overlaid by Last Glacial MIS 2 (dated to 28 to 15 cal ka BP) and Holocene MIS 1 (dated to 7–0 cal ka BP) deposits. The relatively high average chain length (ACL) index of n-alkanes along the cliff profile indicates a predominant contribution of vascular plants to the OM composition. The elevated ratio of iso and anteiso-branched FAs relative to long chain (C ≥ 20) n-FAs in the interstadial MIS 3 and the interglacial MIS 1 deposits, suggests stronger microbial activity and consequently higher input of bacterial biomass during these climatically warmer periods. The overall high carbon preference index (CPI) and higher plant fatty acid (HPFA) values as well as high C / N ratios point to a good quality of the preserved OM and thus to a high potential of the OM for decomposition upon thaw. A decrease of HPFA values downwards along the profile probably indicates a relatively stronger OM decomposition in the oldest (MIS 3) deposits of the cliff.

Description: We present sea ice temperature and salinity data from first-year ice (FYI) and second-year ice (SYI) relevant to the temporal development of sea ice permeability and brine drainage efficiency from the early growth phase in October 2019 to the onset of spring warming in May 2020. Our dataset was collected in the central Arctic Ocean during the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) Expedition in 2019 to 2020. MOSAiC was an international transpolar drift expedition in which the German icebreaker RV Polarstern anchored into an ice floe to gain new insights into Arctic climate over a full annual cycle. In October 2019, RV Polarstern moored to an ice floe in the Siberian sector of the Arctic at 85 degrees north and 137 degrees east to begin the drift towards the North Pole and the Fram Strait via the Transpolar Drift Stream. The data presented here were collected during the first three legs of the expedition, so all the coring activities took place on the same floe. The end dates of legs 1, 2, and 3 were 13 December, 24 February, and 4 June, respectively. The dataset contributed to a baseline study entitled, Deciphering the properties of different Arctic ice types during the growth phase of the MOSAiC floes: Implications for future studies. The study highlights downward directed gas pathways in FYI and SYI by inferring sea ice permeability and potential brine release from several time series of temperature and salinity measurements. The physical properties presented in this paper lay the foundation for subsequent analyses on actual gas contents measured in the ice cores, as well as air-ice and ice-ocean gas fluxes. Sea ice cores were collected with a Kovacs Mark II 9 cm diameter corer. To measure ice temperatures, about 4.5 cm deep holes were drilled into the core (intervals varied by site and leg) . The temperatures were measured by a digital thermometer within minutes after the cores were retrieved. The ice cores were placed into pre-labelled plastic sleeves sealed at the bottom end. The ice cores were transported to RV Polarstern and stored in a -20 degrees Celsius freezer. Each of the cores was sub-sampled, melted at room temperature, and processed for salinity within one or two days. The practical salinity was estimated by measuring the electrical conductivity and temperature of the melted samples using a WTW Cond 3151 salinometer equipped with a Tetra-Con 325 four-electrode conductivity cell. The practical salinity represents the the salinity estimated from the electrical conductivity of the solution. The dataset also contains derived variables, including sea ice density, brine volume fraction, and the Rayleigh number.

Description: These datasets describe sediment samples taken from the Batagay megaslump, located in Yana Uplands in northeastern Siberia. Most sediment samples were taken from the slump headwall (B19-P1) by rapelling down on a rope from the slump surface and taking samples with a hole saw (diameter 55 mm, 40 mm deep) mounted on a handheld power drill. A second profile (B19-02) of the lowest part of the slump headwall was sampled (~100 m south) using a hammer and axe from the slump floor. Two permafrost sediment blocks (B19-03 and B19-04) at the slump bottom that had fallen from the headwall were sampled using a chainsaw. Finally, a baidzherakh (thermokarst mound; B19-05) in the north of the slump was sampled using a hammer and axe. The samples cover 5 stratigraphical units: 1. lower ice complex, 2. lower sand unit, 3. woody layer, 4. upper ice complex, 5. Holocene cover.

Description: Organic carbon (OC) stored in Arctic permafrost represents one of Earth’s largest and most vulnerable terrestrial carbon pools. Amplified climate warming across the Arctic results in widespread permafrost thaw. Permafrost deposits exposed at river cliffs and coasts are particularly susceptible to thawing processes. Accelerating erosion of terrestrial permafrost along shorelines leads to increased transfer of organic matter (OM) to nearshore waters. However, the amount of terrestrial permafrost carbon and nitrogen as well as the OM quality in these deposits are still poorly quantified. Here, we characterise the sources and the quality of OM supplied to the Lena River at a rapidly eroding permafrost river shoreline cliff in the eastern part of the delta (Sobo-Sise Island). Our multi-proxy approach captures bulk elemental, molecular geochemical and carbon isotopic analyses of late Pleistocene Yedoma permafrost and Holocene cover deposits, discontinuously spanning the last ~52 ka. We show that the ancient permafrost exposed in the Sobo-Sise cliff has a high organic carbon content (mean of about 5 wt%).We found that the OM quality, which we define as the intrinsic potential to further transformation, decomposition, and mineralization, is also high as inferred by the lipid biomarker inventory. The oldest sediments stem from Marine Isotope Stage (MIS) 3 interstadial deposits (dated to 52 to 28 cal kyr BP) and is overlaid by Last Glacial MIS 2 (dated to 28 to 15 cal ka BP) and Holocene MIS 1 (dated to 7–0 cal ka BP) deposits. The relatively high average chain length (ACL) index of n-alkanes along the cliff profile indicates a predominant contribution of vascular plants to the OM composition. The elevated ratio of iso and anteiso-branched FAs relative to long chain (C ≥ 20) n-FAs in the interstadial MIS 3 and the interglacial MIS 1 deposits, suggests stronger microbial activity and consequently higher input of bacterial biomass during these climatically warmer periods. The overall high carbon preference index (CPI) and higher plant fatty acid (HPFA) values as well as high C / N ratios point to a good quality of the preserved OM and thus to a high potential of the OM for decomposition upon thaw. A decrease of HPFA values downwards along the profile probably indicates a relatively stronger OM decomposition in the oldest (MIS 3) deposits of the cliff.

Description: We present sea ice temperature and salinity data from first-year ice (FYI) and second-year ice (SYI) relevant to the temporal development of sea ice permeability and brine drainage efficiency from the early growth phase in October 2019 to the onset of spring warming in May 2020. Our dataset was collected in the central Arctic Ocean during the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) Expedition in 2019 to 2020. MOSAiC was an international transpolar drift expedition in which the German icebreaker RV Polarstern anchored into an ice floe to gain new insights into Arctic climate over a full annual cycle. In October 2019, RV Polarstern moored to an ice floe in the Siberian sector of the Arctic at 85 degrees north and 137 degrees east to begin the drift towards the North Pole and the Fram Strait via the Transpolar Drift Stream. The data presented here were collected during the first three legs of the expedition, so all the coring activities took place on the same floe. The end dates of legs 1, 2, and 3 were 13 December, 24 February, and 4 June, respectively. The dataset contributed to a baseline study entitled, Deciphering the properties of different Arctic ice types during the growth phase of the MOSAiC floes: Implications for future studies. The study highlights downward directed gas pathways in FYI and SYI by inferring sea ice permeability and potential brine release from several time series of temperature and salinity measurements. The physical properties presented in this paper lay the foundation for subsequent analyses on actual gas contents measured in the ice cores, as well as air-ice and ice-ocean gas fluxes. Sea ice cores were collected with a Kovacs Mark II 9 cm diameter corer. To measure ice temperatures, about 4.5 cm deep holes were drilled into the core (intervals varied by site and leg) . The temperatures were measured by a digital thermometer within minutes after the cores were retrieved. The ice cores were placed into pre-labelled plastic sleeves sealed at the bottom end. The ice cores were transported to RV Polarstern and stored in a -20 degrees Celsius freezer. Each of the cores was sub-sampled, melted at room temperature, and processed for salinity within one or two days. The practical salinity was estimated by measuring the electrical conductivity and temperature of the melted samples using a WTW Cond 3151 salinometer equipped with a Tetra-Con 325 four-electrode conductivity cell. The practical salinity represents the the salinity estimated from the electrical conductivity of the solution. The dataset also contains derived variables, including sea ice density, brine volume fraction, and the Rayleigh number.

Description: Rock magnetic and paleomagnetic results covering the past 30 ka were constructed from two sediment cores MSM33_856-1 (MSM33-55-1) and MSM33_855-1 (54-3) from the Black Sea. After the Mediterranean Sea water ingression, finely laminated organic-rich sapropelic sediments and coccolith oozes were deposited in the Black Sea since about 8.3 ka. Relict magnetic minerals in the Black Sea sarpoples are ferrous hemoilmenite, Fe-Mn and Fe-Cr spinels, and magnetite inclusions. In sediments deposited between about 14 and 8 ka, greigite and pyrite were formed in sediments because of the seawater penetration from overlying sediments after the seawater ingression. Before ~14 ka, the Black Sea sediments are dominated by detrital (titano-)magnetite minerals and the sporadically formed greigite which has SIRM/kLF ratios 〉 10 kAm-2. By comparison with detrital (titano-)magnetite samples between 20-30 ka, we found that relict magnetic mineral samples between 0-8.3 ka have similar behavior in recording the geomagnetic field. Moreover, the geomagnetic field variations reconstructed from the Black Sea sapropels are comparable with other validated regional datasets for the past 8.3 ka. The natural remanent magnetization (NRM) and the anhysteretic remanent magnetization (ARM) were measured with a 2G Enterprises 755 SRM (cryogenic) long-core magnetometer equipped with a sample holder for eight discrete samples at a separation of 20 cm. The magnetometer's in-line tri-axial alternating field (AF) demagnetizer was used to demagnetize the NRM and ARM of the samples. The NRM was measured after application of AF peak amplitudes of 0, 5, 10, 15, 20, 30, 40, 50, 65, 80, and 100 mT. Directions of the characteristic remanent magnetization (ChRM) were determined by principle component analysis (PCA) according to Kirschvink (1980). The error range of the ChRM is given as the maximum angular deviation (MAD). The ARM was imparted along the samples' z-axis with a static field of 0.05 mT and an AF field of 100 mT. Demagnetization then was performed in steps of 0, 10, 20, 30, 40, 50, 65, and 80 mT. The median destructive field of the ARM (MDFARM) was determined to estimate the coercivity of the sediments. The slope of NRM versus ARM of common demagnetization steps was used to determine the relative paleointensity (RPI). In most cases, demagnetization steps from 20 to 65 mT were used to determine the RPI.

Description: This dataset provides data for four third-degree tidal constituents used in the publication of Sulzbach et al (2022). The tidal constituents provided are the 3M1, 3M3, 3N2 and 3L2 for 134 globally distributed stations. The tide information, such as the nodal modulations of these tides, are taken from Table 1 and Table S2 of Ray (2020). These tidal constants are estimated using the GESLA dataset (Woodworth et al 2014) following the approach presented in Piccioni et al (2019). This record is an add-on to the full TICON dataset (https://doi.org/10.1594/PANGAEA.896587), using exactly the same data format and pre-processing. These steps include using tide gauge data that contains at least ten years of continuous data. Further, the dataset is restricted to only contain open ocean tide gauges by limiting it to a mean surrounding depth of tide gauges to be deeper than 500 meters in a 2-degree radius and excluding stations not native to the ocean domain of the employed tidal model TiME. Duplicate and closely neighbouring tide gauges, found within a 0.2-degree radius, are also removed from the dataset. This resulted in the availability of the four tidal constants for 134 tide gauges. The results are stored in one tab-separated text/ASCII file with 13 columns: 1. Latitude of the tide gauge station 2. Longitude of the tide gauge station 3. Constituent name 4. Amplitude (in cm) 5. Phase (in degrees) 6. Standard deviation of the amplitude (in cm) 7. Standard deviation of the phase (in degrees) 8. Percentage of missing observations 9. Total number of observations analyzed 10. Length of the maximum temporal gap found in the time series in days 11. Date of the first observation 12. Date of the last observation 13. Code that corresponds to the original source of the record TICON is a useful and easy-to-handle data set for tide model validation and allows the users to select the records according to the different criteria most suitable for their purposes. The options span from the choice of a geographical region to the use of single constituents or time periods.

Description: These datasets provide sedimentological data partly at annual resolution and an age model for the lateglacial part of (1) the ICDP sediment core 5017-1-A retrieved from the deep northern Dead Sea basin in 2010/11, and (2) for the Masada outcrop located at the southwestern shore of the Dead Sea sampled in 2018. The here investigated two sediment sections cover the last glacial-interglacial transition (ca. 17-11.5 ka BP) in the hydroclimatically sensitive Levant, when the water level of Lake Lisan – the precursor of the Dead Sea – dropped dramatically from its glacial high-stand to the Holocene low levels. Here, we analyze the interval between the last two gypsum units – the Upper Gypsum Unit (UGU) and the Additional Gypsum Unit (AGU) – which were also used to correlate the two sites. In the ICDP core this section is located between ~101 and 88.5 m sediment depth below lake floor and at Masada it encompasses the uppermost ~3.8 m sediments of the Lisan Formation, which form the terminal deposit at this site. Due to the lake level decline, the complete transition into the Holocene is only recorded in the ICDP core, while sedimentation at Masada terminates earlier. The microfacies was investigated by continuous thin section microscopy, while additional macroscopic information is provided from over- and underlying sediment sections. A revised chronology using age modelling in OxCal (Ramsey 2008; Ramsey 2009; Ramsey and Lee 2013) was developed for the ICDP core and a floating varve chronology was constructed at Masada. Using these new microfacies data from marginal (Masada) and deep-water (ICDP core) sediments, the hydroclimatic variability during the final stage of Lake Lisan can be reconstructed, which could provide important insights into the development of human sedentism in the region at this time.

Description: We present sea ice temperature and salinity data from first-year ice (FYI) and second-year ice (SYI) relevant to the temporal development of sea ice permeability and brine drainage efficiency from the early growth phase in October 2019 to the onset of spring warming in May 2020. Our dataset was collected in the central Arctic Ocean during the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) Expedition in 2019 to 2020. MOSAiC was an international transpolar drift expedition in which the German icebreaker RV Polarstern anchored into an ice floe to gain new insights into Arctic climate over a full annual cycle. In October 2019, RV Polarstern moored to an ice floe in the Siberian sector of the Arctic at 85 degrees north and 137 degrees east to begin the drift towards the North Pole and the Fram Strait via the Transpolar Drift Stream. The data presented here were collected during the first three legs of the expedition, so all the coring activities took place on the same floe. The end dates of legs 1, 2, and 3 were 13 December, 24 February, and 4 June, respectively. The dataset contributed to a baseline study entitled, Deciphering the properties of different Arctic ice types during the growth phase of the MOSAiC floes: Implications for future studies. The study highlights downward directed gas pathways in FYI and SYI by inferring sea ice permeability and potential brine release from several time series of temperature and salinity measurements. The physical properties presented in this paper lay the foundation for subsequent analyses on actual gas contents measured in the ice cores, as well as air-ice and ice-ocean gas fluxes. Sea ice cores were collected with a Kovacs Mark II 9 cm diameter corer. To measure ice temperatures, about 4.5 cm deep holes were drilled into the core (intervals varied by site and leg) . The temperatures were measured by a digital thermometer within minutes after the cores were retrieved. The ice cores were placed into pre-labelled plastic sleeves sealed at the bottom end. The ice cores were transported to RV Polarstern and stored in a -20 degrees Celsius freezer. Each of the cores was sub-sampled, melted at room temperature, and processed for salinity within one or two days. The practical salinity was estimated by measuring the electrical conductivity and temperature of the melted samples using a WTW Cond 3151 salinometer equipped with a Tetra-Con 325 four-electrode conductivity cell. The practical salinity represents the the salinity estimated from the electrical conductivity of the solution. The dataset also contains derived variables, including sea ice density, brine volume fraction, and the Rayleigh number.

Description: This dataset provides annually resolved microfacies data from ICDP core 5017-1-A, retrieved from the deep northern Dead Sea basin in 2010/11, for the last glacial-interglacial transition (ca. 17-11.5 ka BP). Sediments of the Lisan Formation were investigated between ~101 and 88.5 m sediment depth below lake floor by continuous thin section microscopy, while additional macroscopic information is provided from core catchers, as well as from over- and underlying sediment sections. Thin sections were prepared following the standard procedure by Brauer and Casanova (2001) that was adjusted for salty sediments. Thin section analyses were performed on overlapping large-scale thin sections using a Zeiss Axiolab pol microscope at magnifications of 50-400x. Microfacies analyses included varve counting and measurements of varve and sublayer thickness.

Description: Despite the importance of surface energy budgets (SEBs) for land-climate interactions in the Arctic, uncertainties in their prediction persist. In situ observational data of SEB components - useful for research and model validation - are collected at relatively few sites across the terrestrial Arctic, and not all available datasets are readily interoperable. Furthermore, the terrestrial Arctic consists of a diversity of vegetation types, which are generally not well represented in land surface schemes of current Earth system models. This dataset contains metadata information about surface energy budget components measured at 64 tundra and glacier sites 〉60° N across the Arctic. This information was taken from the open-access repositories FLUXNET, Ameriflux, AON, GC-Net and PROMICE. The contained datasets are associated with the publication vegetation type as an important predictor of the Arctic Summer Land Surface Energy Budget by Oehri et al. 2022, and intended to support research of surface energy budgets and their relationship with environmental conditions, in particular vegetation characteristics across the terrestrial Arctic.

Description: This dataset provides the results from Bayesian age depth modelling in OxCal for ICDP core 5017-1-A, retrieved from the deep northern Dead Sea basin in 2010/11, for the last glacial-interglacial transition between ~101 and 88.5 m sediment depth below lake floor (ca. 17-11.5 ka BP). The model was performed in OxCal v.4.4 using a P_Sequence (1,1,C(-2,2)) (Ramsey 2008; Ramsey 2009; Ramsey and Lee 2013) and includes three tephrochronological ages from Neugebauer et al. (2021) and three radiocarbon ages from Kitagawa et al. (2017).

Description: Despite the importance of surface energy budgets (SEBs) for land-climate interactions in the Arctic, uncertainties in their prediction persist. In situ observational data of SEB components - useful for research and model validation - are collected at relatively few sites across the terrestrial Arctic, and not all available datasets are readily interoperable. Furthermore, the terrestrial Arctic consists of a diversity of vegetation types, which are generally not well represented in land surface schemes of current Earth system models. This dataset comprises harmonized, standardized and aggregated in-situ observations of surface energy budget components measured at 64 sites on vegetated and glaciated sites north of 60° latitude, in the time period from 1994 till 2021. The surface energy budget components include net radiation, sensible heat flux, latent heat flux, ground heat flux, net shortwave radiation, net longwave radiation, surface temperature and albedo, which were aggregated to daily mean, minimum and maximum values from hourly and half-hourly measurements. Data were retrieved from the monitoring networks FLUXNET, AmeriFlux, AON, GC-Net and PROMICE.

Description: We present sea ice temperature and salinity data from first-year ice (FYI) and second-year ice (SYI) relevant to the temporal development of sea ice permeability and brine drainage efficiency from the early growth phase in October 2019 to the onset of spring warming in May 2020. Our dataset was collected in the central Arctic Ocean during the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) Expedition in 2019 to 2020. MOSAiC was an international transpolar drift expedition in which the German icebreaker RV Polarstern anchored into an ice floe to gain new insights into Arctic climate over a full annual cycle. In October 2019, RV Polarstern moored to an ice floe in the Siberian sector of the Arctic at 85 degrees north and 137 degrees east to begin the drift towards the North Pole and the Fram Strait via the Transpolar Drift Stream. The data presented here were collected during the first three legs of the expedition, so all the coring activities took place on the same floe. The end dates of legs 1, 2, and 3 were 13 December, 24 February, and 4 June, respectively. The dataset contributed to a baseline study entitled, Deciphering the properties of different Arctic ice types during the growth phase of the MOSAiC floes: Implications for future studies. The study highlights downward directed gas pathways in FYI and SYI by inferring sea ice permeability and potential brine release from several time series of temperature and salinity measurements. The physical properties presented in this paper lay the foundation for subsequent analyses on actual gas contents measured in the ice cores, as well as air-ice and ice-ocean gas fluxes. Sea ice cores were collected with a Kovacs Mark II 9 cm diameter corer. To measure ice temperatures, about 4.5 cm deep holes were drilled into the core (intervals varied by site and leg) . The temperatures were measured by a digital thermometer within minutes after the cores were retrieved. The ice cores were placed into pre-labelled plastic sleeves sealed at the bottom end. The ice cores were transported to RV Polarstern and stored in a -20 degrees Celsius freezer. Each of the cores was sub-sampled, melted at room temperature, and processed for salinity within one or two days. The practical salinity was estimated by measuring the electrical conductivity and temperature of the melted samples using a WTW Cond 3151 salinometer equipped with a Tetra-Con 325 four-electrode conductivity cell. The practical salinity represents the the salinity estimated from the electrical conductivity of the solution. The dataset also contains derived variables, including sea ice density, brine volume fraction, and the Rayleigh number.

Description: We present sea ice temperature and salinity data from first-year ice (FYI) and second-year ice (SYI) relevant to the temporal development of sea ice permeability and brine drainage efficiency from the early growth phase in October 2019 to the onset of spring warming in May 2020. Our dataset was collected in the central Arctic Ocean during the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) Expedition in 2019 to 2020. MOSAiC was an international transpolar drift expedition in which the German icebreaker RV Polarstern anchored into an ice floe to gain new insights into Arctic climate over a full annual cycle. In October 2019, RV Polarstern moored to an ice floe in the Siberian sector of the Arctic at 85 degrees north and 137 degrees east to begin the drift towards the North Pole and the Fram Strait via the Transpolar Drift Stream. The data presented here were collected during the first three legs of the expedition, so all the coring activities took place on the same floe. The end dates of legs 1, 2, and 3 were 13 December, 24 February, and 4 June, respectively. The dataset contributed to a baseline study entitled, Deciphering the properties of different Arctic ice types during the growth phase of the MOSAiC floes: Implications for future studies. The study highlights downward directed gas pathways in FYI and SYI by inferring sea ice permeability and potential brine release from several time series of temperature and salinity measurements. The physical properties presented in this paper lay the foundation for subsequent analyses on actual gas contents measured in the ice cores, as well as air-ice and ice-ocean gas fluxes. Sea ice cores were collected with a Kovacs Mark II 9 cm diameter corer. To measure ice temperatures, about 4.5 cm deep holes were drilled into the core (intervals varied by site and leg) . The temperatures were measured by a digital thermometer within minutes after the cores were retrieved. The ice cores were placed into pre-labelled plastic sleeves sealed at the bottom end. The ice cores were transported to RV Polarstern and stored in a -20 degrees Celsius freezer. Each of the cores was sub-sampled, melted at room temperature, and processed for salinity within one or two days. The practical salinity was estimated by measuring the electrical conductivity and temperature of the melted samples using a WTW Cond 3151 salinometer equipped with a Tetra-Con 325 four-electrode conductivity cell. The practical salinity represents the the salinity estimated from the electrical conductivity of the solution. The dataset also contains derived variables, including sea ice density, brine volume fraction, and the Rayleigh number.

Description: This dataset describes two 17 m long sediment cores taken from beneath two thermokarst lakes in the Yukechi Alas, Central Yakutia, Russia. The first core was taken from below an Alas thermokarst lake (YU-L7; 61.76397°N, 130.46442°E) and the second core below and Yedoma lake (YU-L15; 61.76086°N, 130.47466°E). The dataset presents biogeochemical and biomarker parameters of sediment cores YU-L7 and YU-L15. Biogeochemical analyses include total carbon (TC) content, total organic carbon (TOC) content, total nitrogen (TN) content. Biomarker parameters include the n-alkane concentration, average chain length (ACL), carbon preference index (CPI), brGDGT concentration, archaeol concentration and the isoGDGT-0 concentration. The n-alkanes were measured in the aliphatic fraction by gas chromatography-mass spectromety using a Trace GC Ultra coupled to a DSQ MS. The branched and isoprenoid glycerol dialkyl glycerol tetraethers, as well as the dialkyl glycerol diether lipid (archaeol) were measured in the NSO fraction using a Shimadzu LC-10AD high-performance liquid chromatograph coupled to a Finnigan TSQ 7000 mass spectrometer via an atmospheric pressure chemical ionization interface. The pH soil is the sediment pH which was assessed by adding 6.12 mL of 0.01 M CaCl~2~ to ~2.5 g dried sediment and measuring with a Multilab 540 (WTW) at 20°C.

Description: We present sea ice temperature and salinity data from first-year ice (FYI) and second-year ice (SYI) relevant to the temporal development of sea ice permeability and brine drainage efficiency from the early growth phase in October 2019 to the onset of spring warming in May 2020. Our dataset was collected in the central Arctic Ocean during the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) Expedition in 2019 to 2020. MOSAiC was an international transpolar drift expedition in which the German icebreaker RV Polarstern anchored into an ice floe to gain new insights into Arctic climate over a full annual cycle. In October 2019, RV Polarstern moored to an ice floe in the Siberian sector of the Arctic at 85 degrees north and 137 degrees east to begin the drift towards the North Pole and the Fram Strait via the Transpolar Drift Stream. The data presented here were collected during the first three legs of the expedition, so all the coring activities took place on the same floe. The end dates of legs 1, 2, and 3 were 13 December, 24 February, and 4 June, respectively. The dataset contributed to a baseline study entitled, Deciphering the properties of different Arctic ice types during the growth phase of the MOSAiC floes: Implications for future studies. The study highlights downward directed gas pathways in FYI and SYI by inferring sea ice permeability and potential brine release from several time series of temperature and salinity measurements. The physical properties presented in this paper lay the foundation for subsequent analyses on actual gas contents measured in the ice cores, as well as air-ice and ice-ocean gas fluxes. Sea ice cores were collected with a Kovacs Mark II 9 cm diameter corer. To measure ice temperatures, about 4.5 cm deep holes were drilled into the core (intervals varied by site and leg) . The temperatures were measured by a digital thermometer within minutes after the cores were retrieved. The ice cores were placed into pre-labelled plastic sleeves sealed at the bottom end. The ice cores were transported to RV Polarstern and stored in a -20 degrees Celsius freezer. Each of the cores was sub-sampled, melted at room temperature, and processed for salinity within one or two days. The practical salinity was estimated by measuring the electrical conductivity and temperature of the melted samples using a WTW Cond 3151 salinometer equipped with a Tetra-Con 325 four-electrode conductivity cell. The practical salinity represents the the salinity estimated from the electrical conductivity of the solution. The dataset also contains derived variables, including sea ice density, brine volume fraction, and the Rayleigh number.

Description: Despite the importance of surface energy budgets (SEBs) for land-climate interactions in the Arctic, uncertainties in their prediction persist. In-situ observational data of SEB components - useful for research and model validation - are collected at relatively few sites across the terrestrial Arctic, and not all available datasets are readily interoperable. Furthermore, the terrestrial Arctic consists of a diversity of vegetation types, which are generally not well represented in land surface schemes of current Earth system models. Therefore, we here provide four datasets comprising: 1. Harmonized, standardized and aggregated in situ observations of SEB components at 64 vegetated and glaciated sites north of 60° latitude, in the time period 1994-2021 2. A description of all study sites and associated environmental conditions, including the vegetation types, which correspond to the classification of the Circumpolar Arctic Vegetation Map (CAVM, Raynolds et al. 2019). 3. Data generated in a literature synthesis from 358 study sites on vegetation or glacier (〉=60°N latitude) covered by 148 publications. 4. Metadata, including data contributor information and measurement heights of variables associated with Oehri et al. 2022.

Description: This dataset provides lithological data from ICDP core 5017-1-A, retrieved from the deep northern Dead Sea basin in 2010/11, for the last glacial-interglacial transition (ca. 17-11.5 ka BP). The microfacies of the Lisan Formation was investigated between ~101 and 88.5 m sediment depth below lake floor by continuous thin section microscopy, while additional macroscopic information is provided from core catchers, as well as from over- and underlying sediment sections. Thin sections were prepared following the standard procedure by Brauer and Casanova (2001) that was adjusted for salty sediments. Thin section analyses were performed on overlapping large-scale thin sections using a Zeiss Axiolab pol microscope at magnifications of 50-400x.

Description: This dataset contains observations of water discharge rates and concentrations of dissolved organic carbon (DOC) and dissolved inorganic carbon (DIC) from a polygonal tundra site in the Lena River Delta, Russia. This dataset also contains lateral carbon fluxes of DOC and DIC that were estimated from these observations. Additionally, this dataset contains vertical fluxes of carbon dioxide and methane from the same study site. All observations were recorded on Samoylov Island (N 72.377188, E 126.495144) in the summer of 2014. The abbreviations A1, A2 and B refer to three outflows on the island where the hydrological parameters were observed (A1: N 72.379991, E 126.480886; A2: N 72.380134, E 126.481433; B: N 72.381348, E 126.483482). All outflows were approximately 10 meters. More information can be found in https://doi.org/10.5194/bg-19-3863-2022.

Description: Despite the importance of surface energy budgets (SEBs) for land-climate interactions in the Arctic, uncertainties in their prediction persist. In situ observational data of SEB components - useful for research and model validation - are collected at relatively few sites across the terrestrial Arctic, and not all available datasets are readily interoperable. Furthermore, the terrestrial Arctic consists of a diversity of vegetation types, which are generally not well represented in land surface schemes of current Earth system models. This dataset describes the environmental conditions for 64 tundra and glacier sites (〉=60°N latitude) across the Arctic, for which in situ measurements of surface energy budget components were harmonized (see Oehri et al. 2022). These environmental conditions are (proxies of) potential drivers of SEB-components and could therefore be called SEB-drivers. The associated environmental conditions, include the vegetation types graminoid tundra, prostrate dwarf-shrub tundra, erect-shrub tundra, wetland complexes, barren complexes (≤ 40% horizontal plant cover), boreal peat bogs and glacier. These land surface types (apart from boreal peat bogs) correspond to the main classification units of the Circumpolar Arctic Vegetation Map (CAVM, Raynolds et al. 2019). For each site, additional climatic and biophysical variables are available, including cloud cover, snow cover duration, permafrost characteristics, climatic conditions and topographic conditions.

Description: This dataset contains over 30 marine Electrical Resistivity Tomography (ERT) profiles taken in September 2021 around Tuktoyaktuk Island (NWT / Beaufort Sea, Canada). The measurements were part of the “Mackenzie Delta Permafrost Field Campaign” (mCan2021) within the “Modular Observation solutions for Earth Systems” (MOSES) program. The collected profiles consist of numerous adjacent vertical soundings in a (quasi-symmetric) reciprocal Wenner-Schlumberger array, using a floating cable towed behind a boat. GPS records along the electrode streamer were taken, enabling the improvement of pre- processing by excluding measurements for which the cable was curved and electrode positions deviated too widely. The aim of the study was to determine the depth of the submarine permafrost. Cleaned data is provided in csv format.

Description: This collection contains permafrost related measurements in the Mackenzie Delta, NWT, Canada from the MOSES (Modular Observation Solutions for Earth Systems) field campaign in September 2021. The field campaign was focused on three subaquatic sites: a small thermokarst lake along the ITH just south of Trail Valley Creek, "Lake 3", an elongated lake with known methane occurence in the outer Mackenzie Delta, "Swiss Cheese Lake", and north and south of Tuktoyaktuk Island. At "Swiss Cheese Lake", we measured methane and CO2 concentrations in surface water and in the air above the lake, lake bed temperatures and detailed bathymetry. At "Lake 3" we measured active layer thickness on the lake banks, lake bed temperatures, and detailed bathymetry, as well as an ERT survey to estimate the talik depth below the lake. North and south of Tuktoyaktuk Island, we measured active layer thickness and sea bed temperatures and did an extensive ERT survey to obtain the depth of the subsea permafrost table. An additional passive seismic survey was carried out and the data is available at https://doi.org/10.5880/GIPP.202199.1.

Description: This dataset contains seven Electrical Resistivity Tomography (ERT) profiles taken in September 2021 at “Lake 3”, a thermokarst lake near the Inuvik-Tuktoyaktuk-Highway (ITH), about 50 km north of Inuvik (NWT, Canada). The measurements were part of the “Mackenzie Delta Permafrost Field Campaign” (mCan2021) within the “Modular Observation solutions for Earth Systems” (MOSES) program. The collected profiles consist of numerous adjacent vertical soundings in a (quasi-symmetric) reciprocal Wenner-Schlumberger array. In addition to surveys on the lake, using a floating cable towed behind a boat, two “amphibian” profiles were taken. Starting as purely terrestrial surveys using metal spike electrodes, the cable was then moved towards the lake with some of the electrodes floating on the water surface, and some still on land. The aim of the study was to determine permafrost properties on the land, to detect a possible talik beneath the lake and to especially be able to infer the transition between the two below the shoreline.

Description: Presented are analytical data from lacustrine sediment cores, retrieved from Lake Nam Co (Tibetan Plateau). The sediment core is a composite of one gravity core, taken with a Rumohr-Meischner gravity corer (63 mm diameter) and a piston core, retrieved using an uwitec piston coring system (http://www.uwitec.at; 90 mm diameter). The composite core labelled 〈NC 08/01〉 comprises a total length of 10.378 m. The cores were obtained at N 30.737417, E 090.790333 at a water depth of 93 m on 2008-09-15. The purpose of obtaining this sediment core was to establish a high-resolution record of climate (monsoonal) and environmental change using multiple proxy data. The dataset comprises analytical data based on sedimentological, inorganic geochemical, mineralogical and isotope-geochemical methods. Specifically: sediment water content & density; magnetic susceptibility; particel size data; quantitative inorganic geochemical data (ICP-OES aqua regia and HCL digestions); semi-quantitative XRF elemental data; carbon, nitrogen, sulfur contents; qualitative mineralogical data; bulk sediment stable carbon and oxygen isotope data.

Description: This data set is part of a larger data harmonization effort to make lake sediment core data machine readable and comparable. Here we standardized X-ray fluorescence line scanning (XRF)-based element data of sediment core EN18208, retrieved in 2018 from Lake Ilirney (Chukotka, Russia) at 10.76 m water depth. The glacial lake Ilirney is situated in the forest tundra mountain area and has one outflow, one main inflow and several smaller inflows. It lies at an elevation of ca. 428 m a.s.l. with a surface area of ca. 30 km2 and a maximum lake water depth of estimated 44 m. The 10.76 m sediment core was retrieved by a UWITEC piston corer during the RU-Land_2018_Chukotka expedition of the Alfred Wegener Institute Helmholtz Centre for Polar and Marine Research (AWI, Germany, Potsdam) in cooperation with the North Eastern Federal State University (NEFU, Russia, Yakutsk). The downcore elemental composition was measured using an AVAATECH x-ray fluorescence core scanner at Bundesanstalt für Geowissenschaften und Rohstoffe (BGR) in Berlin, Spandau.

Description: This data set is part of a larger data harmonization effort to make lake sediment core data machine readable and comparable. Here we standardized radiocarbon and OSL age data of sediment core EN18208, retrieved in 2018 from Lake Ilirney (Chukotka, Russia) at 10.76 m water depth. The glacial lake Ilirney is situated in the forest tundra mountain area and has one outflow, one main inflow and several smaller inflows. It lies at an elevation of ca. 428 m a.s.l. with a surface area of ca. 30 km2 and a maximum lake water depth of estimated 44 m. The 10.76 m sediment core was retrieved by a UWITEC piston corer during the RU-Land_2018_Chukotka expedition of the Alfred Wegener Institute Helmholtz Centre for Polar and Marine Research (AWI, Germany, Potsdam) in cooperation with the North Eastern Federal State University (NEFU, Russia, Yakutsk). Radiocarbon data have been analysed from bulk sediment samples in Bremerhaven at the MICADAS laboratory. Optically stimulated luminescence (OSL) dating was performed at the Royal Holloway Luminescence Laboratory using a Risø TL/OSL-DA-15 automated dating system.

Description: A 25-cm long predominantly aragonite stalagmite was collected November 2, 2005 from Dharamjali Cave (29.5°N, 80.2°E) in the central Himalayas. This dataset contains stable isotope, trace element, XRF, U/Th dating, and dripwater data. The age model spans 4.2 to 2.3 ka BP, and the dataset records seasonal shifts in hydroclimate from 4.2 to 3.1 ka BP. Using the DHAR-1A half of the speleothem, 750 samples were milled at 100–300 µm resolution for stable isotope analysis (δ18O and δ13C) and analyzed at GFZ Potsdam. Further high-resolution stable isotope analysis at the University of Cambridge included 876 samples from the bottom 4 cm of the mirroring slab DHAR-1B, covering c. 4.2–3.6 ka BP. The δ44/40Ca measurements were made on 60 aragonite samples of aragonite and 1 calcite sample milled between 4.2 and 2.8 ka BP. The elemental composition of DHAR-1B was determined first with an Avaatech XRF scanner at the University of Cambridge, and later using laser ablation inductively-coupled plasma mass spectrometry (LA-ICP-MS) at the University of Waikato. U-series dating was performed at Caltech on 22 samples. Twelve U-series ages (between 2.55 and 4.14 ka BP) were used to construct the age models, using ensembles of 2000 Monte Carlo simulations for each proxy using the MATLAB-based COPRA script (Breitenbach et al., 2012, https://doi.org/10.5194/cp-8-1765-2012).

Description: This data set is part of a larger data harmonization effort to make lake sediment core data machine readable and comparable. Here we standardized grain size element data of sediment core EN18208, retrieved in 2018 from Lake Ilirney (Chukotka, Russia) at 10.76 m water depth. The glacial lake Ilirney is situated in the forest tundra mountain area and has one outflow, one main inflow and several smaller inflows. It lies at an elevation of ca. 428 m a.s.l. with a surface area of ca. 30 km2 and a maximum lake water depth of estimated 44 m. The 10.76 m sediment core was retrieved by a UWITEC piston corer during the RU-Land_2018_Chukotka expedition of the Alfred Wegener Institute Helmholtz Centre for Polar and Marine Research (AWI, Germany, Potsdam) in cooperation with the North Eastern Federal State University (NEFU, Russia, Yakutsk). Grain-size was measured using a Malvern Mastersizer 3000 laser diffraction particle analyser.

Description: This data set is part of a larger data harmonization effort to make lake sediment core data machine readable and comparable. Here we standardized mineral data of sediment core EN18208, retrieved in 2018 from Lake Ilirney (Chukotka, Russia) at 10.76 m water depth. The glacial lake Ilirney is situated in the forest tundra mountain area and has one outflow, one main inflow and several smaller inflows. It lies at an elevation of ca. 428 m a.s.l. with a surface area of ca. 30 km2 and a maximum lake water depth of estimated 44 m. The 10.76 m sediment core was retrieved by a UWITEC piston corer during the RU-Land_2018_Chukotka expedition of the Alfred Wegener Institute Helmholtz Centre for Polar and Marine Research (AWI, Germany, Potsdam) in cooperation with the North Eastern Federal State University (NEFU, Russia, Yakutsk). Bulk mineralogy was analysed by (x-ray diffractometry (XRD) using a (PHILIPS, Netherlands) PW1820 goniometer.

Description: This data set is part of a larger data harmonization effort to make lake sediment core data machine readable and comparable. Here we standardized radiocarbon and OSL age data of sediment core EN18208, retrieved in 2018 from Lake Ilirney (Chukotka, Russia) at 10.76 m water depth. The glacial lake Ilirney is situated in the forest tundra mountain area and has one outflow, one main inflow and several smaller inflows. It lies at an elevation of ca. 428 m a.s.l. with a surface area of ca. 30 km2 and a maximum lake water depth of estimated 44 m. The 10.76 m sediment core was retrieved by a UWITEC piston corer during the RU-Land_2018_Chukotka expedition of the Alfred Wegener Institute Helmholtz Centre for Polar and Marine Research (AWI, Germany, Potsdam) in cooperation with the North Eastern Federal State University (NEFU, Russia, Yakutsk). Water content and organic matter was analysed at AWI Potsdam. Dried and milled samples were analysed using a Vario EL III carbon-nitrogen-sulphur analyser. Organic carbon content was determined using a Vario MAX C analyser.

Description: Surface air temperature measurements obtained from different sensors are used to construct a unique time series with one minute time-interval. Apart from differences in design and environmental exposition, periods of missing data also exist in the data series of each sensor. A primary data set was selected in terms of quality and temporal extension. A combination of two different techniques is applied to complete this data set: one is based on the autocorrelation of the series and the other on measurements taken from other sensors. The resulting values constitute a complete series of surface air temperature at AGGO.

Description: This dataset contain stable isotope values for water samples collected ~weekly from the Rio Bermejo at the Lavalle bridge (-25.6513, -60.1277) from March 2016 to February 2018. Water samples were filtered to 0.2 micron using a custom filtration device. We measured d2H and d18O on a Picarro L-2140i Cavity Ring-Down Spectrometer at the GFZ Potsdam. Measurements were made in duplicate, normalized to the Vienna Standard Mean Ocean Water (VSMOW), and analytical uncertainty is reported as one standard deviation from the mean. River discharge was measured at the El Colorado gauging station, which is ~100 km down slope from the sampling location.

Description: This dataset provides the geochemistry data for the Holocene sediment sequence retrieved from Lake Uddelermeer (The Netherlands) in 2012. Additionally, alkane concentrations for a set of modern leaf samples are provided. Concentrations of fossil alkanes, GDGTs as well as elemental (C, N, S, H) and compound-specific delta Deuterium measurements are presented against both depth (cm) and age (cal yr. BP). A total of 59 samples were analysed. Modern leaf alkane concentrations are presented as concentrations, 10 samples were analysed. The geochemical data provides information about regional vegetation change as well as changes in effective precipitation. It was produced to inform on the age and duration of major environmental transitions during the middle and late Holocene. Cores were retrieved from the lake using a 3-m long handheld piston corer deployed from a floating coring platform during field work in April and May 2012. Samples were obtained from splits of the core and processed in the laboratory of the University of Amsterdam (the Netherlands) using standard protocols (CNHS, alkane concentrations), the laboratory of Utrecht University (the Netherlands; GDGT concentrations) and at GFZ Potsdam (Germany; delta Deuterium). Name of the Campaign: UDD Event Label: UDD-E Method: Uwitec piston corer Latitude: 52.24652778 Longitude: 5.76097222 Elevation: 24m asl Date/Time of event: 2012-05-01T14:00:00 Further information about event: Lake sediment sequence retrieved using a 60 mm piston corer deployed from a floating platform.

Description: Water samples were filtered to 0.2 micron prior to measurement. Samples for cation analysis were acidified in the field to pH 〈 2 using 6N HNO3. Cation concentrations were measured with a Varian 720 inductively coupled plasma optical emission spectrometer (ICP-OES) at the GFZ Helmholtz Laboratory for the Geochemistry of the Earth Surface (HELGES), using SLRS-5 (Saint-Laurent River Surface, National Research Council - Conseil National de Recherches Canada) and USGS M212 and USGS T187 as external standards. We corrected for instrument drift by measuring an internal standard (GFZ-RW1) every 10 samples and we determined measurement uncertainty using calibration curve uncertainty. Anion concentrations were measured with a Dionex ICS1100 Ion Chromatograph, using USGS standards M206 and M212 as external standards for quality control, with uncertainty determined from triplicate analysis. We corrected cation concentrations for cyclic salt inputs following Bickle et al. (2005, doi:10.1016/j.gca.2004.11.019).

Description: This dataset reports measurements from a laboratory incubation of soils sourced from a boreal peatland and surrounding habitats (Siikaneva Bog, Finland). In August 2021, soil cores were collected from three habitat zones: a well-drained upland forest, an intermediate margin ecotone, and a Sphagnum moss bog. The cores from each habitat were taken from surface to approximately 50cm below surface using an Eijelkamp peat corer and subdivided by soil horizon. The samples were then incubated anaerobically for 140 days in three temperature treatment groups (0, 4, 20°C). Subsamples of the incubations headspace (250 µL) were measured on a gas chromatograph (7890A, Agilent Technologies, USA) with flame ionization detection (FID) for CO2 and CH4 concentrations. The rate of respiration from the samples were calculated per gram carbon and per gram soil as described in the method of Robertson., et al. (1999) and reported here, along with other relevant parameters.

Description: The Global Ocean Data Analysis Project (GLODAP) is a synthesis effort providing regular compilations of surface-to-bottom ocean biogeochemical bottle data, with an emphasis on seawater inorganic carbon chemistry and related variables determined through chemical analysis of seawater samples. GLODAPv2.2022 is an update of the previous version, GLODAPv2.2021 (Lauvset et al., 2021). The major changes are as follows: data from 96 new cruises were added, data coverage was extended until 2021, and for the first time we performed secondary quality control on all sulphur hexafluoride (SF6) data. In addition, a number of changes were made to data included in GLODAPv2.2021. These changes affect specifically the SF6 data, which are now subjected to secondary quality control, and carbon data measured onboard the RV Knorr in the Indian Ocean in 1994–1995 which are now adjusted using CRM measurements made at the time. GLODAPv2.2022 includes measurements from almost 1.4 million water samples from the global oceans collected on 1085 cruises. The data for the now 13 GLODAP core variables (salinity, oxygen, nitrate, silicate, phosphate, dissolved inorganic carbon, total alkalinity, pH, CFC-11, CFC-12, CFC-113, CCl4, and SF6) have undergone extensive quality control with a focus on systematic evaluation of bias. The data are available in two formats: (i) as submitted by the data originator but converted to World Ocean Circulation Experiment (WOCE) exchange format and (ii) as a merged data product with adjustments applied to minimize bias. For the present annual update, adjustments for the 96 new cruises were derived by comparing those data with the data from the 989 quality controlled cruises in the GLODAPv2.2021 data product using crossover analysis. SF6 data from all cruises were evaluated by comparison with CFC-12 data measured on the same cruises. For nutrients and ocean carbon dioxide (CO2) chemistry comparisons to estimates based on empirical algorithms provided additional context for adjustment decisions. The adjustments that we applied are intended to remove potential biases from errors related to measurement, calibration, and data handling practices without removing known or likely time trends or variations in the variables evaluated. The compiled and adjusted data product is believed to be consistent to better than 0.005 in salinity, 1 % in oxygen, 2 % in nitrate, 2 % in silicate, 2 % in phosphate, 4 μmol kg-1 in dissolved inorganic carbon, 4 μmol kg-1 in total alkalinity, 0.01–0.02 in pH (depending on region), and 5 % in the halogenated transient tracers. The other variables included in the compilation, such as isotopic tracers and discrete CO2 fugacity (fCO2), were not subjected to bias comparison or adjustments.

Description: Perturbations in stratospheric aerosol due to explosive volcanic eruptions are a primary contributor to natural climate variability. Observations of stratospheric aerosol are available for the past decades, and information from ice cores has been used to derive estimates of stratospheric sulfur injections and aerosol optical depth over the Holocene (approximately 10 000 BP to present) and into the last glacial period, extending back to 60 000 BP. Tephra records of past volcanism, compared to ice cores, are less complete but extend much further into the past. To support model studies of the potential impacts of explosive volcanism on climate variability across timescales, we present here an ensemble reconstruction of volcanic stratospheric sulfur injection (VSSI) over the last 140 000 years that is based primarily on terrestrial and marine tephra records. VSSI values are computed as a simple function of eruption magnitude based on VSSI estimates from ice cores and satellite observations for identified eruptions. To correct for the incompleteness of the tephra record, we include stochastically generated synthetic eruptions assuming a constant background eruption frequency from the ice core Holocene record. While the reconstruction often differs from ice core estimates for specific eruptions due to uncertainties in the data used and reconstruction method, it shows good agreement with an ice-core-based VSSI reconstruction in terms of millennial-scale cumulative VSSI variations over the Holocene. The PalVol reconstruction provides a new basis to test the contributions of forced vs. unforced natural variability to the spectrum of climate and the mechanisms leading to abrupt transitions in the palaeoclimate record with low- to high-complexity climate models. The PalVol volcanic forcing reconstruction is available at https://doi.org/10.26050/WDCC/PalVolv1 (Toohey and Schindlbeck-Belo, 2023).

Description: A new luminescence erosion meter has huge potential for inferring erosion rates on sub-millennial scales for both steady and transient states of erosion, which is not currently possible with any existing techniques capable of measuring erosion. This study applies new rock luminescence techniques to a well-constrained scenario provided by the Beinn Alligin rock avalanche, NW Scotland. Boulders in this deposit are lithologically consistent and have known cosmogenic nuclide ages and independently derived Holocene erosion rates. We find that luminescence-derived exposure ages for the Beinn Alligin rock avalanche were an order of magnitude younger than existing cosmogenic nuclide exposure ages, suggestive of high erosion rates (as supported by field evidence of quartz grain protrusions on the rock surfaces). Erosion rates determined by luminescence were consistent with independently derived rates measured from boulder edge roundness. Inversion modelling indicates a transient state of erosion reflecting the stochastic nature of erosional processes over the last ∼4.5 kyr in the wet, temperate climate of NW Scotland. Erosion was likely modulated by known fluctuations in moisture availability and to a lesser extent temperature, which controlled the extent of chemical weathering of these highly lithified rocks prior to erosion. The use of a multi-elevated temperature, post-infra-red, infra-red stimulated luminescence (MET-pIRIR) protocol (50, 150 and 225 ∘C) was advantageous as it identified samples with complexities that would not have been observed using only the standard infra-red stimulated luminescence (IRSL) signal measured at 50 ∘C, such as that introduced by within-sample variability (e.g. surficial coatings). This study demonstrates that the luminescence erosion meter can infer accurate erosion rates on sub-millennial scales and identify transient states of erosion (i.e. stochastic processes) in agreement with independently derived erosion rates for the same deposit.

Description: The photo-oxidation of myrcene, a monoterpene species emitted by plants, was investigated at atmospheric conditions in the outdoor simulation chamber SAPHIR (Simulation of Atmospheric PHotochemistry In a Large Reaction Chamber). The chemical structure of myrcene consists of one moiety that is a conjugated π system (similar to isoprene) and another moiety that is a triple-substituted olefinic unit (similar to 2-methyl-2-butene). Hydrogen shift reactions of organic peroxy radicals (RO2) formed in the reaction of isoprene with atmospheric OH radicals are known to be of importance for the regeneration of OH. Structure–activity relationships (SARs) suggest that similar hydrogen shift reactions like in isoprene may apply to the isoprenyl part of RO2 radicals formed during the OH oxidation of myrcene. In addition, SAR predicts further isomerization reactions that would be competitive with bimolecular RO2 reactions for chemical conditions that are typical for forested environments with low concentrations of nitric oxide. Assuming that OH peroxy radicals can rapidly interconvert by addition and elimination of O2 like in isoprene, bulk isomerization rate constants of 0.21 and 0.097 s−1 (T=298 K) for the three isomers resulting from the 3′-OH and 1-OH addition, respectively, can be derived from SAR. Measurements of radicals and trace gases in the experiments allowed us to calculate radical production and destruction rates, which are expected to be balanced. The largest discrepancies between production and destruction rates were found for RO2. Additional loss of organic peroxy radicals due to isomerization reactions could explain the observed discrepancies. The uncertainty of the total radical (ROx=OH+HO2+RO2) production rates was high due to the uncertainty in the yield of radicals from myrcene ozonolysis. However, results indicate that radical production can only be balanced if the reaction rate constant of the reaction between hydroperoxy (HO2) and RO2 radicals derived from myrcene is lower (0.9 to 1.6×10-11 cm3 s−1) than predicted by SAR. Another explanation of the discrepancies would be that a significant fraction of products (yield: 0.3 to 0.6) from these reactions include OH and HO2 radicals instead of radical-terminating organic peroxides. Experiments also allowed us to determine the yields of organic oxidation products acetone (yield: 0.45±0.08) and formaldehyde (yield: 0.35±0.08). Acetone and formaldehyde are produced from different oxidation pathways, so that yields of these compounds reflect the branching ratios of the initial OH addition to myrcene. Yields determined in the experiments are consistent with branching ratios expected from SAR. The yield of organic nitrate was determined from the gas-phase budget analysis of reactive oxidized nitrogen in the chamber, giving a value of 0.13±0.03. In addition, the reaction rate constant for myrcene + OH was determined from the measured myrcene concentration, yielding a value of (2.3±0.3)×10-10 cm3 s−1.

Description: McMurdo Sound sea ice can generally be partitioned into two regimes: (1) a stable fast-ice cover, forming south of approximately 77.6∘ S around March–April and then breaking out the following January–February, and (2) a more dynamic region north of 77.6∘ S that the McMurdo Sound and Ross Sea polynyas regularly impact. In 2019, a stable fast-ice cover formed unusually late due to repeated break-out events. We analyse the 2019 sea-ice conditions and relate them to a modified storm index (MSI), a proxy for southerly wind events. We find there is a strong correlation between the timing of break-out events and several unusually large MSI events.

Description: The source apportionment of aerosol iron (Fe), including natural and combustion Fe, is an important issue because aerosol Fe can enhance oceanic primary production in the surface ocean. Based on our previous finding that combustion Fe emitted by evaporation processes has Fe isotope ratios (δ56Fe) that are approximately 4 ‰ lower than those of natural Fe, this study aimed to distinguish aerosol Fe sources over the northwestern Pacific using two size-fractionated marine aerosols. The δ56Fe values of fine and coarse particles from the eastern or northern Pacific were found to be similar to each other, ranging from 0.0 ‰ to 0.4 ‰. Most of them were close to the crustal average, suggesting the dominance of natural Fe. On the other hand, particles from the direction of East Asia demonstrated lower δ56Fe values in fine particles (−0.5 ‰ to −2.2 ‰) than in coarse particles (on average −0.02 ± 0.12 ‰). The correlations between the δ56Fe values and the enrichment factors of lead and vanadium suggested that the low δ56Fe values obtained were due to the presence of combustion Fe. The δ56Fe values of the soluble component of fine particles in this region were lower than the total, indicating the preferential dissolution of combustion Fe. In addition, we found a negative correlation between the δ56Fe value and the fractional Fe solubility in air masses from the direction of East Asia. These results suggest that the presence of combustion Fe is an important factor in controlling the fractional Fe solubility in air masses from the direction of East Asia, whereas other factors are more important in the other areas. By assuming typical δ56Fe values for combustion and natural Fe, the contribution of combustion Fe to the total (acid-digested) Fe in aerosols was estimated to reach up to 50 % of fine and 21 % of bulk (coarse + fine) particles in air masses from the direction of East Asia, whereas its contribution was small in the other areas. The contribution of combustion Fe to the soluble Fe component estimated for one sample was approximately twice as large as the total, indicating the importance of combustion Fe as a soluble Fe source despite lower emissions than the natural. These isotope-based estimates were compared with those estimated using an atmospheric chemical transport model (IMPACT), in which the fractions of combustion Fe in fine particles, especially in air masses from the direction of East Asia, were consistent with each other. In contrast, the model estimated a relatively large contribution from combustion Fe in coarse particles, probably because of the different characteristics of combustion Fe that are included in the model calculation and the isotope-based estimation. This highlights the importance of observational data on δ56Fe for size-fractionated aerosols to scale the combustion Fe emission by the model. The average deposition fluxes of soluble Fe to the surface ocean were 1.4 and 2.9 nmol m−2 d−1 from combustion and natural aerosols, respectively, in air masses from the direction of East Asia, which suggests that combustion Fe could be an important Fe source to the surface seawater among other Fe sources. Distinguishing Fe sources using the δ56Fe values of marine aerosols and seawater is anticipated to lead to a more quantitative understanding of the Fe cycle in the atmosphere and surface ocean.

Description: Icequakes are the result of processes occurring within the ice mass or between the ice and its environment. Studying icequakes provides a unique view on ice dynamics, specifically on the basal conditions. Changes in conditions due to environmental or climate changes are reflected in icequakes. Counting and characterizing icequakes is thus essential to monitor them. Most of the icequakes recorded by the seismic station at the Belgian Princess Elisabeth Antarctica Station (PE) have small amplitudes corresponding to maximal displacements of a few nanometres. Their detection threshold is highly variable because of the rapid and strong changes in the local seismic noise level. Therefore, we evaluated the influence of katabatic winds on the noise measured by the well-protected PE surface seismometer. Our purpose is to identify whether the lack of icequake detection during some periods could be associated with variations in the processes generating them or simply with a stronger seismic noise linked to stronger wind conditions. We observed that the wind mainly influences seismic noise at frequencies greater than 1 Hz. The seismic noise power exhibits a bilinear correlation with the wind velocity, with two different slopes at a wind velocity lower and greater than 6 m s−1 and with, for example at a period of 0.26 s, a respective variation of 0.4 dB (m −1 s) and 1.4 dB (m −1 s). These results allowed a synthetic frequency and wind-speed-dependent noise model to be presented that explains the behaviour of the wind-induced seismic noise at PE, which shows that seismic noise amplitude increases exponentially with increasing wind speed. This model enables us to study the influence of the wind on the original seismic dataset, which improves the observation of cryoseismic activity near the PE station.

Description: The importance of Antarctic sea ice and Southern Ocean warming has come into the focus of polar research during the last couple of decades. Especially around West Antarctica, where warm water masses approach the continent and where sea ice has declined, the distribution and evolution of sea ice play a critical role in the stability of nearby ice shelves. Organic geochemical analyses of marine seafloor surface sediments from the Antarctic continental margin allow an evaluation of the applicability of biomarker-based sea-ice and ocean temperature reconstructions in these climate-sensitive areas. We analysed highly branched isoprenoids (HBIs), such as the sea-ice proxy IPSO25 and phytoplankton-derived HBI-trienes, as well as phytosterols and isoprenoidal glycerol dialkyl glycerol tetraethers (GDGTs), which are established tools for the assessment of primary productivity and ocean temperatures respectively. The combination of IPSO25 with a phytoplankton marker (i.e. the PIPSO25 index) permits semi-quantitative sea-ice reconstructions and avoids misleading over- or underestimations of sea-ice cover. Comparisons of the PIPSO25-based sea-ice distribution patterns and TEX86L- and RI-OH′-derived ocean temperatures with (1) sea-ice concentrations obtained from satellite observations and (2) instrument measurements of sea surface and subsurface temperatures corroborate the general capability of these proxies to determine oceanic key variables properly. This is further supported by model data. We also highlight specific aspects and limitations that need to be taken into account for the interpretation of such biomarker data and discuss the potential of IPSO25 as an indicator for the former occurrence of platelet ice and/or the export of ice-shelf water.

Description: Aeromagnetic exploration is an important method of geophysical exploration. We study the compensation method of the towed bird system and establish the towed bird interference model. Due to the geomagnetic gradient changing greatly, the geomagnetic gradient is considered in the towed bird interference model. In this paper, we model the geomagnetic field gradient and analyze the influence of the towed bird system on the aeromagnetic compensation results. Finally, we apply the ridge regression method to solve the problem. We verify the feasibility of this compensation method through actual flight tests and further improve the data quality of the towed bird interference.

Description: The earthquake early warning systems (EEWSs) in China have achieved great progress, with warning alerts being successfully delivered to the public in some regions. We examined the performance of the EEWS in China's Sichuan Province during the 2019 Changning earthquake. Although its technical effectiveness was tested with the first alert released 10 s after the earthquake, we found that a big gap existed between the EEWS's message and the public's response. We highlight the importance of EEWS alert effectiveness and public participation for long-term resiliency, such as delivering useful alert messages through appropriate communication channels and training people to understand and properly respond.

Description: Forest fires modify soil organic carbon and suppress soil respiration for many decades after the initial disturbance. The associated changes in soil autotrophic and heterotrophic respiration from the time of the forest fire, however, are less well characterized. The FireResp model predicts soil autotrophic and heterotrophic respiration parameterized with a novel dataset across a fire chronosequence in the Yukon and Northwest Territories of Canada. The dataset consisted of soil incubation experiments and field measurements of soil respiration and soil carbon stocks. The FireResp model contains submodels that consider a Q10 (exponential) model of respiration compared to models of heterotrophic respiration using Michaelis–Menten kinetics parameterized with soil microbial carbon. For model evaluation we applied the Akaike information criterion and compared predicted patterns in components of soil respiration across the chronosequence. Parameters estimated with data from the 5 cm soil depth had better model–data comparisons than parameters estimated with data from the 10 cm soil depth. The model–data fit was improved by including parameters estimated from soil incubation experiments. Models that incorporated microbial carbon with Michaelis–Menten kinetics reproduced patterns in autotrophic and heterotrophic soil respiration components across the chronosequence. Autotrophic respiration was associated with aboveground tree biomass at more recently burned sites, but this association was less robust at older sites in the chronosequence. Our results provide support for more structured soil respiration models than standard Q10 exponential models.

Description: In this paper, we present a new version of the chemistry–climate model SOCOL-AERv2 supplemented by an iodine chemistry module. We perform three 20-year ensemble experiments to assess the validity of the modeled iodine and to quantify the effects of iodine on ozone. The iodine distributions obtained with SOCOL-AERv2-I agree well with AMAX-DOAS observations and with CAM-chem model simulations. For the present-day atmosphere, the model suggests that the iodine-induced chemistry leads to a 3 %–4 % reduction in the ozone column, which is greatest at high latitudes. The model indicates the strongest influence of iodine in the lower stratosphere with 30 ppbv less ozone at low latitudes and up to 100 ppbv less at high latitudes. In the troposphere, the account of the iodine chemistry reduces the tropospheric ozone concentration by 5 %–10 % depending on geographical location. In the lower troposphere, 75 % of the modeled ozone reduction originates from inorganic sources of iodine, 25 % from organic sources of iodine. At 50 hPa, the results show that the impacts of iodine from both sources are comparable. Finally, we determine the sensitivity of ozone to iodine by applying a 2-fold increase in iodine emissions, as it might be representative for iodine by the end of this century. This reduces the ozone column globally by an additional 1.5 %–2.5 %. Our results demonstrate the sensitivity of atmospheric ozone to iodine chemistry for present and future conditions, but uncertainties remain high due to the paucity of observational data of iodine species.

Description: The Arabian Sea (AS) hosts one of the most intense oxygen minimum zones (OMZs) in the world. Observations suggest a decline in O2 in the northern AS over the recent decades accompanied by an intensification of the suboxic conditions there. Over the same period, the local sea surface temperature has risen significantly, particularly over the Arabian Gulf (also known as Persian Gulf, hereafter the Gulf), while summer monsoon winds may have intensified. Here, we simulate the evolution of dissolved oxygen in the AS from 1982 through 2010 and explore its controlling factors, with a focus on changing atmospheric conditions. To this end, we use a set of eddy-resolving hindcast simulations forced with winds and heat and freshwater fluxes from an atmospheric reanalysis. We find a significant deoxygenation in the northern AS, with O2 inventories north of 20∘ N dropping by over 6 % per decade between 100 and 1000 m. These changes cause an expansion of the OMZ volume north of 20∘ N at a rate of 0.6 % per decade as well as an increase in the volume of suboxia and the rate of denitrification by 14 and 15 % per decade, respectively. We also show that strong interannual and decadal variability modulate dissolved oxygen in the northern AS, with most of the O2 decline taking place in the 1980s and 1990s. Using a set of sensitivity simulations we demonstrate that deoxygenation in the northern AS is essentially caused by reduced ventilation induced by the recent fast warming of the sea surface, including in the Gulf, with a contribution from concomitant summer monsoon wind intensification. This is because, on the one hand, surface warming enhances vertical stratification and increases Gulf water buoyancy, thus inhibiting vertical mixing and ventilation of the thermocline. On the other hand, summer monsoon wind intensification causes a rise in the thermocline depth in the northern AS that lowers O2 levels in the upper ocean. Our findings confirm that the AS OMZ is strongly sensitive to upper-ocean warming and concurrent changes in the Indian monsoon winds. Finally, our results also demonstrate that changes in the local climatic forcing play a key role in regional dissolved oxygen changes and hence need to be properly represented in global models to reduce uncertainties in future projections of deoxygenation.

Description: Authorship conflicts are a common occurrence in academic publishing, and they can have serious implications for the careers and well-being of the involved researchers as well as the collective success of research organizations. In addition to not inviting relevant contributors to co-author a paper, the order of authors as well as honorary, gift, and ghost authors are all widely recognized problems related to authorship. Unfair authorship practices disproportionately affect those lower in the power hierarchies – early career researchers, women, researchers from the Global South, and other minoritized groups. Here we propose an approach to preparing author lists based on clear, transparent, and timely communication. This approach aims to minimize the potential for late-stage authorship conflicts during manuscript preparation by facilitating timely and transparent decisions on potential co-authors and their responsibilities. Furthermore, our approach can help avoid imbalances between contributions and credits in published papers by recording planned and executed responsibilities. We present authorship guidelines which also include a novel authorship form along with the documentation of the formulation process for a multidisciplinary and interdisciplinary center with more than 250 researchers. Other research groups, departments, and centers can use or build on this template to design their own authorship guidelines as a practical way to promote fair authorship practices.

Description: This work presents the integration of a gas-phase and particulate atmospheric emission inventory (AEI) for Argentina in high spatial resolution (0.025∘×0.025∘; approx. 2.5 km×2.5 km) considering monthly variability from 1995 to 2020. The new inventory, called GEAA-AEIv3.0M, includes the following activities: energy production, fugitive emissions from oil and gas production, industrial fuel consumption and production, transport (road, maritime, and air), agriculture, livestock production, manufacturing, residential, commercial, and biomass and agricultural waste burning. The following species, grouped by atmospheric reactivity, are considered: (i) greenhouse gases (GHGs) – CO2, CH4, and N2O; (ii) ozone precursors – CO, NOx (NO+NO2), and non-methane volatile organic compounds (NMVOCs); (iii) acidifying gases – NH3 and SO2; and (iv) particulate matter (PM) – PM10, PM2.5, total suspended particles (TSPs), and black carbon (BC). The main objective of the GEAA-AEIv3.0M high-resolution emission inventory is to provide temporally resolved emission maps to support air quality and climate modeling oriented to evaluate pollutant mitigation strategies by local governments. This is of major concern, especially in countries where air quality monitoring networks are scarce, and the development of regional and seasonal emissions inventories would result in remarkable improvements in the time and space chemical prediction achieved by air quality models. Despite distinguishing among different sectoral and activity databases as well as introducing a novel spatial distribution approach based on census radii, our high-resolution GEAA-AEIv3.0M shows equivalent national-wide total emissions compared to the Third National Communication of Argentina (TNCA), which compiles annual GHG emissions from 1990 through 2014 (agreement within ±7.5 %). However, the GEAA-AEIv3.0M includes acidifying gases and PM species not considered in TNCA. Temporal comparisons were also performed against two international databases: Community Emissions Data System (CEDS) and EDGAR HTAPv5.0 for several pollutants; for EDGAR it also includes a spatial comparison. The agreement was acceptable within less than 30 % for most of the pollutants and activities, although a 〉90 % discrepancy was obtained for methane from fuel production and fugitive emissions and 〉120 % for biomass burning. Finally, the updated seasonal series clearly showed the pollution reduction due to the COVID-19 lockdown during the first quarter of year 2020 with respect to same months in previous years. Through an open-access data repository, we present the GEAA-AEIv3.0M inventory as the largest and more detailed spatial resolution dataset for the Argentine Republic, which includes monthly gridded emissions for 12 species and 15 stors between 1995 and 2020. The datasets are available at https://doi.org/10.17632/d6xrhpmzdp.2 (Puliafito et al., 2021), under a CC-BY 4 license.

Description: This article summarises the results of an analysis of solar radio bursts (SRBs) detected by the Compound Astronomical Low-cost Low-frequency Instrument for Spectroscopy and Transportable Observatory (CALLISTO) spectrometer hosted by the University of Rwanda. The data analysed were detected during the first year (2014–2015) of the instrument operation. Using quick plots provided by the e-CALLISTO website, a total of 201 intense and well-separated solar radio bursts detected by the CALLISTO station located in Rwanda, are found consisting of 4 type II, 175 type III and 22 type IV radio bursts. It is found that all analysed type II and ∼ 37 % of type III bursts are associated with impulsive solar flares, while the minority (∼ 13 %) of type IV radio bursts are associated with solar flares. Furthermore, all type II radio bursts are associated with coronal mass ejections (CMEs), ∼ 44 % of type III bursts are associated with CMEs, and the majority (∼ 82 %) of type IV bursts were accompanied by CMEs. With aid of the atmospheric imaging assembly (AIA) images on board the Solar Dynamics Observatory (SDO), the location of open magnetic field lines of non-flare-associated type III radio bursts are shown. The same images are used to show the magnetic loops in the solar corona for type IV radio bursts observed in the absence of solar flares and/or CMEs. Findings from this study indicate that analysis of SRBs that are observed from the ground can provide a significant contribution to the early diagnosis of solar transients phenomena, such as solar flares and CMEs, which are major drivers of potential space weather hazards.

Description: The Arctic is exposed to even faster temperature changes than most other areas on Earth. Constantly increasing temperature will lead to thawing permafrost and changes in the methane (CH4) emissions from wetlands. One of the places exposed to those changes is the Abisko–Stordalen Mire in northern Sweden, where climate and vegetation studies have been conducted since the 1970s. In our study, we analyzed field-scale methane emissions measured by the eddy covariance method at Abisko–Stordalen Mire for 3 years (2014–2016). The site is a subarctic mire mosaic of palsas, thawing palsas, fully thawed fens, and open water bodies. A bimodal wind pattern prevalent at the site provides an ideal opportunity to measure mire patches with different permafrost status with one flux measurement system. The flux footprint for westerly winds was dominated by elevated palsa plateaus, while the footprint was almost equally distributed between palsas and thawing bog-like areas for easterly winds. As these patches are exposed to the same climatic and weather conditions, we analyzed the differences in the responses of their methane emission for environmental parameters. The methane fluxes followed a similar annual cycle over the 3 study years, with a gentle rise during spring and a decrease during autumn, without emission bursts at either end of the ice-free season. The peak emission during the ice-free season differed significantly for the two mire areas with different permafrost status: the palsa mire emitted 19 mg-C m−2 d−1 and the thawing wet sector 40 mg-C m−2 d−1. Factors controlling the methane emission were analyzed using generalized linear models. The main driver for methane fluxes was peat temperature for both wind sectors. Soil water content above the water table emerged as an explanatory variable for the 3 years for western sectors and the year 2016 in the eastern sector. The water table level showed a significant correlation with methane emission for the year 2016 as well. Gross primary production, however, did not show a significant correlation with methane emissions. Annual methane emissions were estimated based on four different gap-filing methods. The different methods generally resulted in very similar annual emissions. The mean annual emission based on all models was 3.1 ± 0.3 g-C m−2 a−1 for the western sector and 5.5 ± 0.5 g-C m−2 a−1 for the eastern sector. The average annual emissions, derived from these data and a footprint climatology, were 2.7 ± 0.5 and 8.2 ± 1.5 g-C m−2 a−1 for the palsa and thawing surfaces, respectively. Winter fluxes were relatively high, contributing 27 %–45 % to the annual emissions.

Description: Two feature-based verification methods, thus far only used for the diagnostic evaluation of atmospheric models, have been applied to compare ∼7 km resolution pre-operational analyses of chlorophyll-a (Chl-a) concentrations to a 1 km gridded satellite-derived Chl-a concentration product. The aim of this study was to assess the value of applying such methods to ocean models. Chl-a bloom objects were identified in both data sets for the 2019 bloom season (1 March to 31 July). These bloom objects were analysed as discrete (2-D) spatial features, but also as space–time (3-D) features, providing the means of defining the onset, duration and demise of distinct bloom episodes and the season as a whole. The new feature-based verification methods help reveal that the model analyses are not able to represent small coastal bloom objects, given the coarser definition of the coastline, also wrongly producing more bloom objects in deeper Atlantic waters. Model analyses' concentrations are somewhat higher overall. The bias manifests itself in the size of the model analysis bloom objects, which tend to be larger than the satellite-derived bloom objects. The onset of the bloom season is delayed by 26 d in the model analyses, but the season also persists for another month beyond the diagnosed end. The season was diagnosed to be 119 d long in the model analyses, compared to 117 d from the satellite product. Geographically, the model analyses and satellite-derived bloom objects do not necessarily exist in a specific location at the same time and only overlap occasionally.

Description: The observing system design of multidisciplinary field measurements involves a variety of considerations on logistics, safety, and science objectives. Typically, this is done based on investigator intuition and designs of prior field measurements. However, there is potential for considerable increases in efficiency, safety, and scientific success by integrating numerical simulations in the design process. Here, we present a novel numerical simulation–environmental response function (NS–ERF) approach to observing system simulation experiments that aids surface–atmosphere synthesis at the interface of mesoscale and microscale meteorology. In a case study we demonstrate application of the NS–ERF approach to optimize the Chequamegon Heterogeneous Ecosystem Energy-balance Study Enabled by a High-density Extensive Array of Detectors 2019 (CHEESEHEAD19). During CHEESEHEAD19 pre-field simulation experiments, we considered the placement of 20 eddy covariance flux towers, operations for 72 h of low-altitude flux aircraft measurements, and integration of various remote sensing data products. A 2 h high-resolution large eddy simulation created a cloud-free virtual atmosphere for surface and meteorological conditions characteristic of the field campaign domain and period. To explore two specific design hypotheses we super-sampled this virtual atmosphere as observed by 13 different yet simultaneous observing system designs consisting of virtual ground, airborne, and satellite observations. We then analyzed these virtual observations through ERFs to yield an optimal aircraft flight strategy for augmenting a stratified random flux tower network in combination with satellite retrievals. We demonstrate how the novel NS–ERF approach doubled CHEESEHEAD19's potential to explore energy balance closure and spatial patterning science objectives while substantially simplifying logistics. Owing to its modular extensibility, NS–ERF lends itself to optimizing observing system designs also for natural climate solutions, emission inventory validation, urban air quality, industry leak detection, and multi-species applications, among other use cases.

Description: The global distribution of cropping intensity (CI) is essential to our understanding of agricultural land use management on Earth. Optical remote sensing has revolutionized our ability to map CI over large areas in a repeated and cost-efficient manner. Previous studies have mainly focused on investigating the spatiotemporal patterns of CI ranging from regions to the entire globe with the use of coarse-resolution data, which are inadequate for characterizing farming practices within heterogeneous landscapes. To fill this knowledge gap, in this study, we utilized multiple satellite data to develop a global, spatially continuous CI map dataset at 30 m resolution (GCI30). Accuracy assessments indicated that GCI30 exhibited high agreement with visually interpreted validation samples and in situ observations from the PhenoCam network. We carried out both statistical and spatial comparisons of GCI30 with six existing global CI estimates. Based on GCI30, we estimated that the global average annual CI during 2016–2018 was 1.05, which is close to the mean (1.09) and median (1.07) CI values of the existing six global CI estimates, although the spatial resolution and temporal coverage vary significantly among products. A spatial comparison with two satellite-based land surface phenology products further suggested that GCI30 was not only capable of capturing the overall pattern of global CI but also provided many spatial details. GCI30 indicated that single cropping was the primary agricultural system on Earth, accounting for 81.57 % (12.28×106 km2) of the world's cropland extent. Multiple-cropping systems, on the other hand, were commonly observed in South America and Asia. We found large variations across countries and agroecological zones, reflecting the joint control of natural and anthropogenic drivers on regulating cropping practices. As the first global-coverage, fine-resolution CI product, GCI30 is expected to fill the data gap for promoting sustainable agriculture by depicting worldwide diversity of agricultural land use intensity. The GCI30 dataset is available on Harvard Dataverse: https://doi.org/10.7910/DVN/86M4PO (Zhang et al., 2020).

Description: Improvements in our capability to reconstruct ancient surface-ocean conditions based on organic-walled dinoflagellate cyst (dinocyst) assemblages from the Southern Ocean provide an opportunity to better establish past position, strength and oceanography of the subtropical front (STF). Here, we aim to reconstruct the late Eocene to early Miocene (37–20 Ma) depositional and palaeoceanographic history of the STF in the context of the evolving Tasmanian Gateway as well as the potential influence of Antarctic circumpolar flow and intense waxing and waning of ice. We approach this by combining information from seismic lines (revisiting existing data and generating new marine palynological data from Ocean Drilling Program (ODP) Hole 1168A) in the western Tasmanian continental slope. We apply improved taxonomic insights and palaeoecological models to reconstruct the sea surface palaeoenvironmental evolution. Late Eocene–early Oligocene (37–30.5 Ma) assemblages show a progressive transition from dominant terrestrial palynomorphs and inner-neritic dinocyst taxa as well as cysts produced by heterotrophic dinoflagellates to predominantly outer-neritic/oceanic autotrophic taxa. This transition reflects the progressive deepening of the western Tasmanian continental margin, an interpretation supported by our new seismic investigations. The dominance of autotrophic species like Spiniferites spp. and Operculodinium spp. reflects relatively oligotrophic conditions, like those of regions north of the modern-day STF. The increased abundance in the earliest Miocene of Nematosphaeropsis labyrinthus, typical for modern subantarctic zone (frontal) conditions, indicates a cooling and/or closer proximity of the STF to the site . The absence of major shifts in dinocyst assemblages contrasts with other records in the region and suggests that small changes in surface oceanographic conditions occurred during the Oligocene. Despite the relatively southerly (63–55∘ S) location of Site 1168, the rather stable oceanographic conditions reflect the continued influence of the proto-Leeuwin Current along the southern Australian coast as Australia continued to drift northward. The relatively “warm” dinocyst assemblages at ODP Site 1168, compared with the cold assemblages at Antarctic Integrated Ocean Drilling Program (IODP) Site U1356, testify to the establishment of a pronounced latitudinal temperature gradient in the Oligocene Southern Ocean.

Description: The mass of the Greenland ice sheet is declining as mass gain from snow accumulation is exceeded by mass loss from surface meltwater runoff, marine-terminating glacier calving and submarine melting, and basal melting. Here we use the input–output (IO) method to estimate mass change from 1840 through next week. Surface mass balance (SMB) gains and losses come from a semi-empirical SMB model from 1840 through 1985 and three regional climate models (RCMs; HIRHAM/HARMONIE, Modèle Atmosphérique Régional – MAR, and RACMO – Regional Atmospheric Climate MOdel) from 1986 through next week. Additional non-SMB losses come from a marine-terminating glacier ice discharge product and a basal mass balance model. From these products we provide an annual estimate of Greenland ice sheet mass balance from 1840 through 1985 and a daily estimate at sector and region scale from 1986 through next week. This product updates daily and is the first IO product to include the basal mass balance which is a source of an additional ∼24 Gt yr−1 of mass loss. Our results demonstrate an accelerating ice-sheet-scale mass loss and general agreement (coefficient of determination, r2, ranges from 0.62 to 0.94) among six other products, including gravitational, volume, and other IO mass balance estimates. Results from this study are available at https://doi.org/10.22008/FK2/OHI23Z (Mankoff et al., 2021).

Description: China has experienced dramatic changes in emissions since 2010, which accelerated following the implementation of the Clean Air Action program in 2013. These changes have resulted in significant air quality improvements that are reflected in observations from both surface networks and satellite observations. Air pollutants, such as PM2.5, surface ozone, and their precursors, have long enough lifetimes in the troposphere to be easily transported downwind. Emission changes in China will thus not only change the domestic air quality but will also affect the air quality in other regions. In this study, we use a global chemistry transport model (CAM-chem) to simulate the influence of Chinese emission changes from 2010 to 2017 on both domestic and foreign air quality. We then quantify the changes in air-pollution-associated (including both PM2.5 and O3) premature mortality burdens at regional and global scales. Within our simulation period, the population-weighted annual PM2.5 concentration in China peaks in 2011 (94.1 µg m−3) and decreases to 69.8 µg m−3 by 2017. These estimated national PM2.5 concentration changes in China are comparable with previous studies using fine-resolution regional models, though our model tends to overestimate PM2.5 from 2013 to 2017 when evaluated with surface observations. Relative to 2010, emission changes in China increased the global PM2.5-associated premature mortality burdens through 2013, among which a majority of the changes (∼ 93 %) occurred in China. The sharp emission decreases after 2013 generated significant benefits for human health. By 2017, emission changes in China reduced premature deaths associated with PM2.5 by 108 800 (92 800–124 800) deaths per year globally, relative to 2010, among which 92 % were realized in China. In contrast, the population-weighted, annually averaged maximum daily 8 h ozone concentration peaked in 2014 and did not reach 2010 levels by 2017. As such, O3 generated nearly 8500 (6500–9900) more premature deaths per year in 2017 compared to 2010. Downwind regions, such as South Korea, Japan, and the United States, generally experienced O3 improvements following 2013 due to the decreased export of ozone and its precursors. Overall, we conclude that the sharp emission reductions in China over the past decade have generated substantial benefits for air quality that have reduced premature deaths associated with air pollution at a global scale.

Description: Recent observational and modeling studies show that variations of stratospheric ozone and the resulting interaction between ozone and the stratospheric circulation play an important role in surface weather and climate. However, in many cases, computationally expensive coupled chemistry models have been used to study these effects. Here, we demonstrate how a much simpler idealized general circulation model (GCM) can be used for studying the impact of interactive stratospheric ozone on the circulation. The model, named Simplified Chemistry-Dynamical Model (SCDM V1.0), is constructed from a preexisting idealized GCM, into which a simplified linear ozone scheme and a parameterization for the shortwave radiative effects of ozone are implemented. The distribution and variability of stratospheric ozone simulated by the new model are in good agreement with the MERRA2 reanalysis, even for extreme circulation events such as Arctic stratospheric sudden warmings. The model thus represents a promising new tool for the study of ozone–circulation interaction in the stratosphere and its associated effects on tropospheric weather and climate.

Description: Frequently occurring mega-droughts under current global climate change have attracted broad social attention. A paleoclimatic perspective is needed to increase our understanding of the causes and effects of droughts. South-western (SW) China has been threatened by severe seasonal droughts. Our current knowledge of millennial-scale dry and wet phases in this region is primarily based on the variability of the Indian summer monsoon. However, water availability over land does not always follow patterns of monsoonal precipitation but also depends on water loss from evaporation and transpiration. Here, we reconstructed precipitation intensity, lake hydrological balance and the soil water stress index (SWSI) for the last 27 000 years. Grain size, geochemical and pollen records from Yilong Lake reveal the long-term relationships and inconsistencies of dry–wet patterns in meteorological, hydrological and soil systems in the central Yunnan region, SW China. Our results show that the long-term trends among precipitation, hydrological balance and soil moisture varied through time. The hydrological balance and soil moisture were primarily controlled by temperature-induced evaporation change during periods of low precipitation such as the Last Glacial Maximum and Younger Dryas. During periods of high precipitation (the early to late Holocene), intensified evaporation from the lake surface offset the effects of increased precipitation on the hydrological balance. However, abundant rainfall and the dense vegetation canopy circumvented a soil moisture deficit that might have resulted from rising temperature. In conclusion, the hydrological balance in the central Yunnan region was more sensitive to temperature change while soil moisture could be further regulated by vegetation changes over millennial timescales. Therefore, under future climate warming, the surface water shortage in the central Yunnan region may become even more serious. Our study suggests that reforestation efforts may provide some relief to soil moisture deficits in this region.

Description: This paper presents a new technique to derive thermospheric temperature from space-based disk observations of far ultraviolet airglow. The technique, guided by findings from principal component analysis of synthetic daytime Lyman–Birge–Hopfield (LBH) disk emissions, uses a ratio of the emissions in two spectral channels that together span the LBH (2,0) band to determine the change in band shape with respect to a change in the rotational temperature of N2. The two-channel-ratio approach limits representativeness and measurement error by only requiring measurement of the relative magnitudes between two spectral channels and not radiometrically calibrated intensities, simplifying the forward model from a full radiative transfer model to a vibrational–rotational band model. It is shown that the derived temperature should be interpreted as a column-integrated property as opposed to a temperature at a specified altitude without utilization of a priori information of the thermospheric temperature profile. The two-channel-ratio approach is demonstrated using NASA GOLD Level 1C disk emission data for the period of 2–8 November 2018 during which a moderate geomagnetic storm has occurred. Due to the lack of independent thermospheric temperature observations, the efficacy of the approach is validated through comparisons of the column-integrated temperature derived from GOLD Level 1C data with the GOLD Level 2 temperature product as well as temperatures from first principle and empirical models. The storm-time thermospheric response manifested in the column-integrated temperature is also shown to corroborate well with hemispherically integrated Joule heating rates, ESA SWARM mass density at 460 km, and GOLD Level 2 column O/N2 ratio.

Description: Non-Gaussian forecast error is a challenge for ensemble-based data assimilation (DA), particularly for more nonlinear convective dynamics. In this study, we investigate the degree of the non-Gaussianity of forecast error distributions at 1 km resolution using a 1000-member ensemble Kalman filter, and how it is affected by the DA update frequency and observation number. Regional numerical weather prediction experiments are performed with the SCALE (Scalable Computing for Advanced Library and Environment) model and the LETKF (local ensemble transform Kalman filter) assimilating phased array radar observations every 30 s. The results show that non-Gaussianity develops rapidly within convective clouds and is sensitive to the DA frequency and the number of assimilated observations. The non-Gaussianity is reduced by up to 40 % when the assimilation window is shortened from 5 min to 30 s, particularly for vertical velocity and radar reflectivity.

Description: Data underlying figures 1, 2, 3 and 5. Figure 1: Monthly mass changes of the Greenland Ice Sheet from GRACE/GRACE-FO and SMB-D (2003-2019) Figure 2: Biennial mass balance and its components from GRACE/GRACE-FO and SMB-D (2003-2018) for the Greenland Ice Sheet, along with regional estimates for 2017-2018 for East and West. Figure 3: Rate of mass change for year 2019 from GRACE/GRACE-FO and SMB-D Figure 5: Annual mass balance and its main components from SMB-D (1948-2019) and GRACE/GRACE-FO (2003-2019)

Description: Here, we provide the raw pollen data archived in three Siberian lake sediment cores spanning the mid-Holocene to the present (7.6-0 cal ka BP), from northern typical tundra to southern open larch forest in the Omoloy region. There are three cores: 1. 14-OM-20B, Lat. / °: 70.53, Lon. / °: 132.91, Ele. / m a.s.l.: 52, Modern vegetation: open larch forest, Lake area / km2: 0.26, Maximal depth / m: 3.4 2. 14-OM-02B, Lat. / °: 70.72, Lon. / °: 132.67, Ele. / m a.s.l.: 58, Modern vegetation: forest tundra, Lake area / km2: 0.08, Maximal depth / m: 3.5 3. 14-OM-12A, Lat. / °: 70.96, Lon. / °: 132.57, Ele. / m a.s.l.: 60, Modern vegetation: tundra, Lake area / km2: 0.09, Maximal depth / m: 4.5 Three lake sediment cores, 14OM12A (33 cm long), 14OM02B (49.5 cm long) and 14OM20B (86 cm long), were recovered from three sites using a UWITEC gravity corer (6 cm internal diameter) equipped with a hammer tool in July 2014. From the three cores, 16 bulk organic carbon samples were selected because of the lack of macrofossil remains and radiocarbon dated using accelerator mass spectrometry (AMS) at Poznań radiocarbon laboratory of Adam Mickiewicz University, Poland. In addition, 30 freeze-dried samples per core at 0.25 or 0.5 cm intervals between 0 and 15 cm were analysed for 210Pb/137Cs at the Liverpool University Environmental Radioactivity Laboratory. In this project, we analyse pollen and sedaDNA (Liu et al., 2020; doi:10.5061/dryad.69p8cz900) from three lake sediment cores from the Omoloy region in north-eastern Siberia (northern Yakutia), which are currently surrounded by different vegetation types ranging from typical tundra to open larch forest. First, our aim is to compare sedaDNA with the pollen data to see whether both methods track the same pattern with respect to compositional changes and diversity changes across the northern Russian treeline zone or are complementary to each other. Second, we reconstruct the mid- to late-Holocene changes of vegetation composition along a north–south transect. Third, we use the sedaDNA data to reconstruct variations in species richness and relate this to vegetation and climate change.

Description: Understanding the resilience of African savannas to global change requires quantitative information on long-term vegetation dynamics. Here we present a reconstruction of past vegetation cover of the northern Namibian savanna obtained after applying the REVEALS model to fossil pollen data from Lake Otjikoto. We also present modern pollen and vegetation data used to calculate pollen productivity estimates for the major Namibian savanna taxa Acacia (Senegalia, Vachellia), Combretaceae, Dichrostachys, Grewia and Poaceae. Data were collected at 10 sites along a rainfall gradient in north central Namibia. Modern pollen was extracted from soil samples collected from plots at the different sites. Vegetation data were extracted from satellite images covering a 1.5 km radius from the plots where pollen was collected. The mean cover of the studied taxa was calculated by 100 m rings.

Description: The Late Holocene is a substantial cultural and economic transition in the eastern Eurasian Steppe and Altai Region, but paleoclimate conditions during this time remain unclear. Therefore, we established a high-resolution paleoclimate record from Lake Khar Nuur in the Mongolian Altai, spanning the last 4200 years. Lake Khar Nuur is a high-altitude lake with a small catchment located at 2,486 m a.s.l. (48°37'22.9"N, 88°56'42.5"E). We recovered the sediment core (that we abbreviate KN18) from the deepest part of the lake (49.4 m) in July 2018 using an Uwitec gravity corer. Within the sediment core KN18, a wide array of lake sediment proxies was measured. While total organic carbon (TOC), total nitrogen, bulk δ^13^C~TOC~, δ^15^N and biogenic silica were measured in 2 cm resolution, the elemental composition (log (Ca/Ti) ratio) was measured in 0.5 cm resolution. Additionally, compound-specific hydrogen isotopic composition of _n_-alkanes was measured in 1 cm resolution.