ALBERT

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: 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: 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: 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 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: 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: 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: Data from the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition allowed us to investigate the temporal dynamics of snowfall, snow accumulation and erosion in great detail for almost the whole accumulation season (November 2019 to May 2020). We computed cumulative snow water equivalent (SWE) over the sea ice based on snow depth and density retrievals from a SnowMicroPen and approximately weekly measured snow depths along fixed transect paths. We used the derived SWE from the snow cover to compare with precipitation sensors installed during MOSAiC. The data were also compared with ERA5 reanalysis snowfall rates for the drift track. We found an accumulated snow mass of 38 m SWE between the end of October 2019 and end of April 2020. The initial SWE over first-year ice relative to second-year ice increased from 50 % to 90 % by end of the investigation period. Further, we found that the Vaisala Present Weather Detector 22, an optical precipitation sensor, and installed on a railing on the top deck of research vessel Polarstern, was least affected by blowing snow and showed good agreements with SWE retrievals along the transect. On the contrary, the OTT Pluvio2 pluviometer and the OTT Parsivel2 laser disdrometer were largely affected by wind and blowing snow, leading to too high measured precipitation rates. These are largely reduced when eliminating drifting snow periods in the comparison. ERA5 reveals good timing of the snowfall events and good agreement with ground measurements with an overestimation tendency. Retrieved snowfall from the ship-based Ka-band ARM zenith radar shows good agreements with SWE of the snow cover and differences comparable to those of ERA5. Based on the results, we suggest the Ka-band radar-derived snowfall as an upper limit and the present weather detector on RV Polarstern as a lower limit of a cumulative snowfall range. Based on these findings, we suggest a cumulative snowfall of 72 to 107 m and a precipitation mass loss of the snow cover due to erosion and sublimation as between 47 % and 68 %, for the time period between 31 October 2019 and 26 April 2020. Extending this period beyond available snow cover measurements, we suggest a cumulative snowfall of 98-114 m.

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: 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: 〈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: 〈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: 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: 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: 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: 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: 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: 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: 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: 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: 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: 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: 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: 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: 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: 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 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: 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 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: 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 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: 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: 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: 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: 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 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: 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: 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: A new 21.3m firn core was drilled in 2015 at a coastal Antarctic high-accumulation site in Adélie Land (66.78◦ S; 139.56◦ E, 602 m a.s.l.), named Terre Adélie 192A (TA192A). The mean isotopic values (−19.3 ‰ ± 3.1 ‰ for δ18O and 5.4 ‰±2.2 ‰ for deuterium excess) are consistent with other coastal Antarctic values. No significant isotope–temperature relationship can be evidenced at any timescale. This rules out a simple interpretation in terms of local temperature. An observed asymmetry in the δ18O seasonal cycle may be explained by the precipitation of air masses coming from the eastern and western sectors in autumn and winter, recorded in the d-excess signal showing outstanding values in austral spring versus autumn. Significant positive trends are observed in the annual d-excess record and local sea ice extent (135–145◦ E) over the period 1998–2014. However, process studies focusing on resulting isotopic compositions and particularly the deuterium excess–δ18O relationship, evidenced as a potential fingerprint of moisture origins, as well as the collection of more isotopic measurements in Adélie Land are needed for an accurate interpretation of our signals.

Description: The timing and intensity of snowmelt processes on sea ice are key drivers determining the seasonal sea-ice energy and mass budgets. In the Arctic, satellite passive microwave and radar observations have revealed a trend towards an earlier snowmelt onset during the last decades, which is an important aspect of Arctic amplification and sea ice decline. Around Antarctica, snowmelt on perennial ice is weak and very different than in the Arctic, with most snow surviving the summer. Here we compile time series of snowmelt-onset dates on seasonal and perennial Antarctic sea ice from 1992 to 2014/15 using active microwave observations from European Remote Sensing Satellite (ERS-1/2), Quick Scatterometer (QSCAT) and Advanced Scatterometer (ASCAT) radar scatterometers. We define two snowmelt transition stages: A weak backscatter rise indicating the initial warming and destructive metamorphism of the snowpack (pre-melt), followed by a rapid backscatter rise indicating the onset of thaw-freeze cycles (snowmelt). Results show large interannual variability with an average pre-melt onset date of 29 November and melt onset of 10 December, respectively, on perennial ice, without any significant trends over the study period, consistent with the small trends of Antarctic sea ice extent. There was a latitudinal gradient from early snowmelt onsets in mid-November in the northern Weddell Sea to late (end-December) or even absent snowmelt conditions in the southern Weddell Sea. We show that QSCAT Ku-band (13.4 GHz signal frequency) derived pre-melt and snowmelt onset dates are earlier by 20 and 18 days, respectively, than ERS and ASCAT C-band (5.6 GHz) derived dates. This offset has been considered when constructing the time series. Snowmelt onset dates from passive microwave observations (37 GHz) are later by 14 and 6 days than those from the scatterometers, respectively. Based on these characteristic differences between melt onset dates observed by different microwave wavelengths, we developed a conceptual model which illustrates how the seasonal evolution of snow temperature profiles may affect different microwave bands with different penetration depths. These suggest that future multi-frequency active/passive microwave satellite missions could be used to resolve melt processes throughout the vertical snow column of thick snow on perennial Antarctic sea ice.

Description: Although quantitative isotope data from speleothems has been used to evaluate isotope-enabled model simulations, currently no consensus exists regarding the most appropriate methodology through which to achieve this. A number of modelling groups will be running isotope-enabled palaeoclimate simulations in the framework of the Coupled Model Intercomparison Project Phase 6, so it is timely to evaluate different approaches to using the speleothem data for data–model comparisons. Here, we illustrate this using 456 globally distributed speleothem δ18O records from an updated version of the Speleothem Isotopes Synthesis and Analysis (SISAL) database and palaeoclimate simulations generated using the ECHAM5-wiso isotope-enabled atmospheric circulation model. We show that the SISAL records reproduce the first-order spatial patterns of isotopic variability in the modern day, strongly supporting the application of this dataset for evaluating model-derived isotope variability into the past. However, the discontinuous nature of many speleothem records complicates the process of procuring large numbers of records if data–model comparisons are made using the traditional approach of comparing anomalies between a control period and a given palaeoclimate experiment. To circumvent this issue, we illustrate techniques through which the absolute isotope values during any time period could be used for model evaluation. Specifically, we show that speleothem isotope records allow an assessment of a model’s ability to simulate spatial isotopic trends. Our analyses provide a protocol for using speleothem isotope data for model evaluation, including screening the observations to take into account the impact of speleothem mineralogy on δ18O values, the optimum period for the modern observational baseline and the selection of an appropriate time window for creating means of the isotope data for palaeo-time-slices.

Description: We present here the first results, for the preindustrial and mid-Holocene climatological periods, of the newly developed isotope-enhanced version of the fully coupled Earth system model MPI-ESM, called hereafter MPI-ESM-wiso. The water stable isotopes H16O, H18O and HDO have been implemented into all components of the coupled model setup. The mid-Holocene provides the opportunity to evaluate the model response to changes in the seasonal and latitudinal distribution of insolation induced by different orbital forcing conditions. The results of our equilibrium simulations allow us to evaluate the performance of the isotopic model in simulating the spatial and temporal variations of water isotopes in the different compartments of the hydrological system for warm climates. For the preindustrial climate, MPI-ESM-wiso reproduces very well the observed spatial distribution of the isotopic content in precipitation linked to the spatial variations in temperature and precipitation rate. We also find a good model–data agreement with the observed distribution of isotopic composition in surface seawater but a bias with the presence of surface seawater that is too 18O-depleted in the Arctic Ocean. All these results are improved compared to the previous model version ECHAM5/MPIOM. The spatial relationships of water isotopic composition with temperature, precipitation rate and salinity are consistent with observational data. For the preindustrial climate, the interannual relationships of water isotopes with temperature and salinity are globally lower than the spatial ones, consistent with previous studies. Simulated results under mid-Holocene conditions are in fair agreement with the isotopic measurements from ice cores and continental speleothems. MPI-ESM-wiso simulates a decrease in the isotopic composition of precipitation from North Africa to the Tibetan Plateau via India due to the enhanced monsoons during the mid-Holocene. Over Greenland, our simulation indicates a higher isotopic composition of precipitation linked to higher summer temperature and a reduction in sea ice, shown by positive isotope–temperature gradient. For the Antarctic continent, the model simulates lower isotopic values over the East Antarctic plateau, linked to the lower temperatures during the mid-Holocene period, while similar or higher isotopic values are modeled over the rest of the continent. While variations of isotopic contents in precipitation over West Antarctica between mid-Holocene and preindustrial periods are partly controlled by changes in temperature, the transport of relatively 18O-rich water vapor near the coast to the western ice core sites could play a role in the final isotopic composition. So, more caution has to be taken about the reconstruction of past temperature variations during warm periods over this area. The coupling of such a model with an ice sheet model or the use of a zoomed grid centered on this region could help to better describe the role of the water vapor transport and sea ice around West Antarctica. The reconstruction of past salinity through isotopic content in sea surface waters can be complicated for regions with strong ocean dynamics, variations in sea ice regimes or significant changes in freshwater budget, giving an extremely variable relationship between the isotopic content and salinity of ocean surface waters over small spatial scales. These complicating factors demonstrate the complexity of interpreting water isotopes as past climate signals of warm periods like the mid-Holocene. A systematic isotope model intercomparison study for further insights on the model dependency of these results would be beneficial.

Description: In this study, we present results obtained from modelling the mid-Pliocene warm period using the Community Earth System Models (COSMOS, version: COSMOS-landveg r2413, 2009) with the two different sets of boundary conditions prescribed for the two phases of the Pliocene Model Intercomparison Project (PlioMIP). Boundary conditions, model forcing, and modelling methodology of the two phases of PlioMIP, tagged PlioMIP1 and PlioMIP2,differ considerably in palaeogeography, in particular with regards to the state of ocean gateways, ice-masks, vegetation and topography. Further differences between model setups as suggested for PlioMIP1 and PlioMIP2 consider updates to the concentration of atmospheric carbon dioxide (CO2), that is specified as 405 and 400 parts per million by volume (ppmv) for PlioMIP1 and PlioMIP2, respectively, as well as minor differences in the concentrations of methane (CH4) and nitrous oxide (N2O) due to changes in the protocol of the Paleoclimate Model Intercomparison Project (PMIP) from phase 3 to phase 4. With this manuscript, we bridge the gap between our contributions to PlioMIP1 (Stepanek and Lohmann, 2012) and PlioMIP2 (Stepanek et al., 2019). We highlight some of the effects that differences in the chosen Mid-Pliocene model setup (PlioMIP2 vs. PlioMIP1) have on the climate state as derived with the COSMOS, as this information will be valuable in the framework of the model-model and model-data-comparison within PlioMIP2. We evaluate the model sensitivity to improved mid-Pliocene boundary conditions using PlioMIP’s core mid-Pliocene experiments for PlioMIP1 and PlioMIP2, and present further simulations where we test model sensitivity to variations in palaeogeography, orbit and concentration of CO2. Firstly,we highlight major changes in boundary conditions from PlioMIP1 to PlioMIP2 and also the limitations recorded from the initial effort. The results derived from of our simulations show that COSMOS simulates a mid-Pliocene climate state that is 0.08 K colder in PlioMIP2, if compared to PlioMIP1. On one hand, high-latitude warming,which is supported by proxy evidence of the mid-Pliocene, is underestimated in simulations of both PlioMIP1 andPlioMIP2. On the other hand, spatial variations in surface air temperature (SAT), sea surface temperature (SST) as well as the distribution of sea ice suggest improvement of simulated SAT and SST in PlioMIP2 if employing the updated palaeogeography. The PlioMIP2 Mid-Pliocene simulation produces warmer SSTs in the Arctic and North Atlantic Ocean than derived from the respective PlioMIP1 climate state. The difference in prescribed CO2accountsfor 1.1 K of warming in the Arctic, leading to an ice-free summer in the PlioMIP1 simulation, and a quasi-ice-free summer in PlioMIP2. Furthermore, employing different orbital forcings in simulating the Mid-Pliocene lead to pronounced annual and seasonal variations, which is not accounted for by marine and terrestrial reconstruction of the time-slice.

Description: The Antarctic temperature changes over the past millennia remain more uncertain than in many other continental regions. This has several origins: (1) the number of high-resolution ice cores is small, in particular on the East Antarctic plateau and in some coastal areas in East Antarctica; (2) the short and spatially sparse instrumental records limit the calibration period for reconstructions and the assessment of the methodologies; (3) the link between isotope records from ice cores and local climate is usually complex and dependent on the spatial scales and timescales investigated. Here, we use climate model results, pseudo-proxy experiments and data assimilation experiments to assess the potential for reconstructing the Antarctic temperature over the last 2 millennia based on a new database of stable oxygen isotopes in ice cores compiled in the frame- work of Antarctica2k (Stenni et al., 2017). The well-known covariance between δ18O and temperature is reproduced in the two isotope-enabled models used (ECHAM5/MPI-OM and ECHAM5-wiso), but is generally weak over the different Antarctic regions, limiting the skill of the reconstructions. Furthermore, the strength of the link displays large variations over the past millennium, further affecting the potential skill of temperature reconstructions based on statistical methods which rely on the assumption that the last decades are a good estimate for longer temperature reconstructions. Using a data assimilation technique allows, in theory, for changes in the δ18O–temperature link through time and space to be taken into account. Pseudoproxy experiments confirm the benefits of using data assimilation methods instead of statistical methods that provide reconstructions with unrealistic variances in some Antarctic subregions. They also confirm that the relatively weak link between both variables leads to a limited potential for reconstructing temperature based on δ18O. However, the reconstruction skill is higher and more uniform among reconstruction methods when the reconstruction target is the Antarctic as a whole rather than smaller Antarctic subregions. This consistency between the methods at the large scale is also observed when reconstructing temperature based on the real δ18O regional composites of Stenni et al. (2017). In this case, temperature reconstructions based on data assimilation confirm the long-term cooling over Antarctica during the last millennium, and the later onset of anthropogenic warming compared with the simulations without data assimilation, which is especially visible in West Antarctica. Data assimilation also allows for models and direct observations to be reconciled by reproducing the east–west contrast in the recent temperature trends. This recent warming pattern is likely mostly driven by internal variability given the large spread of individual Paleoclimate Modelling Intercomparison Project (PMIP)/Coupled Model Intercomparison Project (CMIP) model realizations in simulating it. As in the pseudoproxy framework, the reconstruction methods perform differently at the subregional scale, especially in terms of the variance of the time series produced. While the potential benefits of using a data assimilation method instead of a statistical method have been highlighted in a pseudoproxy framework, the instrumental series are too short to confirm this in a realistic setup.

Description: Ice-wedge polygons are common features of lowland tundra in the continuous permafrost zone and prone to rapid degradation through melting of ground ice. There are many interrelated processes involved in ice-wedge thermokarst and it is a major challenge to quantify their influence on the stability of the permafrost underlying the landscape. In this study we used a numerical modelling approach to investigate the degradation of ice wedges with a focus on the influence of hydrological conditions. Our study area was Samoylov Island in the Lena River delta of northern Siberia, for which we had in situ measurements to evaluate the model. The tailored version of the CryoGrid 3 land surface model was capable of simulating the changing microtopography of polygonal tundra and also regarded lateral fluxes of heat, water, and snow. We demonstrated that the approach is capable of simulating ice-wedge degradation and the associated transition from a low-centred to a high-centred polygonal microtopography. The model simulations showed ice-wedge degradation under recent climatic conditions of the study area, irrespective of hydrological conditions. However, we found that wetter conditions lead to an earlier onset of degradation and cause more rapid ground subsidence. We set our findings in correspondence to observed types of ice-wedge polygons in the study area and hypothesized on remaining discrepancies between modelled and observed ice-wedge thermokarst activity. Our quantitative approach provides a valuable complement to previous, more qualitative and conceptual, descriptions of the possible pathways of ice-wedge polygon evolution. We concluded that our study is a blueprint for investigating thermokarst landforms and marks a step forward in understanding the complex interrelationships between various processes shaping ice-rich permafrost landscapes.

Description: Warming of the Arctic led to an increase in permafrost temperatures by about 0.3 �C during the last decade. Permafrost warming is associated with increasing sediment water content, permeability, and diffusivity and could in the long term alter microbial community composition and abundance even before permafrost thaws. We studied the long-term effect (up to 2500 years) of submarine permafrost warming on microbial communities along an onshore–offshore transect on the Siberian Arctic Shelf displaying a natural temperature gradient of more than 10 �C. We analysed the in situ development of bacterial abundance and community composition through total cell counts (TCCs), quantitative PCR of bacterial gene abundance, and amplicon sequencing and correlated the microbial community data with temperature, pore water chemistry, and sediment physicochemical parameters. On timescales of centuries, permafrost warming coincided with an overall decreasing microbial abundance, whereas millennia after warming microbial abundance was similar to cold onshore permafrost. In addition, the dissolved organic carbon content of all cores was lowest in submarine permafrost after millennial-scale warming. Based on correlation analysis, TCC, unlike bacterial gene abundance, showed a significant rank-based negative correlation with increasing temperature, while bacterial gene copy numbers showed a strong negative correlation with salinity. Bacterial community composition correlated only weakly with temperature but strongly with the pore water stable isotopes �18O and �D, as well as with depth. The bacterial community showed substantial spatial variation and an overall dominance of Actinobacteria, Chloroflexi, Firmicutes, Gemmatimonadetes, and Proteobacteria, which are amongst the microbial taxa that were also found to be active in other frozen permafrost environments. We suggest that, millennia after permafrost warming by over 10 �C, microbial community composition and abundance show some indications for proliferation but mainly reflect the sedimentation history and paleoenvironment and not a direct effect through warming.

Description: The Global Ocean Data Analysis Project (GLODAP) is a synthesis effort providing regular compilations of surface to bottom ocean biogeochemical data, with an emphasis on seawater inorganic carbon chemistry and related variables determined through chemical analysis of water samples. This update of GLODAPv2, v2.2019, adds data from 116 cruises to the previous version, extending its coverage in time from 2013 to 2017, while also adding some data from prior years. GLODAPv2.2019 includes measurements from more than 1.1 million water samples from the global oceans collected on 840 cruises. The data for the 12 GLODAP core variables (salinity, oxygen, nitrate, silicate, phosphate, dissolved inorganic carbon, total alkalinity, pH, CFC-11, CFC-12, CFC-113, and CCl4) have undergone extensive quality control, especially systematic evaluation of bias. The data are available in two formats: (i) as submitted by the data originator but updated to WOCE exchange format and (ii) as a merged data product with adjustments applied to minimize bias. These adjustments were derived by comparing the data from the 116 new cruises with the data from the 724 quality-controlled cruises of the GLODAPv2 data product. They correct for errors related to measurement, calibration, and data handling practices, taking into account any known or likely time trends or variations. 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, and 5 % in the halogenated transient tracers. The compilation also includes data for several other variables, such as isotopic tracers. These were not subjected to bias comparison or adjustments. The original data, their documentation and DOI codes are available in the Ocean Carbon Data System of NOAA NCEI (https://www.nodc.noaa.gov/ocads/oceans/GLODAPv2_2019/, last access: 17 September 2019). This site also provides access to the merged data product, which is provided as a single global file and as four regional ones – the Arctic, Atlantic, Indian, and Pacific oceans – under https://doi.org/10.25921/xnme-wr20 (Olsen et al., 2019). The product files also include significant ancillary and approximated data. These were obtained by interpolation of, or calculation from, measured data. This paper documents the GLODAPv2.2019 methods and provides a broad overview of the secondary quality control procedures and results.

Description: We present a Lagrangian convective transport scheme developed for global chemistry and transport models, which considers the variable residence time that an air parcel spends in convection. This is particularly important for accurately simulating the tropospheric chemistry of short-lived species, e.g., for determining the time available for heterogeneous chemical processes on the surface of cloud droplets. In current Lagrangian convective transport schemes air parcels are stochastically redistributed within a fixed time step according to estimated probabilities for convective entrainment as well as the altitude of detrainment. We introduce a new scheme that extends this approach by modeling the variable time that an air parcel spends in convection by estimating vertical updraft velocities. Vertical updraft velocities are obtained by combining convective mass fluxes from meteorological analysis data with a parameterization of convective area fraction profiles. We implement two different parameterizations: a parameterization using an observed constant convective area fraction profile and a parameterization that uses randomly drawn profiles to allow for variability. Our scheme is driven by convective mass fluxes and detrainment rates that originate from an external convective parameterization, which can be obtained from meteorological analysis data or from general circulation models. We study the effect of allowing for a variable time that an air parcel spends in convection by performing simulations in which our scheme is implemented into the trajectory module of the ATLAS chemistry and transport model and is driven by the ECMWF ERA-Interim reanalysis data. In particular, we show that the redistribution of air parcels in our scheme conserves the vertical mass distribution and that the scheme is able to reproduce the convective mass fluxes and detrainment rates of ERA-Interim. We further show that the estimated vertical updraft velocities of our scheme are able to reproduce wind profiler measurements performed in Darwin, Australia, for velocities larger than 0.6 m s−1. SO2 is used as an example to show that there is a significant effect on species mixing ratios when modeling the time spent in convective updrafts compared to a redistribution of air parcels in a fixed time step. Furthermore, we perform long-time global trajectory simulations of radon-222 and compare with aircraft measurements of radon activity.

Description: Holocene permafrost from ice wedge polygons in the vicinity of large seabird breeding colonies in the Thule District, NW Greenland, was drilled to explore the relation between permafrost aggradation and seabird presence. The latter is reliant on the presence of the North Water Polynya (NOW) in the northern Baffin Bay. The onset of peat accumulation associated with the arrival of little auks (Alle alle) in a breeding colony at Annikitisoq, north of Cape York, is radiocarbon-dated to 4400 cal BP. A thick-billed murre (Uria lomvia) colony on Appat (Saunders Island) in the mouth of the Wolstenholme Fjord started 5650 cal BP. Both species provide marine-derived nutrients (MDNs) that fertilize vegetation and promote peat growth. The geochemical signature of organic matter left by the birds is traceable in the frozen Holocene peat. The peat accumulation rates at both sites are highest after the onset, decrease over time, and were about 2-times faster at the little auk site than at the thick-billed murre site. High accumulation rates induce shorter periods of organic matter (OM) decomposition before it enters the perennially frozen state. This is seen in comparably high C=N ratios and less depleted 13C, pointing to a lower degree of OM decomposition at the little auk site, while the opposite pattern can be discerned at the thick-billed murre site. Peat accumulation rates correspond to 15N trends, where decreasing accumulation led to increasing depletion in 15N as seen in the little-auk-related data. In contrast, the more decomposed OM of the thick-billed murre site shows almost stable 15N. Late Holocene wedge ice fed by cold season precipitation was studied at the little auk site and provides the first stable-water isotopic record from Greenland with mean 18O of 8:00:8, mean D of 36:25:7, mean d excess of 7:70:7, and a 18O-D slope of 7.27, which is close to those of the modern Thule meteoric water line. The syngenetic ice wedge polygon development is mirrored in testacean records of the little auk site and delineates polygon low-center, dry-out, and polygon-high-center stages. The syngenetic permafrost formation directly depending on peat growth (controlled by bird activity) falls within the period of neoglacial cooling and the establishment of the NOW, thus indirectly following the Holocene climate trends.

Description: 〈b〉Ensemble forecasts of air quality in eastern China – Part 1: Model description and implementation of the MarcoPolo–Panda prediction system, version 1〈/b〉〈br〉 Guy P. Brasseur, Ying Xie, Anna Katinka Petersen, Idir Bouarar, Johannes Flemming, Michael Gauss, Fei Jiang, Rostislav Kouznetsov, Richard Kranenburg, Bas Mijling, Vincent-Henri Peuch, Matthieu Pommier, Arjo Segers, Mikhail Sofiev, Renske Timmermans, Ronald van der A, Stacy Walters, Jianming Xu, and Guangqiang Zhou〈br〉 Geosci. Model Dev., 12, 33-67, https://doi.org/10.5194/gmd-12-33-2019, 2019〈br〉 An operational multi-model forecasting system for air quality provides daily forecasts of ozone, nitrogen oxides, and particulate matter for 37 urban areas in China. The paper presents an intercomparison of the different forecasts performed during a specific period of time and highlights recurrent differences between the model output. Pathways to improve the forecasts by the multi-model system are suggested.

Description: 〈b〉RandomFront 2.3: a physical parameterisation of fire spotting for operational fire spread models – implementation in WRF-SFIRE and response analysis with LSFire+〈/b〉〈br〉 Andrea Trucchia, Vera Egorova, Anton Butenko, Inderpreet Kaur, and Gianni Pagnini〈br〉 Geosci. Model Dev., 12, 69-87, https://doi.org/10.5194/gmd-12-69-2019, 2019〈br〉 Wildfires are a concrete problem and impact on human life, property and the environment. An extremely dangerous phenomenon is so-called 〈q〉fire spotting〈/q〉, i.e., the generation of secondary ignitions responsible for dangerous flare-ups during wildfires. The aim of this research was to improve the tools used for risk management through the inclusion of fire spotting in operational wildfire simulators used by forest service agencies.

Description: 〈b〉Vertical characterization of aerosol optical properties and brown carbon in winter in urban Beijing, China〈/b〉〈br〉 Conghui Xie, Weiqi Xu, Junfeng Wang, Qingqing Wang, Dantong Liu, Guiqian Tang, Ping Chen, Wei Du, Jian Zhao, Yingjie Zhang, Wei Zhou, Tingting Han, Qingyun Bian, Jie Li, Pingqing Fu, Zifa Wang, Xinlei Ge, James Allan, Hugh Coe, and Yele Sun〈br〉 Atmos. Chem. Phys., 19, 165-179, https://doi.org/10.5194/acp-19-165-2019, 2019〈br〉 We present the first simultaneous real-time online measurements of aerosol optical properties at ground level and at 260 m on a meteorological tower in urban Beijing in winter. The vertical similarities and differences in scattering and absorption coefficients were characterized. The increases in MAC of BC were mainly associated with the coating materials on rBC. Coal combustion was the dominant source contribution of brown carbon followed by biomass burning and SOA in winter in Beijing.

Description: 〈b〉Dynamic changes in optical and chemical properties of tar ball aerosols by atmospheric photochemical aging〈/b〉〈br〉 Chunlin Li, Quanfu He, Julian Schade, Johannes Passig, Ralf Zimmermann, Daphne Meidan, Alexander Laskin, and Yinon Rudich〈br〉 Atmos. Chem. Phys., 19, 139-163, https://doi.org/10.5194/acp-19-139-2019, 2019〈br〉 〈p〉Following wood pyrolysis, tar ball aerosols were laboratory generated from wood tar separated into polar and nonpolar phases. Chemical information of fresh tar balls was obtained from a high-resolution time-of-flight aerosol mass spectrometer (HR-ToF-AMS) and single-particle laser desorption/resonance enhanced multiphoton ionization mass spectrometry (SP-LD-REMPI-MS). Their continuous refractive index (RI) between 365 and 425 nm was retrieved using a broadband cavity enhanced spectroscopy (BBCES). Dynamic changes in the optical and chemical properties for the nonpolar tar ball aerosols in 〈span〉NO〈sub〉〈i〉x〈/i〉〈/sub〉〈/span〉-dependent photochemical process were investigated in an oxidation flow reactor (OFR). Distinct differences in the chemical composition of the fresh polar and nonpolar tar aerosols were identified. Nonpolar tar aerosols contain predominantly high-molecular weight unsubstituted and alkyl-substituted polycylic aromatic hydrocarbons (PAHs), while polar tar aerosols consist of a high number of oxidized aromatic substances (e.g., methoxy-phenols, benzenediol) with higher O : C ratios and carbon oxidation states. Fresh tar balls have light absorption characteristics similar to atmospheric brown carbon (BrC) aerosol with higher absorption efficiency towards the UV wavelengths. The average retrieved RI is 〈span〉1.661+0.020〈i〉i〈/i〉〈/span〉 and 〈span〉1.635+0.003〈i〉i〈/i〉〈/span〉 for the nonpolar and polar tar aerosols, respectively, with an absorption Ångström exponent (AAE) between 5.7 and 7.8 in the detected wavelength range. The RI fits a volume mixing rule for internally mixed nonpolar/polar tar balls. The RI of the tar ball aerosols decreased with increasing wavelength under photochemical oxidation. Photolysis by UV light (254 nm), without strong oxidants in the system, slightly decreased the RI and increased the oxidation state of the tar balls. Oxidation under varying OH exposure levels and in the absence of 〈span〉NO〈sub〉〈i〉x〈/i〉〈/sub〉〈/span〉 diminished the absorption (bleaching) and increased the O : C ratio of the tar balls. The photobleaching via OH radical initiated oxidation is mainly attributed to decomposition of chromophoric aromatics, nitrogen-containing organics, and high-molecular weight components in the aged particles. Photolysis of nitrous oxide (〈span〉N〈sub〉2〈/sub〉O〈/span〉) was used to simulate 〈span〉NO〈sub〉〈i〉x〈/i〉〈/sub〉〈/span〉-dependent photochemical aging of tar balls in the OFR. Under high-〈span〉NO〈sub〉〈i〉x〈/i〉〈/sub〉〈/span〉 conditions with similar OH exposure, photochemical aging led to the formation of organic nitrates, and increased both oxidation degree and light absorption for the aged tar ball aerosols. These observations suggest that secondary organic nitrate formation counteracts the bleaching by OH radical photooxidation to eventually regain some absorption of the aged tar ball aerosols. The atmospheric implication and climate effects from tar balls upon various oxidation processes are briefly discussed.〈/p〉

Description: 〈b〉Atmospheric new particle formation in China〈/b〉〈br〉 Biwu Chu, Veli-Matti Kerminen, Federico Bianchi, Chao Yan, Tuukka Petäjä, and Markku Kulmala〈br〉 Atmos. Chem. Phys., 19, 115-138, https://doi.org/10.5194/acp-19-115-2019, 2019〈br〉 The characteristics of new particle formation (NPF) in China, including frequency, formation rate, and particle growth rate, were summarized comprehensively and were compared among observations in different environments. The interactions between air pollution and NPF are discussed, as well as the possible reasons for more frequent NPF under heavy pollution conditions than in our current understanding. Significant and future research directions for NPF studies in China are also summarized.

Description: 〈b〉Dynamic changes in optical and chemical properties of tar ball aerosols by atmospheric photochemical aging〈/b〉〈br〉 Chunlin Li, Quanfu He, Julian Schade, Johannes Passig, Ralf Zimmermann, Daphne Meidan, Alexander Laskin, and Yinon Rudich〈br〉 Atmos. Chem. Phys., 19, 139-163, https://doi.org/10.5194/acp-19-139-2019, 2019〈br〉 〈p〉Following wood pyrolysis, tar ball aerosols were laboratory generated from wood tar separated into polar and nonpolar phases. Chemical information of fresh tar balls was obtained from a high-resolution time-of-flight aerosol mass spectrometer (HR-ToF-AMS) and single-particle laser desorption/resonance enhanced multiphoton ionization mass spectrometry (SP-LD-REMPI-MS). Their continuous refractive index (RI) between 365 and 425 nm was retrieved using a broadband cavity enhanced spectroscopy (BBCES). Dynamic changes in the optical and chemical properties for the nonpolar tar ball aerosols in 〈span〉NO〈sub〉〈i〉x〈/i〉〈/sub〉〈/span〉-dependent photochemical process were investigated in an oxidation flow reactor (OFR). Distinct differences in the chemical composition of the fresh polar and nonpolar tar aerosols were identified. Nonpolar tar aerosols contain predominantly high-molecular weight unsubstituted and alkyl-substituted polycylic aromatic hydrocarbons (PAHs), while polar tar aerosols consist of a high number of oxidized aromatic substances (e.g., methoxy-phenols, benzenediol) with higher O : C ratios and carbon oxidation states. Fresh tar balls have light absorption characteristics similar to atmospheric brown carbon (BrC) aerosol with higher absorption efficiency towards the UV wavelengths. The average retrieved RI is 〈span〉1.661+0.020〈i〉i〈/i〉〈/span〉 and 〈span〉1.635+0.003〈i〉i〈/i〉〈/span〉 for the nonpolar and polar tar aerosols, respectively, with an absorption Ångström exponent (AAE) between 5.7 and 7.8 in the detected wavelength range. The RI fits a volume mixing rule for internally mixed nonpolar/polar tar balls. The RI of the tar ball aerosols decreased with increasing wavelength under photochemical oxidation. Photolysis by UV light (254 nm), without strong oxidants in the system, slightly decreased the RI and increased the oxidation state of the tar balls. Oxidation under varying OH exposure levels and in the absence of 〈span〉NO〈sub〉〈i〉x〈/i〉〈/sub〉〈/span〉 diminished the absorption (bleaching) and increased the O : C ratio of the tar balls. The photobleaching via OH radical initiated oxidation is mainly attributed to decomposition of chromophoric aromatics, nitrogen-containing organics, and high-molecular weight components in the aged particles. Photolysis of nitrous oxide (〈span〉N〈sub〉2〈/sub〉O〈/span〉) was used to simulate 〈span〉NO〈sub〉〈i〉x〈/i〉〈/sub〉〈/span〉-dependent photochemical aging of tar balls in the OFR. Under high-〈span〉NO〈sub〉〈i〉x〈/i〉〈/sub〉〈/span〉 conditions with similar OH exposure, photochemical aging led to the formation of organic nitrates, and increased both oxidation degree and light absorption for the aged tar ball aerosols. These observations suggest that secondary organic nitrate formation counteracts the bleaching by OH radical photooxidation to eventually regain some absorption of the aged tar ball aerosols. The atmospheric implication and climate effects from tar balls upon various oxidation processes are briefly discussed.〈/p〉

Description: 〈b〉Watershed classification for the Canadian prairie〈/b〉〈br〉 Jared D. Wolfe, Kevin R. Shook, Chris Spence, and Colin J. Whitfield〈br〉 Hydrol. Earth Syst. Sci. Discuss., https//doi.org/10.5194/hess-2018-625,2019〈br〉 〈b〉Manuscript under review for HESS〈/b〉 (discussion: open, 0 comments)〈br〉 Watershed classification can identify regions expected to respond similarly to disturbance. Methods should extend beyond hydrology to include other environmental questions, such as ecology and water quality. We developed a classification for the Canadian Prairie, and identified seven classes defined by watershed characteristics, including elevation, climate, wetland density, and surficial geology. Results provide a basis for evaluating watershed response to land management and climate condition.

Description: 〈b〉Implications of water management representations for watershed hydrologic modeling in the Yakima River basin〈/b〉〈br〉 Jiali Qiu, Qichun Yang, Xuesong Zhang, Maoyi Huang, Jennifer C. Adam, and Keyvan Malek〈br〉 Hydrol. Earth Syst. Sci., 23, 35-49, https://doi.org/10.5194/hess-23-35-2019, 2019〈br〉 Complex water management activities challenge hydrologic modeling. We evaluated how different representations of reservoir operation and agricultural irrigation affect streamflow simulations in the Yakima River basin. Results highlight the importance of the inclusion of reliable reservoir and irrigation information in watershed models for improving watershed hydrology modeling. Models used here are public and hold the promise to benefit water assessment and management in other basins.

Description: 〈b〉Aircraft-based measurements of High Arctic springtime aerosol show evidence for vertically varying sources, transport and composition〈/b〉〈br〉 Megan D. Willis, Heiko Bozem, Daniel Kunkel, Alex K. Y. Lee, Hannes Schulz, Julia Burkart, Amir A. Aliabadi, Andreas B. Herber, W. Richard Leaitch, and Jonathan P. D. Abbatt〈br〉 Atmos. Chem. Phys., 19, 57-76, https://doi.org/10.5194/acp-19-57-2019, 2019〈br〉 The vertical distribution of Arctic aerosol is an important driver of its climate impacts. We present vertically resolved measurements of aerosol composition and properties made in the High Arctic during spring on an aircraft platform. We explore how aerosol properties are related to transport history and show evidence of vertical trends in aerosol sources, transport mechanisms and composition. These results will help us to better understand aerosol–climate interactions in the Arctic.

Description: 〈b〉Estimating the soil N〈sub〉2〈/sub〉O emission intensity of croplands in northwest Europe〈/b〉〈br〉 Vasileios Myrgiotis, Mathew Williams, Robert M. Rees, and Cairistiona F. E. Topp〈br〉 Biogeosciences Discuss., https//doi.org/10.5194/bg-2018-490,2019〈br〉 〈b〉Manuscript under review for BG〈/b〉 (discussion: open, 0 comments)〈br〉 This study focuses on a northwest European cropland region and shows that the type of crop growing on a soil has notable effects on the emission of nitrous oxide (N〈sub〉2〈/sub〉O – a greenhouse gas) from that soil. It was found that N〈sub〉2〈/sub〉O emissions from soils under oilseed cultivation are significantly higher than soils under cereal cultivation. This variation is mostly explained by the fact that oilseeds require more nitrogen (fertiliser) than cereals especially at early crop growth stages.

Description: 〈b〉A Double ITCZ Phenomenology of Wind Errors in the Equatorial Atlantic in Seasonal Forecasts with ECMWF Models〈/b〉〈br〉 Jonathan K. P. Shonk, Teferi D. Demissie, and Thomas Toniazzo〈br〉 Atmos. Chem. Phys. Discuss., https//doi.org/10.5194/acp-2018-1316,2019〈br〉 〈b〉Manuscript under review for ACP〈/b〉 (discussion: open, 0 comments)〈br〉 Modern climate models are affected by systematic biases that harm their ability to produce reliable seasonal forecasts and climate projections. In this study, we investigate causes of biases in wind patterns over the tropical Atlantic during northern spring in three related models. We find that the wind biases are associated with an increase in excess rainfall and convergence in the tropical western Atlantic at the start of April, leading to the redirection of trade winds away from the equator.

Description: 〈b〉Spiky Fluctuations and Scaling in High-Resolution EPICA Ice Core Dust Fluxes〈/b〉〈br〉 Shaun Lovejoy and Fabrice Lambert〈br〉 Clim. Past Discuss., https//doi.org/10.5194/cp-2018-171,2019〈br〉 〈b〉Manuscript under review for CP〈/b〉 (discussion: open, 0 comments)〈br〉 We analyze the statistical properties of the eight past glacial-interglacial cycles as well as subsections of a generic glacial cycle using the high-resolution dust flux dataset from the Antarctic EPICA Dome C ice core. We show that the high southern latitude climate during glacial maxima, interglacial, and glacial inception is generally more stable but more drought-prone than during mid-glacial conditions.

Description: 〈b〉Nutrient distribution and nitrogen and oxygen isotopic composition of nitrate in water masses of the subtropical South Indian Ocean〈/b〉〈br〉 Natalie C. Harms, Niko Lahajnar, Birgit Gaye, Tim Rixen, Kirstin Dähnke, Markus Ankele, Ulrich Schwarz-Schampera, and Kay-Christian Emeis〈br〉 Biogeosciences Discuss., https//doi.org/10.5194/bg-2018-511,2019〈br〉 〈b〉Manuscript under review for BG〈/b〉 (discussion: open, 0 comments)〈br〉 The Indian Ocean subtropical gyre is a large oligotrophic area that is likely to adjust to continued warming by increasing stratification, reduced nutrient supply, and decreasing biological production. In this study, we investigated concentrations of nutrients and stable isotopes of nitrate. We determine the lateral influence of water masses entering the gyre from the northern Indian Ocean and from the Southern Ocean and quantify the input of nitrogen by N〈sub〉2〈/sub〉-fixation into the surface layer.

Description: 〈b〉New particle formation in the marine atmosphere during seven cruise campaigns〈/b〉〈br〉 Yujiao Zhu, Kai Li, Yanjie Shen, Yang Gao, Xiaohuan Liu, Yang Yu, Huiwang Gao, and Xiaohong Yao〈br〉 Atmos. Chem. Phys., 19, 89-113, https://doi.org/10.5194/acp-19-89-2019, 2019〈br〉 In this paper, we investigate new particle formation (NPF) events during seven cruises. NPF events were observed on 25 days and were most likely associated with the long-range transport of anthropogenic air pollutants. The relationship between the net generated amount of new particles and their apparent formation rate is established and explained in terms of the roles of different vapor precursors. The survival probability of new particles to CCN size is also discussed.

Description: 〈b〉Aircraft-based measurements of High Arctic springtime aerosol show evidence for vertically varying sources, transport and composition〈/b〉〈br〉 Megan D. Willis, Heiko Bozem, Daniel Kunkel, Alex K. Y. Lee, Hannes Schulz, Julia Burkart, Amir A. Aliabadi, Andreas B. Herber, W. Richard Leaitch, and Jonathan P. D. Abbatt〈br〉 Atmos. Chem. Phys., 19, 57-76, https://doi.org/10.5194/acp-19-57-2019, 2019〈br〉 The vertical distribution of Arctic aerosol is an important driver of its climate impacts. We present vertically resolved measurements of aerosol composition and properties made in the High Arctic during spring on an aircraft platform. We explore how aerosol properties are related to transport history and show evidence of vertical trends in aerosol sources, transport mechanisms and composition. These results will help us to better understand aerosol–climate interactions in the Arctic.

Description: 〈b〉Simulated single-layer forest canopies delay Northern Hemisphere snowmelt〈/b〉〈br〉 Markus Todt, Nick Rutter, Christopher G. Fletcher, and Leanne M. Wake〈br〉 The Cryosphere Discuss., https//doi.org/10.5194/tc-2018-270,2019〈br〉 〈b〉Manuscript under review for TC〈/b〉 (discussion: open, 0 comments)〈br〉 Vegetation is often represented by a single layer in global land models. Studies have found deficient simulation of thermal radiation beneath forest canopies when represented by single-layer vegetation. This study corrects thermal radiation in forests for a global land model using single-layer vegetation in order to assess the effect of deficient thermal radiation on snow cover and snowmelt. Results indicate that single-layer vegetation causes snow in forests being too cold and melting too late.

Description: 〈b〉Brief Communication: Early season snowpack loss and implications for oversnow vehicle recreation travel planning〈/b〉〈br〉 Benjamin J. Hatchett and Hilary G. Eisen〈br〉 The Cryosphere, 13, 21-28, https://doi.org/10.5194/tc-13-21-2019, 2019〈br〉 We examine the timing of early season snowpack relevant to oversnow vehicle (OSV) recreation over the past 3 decades in the Lake Tahoe region (USA). Data from two independent data sources suggest that the timing of achieving sufficient snowpack has shifted later by 2 weeks. Increasing rainfall and more dry days play a role in the later onset. Adaptation strategies are provided for winter travel management planning to address negative impacts of loss of early season snowpack for OSV usage.

Description: 〈b〉Aerosol radiative effects with MACv2〈/b〉〈br〉 Stefan Kinne〈br〉 Atmos. Chem. Phys. Discuss., https//doi.org/10.5194/acp-2018-949,2019〈br〉 〈b〉Manuscript under review for ACP〈/b〉 (discussion: open, 0 comments)〈br〉 〈p〉onthly global maps for aerosol properties of the MACv2 climatology are applied in an off-line radiative transfer model to determine aerosol radiative effects. For details beyond global averages in most cases global maps are presented to visualize regional and seasonal details. Aside from the direct radiative (aerosol presence) effect, including those for aerosol components as extracted from MACv2 aerosol optics, also the major aerosol indirect radiative effect is covered. Hereby, the impact of smaller drops in water clouds due to added anthropogenic aerosol was simulated by applying a satellite retrieval based fit from locally associations between aerosol and drop concentrations over oceans. Present-day anthropogenic aerosols of MACv2 – on a global average basis – reduce the radiative net-fluxes at the top of the atmosphere (TOA) by −1.0 W/m〈sup〉2〈/sup〉 and at the surface by −2.1 W/m〈sup〉2〈/sup〉. Direct cooling contributions are only about half of indirect contributions (−.35 vs −.65) at TOA, but about twice at the surface (−1.45 vs −.65), as solar absorption of the direct effect warms the atmosphere by +1.1 W/m〈sup〉2〈/sup〉. Natural aerosols are on average less absorbing (for a relatively larger solar TOA cooling) and larger in size (now contributing with IR greenhouse warming). Thus, average TOA direct forcing efficiencies for total and anthropogenic aerosol happen to be similar: −11 W/m〈sup〉2〈/sup〉/AOD at all-sky and −24 W/m〈sup〉2〈/sup〉/AOD at clear-sky conditions. The present-day direct impact by all soot (BC) is globally averaged +0.55W/m〈sup〉2〈/sup〉 and at least half of it should be attributed to anthropogenic sources. Hereby any accuracy of anthropogenic impacts, not just for soot, suffers from the limited access to a pre-industrial reference. Anthropogenic uncertainty has a particular strong impact on aerosol indirect effects, which dominate the (TOA) forcing. Accounting for uncertainties in the anthropogenic definition, present-day aerosol forcing is estimated to stay within the −0.7 to −1.6 W/m〈sup〉2〈/sup〉 range, with a best estimate at −1 W/m〈sup〉2〈/sup〉. Calculations with model predicted temporal changes to anthropogenic AOD indicate that qualitatively the anthropogenic aerosol forcing has not changed much over the last decades and is not likely to increase over the next decades, despite strong regional shifts. These regional shifts explain most solar insolation (brightening or dimming) trends that have been observed by ground-based radiation data.〈/p〉

Description: 〈b〉Thermal structure of the mesopause region during the WADIS-2 rocket campaign〈/b〉〈br〉 Raimund Wörl, Boris Strelnikov, Timo P. Viehl, Josef Höffner, Pierre-Dominique Pautet, Michael J. Taylor, Yucheng Zhao, and Franz-Josef Lübken〈br〉 Atmos. Chem. Phys., 19, 77-88, https://doi.org/10.5194/acp-19-77-2019, 2019〈br〉 Simultaneous temperature measurements during the WADIS-2 rocket campaign are used to investigate the thermal structure of the mesopause region. Vertically and horizontally resolved in situ and remote measurements are in good agreement and show dominating long-term and large-scale waves with periods of 24 h and higher tidal harmonics. Only a few gravity waves with periods shorter than 6 h and small amplitudes are there.

Description: 〈b〉Verification of anthropogenic VOC emission inventory through ambient measurements and satellite retrievals〈/b〉〈br〉 Jing Li, Yufang Hao, Maimaiti Simayi, Yuqi Shi, Ziyan Xi, and Shaodong Xie〈br〉 Atmos. Chem. Phys. Discuss., https//doi.org/10.5194/acp-2018-1133,2019〈br〉 〈b〉Manuscript under review for ACP〈/b〉 (discussion: open, 0 comments)〈br〉 We established an emission inventory of anthropogenic VOCs in the Beijing–Tianjin–Hebei region of China. The developed emission inventory was evaluated through ambient measurements and satellite retrieval results. To obtain a more accurate emission inventory, we propose the investigation of the household coal consumption, the adjustment of EFs based on the latest pollution control policies, and the verification of the source profiles of OVOCs and halocarbons.

Description: 〈b〉Classification of aerosol population type and cloud condensation nuclei properties in a coastal California littoral environment using an unsupervised cluster model〈/b〉〈br〉 Samuel A. Atwood, Sonia M. Kreidenweis, Paul J. DeMott, Markus D. Petters, Gavin C. Cornwell, Andrew C. Martin, and Kathryn A. Moore〈br〉 Atmos. Chem. Phys. Discuss., https//doi.org/10.5194/acp-2018-1297,2019〈br〉 〈b〉Manuscript under review for ACP〈/b〉 (discussion: open, 0 comments)〈br〉 This paper presents measurements of aerosol particles at a coastal location. The particles were classified into distinct aerosol types using both microphysical measurements and meteorological information, allowing rapid changes between the aerosol types to be reliably identified. These particles can alter cloud and precipitation processes, and inclusion of the differences between types can improve atmospheric models and remote sensing retrievals in littoral zones.

Description: 〈b〉Development of a balloon-borne instrument for CO〈sub〉2〈/sub〉 vertical profile observations in the troposphere〈/b〉〈br〉 Mai Ouchi, Yutaka Matsumi, Tomoki Nakayama, Kensaku Shimizu, Takehiko Sawada, Toshinobu Machida, Hidekazu Matsueda, Yousuke Sawa, Isamu Morino, Osamu Uchino, Tomoaki Tanaka, and Ryoichi Imasu〈br〉 Atmos. Meas. Tech. Discuss., https//doi.org/10.5194/amt-2018-376,2019〈br〉 〈b〉Manuscript under review for AMT〈/b〉 (discussion: open, 0 comments)〈br〉 A novel, practical observation system for measuring tropospheric carbon dioxide (CO〈sub〉2〈/sub〉) concentrations carried by a small helium-filled balloon (CO〈sub〉2〈/sub〉 sonde), has been developed for the first time. The low-cost CO〈sub〉2〈/sub〉 sondes can potentially be used for frequently measurements of vertical profiles of CO〈sub〉2〈/sub〉 in any parts of the world providing as useful information to understand the global and regional carbon budgets by replenishing the present sparse observation coverage.

Description: 〈b〉Comparison of Different Sequential Assimilation Algorithms for Satellite-derived Leaf Area Index Using the Data Assimilation Research Testbed (lanai)〈/b〉〈br〉 Xiao-Lu Ling, Cong-Bin Fu, Zong-Liang Yang, and Wei-Dong Guo〈br〉 Geosci. Model Dev. Discuss., https//doi.org/10.5194/gmd-2018-232,2019〈br〉 〈b〉Manuscript under review for GMD〈/b〉 (discussion: open, 0 comments)〈br〉 Both observation and simulation can provide the temporal and spatial variation of vegetation characteristic, while they are not satisfactory for understanding the mechanism of the exchange between ecosystems and atmosphere. Data assimilation (DA) can combine the observation and models via mathematical statistical analysis.The results show that the Ensemble Adjust Kalman Filter (EAKF) is the optical algorithm. In addition, models perform better when the DA accept more proportion of observation.

Description: 〈b〉The importance of crystalline phases in ice nucleation by volcanic ash〈/b〉〈br〉 Elena C. Maters, Donald B. Dingwell, Corrado Cimarelli, Dirk Müller, Thomas F. Whale, and Benjamin J. Murray〈br〉 Atmos. Chem. Phys. Discuss., https//doi.org/10.5194/acp-2018-1326,2019〈br〉 〈b〉Manuscript under review for ACP〈/b〉 (discussion: open, 0 comments)〈br〉 This experimental study investigates the influence of volcanic ash chemical composition, crystallinity, and mineralogy on its ability to promote freezing of supercooled water. The results indicate that crystals in ash play a key role in this process, and suggest that alkali and plagioclase feldspars in ash may be highly ice-active. These findings contribute to improving understanding of the potential of ash emissions from different explosive eruptions to impact ice formation in the atmosphere.

Description: 〈b〉Large increases in N〈sub〉cn〈/sub〉 and N〈sub〉ccn〈/sub〉 together with a nucleation-modeparticle pool over the northwestern Pacific Ocean in the spring of 2014〈/b〉〈br〉 Juntao Wang, Yanjie Shen, Kai Li, Yang Gao, Huiwang Gao, and Xiaohong Yao〈br〉 Atmos. Chem. Phys. Discuss., https//doi.org/10.5194/acp-2018-1089,2019〈br〉 〈b〉Manuscript under review for ACP〈/b〉 (discussion: open, 0 comments)〈br〉 In this paper, we studied the spatiotemporal variability of N〈sub〉cn〈/sub〉 and particle number size distributions, as well as N〈sub〉ccn〈/sub〉 and CCN activities over the NWPO in the spring of 2014. Several new findings have been revealed and discussed, e.g., large increases of N〈sub〉cn〈/sub〉 and N〈sub〉ccn〈/sub〉, against historical data and small contributions of dust and BB aerosols to N〈sub〉cn〈/sub〉 and N〈sub〉ccn〈/sub〉 on the monthly time scale, etc. This work may help understand the influence of the current outflow from the Asian continent on the climate.

Description: 〈b〉Analyses of temperature and precipitation in the Indian Jammu and Kashmir region for the 1980–2016 period: implications for remote influence and extreme events〈/b〉〈br〉 Sumira Nazir Zaz, Shakil Ahmad Romshoo, Ramkumar Thokuluwa Krishnamoorthy, and Yesubabu Viswanadhapalli〈br〉 Atmos. Chem. Phys., 19, 15-37, https://doi.org/10.5194/acp-19-15-2019, 2019〈br〉 This paper is of first of its kind for the Jammu and Kashmir (western Himalayas) region, India. It shows the clear relation between the upper tropospheric Rossby wave activity (potential vorticity at the 350 K potential temperature and 200 mb level surface pressure) and the surface weather parameters (e.g., precipitation) over a period of 3 decades during 1980–2016. This indicates that the climatic weather pattern over the Kashmir region is influenced mostly by global climate change processes.

Description: 〈b〉Secondary organic aerosol formation from smoldering and flaming combustion of biomass: a box model parametrization based on volatility basis set〈/b〉〈br〉 Giulia Stefenelli, Jianhui Jiang, Amelie Bertrand, Emily A. Bruns, Simone M. Pieber, Urs Baltensperger, Nicolas Marchand, Sebnem Aksoyoglu, André S. H. Prévôt, Jay G. Slowik, and Imad El Haddad〈br〉 Atmos. Chem. Phys. Discuss., https//doi.org/10.5194/acp-2018-1308,2019〈br〉 〈b〉Manuscript under review for ACP〈/b〉 (discussion: open, 0 comments)〈br〉 Box model simulations, based on the volatility basis set (VBS) approach, of smog chamber wood combustion experiments conducted at different temperatures (−10 °C, 2 °C, 15 °C), emission loads, combustion conditions (flaming and smoldering) and residential stoves fabricated in the last two decades. Novel parameterization methods based on a genetic algorithm approach allowed estimation of precursor class contributions to SOA and evaluation of the effect of emission variability on SOA yield predictions.

Description: 〈b〉Impacts of Horizontal Resolution and Air–Sea Flux Parameterization on the Intensity and Structure of simulated Typhoon Haiyan (2013)〈/b〉〈br〉 Mien-Tze Kueh, Wen-Mei Chen, Yang-Fan Sheng, Simon C. Lin, Tso-Ren Wu, Eric Yen, Yu-Lin Tsai, and Chuan-Yao Lin〈br〉 Nat. Hazards Earth Syst. Sci. Discuss., https//doi.org/10.5194/nhess-2018-333,2019〈br〉 〈b〉Manuscript under review for NHESS〈/b〉 (discussion: open, 0 comments)〈br〉 In this study, we show that both the model horizontal resolution and air-sea flux parameterization can exert large influence on tropical cyclone intensity simulation but with different impacts on wind structures. We highlight the intensification and contraction of TC eyewall in response to the reduction of grid spacing. We also suggest that a well-developed eyewall is more conducive to the positive effect of flux formulas on TC development.

Description: 〈b〉Primary and secondary sources of ambient formaldehyde in the Yangtze River Delta based on OMPS observation〈/b〉〈br〉 Wenjing Su, Cheng Liu, Qihou Hu, Shaohua Zhao, Youwen Sun, Wei Wang, Yizhi Zhu, Jianguo Liu, and Jhoon Kim〈br〉 Atmos. Chem. Phys. Discuss., https//doi.org/10.5194/acp-2018-1192,2019〈br〉 〈b〉Manuscript under review for ACP〈/b〉 (discussion: open, 0 comments)〈br〉 For a better understanding of the pollution and atmospheric chemistry of HCHO, we evaluated its primary and secondary sources in the Yangtze River Delta based on HCHO column density from Ozone Mapping and Profiler Suite (OMPS) and combining with surface in-situ measurements. We found that secondary formation contributed most to ambient HCHO for a long time, but primary emission could be dominant in the winter. The usability of total HCHO as the proxy of VOCs reactivity depends on time scale.

Description: 〈b〉Local and regional contributions to fine particulate matter in the 18 cities of Sichuan Basin, southwestern China〈/b〉〈br〉 Xue Qiao, Hao Guo, Ya Tang, Pengfei Wang, Wenye Deng, Xing Zhao, Jianlin Hu, Qi Ying, and Hongliang Zhang〈br〉 Atmos. Chem. Phys. Discuss., https//doi.org/10.5194/acp-2018-1214,2019〈br〉 〈b〉Manuscript under review for ACP〈/b〉 (discussion: open, 0 comments)〈br〉 A source-oriented version of the CMAQ model was used to quantify contributions from nine regions to PM〈sub〉2.5〈/sub〉 and its components in the 18 cities within the Sichuan Basin. Non-local emissions contribute to 39~66 % and 25~52 % to the citywide average PM〈sub〉2.5〈/sub〉 concentrations of 45~126 and 14~31 μg m〈sup〉−3〈/sup〉 in the winter and summer, respectively. This study demonstrates the importance of joint emission control efforts among cities within the SCB and neighbor regions to the east.

Description: 〈b〉MAX-DOAS measurements of tropospheric NO〈sub〉2〈/sub〉 and HCHO in Nanjing and the comparison to OMI observations〈/b〉〈br〉 Ka Lok Chan, Zhuoru Wang, Aijun Ding, Klaus-Peter Heue, Yicheng Shen, Jing Wang, Feng Zhang, Nan Hao, and Mark Wenig〈br〉 Atmos. Chem. Phys. Discuss., https//doi.org/10.5194/acp-2018-1266,2019〈br〉 〈b〉Manuscript under review for ACP〈/b〉 (discussion: open, 0 comments)〈br〉 The manuscript presents long term observations of atmospheric nitrogen dioxide (NO〈sub〉2〈/sub〉) and formaldehyde (HCHO) in Nanjing using a MAX-DOAS instrument. The measurements were performed from April 2013 to February 2017. The MAX-DOAS measurements of NO〈sub〉2〈/sub〉 and HCHO are used to validate OMI satellite observations and to investigate the influences of region transport of air pollutants on the air quality in Nanjing.

Description: 〈b〉Detection of Outflow of Formaldehyde and Glyoxal from the African continent to the Atlantic Ocean with a MAX-DOAS Instrument〈/b〉〈br〉 Lisa K. Behrens, Andreas Hilboll, Andreas Richter, Enno Peters, Leonardo M. A. Alvarado, Anna B. Kalisz Hedegaard, Folkard Wittrock, John P. Burrows, and Mihalis Vrekoussis〈br〉 Atmos. Chem. Phys. Discuss., https//doi.org/10.5194/acp-2018-1286,2019〈br〉 〈b〉Manuscript under review for ACP〈/b〉 (discussion: open, 0 comments)〈br〉 MAX-DOAS measurements were conducted on the research vessel Maria S. Merian during a cruise from the Azores to South Africa in October 2016. The measurements indicate enhanced levels of HCHO and CHOCHO over the remote Atlantic Ocean, which is unexpected due to their short lifetime. Precursors of these gases or gas-aerosol combinations might be transported. Model simulations indicate potential source regions over the African continent, probably related to biomass burning or biogenic emissions.

Description: 〈b〉Develop a coupled agent-based modeling approach for uncertain water management decisions〈/b〉〈br〉 Jin-Young Hyun, Shih-Yu Huang, Y. C. Ethan Yang, Vincent Tidwell, and Jordan Macknick〈br〉 Hydrol. Earth Syst. Sci. Discuss., https//doi.org/10.5194/hess-2018-555,2019〈br〉 〈b〉Manuscript under review for HESS〈/b〉 (discussion: open, 0 comments)〈br〉 This study applies a 〈q〉two-way〈/q〉 coupled agent-based model (ABM) with a River-Reservoir management model (RiverWare) to address the interaction between human and natural systems using Bayesian Inference (BI) mapping joined with Cost-Loss (CL). The calibration results show that this methodology can capture the historical pattern of both human activities and natural dynamics and outperforms those without using BI and CL.

Description: 〈b〉Potentials and pitfalls of permafrost active layer monitoring using the HVSR method: a case study in Svalbard〈/b〉〈br〉 Andreas Köhler and Christian Weidle〈br〉 Earth Surf. Dynam., 7, 1-16, https://doi.org/10.5194/esurf-7-1-2019, 2019〈br〉 The uppermost part of permanently frozen ground can thaw during summer and refreeze during winter. We use a method based on naturally generated seismic waves to continuously monitor these changes close to the research settlement of Ny-Ålesund in Svalbard between April and August 2016. Our results reveal some potential pitfalls when interpreting temporal variations in the data. However, we show that a careful data analysis makes this method a very useful tool for long-term permafrost monitoring.

Description: 〈b〉Long-profile evolution of transport-limited gravel-bed rivers〈/b〉〈br〉 Andrew D. Wickert and Taylor F. Schildgen〈br〉 Earth Surf. Dynam., 7, 17-43, https://doi.org/10.5194/esurf-7-17-2019, 2019〈br〉 Rivers can raise or lower their beds by depositing or eroding sediments. We combine equations for flow, channel/valley geometry, and gravel transport to learn how climate and tectonics shape down-valley profiles of river-bed elevation. Rivers steepen when they receive more sediment (relative to water) and become straighter with tectonic uplift. Weathering and breakdown of gravel is needed to produce gradually widening river channels with concave-up profiles that are often observed in the field.

Description: 〈b〉Definition differences and internal variability affect the simulated Arctic sea ice melt season〈/b〉〈br〉 Abigail Smith and Alexandra Jahn〈br〉 The Cryosphere, 13, 1-20, https://doi.org/10.5194/tc-13-1-2019, 2019〈br〉 Here we assessed how natural climate variations and different definitions impact the diagnosed and projected Arctic sea ice melt season length using model simulations. Irrespective of the definition or natural variability, the sea ice melt season is projected to lengthen, potentially by as much as 4–5 months by 2100 under the business as usual scenario. We also find that different definitions have a bigger impact on melt onset, while natural variations have a bigger impact on freeze onset.

Description: 〈b〉Three decades of simulated global terrestrial carbon fluxes from a data assimilation system confronted to different periods of observations〈/b〉〈br〉 Karel Castro-Morales, Gregor Schürmann, Christoph Köstler, Christian Rödenbeck, Martin Heimann, and Sönke Zaehle〈br〉 Biogeosciences Discuss., doi:10.5194/bg-2018-517,2019〈br〉 〈b〉Manuscript under review for BG〈/b〉 (discussion: open, 0 comments)〈br〉 To obtain nearly thirty years of global terrestrial carbon fluxes, we simultaneously incorporated in a land surface model three different time periods of two observational data sets: absorbed photosynthetic active radiation and atmospheric CO〈sub〉2〈/sub〉 concentrations. One decade of data is enough to improve the modeled long-term trends and seasonal amplitudes of the assimilated variables, particularly in boreal regions. This model has the potential to provide short-term predictions of land carbon fluxes.

Description: 〈b〉Ostracods as ecological and isotopic indicators of lake water salinity changes: The Lake Van example〈/b〉〈br〉 Jeremy McCormack, Finn Viehberg, Derya Akdemir, Adrian Immenhauser, and Ola Kwiecien〈br〉 Biogeosciences Discuss., doi:10.5194/bg-2018-476,2019〈br〉 〈b〉Manuscript under review for BG〈/b〉 (discussion: open, 0 comments)〈br〉 We juxtapose changes in ostracod taxonomy, morphology (noding) and oxygen (δ〈sup〉18〈/sup〉O) and carbon (δ〈sup〉13〈/sup〉C) isotopic composition for the last 150 kyr with independent low-resolution salinity proxies. We demonstrate that for Lake Van, salinity is the most important factor influencing the composition of the ostracod assemblage and the formation of nodes on the valves of limnocytherinae species. Ostracod δ〈sup〉18〈/sup〉O shows a higher sensibility towards climatic and hydrological variations than the bulk isotopy.

Description: 〈b〉Organic signatures in Pleistocene cherts from Lake Magadi (Kenya), analogs for early Earth hydrothermal deposits〈/b〉〈br〉 Manuel Reinhardt, Walter Goetz, Jan-Peter Duda, Christine Heim, Joachim Reitner, and Volker Thiel〈br〉 Biogeosciences Discuss., https//doi.org/10.5194/bg-2018-513,2019〈br〉 〈b〉Manuscript under review for BG〈/b〉 (discussion: open, 0 comments)〈br〉 Organic matter in Archean hydrothermal cherts may contain molecular traces of early life. Alteration processes during and after deposition, however, may have obliterated potential biosignatures. Our results from modern analog samples (Pleistocene cherts from Lake Magadi, Kenya) show that biomolecules can survive early hydrothermal destruction in the macromolecular fraction of the organic matter. A conservation of molecular biosignatures in Archean hydrothermal cherts therefore seems possible.